[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^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int inf=0x7fffffff;
const int N=2000,M=2001,V=301;
int c[N],d[N],dis[V][V];
double k[N],f[2][M][2];
int main() {
    int n=getint(),m=getint(),v=getint(),e=getint();
    for(register int i=0;i<n;i++) c[i]=getint();
    for(register int i=0;i<n;i++) d[i]=getint();
    for(register int i=0;i<n;i++) scanf("%lf",&k[i]);
    for(register int i=1;i<=v;i++) {
        for(register int j=1;j<=v;j++) {
            if(i!=j) dis[i][j]=inf;
        }
    }
    for(register int i=0;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=1;k<=v;k++) {
        for(register int i=1;i<=v;i++) {
            if(dis[i][k]==inf) continue;
            for(register int j=1;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=0;j<=m;j++) {
        f[0][j][0]=f[0][j][1]=1e9;
    }
    f[0][0][0]=f[0][1][1]=0;
    for(register int i=1;i<n;i++) {
        for(register int j=0;j<=m;j++) {
            f[i&1][j][0]=f[i&1][j][1]=1e9;
        }
        f[i&1][0][0]=std::min(f[!(i&1)][0][0]+dis[c[i-1]][c[i]],f[!(i&1)][0][1]+dis[c[i-1]][c[i]]*(1-k[i-1])+dis[d[i-1]][c[i]]*k[i-1]);
        for(register int j=1;j<=m;j++) {
            f[i&1][j][0]=std::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&1][j][1]=std::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]));
        }
    }
    double ans=1e9;
    for(register int i=0;i<=m;i++) {
        ans=std::min(ans,std::min(f[!(n&1)][i][0],f[!(n&1)][i][1]));
    }
    printf("%.2f
",ans);
    return 0;
}