传送门!

Simple Description:

牛牛有$n$节课要上,每个时段都有$2$个教室在上课(具体来说是第$c[i]$与第$d[i]$个教室),牛牛在其中的一个教室上课,牛牛起初每节课分别被安排在$c[i](i=1..n)$节教室上课,他可以提出申请到教室$d[i]$上课,申请通过的概率为$pi$。牛牛最多能申请m次。题目给出一些联通教室的路径和通过路径要花费的体力值。求牛牛耗费体力值的总和的最小期望值。

(希望我说清楚了)

Solution:

数学期望+DP

求最值啊就会想到DP的啊,而且这题看起来很DP。

状态

也很容易想,$f[i][j]$表示现在到第i节课,已经申请了m次的最小期望值。然后就会发现不好转移...一般套路就是再加一维表示上一节课是否申请了。于是,最终是:

$f[i][j][0/1]$表示现在到第i节课,已经申请了m次,上一次课没申请/申请了的最小期望值

转移方程

设$i,j$之间的最短路为$dis[i][j]$(这里初始化一下)

这节没申请($f[i][j][0]$)

上节也没申请,这种情况乘的概率就是百分百啊:  

f[i][j][0]=f[i-1][j][0]+dis[c[i-1]][c[i]]

上节申请了,分申请通过和没通过两种情况,分别乘上自己的概率:

f[i][j][0]=f[i-1][j] [1]+dis[d[i-1] ][c[i]]*p[i-1]+ dis[c[i-1]][c[i]]*(1-p[i-1])

综合起来取小:

f[i][j][0]=min(f[i-1][j][0]+dis[c[i-1]][c[i]],f[i-1][j][1]+dis[d[i-1]][c[i]]*p[i-1]+dis[c[i-1]][c[i]]*(1-p[i-1]))

这节申请了($f[i][j][1]$)

上节没申请:

f[i][j][1]=f[i-1][j-1][0]+dis[c[i-1]][d[i]]*p[i]+dis[c[i-1]][c[i]]*(1-p[i])

上节申请了:

f[i][j][1]=f[i-1][j-1][1]+dis[d[i-1]][d[i]]*p[i-1]*p[i]+dis[c[i-1]][d[i]]*(1-p[i-1])*p[i]+dis[d[i-1]][c[i]]*p[i-1]*(1-p[i])+dis[c[i-1]][c[i]]*(1-p[i-1]*(1-p[i]))

综合起来取小:

f[i][j][1]=min(f[i-1][j-1][0]+dis[c[i-1]][d[i]]*p[i]+dis[c[i-1]][c[i]]*(1-p[i]),f[i-1][j-1][1]+dis[d[i-1]][d[i]]*p[i-1]*p[i]+dis[c[i-1]][d[i]]*(1-p[i-1])*p[i]+dis[d[i-1]][c[i]]*p[i-1]*(1-p[i])+dis[c[i-1]][c[i]]*(1-p[i-1]*(1-p[i])))

Code!

#include<iostream>
#include<cstdio>
#define Rg register
#define il inline
#define db double
#define inf 200000001
#define int long long
#define go(i,a,b) for(Rg int i=a;i<=b;i++)
#define yes(i,a,b) for(Rg int i=a;i>=b;i--)
using namespace std;
il int read()
{int x=0,y=1;char c=getchar();while(c<'0'||c>'9'){if(c=='-')y=-1;c=getchar();}while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}return x*y;
}
const int N=2001,M=301;
int m,n,v,e,c[N],d[N],dis[M][M];
db f[N][N][2],p[N],ans=inf;
void Read()
{n=read(),m=read(),v=read(),e=read();go(i,1,v) go(j,1,v)if(i!=j)dis[i][j]=inf;go(i,1,n)c[i]=read(); go(i,1,n)d[i]=read(); go(i,1,n)scanf("%lf",&p[i]);go(i,1,e){int x=read(),y=read();dis[x][y]=dis[y][x]=min(dis[x][y],read());}
}
void Floyed()
{go(k,1,v) go(j,1,v) go(i,1,j)if(dis[i][k]+dis[k][j]<dis[i][j])dis[i][j]=dis[j][i]=dis[i][k]+dis[k][j];
}
main()
{Read();Floyed();go(i,0,n) go(j,0,n)f[i][j][0]=f[i][j][1]=f[j][i][0]=f[j][i][1]=inf;f[1][0][0]=f[1][1][1]=0;go(i,2,n){int maxn=min(i,m);go(j,0,maxn){f[i][j][0]=min((db)(f[i-1][j][0]+dis[c[i-1]][c[i]]),(db)(f[i-1][j][1]+dis[d[i-1]][c[i]]*p[i-1]+dis[c[i-1]][c[i]]*(1-p[i-1])));if(j>0)f[i][j][1]=min((db)(f[i-1][j-1][0]+dis[c[i-1]][d[i]]*p[i]+dis[c[i-1]][c[i]]*(1-p[i])),(db)(f[i-1][j-1][1]+dis[c[i-1]][c[i]]*(1-p[i-1])*(1-p[i])+dis[c[i-1]][d[i]]*(1-p[i-1])*p[i]+dis[d[i-1]][c[i]]*(1-p[i])*p[i-1]+dis[d[i-1]][d[i]]*p[i-1]*p[i]));}}go(i,0,m)ans=min(min(f[n][i][0],f[n][i][1]),ans);printf("%.2lf",ans);return 0;
}

View Code

转载于:https://www.cnblogs.com/forward777/p/10827409.html

Noip2016换教室相关推荐

  1. BZOJ 4720 [Noip2016]换教室

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

  2. [NOIP2016]换教室(概率期望$DP$)

    其实吧我老早就把这题切了--因为说实话,这道题确实不难啊--李云龙:比他娘的状压DP简单多了 今天我翻以前在Luogu上写的题解时,突然发现放错代码了,然后被一堆人\(hack\)--蓝瘦啊\(ORZ ...

  3. noip2016 换教室

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

  4. BZOJ 4720: [Noip2016]换教室

    错是被学长挑出来的...QWQ@Monster_Yi f[i][j][0/1] 为 第i天,换j次,换1/没换0 然后大力DP #include<cstdio> #include<i ...

  5. luogu1850 [NOIp2016]换教室 (floyd+dp)

    首先floyd求出每两点间的距离(注意自己到自己的距离要设成0) 然后就是dp了 一开始照着Lifeguards的样子,钦定了一下i这个点一定要选,然后发现复杂度不对,还想了好长时间优化 然后一翻题解 ...

  6. 洛谷1850(NOIp2016) 换教室——期望dp

    题目:https://www.luogu.org/problemnew/show/P1850 状态里记录的是"上一回有没有申请",而不是"上一回申请成功否",不 ...

  7. P1850 [NOIP2016 提高组] 换教室

    P1850 [NOIP2016 提高组] 换教室 题意: 有2n个课安排在n个时间段上,每个时间段上都有两个一样的课同时在不同地方上,起初牛牛被所有课都被安排在Ci上课,另一节课在Di上课.牛牛现在想 ...

  8. 【NOIp2016 day1t3】换教室

    NOIP第一次考期望,着实吓一跳... 读入之后, 由于 n<=200 n <= 200 n,给我们建立了天然的 floyd f l o y d floyd的机会, 建完图之后,能够顺利的 ...

  9. 【DP】【期望】$P1850$换教室

    [DP][期望]\(P1850\)换教室 链接 题目描述 有 \(2n\) 节课程安排在$ n$ 个时间段上.在第 \(i\)(\(1 \leq i \leq n\))个时间段上,两节内容相同的课程同 ...

最新文章

  1. 《Linux嵌入式实时应用开发实战(原书第3版)》——1.6 资源
  2. Rumor CodeForces - 893C(并查集)
  3. stm32 IAP + APP ==双剑合一
  4. 2017.9.28 CF #R2 B 思考记录
  5. ERP实施--常见问题
  6. 关于CocosStudioV3.10的场景导入Cocos2d-xV3.17.2场景未发生变化的解决方案
  7. 自增字段不连续_MySQL中自增主键不连续之解决方案。(20131109)
  8. c++游戏编程100事列_游戏编程入门(c++小游戏编程100例)
  9. android超级终端使用,android系统超级终端怎么用?
  10. 在虚拟机中安装Linux系统
  11. 计算机串口是232吗,RS232串口简介
  12. java.lang.IllegalStateException: focus search returned a view that wasn‘t able to take focus!
  13. 揭密备份恢复的原理!
  14. 在C中将二进制转换为十进制
  15. 阿里菜鸟愚人节最新力作,高大上宠物递送萌翻铲屎官
  16. 从校训、企业文化到团队文化
  17. 任意函数展开为各阶Taylor多项式的matlab程序
  18. 线性规划以及二次规划
  19. 迅雷之CDN加速下载/离线下载
  20. SMART 200PLC指针寻址

热门文章

  1. 用jquery的ajax功能获取网站alexa的方法
  2. 下列不是python对文件的读操作方法是-小猿圈python之python期末考试测试题(一)...
  3. 下列关于python的描述正确的是-以下关于Python循环结构的描述中,错误的是( )...
  4. python爬虫入门代码-Python爬虫入门
  5. python自动化测试看什么书-《Python自动化测试实战》终于出版!!!
  6. python如何爬虫网页数据-python爬虫——爬取网页数据和解析数据
  7. python while循环语句-Python中的while循环语句怎么用?
  8. python读文件路径-python获取程序执行文件路径的方法(推荐)
  9. python小游戏编程实例-10分钟教你用Python写一个贪吃蛇小游戏,适合练手项目
  10. python有趣代码-你都知道哪些有趣的Python代码?