NOIP2016换教室
题目
一道毒瘤概率期望DP。点这里
首先
对于每个时间段(就是每节课),我们有申请和不申请两种情况。如果不申请的话,一定在\(c[ i ]\),否则,可能在\(c[ i ]\),也可能在\(d[ i ]\)。所以我们先预处理一些东西。 先用弗洛伊德处理出最短路(这就不用详细介绍了吧)。我们先考虑相邻的两间教室,如果这两间教室都申请的话,就会有四种情况(①都通过,②都没通过,③第一个通过,④第二个通过)感觉又长又臭。所以我们先处理出两间教室都申请的情况。 用\(a[ i ]\)表示在都申请的前提下从第i个教室走到第i+1个教室的期望长度 \(dis[ i ] [ j ]\)表示从教室i到j的最短路。
情况①的概率为 $ k[ i ] k[ i+1 ] $ ,再乘上路径的长度 \(dis[ c[ i ] ][ c[ i+1 ] ]\)
同理 情况② \(1 - k[ i ]) (1 - k[ i+1 ]) dis[ d[ i ] ][ d[ i+1 ] ]\)
情况③ \(k[ i ] (1 - k[ i+1 ]) dis[ c[ i ] ][ d[ i+1 ] ]\)
情况④$(1 - k[ i ]) k[ i+1 ] * dis[ d[ i ] ][ c[ i+1 ] ] $
所以:
for(int i=1;i<=n;i++)
a[i]=(1-k[i])*(1-k[i+1])*dis[c[i]][c[i+1]]+
k[i]*(1-k[i+1])*dis[d[i]][c[i+1]]+k[i]*k[i+1]*dis[d[i]][d[i+1]]+
(1-k[i])*k[i+1]*dis[c[i]][d[i+1]]
太丑了
然后,继续预处理......
\(dp[ i ][ 0/1 ][ 0/1 ]\) 表示从第i节课的教室到第i+1节课的教室,第i节课的教室有(1)/没有(0)申请,第i+1个教室有(1)/没有(0)申请的期望长度
然后\(dp[ i ][ 1 ][ 1 ]\)的情况就是刚才的\(a[ i ]\)数组啦
其余情况同理可得,于是就:
for(int i=1;i<n;i++)
{
dp[i][0][0]=dis[c[i]][c[i+1]];
dp[i][1][1]=a[i];
dp[i][0][1]=k[i+1]*dis[c[i]][d[i+1]]+(1-k[i+1])*dis[c[i]][c[i+1]];
dp[i][1][0]=k[i]*dis[d[i]][c[i+1]]+(1-k[i])*dis[c[i]][c[i+1]];
}
然后,终于开始DP啦
用\(f[ i ][ j ][ 0/1 ]\) 示到达前i个教室,申请了j次,第i间教室有没有申请的期望最短长度。
!!!注意,j要从0开始,因为可以一次也不申请!!!
\(f[ i ][ j ][ 0 ]\)可以由\(f[ i - 1 ][ j ][ 0 ]\)和\(f[ i - 1 ][ j ][ 1 ]\)转移过来
\(f[ i ][ j ][ 1 ]\)可以由\(f[ i - 1][ j - 1 ][ 0 ]\)和\(f[ i - 1][ j - 1 ][ 1 ]\)转移过来
先初始化一下 \(f[ i ][ 0 ][ 0 ]\) 的情况
于是就得到了:
f[1][1][1]=0; f[1][0][0]=0;
for(int i=2;i<=n;i++)
f[i][0][0]=f[i-1][0][0]+dis[c[i-1]][c[i]];
for(int i=2;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
f[i][j][0]=min(f[i-1][j][1]+dp[i-1][1][0],f[i-1][j][0]+dp[i-1][0][0]);
f[i][j][1]=min(f[i-1][j-1][1]+dp[i-1][1][1],f[i-1][j-1][0]+dp[i-1][0][1]);
}
}
最后统计答案就行啦,要注意可以申请 0到m次。
for(int i=1;i<=m;i++)
ans=min(ans,min(f[n][i][0],f[n][i][1]));
于是,这道题就做完了!
最后贴上完整版代码(大佬勿喷):
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 2005;
const int INF = 27654134;
int n,m,v,e,c[N],d[N],dis[N][N];
double a[N],k[N],f[N][N][2],dp[N][2][2],ans=99999999.9;
//a[i]表示在都申请的前提下从第i个教室走到第i+1个教室的期望长度
// dp[i][0/1][0/1]从第i个教室到第i+1个教室,第i个教室有没有申请,第i+1个教室有没有申请的期望长度
//f[i][j][0/1]表示到达前i个教室,申请了j次,第i间教室有没有申请的期望最短长度
int main()
{
memset(f,0x7f,sizeof(f));
// memset(dis,0x7f,sizeof(dis));
scanf("%d%d%d%d",&n,&m,&v,&e);
for(int i=1;i<=v;i++)
for(int j=1;j<=v;j++)
dis[i][j]=INF;
for(int i=1;i<=n;i++) scanf("%d",&c[i]);
for(int i=1;i<=n;i++) scanf("%d",&d[i]);
for(int i=1;i<=n;i++) cin>>k[i];
int x,y,z;
for(int i=1;i<=e;i++)
{
scanf("%d%d%d",&x,&y,&z);
if(dis[x][y]) dis[x][y]=min(dis[x][y],z),dis[y][x]=min(dis[y][x],z);
else dis[x][y]=z,dis[y][x]=z;
}
for(int i=1;i<=v;i++) dis[i][i]=0;
for(int kk=1;kk<=v;kk++)
{
for(int i=1;i<=v;i++)
{
if(kk==i) continue;
for(int j=1;j<=v;j++)
{
if(kk==j||i==j) continue;
if(dis[i][kk]+dis[kk][j]<dis[i][j])
dis[i][j]=dis[i][kk]+dis[kk][j];
}
}
}
for(int i=1;i<=n;i++)
a[i]=(1-k[i])*(1-k[i+1])*dis[c[i]][c[i+1]]+
k[i]*(1-k[i+1])*dis[d[i]][c[i+1]]+k[i]*k[i+1]*dis[d[i]][d[i+1]]+
(1-k[i])*k[i+1]*dis[c[i]][d[i+1]];
for(int i=1;i<n;i++)
{
dp[i][0][0]=dis[c[i]][c[i+1]];
dp[i][1][1]=a[i];
dp[i][0][1]=k[i+1]*dis[c[i]][d[i+1]]+(1-k[i+1])*dis[c[i]][c[i+1]];
dp[i][1][0]=k[i]*dis[d[i]][c[i+1]]+(1-k[i])*dis[c[i]][c[i+1]];
}
f[1][1][1]=0; f[1][0][0]=0;
for(int i=2;i<=n;i++)
f[i][0][0]=f[i-1][0][0]+dis[c[i-1]][c[i]];
for(int i=2;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
f[i][j][0]=min(f[i-1][j][1]+dp[i-1][1][0],f[i-1][j][0]+dp[i-1][0][0]);
f[i][j][1]=min(f[i-1][j-1][1]+dp[i-1][1][1],f[i-1][j-1][0]+dp[i-1][0][1]);
}
}
for(int i=0;i<=m;i++)
ans=min(ans,min(f[n][i][0],f[n][i][1]));
printf("%.2lf\n",ans);
return 0;
}
NOIP2016换教室的更多相关文章
- [NOIP2016]换教室 D1 T3 Floyed+期望DP
[NOIP2016]换教室 D1 T3 Description 对于刚上大学的牛牛来说, 他面临的第一个问题是如何根据实际情况中情合适的课程. 在可以选择的课程中,有2n节课程安排在n个时间段上.在第 ...
- BZOJ 4720 [Noip2016]换教室
4720: [Noip2016]换教室 Description 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节课程安排在n个时间段上.在第i( ...
- 【BZOJ】4720: [Noip2016]换教室
4720: [Noip2016]换教室 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1690 Solved: 979[Submit][Status ...
- bzoj4720: [Noip2016]换教室(期望dp)
4720: [Noip2016]换教室 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1294 Solved: 698[Submit][Status ...
- [NOIP2016]换教室 题解(奇怪的三种状态)
2558. [NOIP2016]换教室 [题目描述] 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程. 在可以选择的课程中,有2n节课程安排在n个时间段上.在第i(1< ...
- 【bzoj4720】[NOIP2016]换教室
题目描述 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节课程安排在n个时间段上.在第i(1≤i≤n)个时间段上,两节内容相同的课程同时在不同的 ...
- [NOIp2016] 换教室
题目类型:期望\(DP\) 传送门:>Here< 题意:现有\(N\)个时间段,每个时间段上一节课.如果不申请换教室,那么时间段\(i\)必须去教室\(c[i]\)上课,如果申请换课成功, ...
- [NOIP2016]换教室(概率期望$DP$)
其实吧我老早就把这题切了--因为说实话,这道题确实不难啊--李云龙:比他娘的状压DP简单多了 今天我翻以前在Luogu上写的题解时,突然发现放错代码了,然后被一堆人\(hack\)--蓝瘦啊\(ORZ ...
- 【bzoj4720】[NOIP2016]换教室 期望dp
题目描述 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节课程安排在n个时间段上.在第i(1≤i≤n)个时间段上,两节内容相同的课程同时在不同的 ...
- NOIP2016换教室 BZOJ 4720
BZOJ 4720 换教室 题目描述: 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程.在可以选择的课程中,有2n节 课程安排在n个时间段上.在第i(1≤i≤n)个时间段上 ...
随机推荐
- 【转】shell脚本中如何传入参数
(1)直接用$1,$2取传入的参数vim /root/test.sh#!/bin/bashif [ $1 == "start" ] then echo "do ...
- JAVA ==号和equals()的区别
==号和equals()方法都是比较是否相等的方法,那它们有什么区别和联系呢? 首先,==号在比较基本数据类型时比较的是值,而用==号比较两个对象时比较的是两个对象的地址值: int x = 10; ...
- ABP中的拦截器之ValidationInterceptor(下)
在上篇我分析了整个ABP中ValitationInterceptor的整个过程,就其中涉及到的Validator过程没有详细的论述,这篇文章就这个过程进行详细的论述,另外任何一个重要的特性如何应用是最 ...
- springboot 静态注入 单例
package com.b2q.web_push.util; import io.goeasy.GoEasy; import org.springframework.beans.factory.ann ...
- Batch Normalization原理
Batch Normalization导读 博客转载自:https://blog.csdn.net/malefactor/article/details/51476961 作者: 张俊林 为什么深度神 ...
- python之异常处理和re模块补充
一.re模块的补充 1.从一个字符串中获取要匹配的内容 findall:返回一个列表 2.search ***** 验证用户输入内容 '^正则规则$':返回一个对象,用group()取值 3.matc ...
- shell实战之日志备份
util.sh #!/bin/bash - # Read config # kay # Version: 1.0 # // # configuration file path config=${log ...
- 微信小程序 mpvue vant
Mpvue中使用Vant Weapp组件库 https://segmentfault.com/a/1190000016228410?utm_source=tag-newest 小程序采坑记 mpvue ...
- ue4 材质表达式分类
绿色节点 颜色 Color Desaturation 数学 Math GO 字体 Font FontSample,FontSampleParameter 实用程序 Utility 常用: Desatu ...
- How to avoid the 0-byte file on Webclient download error
_client.DownloadDataAsync(new Uri(url)); _client.DownloadDataCompleted += (sender, e) => { try { ...