[luogu-1850]noip2016day1-T3 换教室 题解

题目传送门
题意解析：题目就是给了你相对应的2n个教室，m次换的机会。对于每一次去教师上课，你都可以申请或者不申请换教室，换的话就会浪费一次机会，当然每次换都有成功或者不成功的概率。然后题目一开始就告诉你了，每次你在哪个教室上课，还有换的话去哪个教师，以及每次的成功概率，还有每两个教室之间的距离。最后求最小的期望，每次期望就是路径长度乘以概率。

---

My opinion：一开始在比赛看到这题………………只会暴力。暴力还有80分，所以正解无用。现在又来做这题，可以明显得看出来每个教室之间距离一定要是最小的，两个教室之间到达的概率是不变的，为了期望最小，所以走最短路一定是最优的。那么最短路就十分明显了，而且v才300，完全可以用floyd。然后像我这种菜鸡就只会暴力了就有dalao点醒于我。瞬间我就燃了题目瞬间变简单了，每次上课，到达这个教室的期望只跟前面那个上课的教室有关，所以完全可以用上dp，状态也十分明显（至少比day2-T3明显），dp[i][j][0/1]表示第i次课，使用了j次换的机会，这次换成功/不成功的最小期望（这里我用0表示失败，1表示成功）。
总结：
1、输入。
2、预处理出每两间教室之间的距离。
3、dp一次。
4、输出。
dp方程式及其恶心长

---

有一行特别长的代码：

#include<iostream> 
#include<cstdio> 
#include<cmath> 
#include<cstring> 
#include<algorithm> 
#include<queue> 
#define rep(i,a,n) for (int i=a;i<=n;i++) 
#define per(i,a,n) for (int i=a;i>=n;i--) 
#define Clear(a,x) memset(a,x,sizeof(a)) 
#define ll long long 
#define INF 1000000000 
#define eps 1e-8 
#define db double 
using namespace std; 
ll read(){ 
    ll x=0,f=1; 
    char ch=getchar(); 
    while (ch<'0'||ch>'9') f=ch=='-'?-1:f,ch=getchar(); 
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar(); 
    return x*f; 
} 
const int maxn=2005,maxv=305,maxe=90005; 
int n,m,v,e; 
int c[maxn],d[maxn]; 
db g[maxn],f[maxn][maxn][2]; 
int map[maxv][maxv]; 
int main(){ 
    n=read(),m=read(),v=read(),e=read(); 
    rep(i,1,n) c[i]=read(); 
    rep(i,1,n) d[i]=read(); 
    rep(i,1,n) scanf("%lf",&g[i]); 
    rep(i,1,v) 
        rep(j,1,v)  
            if (i!=j) map[i][j]=INF; 
    rep(i,1,e){ 
        int u=read(),v=read(),w=read(); 
        map[u][v]=min(map[u][v],w); 
        map[v][u]=min(map[v][u],w); 
    } 
    rep(k,1,v) 
        rep(i,1,v) 
            rep(j,1,v) 
                if (map[i][j]>map[i][k]+map[k][j]) 
                    map[i][j]=map[i][k]+map[k][j]; 
    rep(i,1,n) 
        rep(j,0,m) 
            rep(k,0,1) 
                f[i][j][k]=INF; 
    f[1][0][0]=0,f[1][1][1]=0; 
    rep(i,2,n){ 
        int u=c[i],v=c[i-1],U=d[i],V=d[i-1]; 
        f[i][0][0]=f[i-1][0][0]+map[u][v]; 
        rep(j,1,min(i,m)){ 
            f[i][j][0]=min(f[i-1][j][0]+map[u][v],f[i-1][j][1]+g[i-1]*map[u][V]+(1-g[i-1])*map[u][v]); 
            f[i][j][1]=min(f[i-1][j-1][0]+g[i]*map[U][v]+(1-g[i])*map[u][v],f[i-1][j-1][1]+g[i]*g[i-1]*map[U][V]+g[i]*(1-g[i-1])*map[U][v]+g[i-1]*(1-g[i])*map[u][V]+(1-g[i])*(1-g[i-1])*map[u][v]); 
        } 
    } 
    db ans=f[n][0][0]; 
    rep(i,1,m) 
        ans=min(ans,min(f[n][i][0],f[n][i][1])); 
    printf("%.2lf\n",ans); 
    return 0; 
}