Noip2016 Day1 T3 换教室（Floyd + Dp）

题目描述

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

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

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

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

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

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

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

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

输入输出格式

输入格式

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

第二行n个正整数，第i (1 <= i<= n)个正整数表示ci，即第i个时间段牛牛被安排上课的教室；保证1 <= Ci <= v。

第三行n个正整数，第i (1<=i<=n)个正整数表示di，即第i个时间段另一间上同样课程的教室；保证1<=di<=v。

第四行n个实数，第i (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 x 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：图上任意两点ai , bi , ai ≠ bi间，存在一条耗费体力最少的路径只包含一条道路。
特殊性质2：对于所有的1 <= i <= n，ki = 1。

提示

【样例1说明】
所有可行的申请方案和期望收益如下表：

【提示】
道路中可能会有多条双向道路连接相同的两间教室。也有可能有道路两端连接的是 同一间 教室。

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

思路

讲真，看到这么长的题面我一开始是拒绝的，而且还是浮点类型…

而且确实也没有思路，只能先写个最短路骗一下n=1和m=0,1时的分（也是有32分的嘛）

然后我偷偷去看了一下别人的题解 实际上这道题的正解就是最短路加Dp

最短路还是很好写，因为最多只有300个教室我们完全可以用Floyd来写，四行代码解决

而动规方程则有点难想到也难实现，我们要建立f [N] [N] [2] （N>300）这样一个数组

第一维表示当前是第几节课，第二维表示目前换了几间教室，最后一维表示当前教室换不换

设原来的教室是 now[i]，可以去的教室是 to[i]，可能性是 k[i] 那么我们就可以得到方程

f[i][j][0] = min(f[i-1][j][0] + d[now[i-1]][now[i]], f[i-1][j][1] + d[now[i-1]][now[i]] (1 - k[i-1]) + d[to[i-1]][now[i]] k[i-1]);

f[i][j][1] = min(f[i-1][j-1][0] + d[now[i-1]][now[i]] (1.0-k[i]) + d[now[i-1]][to[i]] k[i], f[i-1][j-1][1] + d[now[i-1]][now[i]] (1.0-k[i-1]) (1.0-k[i]) + d[to[i-1]][now[i]] k[i-1] (1.0-k[i]) + d[now[i-1]][to[i]] (1.0-k[i-1]) k[i] + d[to[i-1]][to[i]] k[i-1] k[i])

直接看有点辣眼睛，建议复制到编译器上去理解（也可以看下面我的代码）

如果可以成功理解上面的方程那么这道题也就可以迎刃而解了

（其实初始化和Floyd那几行可以写一个函数增加代码可读性）

（但是会占用时间我就直接放在了主函数里面去了）

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=2005,V=305;
int n,m,v,e；
int now[N],to[N]; //now[]代表本来要去的教室，to[]代表上同一内容的另一教室
double dis[V][V]; //存两间教室之间的距离
double k[N],ans;  //k[]是更换教室的可能性，ans是最后的答案
double f[N][N][2];//f[][][]是期望值之和
int main() 
{
    int x,y,z;
    scanf("%d%d%d%d",&n,&m,&v,&e);

    //初始化 dis 和 f 数组  
    for(int i=0;i<V;i++)  
     for(int j=0;j<V;j++)
      if(i!=j) dis[i][j]=1e9;
    //一定要记得是从0到大写V啊，不是小写v
    for(int i=0;i<N;i++)
     for(int j=0;j<N;j++)
      f[i][j][1]=f[i][j][0]=1e9;
    f[1][0][0]=0,f[1][1][1]=0;  

    for(int i=1;i<=n;i++) scanf("%d",&now[i]);
    for(int i=1;i<=n;i++) scanf("%d",&to[i]);
    for(int i=1;i<=n;i++) scanf("%lf",&k[i]);

    for(int i=1;i<=e;i++) 
    {
        scanf("%d%d",&x,&y),scanf("%d",&z);
        if(x==y) continue;
        if(z<dis[x][y]) dis[x][y]=z; //去重边 
        dis[y][x]=dis[x][y];
    }

    //Floyd求最短路 
    for(int l=1;l<=v;l++)
     for(int i=1;i<=v;i++)
      for(int j=1;j<=v;j++)
       dis[i][j]=min(dis[i][j],dis[i][l]+dis[l][j]);

    for(int i=2;i<=n;i++)
    {
        f[i][0][0]=f[i-1][0][0]+dis[now[i-1]][now[i]];
        for(int j=1;j<=min(m,i);j++)
        {
            f[i][j][0]=min(f[i-1][j][0]+dis[now[i-1]][now[i]],f[i-1][j][1]+dis[now[i-1]][now[i]]*(1.0-k[i-1])+dis[to[i-1]][now[i]]*k[i-1]);
            f[i][j][1]=f[i-1][j-1][0]+dis[now[i-1]][now[i]]*(1.0-k[i])+dis[now[i-1]][to[i]]*k[i];
            double mid=f[i-1][j-1][1]+dis[now[i-1]][now[i]]*(1.0-k[i-1])*(1.0-k[i]);
            mid+=dis[to[i-1]][now[i]]*k[i-1]*(1.0-k[i]);
            mid+=dis[now[i-1]][to[i]]*(1.0-k[i-1])*k[i];
            mid+=dis[to[i-1]][to[i]]*k[i-1]*k[i];
            //实在是太长了所以分开写了 
            f[i][j][1]=min(f[i][j][1],mid);
        }
    } 

    ans=f[n][0][0];
    for(int i=1;i<=m;i++)
     ans=min(ans,min(f[n][i][0],f[n][i][1]));
    //找答案 

    printf("%.2lf",ans);
    return 0;
}