P1850 换教室

题目描述

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

在可以选择的课程中，有 2n2n 节课程安排在 n 个时间段上。在第 i（\(1 \leq i \leq n\)）个时间段上，两节内容相同的课程同时在不同的地点进行，其中，牛牛预先被安排在教室 ci 上课，而另一节课程在教室 di 进行。

在不提交任何申请的情况下，学生们需要按时间段的顺序依次完成所有的 n 节安排好的课程。如果学生想更换第 i 节课程的教室，则需要提出申请。若申请通过，学生就可以在第 i 个时间段去教室 di 上课，否则仍然在教室 ci 上课。

由于更换教室的需求太多，申请不一定能获得通过。通过计算，牛牛发现申请更换第 i 节课程的教室时，申请被通过的概率是一个已知的实数 ki，并且对于不同课程的申请，被通过的概率是互相独立的。

学校规定，所有的申请只能在学期开始前一次性提交，并且每个人只能选择至多 m 节课程进行申请。这意味着牛牛必须一次性决定是否申请更换每节课的教室，而不能根据某些课程的申请结果来决定其他课程是否申请；牛牛可以申请自己最希望更换教室的 m 门课程，也可以不用完这 m 个申请的机会，甚至可以一门课程都不申请。

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

牛牛所在的大学有 v 个教室，有 e 条道路。每条道路连接两间教室，并且是可以双向通行的。由于道路的长度和拥堵程度不同，通过不同的道路耗费的体力可能会有所不同。 当第 i（\(1 \leq i \leq n-1\)）节课结束后，牛牛就会从这节课的教室出发，选择一条耗费体力最少的路径前往下一节课的教室。

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

输入格式

​ 第一行四个整数 n,m,v,e。n 表示这个学期内的时间段的数量；m 表示牛牛最多可以申请更换多少节课程的教室；v 表示牛牛学校里教室的数量；e 表示牛牛的学校里道路的数量。

​ 第二行 n 个正整数，第 iii（1≤i≤n）个正整数表示 ci，即第 i 个时间段牛牛被安排上课的教室；保证 1≤ci≤v。

​ 第三行 n 个正整数，第 iii（1≤i≤n）个正整数表示 di，即第 i 个时间段另一间上同样课程的教室；保证 1≤di≤v。

​ 第四行 n 个实数，第 iii（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 位小数。

输出格式

输出一行，包含一个实数，四舍五入精确到小数点后恰好2位，表示答案。你的输出必须和标准输出完全一样才算正确。

测试数据保证四舍五入后的答案和准确答案的差的绝对值不大于 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

​ 首先用\(Floyd\)求出任意两点之间的距离。

​ 用\(f[i][j][0/1]\)表示前i节课，用了j次交换机会，第i节课是/否要换的期望总路程。

​ 那么转移方程就是：

for(int i = 2;i <= n; i++) {
        for(int j = 0;j <= min(m, i); j++) {
            f[i][j][0] = min(f[i - 1][j][0] + dis[c[i - 1]][c[i]], 
                        f[i - 1][j][1] + (1 - k[i - 1]) * dis[c[i - 1]][c[i]] + k[i - 1] * dis[d[i - 1]][c[i]]);
            if(j != 0)  
            f[i][j][1] = min(f[i - 1][j - 1][0] + (1 - k[i]) * dis[c[i - 1]][c[i]] + k[i] * dis[c[i - 1]][d[i]], 
                        f[i - 1][j - 1][1] + (1 - k[i]) * (1 - k[i - 1]) * (dis[c[i - 1]][c[i]]) +
                        (1 - k[i]) * k[i - 1] * dis[d[i - 1]][c[i]] + k[i] * (1 - k[i - 1]) * dis[c[i - 1]][d[i]] +
                        k[i] * k[i - 1] * dis[d[i - 1]][d[i]]);
        }
}  

​ 第i个教室不换的期望是：min（上次也不换的期望+这两次之间的路程，上次换了之后的期望+不成功换的概率\(*\)上次不换和这次的路程+成功的概率\(*\)上次换和这次的路程）

​ 第i个教室换也是同理，只不过稍微麻烦一些。

#include <iostream>
#include <cstdio>
#include <cctype>
#include <cstring>

using namespace std;

inline long long read() {
    long long s = 0, f = 1; char ch;
    while(!isdigit(ch = getchar())) (ch == '-') && (f = -f);
    for(s = ch ^ 48;isdigit(ch = getchar()); s = (s << 1) + (s << 3) + (ch ^ 48));
    return s * f;
}

const int N = 2005;
int n, m, v, e;
int c[N], d[N];
double k[N], f[N][N][2], dis[N][N];

void init() {
    n = read(); m = read(); v = read(); e = read();
    for(int i = 1;i <= n; i++) c[i] = read();
    for(int i = 1;i <= n; i++) d[i] = read();
    for(int i = 1;i <= n; i++) scanf("%lf", &k[i]);
    for(int i = 1;i <= v; i++) 
        for(int j = 1;j < i; j++) 
            dis[j][i] = dis[i][j] = 9e9;
    for(int i = 1, x, y;i <= e; i++) {
        x = read(); y = read(); 
        double z = read(); 
        dis[x][y] = dis[y][x] = min(dis[x][y], z);
    }
    for(int i = 1;i <= n; i++) dis[i][i] = 0;
    for(int k = 1;k <= v; k++) 
        for(int i = 1;i <= v; i++) 
            for(int j = 1;j < i; j++) 
                dis[j][i] = dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
}

void work() {
    for(int i = 0;i <= n; i++) 
        for(int j = 0;j <= m; j++) 
            f[i][j][0] = f[i][j][1] = 9e9;
    f[1][0][0] = f[1][1][1] = 0;
    for(int i = 2;i <= n; i++) {
        for(int j = 0;j <= min(m, i); j++) {
            f[i][j][0] = min(f[i - 1][j][0] + dis[c[i - 1]][c[i]], 
                        f[i - 1][j][1] + (1 - k[i - 1]) * dis[c[i - 1]][c[i]] + k[i - 1] * dis[d[i - 1]][c[i]]);
            if(j != 0)  
            f[i][j][1] = min(f[i - 1][j - 1][0] + (1 - k[i]) * dis[c[i - 1]][c[i]] + k[i] * dis[c[i - 1]][d[i]], 
                        f[i - 1][j - 1][1] + (1 - k[i]) * (1 - k[i - 1]) * (dis[c[i - 1]][c[i]]) +
                        (1 - k[i]) * k[i - 1] * dis[d[i - 1]][c[i]] + k[i] * (1 - k[i - 1]) * dis[c[i - 1]][d[i]] +
                        k[i] * k[i - 1] * dis[d[i - 1]][d[i]]);
        }
    }        
    double ans = 9e9;
    for(int i = 0;i <= m; i++)
        ans = min(ans, min(f[n][i][0], f[n][i][1]));
    printf("%.2lf", ans);
}

int main() {

    init();
    work();

    return 0;
}