换教室

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

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

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

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

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

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

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

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

 

一定要记住期望的线性性，就是一条道路的期望就是每段路的期望之和

代码有许多的细节，要注意

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=2e3+7;
const int INF=2e8+7;//double,0x7f7f7f7f太大了，最大可能2e8
int n,m,v,e,c[maxn],d[maxn],map[maxn][maxn];
double p[maxn],f[maxn][maxn][7],ans;
int main(){
  cin>>n>>m>>v>>e;
  for(int i=1;i<=n;i++)
    for(int j=0;j<=m;j++)
      f[i][j][0]=f[i][j][1]=INF;//这样不全是1或全是0的数不能用memset 
  f[1][0][0]=f[1][1][1]=0;
  for(int i=1;i<=v;i++)
    for(int j=1;j<i;j++)
      map[i][j]=map[j][i]=INF;//注意，这样就忽略了map[i][i]的情况 
  for(int i=1;i<=n;i++) cin>>c[i];
  for(int i=1;i<=n;i++) cin>>d[i];
  for(int i=1;i<=n;i++) cin>>p[i];
  for(int i=1;i<=e;i++){
    int u,v,w;cin>>u>>v>>w;
    map[u][v]=map[v][u]=min(map[u][v],w); 
  }
  for(int k=1;k<=v;k++){
    for(int i=1;i<=v;i++){
      for(int j=1;j<=v;j++){
        if(map[i][j]>map[i][k]+map[k][j])
          map[i][j]=map[i][k]+map[k][j];
      }
    }
  }
  for(int i=2;i<=n;i++){
    for(int j=0;j<=min(i,m);j++){
      f[i][j][0]=min(f[i][j][0],f[i-1][j][0]+(double)map[c[i-1]][c[i]]);
      f[i][j][0]=min(f[i][j][0],f[i-1][j][1]+(double)map[c[i-1]][c[i]]*(1-p[i-1])+(double)map[d[i-1]][c[i]]*p[i-1]);
      if(!j) continue;
      f[i][j][1]=min(f[i][j][1],f[i-1][j-1][0]+(double)map[c[i-1]][c[i]]*(1-p[i])+(double)map[c[i-1]][d[i]]*p[i]);
      f[i][j][1]=min(f[i][j][1],f[i-1][j-1][1]+(double)map[c[i-1]][c[i]]*(1-p[i-1])*(1-p[i])+(double)map[c[i-1]][d[i]]*(1-p[i-1])*p[i]
                    +(double)map[d[i-1]][c[i]]*p[i-1]*(1-p[i])+(double)map[d[i-1]][d[i]]*p[i-1]*p[i]);
    }
  }
  ans=(double)INF;
  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;
}