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. 【php增删改查实例】 第三节 - mysql 建表

    这一节,来给数据库中添加一些测试数据. 登陆mysql: 找到%xampp%\mysql\bin 目录, 在此处打开命令窗口,用root用户登陆mysql 用户表建表sql: CREATE TABLE ...

  2. 微信小程序之路由

    1. 路由方式 路由方式 触发时机 路由前页面 路由后页面 初始化 小程序打开的第一个页面 onLoad, onShow 打开新页面 调用 API wx.navigateTo 或使用组件 onHide ...

  3. Flask学习-Flask app启动过程

    因为0.1版本整体代码大概只有350行,比较简单.所以本篇文章会以Flask 0.1版本源码为基础进行剖析Flask应用的启动过程. Flask参考资料flask,官网有一个最简单app: from ...

  4. (2)学习笔记 ) ASP.NET CORE微服务 Micro-Service ---- .NetCore启动配置 和 .NetCoreWebApi

    什么是.Net Core?.Net Core是微软开发的另外一个可以跨Linux.Windows.mac等平台的.Net.Net Core相关知识看文章地步dotnet dllname.dll 运行P ...

  5. 机器视觉及图像处理系列之二(C++,VS2015)——图像级的人脸识别(1)

    接上一篇,一切顺利的话,你从github上clone下来的整个工程应该已经成功编译并生成dll和exe文件了:同时,ImageMagic程序亦能够打开并编辑图像了,如此,证明接下来的操练你不会有任何障 ...

  6. JavaScript中执行环境和栈

    在这篇文章中,我会深入理解JavaScript最根本的组成之一 : "执行环境(执行上下文)".文章结束后,你应该对解释器试图做什么,为什么一些函数/变量在未声明时就可以调用并且他 ...

  7. js执行问题

    金三银四搞事季,前端这个近年的热门领域,搞事气氛特别强烈,我朋友小伟最近就在疯狂面试,遇到了许多有趣的面试官,有趣的面试题,我来帮这个搞事 boy 转述一下. 以下是我一个朋友的故事,真的不是我. f ...

  8. 20135202闫佳歆--week1 计算机是如何工作的

    计算机是如何工作的 这一周我学习了计算机工作的相关知识. 最基础的,就是冯诺依曼体系结构结构,它最核心的思想是存储程序计算机,要点是:数字计算机的数制采用二进制:计算机应该按照程序顺序执行. 除了思想 ...

  9. 小学生二元四则运算(F)

      整体功能简介: 1.题目不重复: 2.可以定制数量: 3.可以自己选择输入范围: 4.可以选择是否添加乘除法: 5.可以选择除法结果是取整或商加余数形式表示或小数方式(默认小数点后两位)表示: 6 ...

  10. 团队作业一 庆祝"十五万的总冠军"成立

    很荣幸能够撰写我们团队的第一篇博客. 我们这些同学能组成一个新的团队真的很高兴,团队中的每一个人都有自己的优点的长处.希望在工作中我们能竭尽 所能,充分发挥我们的本事,让我们大家各自发挥自己的才能.. ...