[NOIp2016提高组]换教室

题目大意：
　　有n节课，第i节课在c[i]上课，同时d[i]也有一节课d[i]。
　　你有权利向教务处发出m次申请把自己的教室改到d[i]，相应的批准概率是k[i]。
　　教室是图上的一些点，其中每条边都有边权。
　　问你上完所有课走的总路程的期望。

思路：
　　动态规划。
　　f[i][j][0/1]表示上了i个课，换了j次教室，这门课有没有换教室。
　　转移方程为：
　　f[i][j][0]=min(f[i-1][j][0]+dis[c[i-1]][c[i]],f[i-1][j][1]+dis[c[i-1]][c[i]]*(1-k[i-1])+dis[d[i-1]][c[i]]*k[i-1]);
　　f[i][j][1]=min(f[i-1][j-1][0]+dis[c[i-1]][d[i]]*k[i]+dis[c[i-1]][c[i]]*(1-k[i]),f[i-1][j-1][1]+dis[d[i-1]][d[i]]*k[i-1]*k[i]+dis[d[i-1]][c[i]]*k[i-1]*(1-k[i])+dis[c[i-1]][d[i]]*(1-k[i-1])*k[i]+dis[c[i-1]][c[i]]*(1-k[i-1])*(1-k[i]));

 #include<cstdio>
 #include<cctype>
 #include<algorithm>
 inline int getint() {
     register char ch;
     while(!isdigit(ch=getchar()));
     register int x=ch^'';
     while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int inf=0x7fffffff;
 const int N=,M=,V=;
 int c[N],d[N],dis[V][V];
 double k[N],f[][M][];
 int main() {
     int n=getint(),m=getint(),v=getint(),e=getint();
     for(register int i=;i<n;i++) c[i]=getint();
     for(register int i=;i<n;i++) d[i]=getint();
     for(register int i=;i<n;i++) scanf("%lf",&k[i]);
     for(register int i=;i<=v;i++) {
         for(register int j=;j<=v;j++) {
             if(i!=j) dis[i][j]=inf;
         }
     }
     for(register int i=;i<e;i++) {
         const int u=getint(),v=getint(),w=getint();
         dis[u][v]=dis[v][u]=std::min(dis[u][v],w);
     }
     for(register int k=;k<=v;k++) {
         for(register int i=;i<=v;i++) {
             if(dis[i][k]==inf) continue;
             for(register int j=;j<=v;j++) {
                 if(dis[k][j]==inf) continue;
                 dis[i][j]=std::min(dis[i][j],dis[i][k]+dis[k][j]);
             }
         }
     }
     for(register int j=;j<=m;j++) {
         f[][j][]=f[][j][]=1e9;
     }
     f[][][]=f[][][]=;
     for(register int i=;i<n;i++) {
         for(register int j=;j<=m;j++) {
             f[i&][j][]=f[i&][j][]=1e9;
         }
         f[i&][][]=std::min(f[!(i&)][][]+dis[c[i-]][c[i]],f[!(i&)][][]+dis[c[i-]][c[i]]*(-k[i-])+dis[d[i-]][c[i]]*k[i-]);
         for(register int j=;j<=m;j++) {
             f[i&][j][]=std::min(f[!(i&)][j][]+dis[c[i-]][c[i]],f[!(i&)][j][]+dis[c[i-]][c[i]]*(-k[i-])+dis[d[i-]][c[i]]*k[i-]);
             f[i&][j][]=std::min(f[!(i&)][j-][]+dis[c[i-]][d[i]]*k[i]+dis[c[i-]][c[i]]*(-k[i]),f[!(i&)][j-][]+dis[d[i-]][d[i]]*k[i-]*k[i]+dis[d[i-]][c[i]]*k[i-]*(-k[i])+dis[c[i-]][d[i]]*(-k[i-])*k[i]+dis[c[i-]][c[i]]*(-k[i-])*(-k[i]));
         }
     }
     double ans=1e9;
     for(register int i=;i<=m;i++) {
         ans=std::min(ans,std::min(f[!(n&)][i][],f[!(n&)][i][]));
     }
     printf("%.2f\n",ans);
     return ;
 }