4720: [Noip2016]换教室

Description

对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程。在可以选择的课程中,有2n节课程安排在n个时间段上。在第i(1≤i≤n)个时间段上,两节内容相同的课程同时在不同的地点进行,其中,牛牛预先被安排在教室ci上课,而另一节课程在教室di进行。在不提交任何申请的情况下,学生们需要按时间段的顺序依次完成所有的n节安排好的课程。如果学生想更换第i节课程的教室,则需要提出申请。若申请通过,学生就可以在第i个时间段去教室di上课,否则仍然在教室ci上课。由于更换教室的需求太多,申请不一定能获得通过。通过计算,牛牛发现申请更换第i节课程的教室时,申请被通过的概率是一个已知的实数ki,并且对于不同课程的申请,被通过的概率是互相独立的。学校规定,所有的申请只能在学期开始前一次性提交,并且每个人只能选择至多m节课程进行申请。这意味着牛牛必须一次性决定是否申请更换每节课的教室,而不能根据某些课程的申请结果来决定其他课程是否申请;牛牛可以申请自己最希望更换教室的m门课程,也可以不用完这m个申请的机会,甚至可以一门课程都不申请。因为不同的课程可能会被安排在不同的教室进行,所以牛牛需要利用课间时间从一间教室赶到另一间教室。牛牛所在的大学有v个教室,有e条道路。每条道路连接两间教室,并且是可以双向通行的。由于道路的长度和拥堵程度不同,通过不同的道路耗费的体力可能会有所不同。当第i(1≤i≤n-1)节课结束后,牛牛就会从这节课的教室出发,选择一条耗费体力最少的路径前往下一节课的教室。现在牛牛想知道,申请哪几门课程可以使他因在教室间移动耗费的体力值的总和的期望值最小,请你帮他求出这个最小值。

Input

第一行四个整数n,m,v,e。n表示这个学期内的时间段的数量;m表示牛牛最多可以申请更换多少节课程的教室;
v表示牛牛学校里教室的数量;e表示牛牛的学校里道路的数量。
第二行n个正整数,第i(1≤i≤n)个正整数表示c,,即第i个时间段牛牛被安排上课的教室;保证1≤ci≤v。
第三行n个正整数,第i(1≤i≤n)个正整数表示di,即第i个时间段另一间上同样课程的教室;保证1≤di≤v。
第四行n个实数,第i(1≤i≤n)个实数表示ki,即牛牛申请在第i个时间段更换教室获得通过的概率。保证0≤ki≤1。
接下来e行,每行三个正整数aj,bj,wj,表示有一条双向道路连接教室aj,bj,通过这条道路需要耗费的体力值是Wj;
保证1≤aj,bj≤v,1≤wj≤100。
保证1≤n≤2000,0≤m≤2000,1≤v≤300,0≤e≤90000。
保证通过学校里的道路,从任何一间教室出发,都能到达其他所有的教室。
保证输入的实数最多包含3位小数。

Output

输出一行,包含一个实数,四舎五入精确到小数点后恰好2位,表示答案。你的
输出必须和标准输出完全一样才算正确。
测试数据保证四舎五入后的答案和准确答案的差的绝对值不大于4*10^-3。(如果你不知道什么是浮点误差,这段话可以理解为:对于大多数的算法,你可以正常地使用浮点数类型而不用对它进行特殊的处理)

Sample Input

3 2 3 3
2 1 2
1 2 1
0.8 0.2 0.5
1 2 5
1 3 3
2 3 1

Sample Output

2.80

  这道题也可以算得上是经典了。当时考NOIP碰到这道题,当时还壮志满怀,很长的一道题(至少当时看来是这样的)很快就被我读完了,发现可以分成三份:课程、教室,申请、概率,道路、耗费体力最少。
  当时,我的很多同学见此题太长,就根本没有看。此题虽然沾期望,但其实很简单。
----------因为后面扯的东西比较多,就在这里讲一下思路。----------
  我们这样开数组:f[2][2][M],第一维枚举课程但可以“滚动”(cur),第二位判断在第i门课程是否提交申请,第三维给现在已申请的课程计一下数。空间是O(M)的,时间是O(NM)的,最后输出f[cur][j][p](j=0,1; p=0,1,...,m)
  这样,我们可以考虑转移。一个课程如果没有申请,那就是100%的“纯爷们儿”教室c。但如果申请了呢?可以YY一下薛定谔的猫,有概率的事件可以当做对立统一的个体。也就是说一个课程如果申请了,就可以看做k[i]个d教室和1-k[i]个c教室,实际转移时思考一下即可。
  关于方程可以看看我的代码(本文末尾)。
----------好的,思路讲完了。我可以开始扯了。----------
  当时我高兴极了,发现可以把课程当做i,把申请数当做j,把此时是否申请当0/1。当时才学了DP,只会背包(还是最简单的)和floyd,什么LIS、数位、状压(NOIPD2T3ANGRYBIRD)、树形都不知道,更不用说什么概率和期望了……然而就是考了。YJQ表示,不要再逼我吃键盘了。
  关于吃键盘一事,都是因为——YJQ说NOIP绝对不考期望DPNOIP这么水的考试怎么可能考这种东西写个暴力就A题的比赛要是考了这种东西我就直播吃键盘……YJQ其实是一个很谨慎的人,不像lemonoil每天立下flag要直播吃掉∞个键盘。
兜售巧克力味键盘了喂!专门为立flag死掉的人设计的喂~~~”
  所以就是这样的。然而我当时就是想到了如何去做,虽然现在看完题就可以切过去。因为当时我既会背包,也会floyd。不过,天有不测风云,我当时的代码实现能力实在是太弱了……
  弱到什么地步?也就是说,我发现当时k数组定义的是int类型……也许,质问当时的我,我还会义愤填膺地说,int不就是数么整数是数小数难道不是数么?!
  呃呃呃呃呃。
  然而,当我现在做到这道题时,当我已经“掌握”了概率与期望的时候,还是WA。
  恩。请来客忽略我的强制调剂(eyeprotect)。一边调,一会儿AC,一会儿又WA的体验真的是第一回碰到。更令人惊奇的是,中间居然还有CE,让人颇为尴尬。要是NOIP也这样,那恐怕就只有再见OI了。
  最后突然发现,卡了这么久竟是这两个原因。
  1.memset(dis,65,sizeof(dis))很明显是错的,因为在floyd内部,dis+dis会跳负。应将65改成63或127/2。
  2.这个地方才真正调了我很久,这个地方是让人非常尴尬的。DP中常常有x=std::min(x,y);的语句出现,但我在修改时只把第一个x改对了……所以说需要有#define smin(x,y) x=std::min(x,y)啊(YYF、mcfx等大佬强力推荐)!
  最后挂一波代码吧。
 /**************************************************************
Problem: 4720
User: Doggu
Language: C++
Result: Accepted
Time:1960 ms
Memory:1292 kb
****************************************************************/ #define PN "classroom"
#include <cstdio>
#include <cstring>
#include <algorithm>
#define smin(x,y) x=std::min(x,y)
const int N = +;
const int M = +;
const int V = +;
int n, m, v, e, c[N], d[N], dis[V][V], cur=;
double k[N], f[][][M]; int main() {
scanf("%d%d%d%d",&n,&m,&v,&e);
for( int i = ; i <= n; i++ ) scanf("%d",&c[i]);
for( int i = ; i <= n; i++ ) scanf("%d",&d[i]);
for( int i = ; i <= n; i++ ) scanf("%lf",&k[i]);
memset(dis,,sizeof(dis));
for( int i = ,a,b,w; i <= e; i++ ) {
scanf("%d%d%d",&a,&b,&w);
if(dis[a][b]>w) dis[a][b]=dis[b][a]=w;
}
for( int i = ; i <= v; i++ ) dis[i][i]=;
for( int p = ; p <= v; p++ )
for( int i = ; i <= v; i++ )
for( int j = ; j <= v; j++ )
if (dis[i][j]>dis[i][p]+dis[p][j])
dis[i][j]=dis[i][p]+dis[p][j],
dis[j][i]=dis[i][j];
memset(f[cur],,sizeof(f[cur]));f[cur][][]=f[cur][][]=;
for( int i = ; i <= n; i++ ) {
cur^=;memset(f[cur],,sizeof(f[cur]));
for( int p = ; p <= m; p++ ) {
smin(f[cur][][p],f[cur^][][p]+k[i-]*dis[c[i]][d[i-]]+(-k[i-])*dis[c[i]][c[i-]]);
smin(f[cur][][p],f[cur^][][p]+dis[c[i]][c[i-]]);
smin(f[cur][][p+],f[cur^][][p]+k[i]*k[i-]*dis[d[i]][d[i-]]+(-k[i])*k[i-]*dis[c[i]][d[i-]]+k[i]*(-k[i-])*dis[d[i]][c[i-]]+(-k[i])*(-k[i-])*dis[c[i]][c[i-]]);
smin(f[cur][][p+],f[cur^][][p]+k[i]*dis[d[i]][c[i-]]+(-k[i])*dis[c[i]][c[i-]]);
}
}
double ans=;
for( int j = ; j <= ; j++ )
for( int p = ; p <= m; p++ )
ans=std::min(ans,f[cur][j][p]);
printf("%.2lf\n",ans);
return ;
}

NOIP第一道期望DP

BZOJ 4720 [Noip2016]换教室的更多相关文章

  1. bzoj 4720: [Noip2016]换教室【期望dp】

    状压dp,设f[i][j][0/1]为前i个时间段换了j间教室的期望体力消耗,转移很好想(但是写起来好长= =) #include<iostream> #include<cstdio ...

  2. 【BZOJ】4720: [Noip2016]换教室

    4720: [Noip2016]换教室 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1690  Solved: 979[Submit][Status ...

  3. BZOJ P4720[Noip2016]换教室____solution

    题目太长不表 <--无形传送,最为致命 学习一点数学期望的基础,预处理最短路,然后加上DP即可.(废话) 理解决策和结果的差别: 在这里每阶段的决策有两个:申请|不申请 结果有两个:换|不换 然 ...

  4. bzoj4720: [Noip2016]换教室(期望dp)

    4720: [Noip2016]换教室 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1294  Solved: 698[Submit][Status ...

  5. [NOIP2016]换教室 D1 T3 Floyed+期望DP

    [NOIP2016]换教室 D1 T3 Description 对于刚上大学的牛牛来说, 他面临的第一个问题是如何根据实际情况中情合适的课程. 在可以选择的课程中,有2n节课程安排在n个时间段上.在第 ...

  6. [NOIP2016]换教室 题解(奇怪的三种状态)

    2558. [NOIP2016]换教室 [题目描述] 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程. 在可以选择的课程中,有2n节课程安排在n个时间段上.在第i(1< ...

  7. NOIP2016换教室 BZOJ 4720

    BZOJ 4720 换教室 题目描述: 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节 课程安排在n个时间段上.在第i(1≤i≤n)个时间段上 ...

  8. 【bzoj4720】[NOIP2016]换教室

    题目描述 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节课程安排在n个时间段上.在第i(1≤i≤n)个时间段上,两节内容相同的课程同时在不同的 ...

  9. [NOIp2016] 换教室

    题目类型:期望\(DP\) 传送门:>Here< 题意:现有\(N\)个时间段,每个时间段上一节课.如果不申请换教室,那么时间段\(i\)必须去教室\(c[i]\)上课,如果申请换课成功, ...

随机推荐

  1. 解决debug到jdk源码时不能查看变量值的问题

    目录 如何跟踪jdk源码 1. 编译源码 2. 关联源码 3. 大功告成 如何跟踪jdk源码 看到这个标题大概大家都会在心里想谁还跟踪个源码呀,在eclipse中打个断点,以debug的方式运行,然后 ...

  2. Kubernetes学习之路(二十四)之Prometheus监控

    目录 1.Prometheus概述 2.Prometheus部署 2.1.创建名称空间prom 2.2.部署node_exporter 2.3.部署prometheus-server 2.4.部署ku ...

  3. opencv配置(转)

    1. 下载安装Opencv,去官网http://opencv.org/即可下载最新版本的Opencv,此处用的是Opencv 2.4.10 安装时傻瓜式的,最新版本的安装就是相当于解压到你指定的安装目 ...

  4. 继承类中static数据值

    class A{ static int num = 1; public static void Display(){ System.out.println( num ); } } class B ex ...

  5. JS跨浏览器的事件处理

    1. 事件流 事件:用户或浏览器自身执行的某种动作.如click(点击事件).mouse***(鼠标事件). 事件流:页面中接收事件的顺序,也可理解为事件在页面中传播的顺序. DOM事件流包括三个阶段 ...

  6. CS50.2

    1,ssd硬盘(solid state disk),固态硬盘,固盘. 2,把电磁信号转化为0或者1 ps:记得吧图给加上 反向即从磁盘中得到数据 3,软盘,floppy disk.早期使用的一种存储. ...

  7. 一步一步来熟悉Akka.Net(一)

    一步一步来熟悉Akka.Net(一) 标签(空格分隔): .netcore 分布式 一.不利flag   好久没写过文章了,翻开前几年写的博客,看到有两个目标"代码生成器"和&qu ...

  8. 【GAN与NLP】GAN的原理 —— 与VAE对比及JS散度出发

    0. introduction GAN模型最早由Ian Goodfellow et al于2014年提出,之后主要用于signal processing和natural document proces ...

  9. 一个http请求发送到后端的详细过程

    我们来看当我们在浏览器输入http://www.mycompany.com:8080/mydir/index.html,幕后所发生的一切. 首先http是一个应用层的协议,在这个层的协议,只是一种通讯 ...

  10. jmeter学习(1)基础支持+安装部署

    1. Apache jmeter 是100%的java桌面应用程序 支持的协议有:WEB-HTTP/HTTPS   , SOAP,   FTP,  JDBC,  LDAP, MAIL, MongoDB ...