NOIP 2016 换教室 (luogu 1850 & uoj 262) - 概率与期望 - 动态规划

题目描述

对于刚上大学的牛牛来说，他面临的第一个问题是如何根据实际情况申请合适的课程。

在可以选择的课程中，有 $\leq i \leq n1≤i≤n）个时间段上，两节内容相同的课程同时在不同的地点进行，其中，牛牛预先被安排在教室 c_ici 上课，而另一节课程在教室 d_idi 进行。$

在不提交任何申请的情况下，学生们需要按时间段的顺序依次完成所有的 $d_idi 上课，否则仍然在教室 c_ici 上课。$

由于更换教室的需求太多，申请不一定能获得通过。通过计算，牛牛发现申请更换第 $k_iki，并且对于不同课程的申请，被通过的概率是互相独立的。$

学校规定，所有的申请只能在学期开始前一次性提交，并且每个人只能选择至多

因为不同的课程可能会被安排在不同的教室进行，所以牛牛需要利用课间时间从一间教室赶到另一间教室。

牛牛所在的大学有 $\leq i \leq n-11≤i≤n−1）节课结束后，牛牛就会从这节课的教室出发，选择一条耗费体力最少的路径前往下一节课的教室。$

现在牛牛想知道，申请哪几门课程可以使他因在教室间移动耗费的体力值的总和的期望值最小，请你帮他求出这个最小值。

输入输出格式

输入格式：

第一行四个整数

第二行 $\leq i \leq n1≤i≤n）个正整数表示 c_ici，即第 ii 个时间段牛牛被安排上课的教室；保证 1 \le c_i \le v1≤ci≤v。$

第三行 $\leq i \leq n1≤i≤n）个正整数表示 d_idi，即第 ii 个时间段另一间上同样课程的教室；保证 1 \le d_i \le v1≤di≤v。$

第四行 $\leq i \leq n1≤i≤n）个实数表示 k_iki，即牛牛申请在第 ii 个时间段更换教室获得通过的概率。保证 0 \le k_i \le 10≤ki≤1。$

接下来 $a_j, b_j, w_jaj,bj,wj，表示有一条双向道路连接教室 a_j, b_jaj,bj，通过这条道路需要耗费的体力值是 w_jwj；保证 1 \le a_j, b_j \le v1≤aj,bj≤v， 1 \le w_j \le 1001≤wj≤100。$

保证 $\leq n \leq 20001≤n≤2000，0 \leq m \leq 20000≤m≤2000，1 \leq v \leq 3001≤v≤300，0 \leq e \leq 900000≤e≤90000。$

保证通过学校里的道路，从任何一间教室出发，都能到达其他所有的教室。

保证输入的实数最多包含

输出格式：

输出一行，包含一个实数，四舍五入精确到小数点后恰好

测试数据保证四舍五入后的答案和准确答案的差的绝对值不大于 $\times 10^{-3}4×10−3。（如果你不知道什么是浮点误差，这段话可以理解为：对于大多数的算法，你可以正常地使用浮点数类型而不用对它进行特殊的处理）$

输入输出样例

输入样例#1：复制

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

输出样例#1：复制

2.80

说明

【样例1说明】

所有可行的申请方案和期望收益如下表:

【提示】

道路中可能会有多条双向道路连接相同的两间教室。也有可能有道路两端连接

的是同一间教室。

2.请注意区分n,m,v,e的意义, n不是教室的数量, m不是道路的数量。

特殊性质1:图上任意两点 $a_iai, b_ibi, a_iai≠ b_ibi间,存在一条耗费体力最少的路径只包含一条道路。$

特殊性质2:对于所有的 $1≤ i≤ n$， $k_i= 1ki=1 。$

　　题目大意 (题目太"简洁"，不需要大意)

　　显然动态规划。

　　很轻松就能想出 f[i][j][0 / 1] 表示上到第i节课，已经申请了j次，第i节课有没有申请(这个会影响结果)

　　然后考虑转移。

　　对于状态 f[i][j][0] ,它可以通过上次申请了或者没有申请转移过来，即 f[i][j][0] = min(f[i - 1][j][0] + dcc, f[i - 1][j][1] + dcc * (1 - ip) + ddc * ip)

　　对于状态 f[i][j][1] ，同理，然后加上对应距离乘对应概率就好了。(真的觉得没什么好说的。。)

　　记得去年考NOIP的时候做这道题，都不知道什么是期望，一脸懵逼。

　　其实这道题并不难。然而我把Floyd的for k写到了最里面，于是我就笑笑（1个小时就这么飞走了，我是不是学了假的Floyd?）。。

Code

 1 /**
 2  * luogu
 3  * Problem#1850
 4  * Accepted
 5  * Time: 2724ms
 6  * Memory: 63597k
 7  */
 8 #include <bits/stdc++.h>
 9 using namespace std;
10 #define smin(_a, _b) _a = min(_a, _b)
11 #define smax(_a, _b) _a = max(_a, _b) 
12 
13 const int V = 302;
14 const int N = 2002, M = 2002;
15 
16 int n, m, v, e;
17 int dis[V][V];
18 double f[N][M][2];
19 int *c, *d;
20 double *p;
21 
22 inline void init() {
23     scanf("%d%d%d%d", &n, &m, &v, &e);
24     c = new int[(n + 1)];
25     d = new int[(n + 1)];
26     p = new double[(n + 1)];
27     for(int i = 1; i <= n; i++)
28         scanf("%d", c + i);
29     for(int i = 1; i <= n; i++)
30         scanf("%d", d + i);
31     for(int i = 1; i <= n; i++)
32         scanf("%lf", p + i);
33     for(int i = 1; i <= v; i++)
34         for(int j = 1; j <= v; j++)
35             dis[i][j] = 1 << 29;
36     for(int i = 1, u, v, w; i <= e; i++) {
37         scanf("%d%d%d", &u, &v, &w);
38         smin(dis[u][v], w);
39         dis[v][u] = dis[u][v];
40     }
41 }
42 
43 inline void floyed() {
44     for(int i = 1; i <= v; i++)
45         dis[i][i] = 0;
46     for(int k = 1; k <= v; k++)
47         for(int i = 1; i <= v; i++)
48             for(int j = 1; j <= v; j++)
49                 smin(dis[i][j], dis[i][k] + dis[k][j]);
50 }
51 
52 inline void solve() {
53     for(int i = 1; i <= n; i++)
54         for(int j = 0; j <= m; j++)
55             f[i][j][0] = f[i][j][1] = 1e100;
56     f[1][0][0] = f[1][1][1] = 0;
57     for(int i = 2; i <= n; i++)
58         for(int j = 0; j <= m && j <= i; j++) {
59             int &dcc = dis[c[i - 1]][c[i]], &dcd = dis[c[i - 1]][d[i]];
60             int &ddc = dis[d[i - 1]][c[i]], &ddd = dis[d[i - 1]][d[i]];
61             double &ip = p[i - 1], &np = p[i];
62             if(j <= i - 1)
63                 f[i][j][0] = min(f[i - 1][j][0] + dcc, f[i - 1][j][1] + dcc * (1 - ip) + ddc * ip);
64             if(j) {
65                 smin(f[i][j][1], f[i - 1][j - 1][0] + dcc * (1 - np) + dcd * np);
66                 smin(f[i][j][1], f[i - 1][j - 1][1] + (1 - ip) * (dcc * (1 - np) + dcd * np) + ip * (ddc * (1 - np) + ddd * np));
67             }
68         }
69     double ans = 1e100;
70     for(int i = 0; i <= m && i <= n; i++)
71         for(int k = 0; k < 2; k++)
72             smin(ans, f[n][i][k]);
73     printf("%.2lf", ans); 
74 }
75 
76 int main() {
77     init();
78     floyed();
79     solve();
80     return 0;
81 }

posted @ 2017-10-29 16:40 阿波罗2003 阅读(250) 评论(0) 编辑收藏举报

会员力量，点亮园子希望

刷新页面返回顶部

NOIP 2016 换教室 (luogu 1850 & uoj 262) - 概率与期望 - 动态规划

题目描述

输入输出格式

输入输出样例

说明

Code

公告