题目链接:点击这里献出你宝贵的时间(是用来做题不是捐赠Emmmm)..

Emmmm我太弱了= =

做完这题我觉得我应该去打星际..这题怎么就有重边了呢..

这题就是一道期望= =当时考场上好像完全不会期望,dp一直也都不会推式子= =
然后就用floyd水掉了没有期望的点,分数还是太少了,最后就省二滚粗了..

现在觉得应该在NOIP2017之前学习一下期望,然后就想起了这道题,lhr大佬之前也做了这道题但好像没做完??然后就和我一起做= =他代码跑飞快我也不知道为啥= =

好吧言归正传,这题就是floyd+dp,当然了这题当然和期望有关,你要接下去看我就默认你对期望有所了解了哈= =不然去找度娘问问去= = 数学期望

由于我们要求出多源最短路,所以floyd.. 而且v<=300不是嘛..
n和v真的挺坑的,刚开始还RE了不少遍啊= =

Emmmm,我们发现由于只有上一节课的选择会影响到即将加一条边的期望值,所以我们可以dp(其实我不知道为什么嘛~)
我们设f[i][j][0/1]表示处理第i节课,已经申请了j次,这次选择/不选择申请的期望…
这样的话状态转移方程就呼之欲出了(不存在的,年轻人,要多想):
然而要做很复杂的分类讨论:

*我们先做一下符号的声明:
*- pre[i]表示第i节本来的安排
*- now[i]表示第i次申请能换到的教室
*- p[i]表示第i次申请通过的概率
*- dis[i][j]表示i,j间的距离

  • 这次不申请,我们更新f[i][j][0]的值,那么新加的边的终点就是pre[i]..
    1.如果上次不申请,那么新加的边的起点就是pre[i-1],不需要考虑期望..
    2.如果上次申请了,那么成功了起点就是now[i-1],不成功起点就是pre[i-1]
    所以:
f[i][j][0]=min(
f[i-1][j][0]+dis[pre[i-1]][pre[i]],
f[i-1][j][1]+dis[now[i-1]][pre[i]]*p[i-1] //若上一次申请成功
+dis[pre[i-1]][pre[i]]*(1-p[i-1]) //若上一次申请不成功
//根据期望的性质易得,申请成功时的期望*成功率+不成功的期望*不成功率
);
  • 而这次申请的话,我们更新f[i][j][1]的值,
    1.如果上次不申请,起点就是pre[i-1]…终点还是根据期望的性质推:
f[i-1][j-1][0]+ //这次申请所以之前申请j-1次
dis[pre[i-1]][now[i]]*p[i] //这次申请成功
+dis[pre[i-1]][pre[i]]*(1-p[i]) //这次申请不成功

2.如果上次申请,起点和终点我们就要根据期望的性质枚举一下(太麻烦了= =)

f[i-1][j-1][1]+
dis[now[i-1]][now[i]]*p[i-1]*p[i] //两次申请都成功
+dis[now[i-1]][pre[i]]*p[i-1]*(1-p[i]) //上次申请成功,这次不成功
+dis[pre[i-1]][now[i]]*(1-p[i-1])*p[i] //上次申请不成功,这次成功
+dis[pre[i-1]][pre[i]]*(1-p[i-1])*(1-p[i]) //两次申请都不成功

然后我们就对以上两个式子取一下min:

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

Emmmm这个状态转移方程看上去很简单不是么(大雾)
边界条件的话,

f[1][0][0]=f[1][1][1]=0;

然后一路推过去就行了~~
代码:

#include <cstdio>
const int V = 303;
const int N = 2002;
template<class T>
inline T min(const T &a, const T &b) {
if (a<b) return a; return b;
}
int d[V][V], n, m, v, e;
double f[N][N][2];
struct apply {
double p; int pre, now;
}a[N]; void floyd() {
for (int k = 1; k <= v; k++)
for (int i = 1; i <= v; i++)
for (int j = 1; j <= v; j++)
if (i != j&&i != k&&k != j)
d[i][j] = min(d[i][j], d[i][k] + d[k][j]);
}
inline int getint() {
int a = 0; char c = getchar(); bool f = 0;
for (; (c<'0' || c>'9') && c != '-'; c = getchar());
if (c == '-') f = 1, c = getchar();
for (; c >= '0'&&c <= '9'; c = getchar()) a = (a << 1) + (a << 3) + c - '0';
if (f) return -a; return a;
}
int main() {
n = getint(), m = getint(), v = getint(), e = getint();
for (int i = 1; i <= n; i++) a[i].pre = getint();
for (int i = 1; i <= n; i++) a[i].now = getint();
for (int i = 1; i <= n; i++) scanf("%lf", &a[i].p);
for (int i = 0; i <= v; i++) for (int j = 0; j <= v; j++) d[i][j] = 0x3fffffff;
for (int i = 0; i <= v; i++) d[i][i] = 0;
for (int i = 1; i <= e; i++) { int x = getint(), y = getint(); d[x][y] = d[y][x] = min(d[x][y], getint()); }floyd(); for (int i = 1; i <= n; i++) for (int j = 0; j <= m; j++) f[i][j][0] = f[i][j][1] = 0x3fffffff;
f[1][0][0] = f[1][1][1] = 0;
for (int i = 2; i <= n; i++)
for (int j = 0; j <= m; j++) {
f[i][j][0] = min(
f[i - 1][j][0] + d[a[i - 1].pre][a[i].pre],
f[i - 1][j][1] + d[a[i - 1].now][a[i].pre] * a[i - 1].p + d[a[i - 1].pre][a[i].pre] * (1 - a[i - 1].p)
);
if(j) f[i][j][1] = min(
f[i - 1][j - 1][0] + d[a[i - 1].pre][a[i].now] * a[i].p + d[a[i - 1].pre][a[i].pre] * (1 - a[i].p),
f[i - 1][j - 1][1] + d[a[i - 1].now][a[i].now] * a[i - 1].p * a[i].p
+ d[a[i - 1].now][a[i].pre] * a[i - 1].p * (1 - a[i].p)
+ d[a[i - 1].pre][a[i].now] * (1 - a[i - 1].p) * a[i].p
+ d[a[i - 1].pre][a[i].pre] * (1 - a[i - 1].p) * (1 - a[i].p)
);
}
double ans = 1e10;
for (int i = 0; i <= m; i++) ans = min(ans, min(f[n][i][0], f[n][i][1])); printf("%.2lf", ans);
}

跑得还是挺慢的……Emmmm……

就这样吧= =

【学术篇】NOIP2016 D1T3 luogu1850换教室的更多相关文章

  1. [NOIp2016提高组]换教室

    题目大意: 有n节课,第i节课在c[i]上课,同时d[i]也有一节课d[i]. 你有权利向教务处发出m次申请把自己的教室改到d[i],相应的批准概率是k[i]. 教室是图上的一些点,其中每条边都有边权 ...

  2. NOIP2016 DAY1 T3 换教室

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

  3. Luogu P1850 [NOIp2016提高组]换教室 | 期望dp

    题目链接 思路: <1>概率与期望期望=情况①的值*情况①的概率+情况②的值*情况②的概率+--+情况n的值*情况n的概率举个例子,抛一个骰子,每一面朝上的概率都是1/6,则这一个骰子落地 ...

  4. 洛谷 1850 NOIP2016提高组 换教室

    [题解] 先用floyed处理出两点间的最短路. 设f[i][j][k]表示走到第i个教室,总共换了j次,当前换或者不换,期望的最小移动距离. 分情况讨论来转移即可. #include<cstd ...

  5. 题解 P1850 [NOIP2016 提高组] 换教室

    做完这道题才略微感觉自己懂了一点关于概率与期望的知识QAQ... 一:关于概率与期望的定义 转载节选于blog 1.什么是数学期望? 数学期望亦称期望.期望值等.在概率论和统计学中,一个离散型随机变量 ...

  6. [NOIP2016][luogu]换教室[DP]

    [NOIP2016] Day1 T3 换教室 ——!x^n+y^n=z^n 题目描述 对于刚上大学的牛牛来说,他面临的第一个问题是如何根据实际情况申请合适的课程. 在可以选择的课程中,有 2n 节课程 ...

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

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

  8. BZOJ 4720 [Noip2016]换教室

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

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

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

随机推荐

  1. NX二次开发-UFUN关闭STL文件函数UF_STD_close_stl_file

    NX9+VS2012 #include <uf.h> #include <uf_obj.h> #include <uf_modl.h> #include <u ...

  2. NX二次开发-获得制图中对象的坐标点UF_DRF_ask_origin

    #include <uf.h> #include <uf_ui.h> #include <uf_drf.h> #include <uf_obj.h> # ...

  3. hdu多校第五场1002 (hdu6625) three arrays 字典树/dfs

    题意: 给你两个序列a,b,序列c的某位是由序列a,b的此位异或得来,让你重排序列ab,找出字典序最小的序列c. 题解: 如果能找到a,b序列中完全一样的值当然最好,要是找不到,那也尽量让低位不一样. ...

  4. Java设计模式简单总结

    1.单例模式:在应用程序整个生命周期中,单例类的实例只有一个,并且会自动实例化.单例类的构造方法必须为私有,并且提供一个全局访问点 public class Test { private Test() ...

  5. SSM 整合 Shiro

    1. 导包 <!-- spring --> <dependency> <groupId>org.springframework</groupId> &l ...

  6. Python3入门机器学习经典算法与应用✍✍✍

    Python3入门机器学习经典算法与应用 整个课程都看完了,这个课程的分享可以往下看,下面有链接,之前做java开发也做了一些年头,也分享下自己看这个视频的感受,单论单个知识点课程本身没问题,大家看的 ...

  7. Python匹马天下行之python基本语法

    输出 输出皆是在控制台打印出要打印的内容,在python中用print("内容") 在括号中加入要输出的内容就可以打印出来了. 例: print("Hello Pytho ...

  8. RDLC报表问题:尚未指定报表“Report1”的报表定义

    原文:尚未指定报表“Report1”的报表定义 在做RDLC项目中遇到这样的错误 本地报表处理期间出错. 尚未指定报表“Report1”的报表定义 未将对象引用设置到对象的实例. 解决方案: 打开re ...

  9. Django ORM 之基于对象、双下划线查询

    返回ORM目录 Django ORM 内容目录: 一. 基于对象的表查询 二. 基于双下划线的查询 三. 聚合查询 aggregate 四. 分组查询 annotate 一. 基于对象的表查询 1.正 ...

  10. 笔记:简单的面向对象-web服务器

    import socket import re import multiprocessing import time import mini_frame class WSGIServer(object ...