BZOJ4720-换教室

题目很长，是一道概率dp题，一般需要逆推，但这题结局不确定所以要顺推。

用f[i][j][k]，i表示第i段时间，j表示用了j次申请，k就表示这轮是否用申请。

那么要求min(f[n][0~m][0],f[n][0~m][1])

状态转移方程有点麻烦，需要细心。

①前一个不需要申请，当前的时间也不需要申请。这个最容易推出来f[i][j][0]=min(f[i][j][0],（f[i-1][j][0]+dis[c[i-1]][c[i]])）×1.0)

②当前j>=1,那么至少可以申请一次。也就是可以前一个时间没申请，当前申请。或者前一个时间申请，当前没申请。

③当前j>=2，那么就可以考虑前一段时间也申请了。因为申请分成功和失败，所以这里有4个不同的状况要考虑。分别是：

申请成功->申请成功

申请成功->申请失败

申请失败->申请成功

申请失败->申请失败

代码中的_11s表示前一段时间申请了，这一次也申请了，并申请成功。1表示申请了，0表示没申请。s表示成功，f表示失败。

代码如下：

#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
typedef long long ll;

int n,m,v,e;
int c[],d[];
double k[];
int dis[][],a,b,cc;
double f[][][];
int  main(){
    scanf("%d%d%d%d",&n,&m,&v,&e);
    for(int i=;i<=n;++i)scanf("%d",&c[i]);
    for(int i=;i<=n;++i)scanf("%d",&d[i]);
    for(int i=;i<=n;++i)scanf("%lf",&k[i]);
    memset(dis,0x3f,sizeof dis);
    for(int i=;i<=v;++i)dis[i][i]=;
    for(int i=;i<=e;++i){
        scanf("%d%d%d",&a,&b,&cc);
        dis[a][b]=dis[b][a]=min(dis[a][b],cc);
    }
    for(int kk=;kk<=v;++kk)
        for(int i=;i<=v;++i)
            for(int j=;j<=v;++j)
                dis[i][j]=dis[j][i]=min(dis[i][j],dis[i][kk]+dis[kk][j]);
    for(int i=;i<=n;++i)for(int j=;j<=m;++j)
        f[i][j][]=f[i][j][]=0x3f3f3f3f;
    f[][][]=f[][][]=0.0;
    for(int i=;i<=n;++i){
        for(int j=;j<=min(i,m);++j){
            f[i][j][]=min(f[i][j][],f[i-][j][]+dis[c[i-]][c[i]]);
            if(j>=){
                double _10s=(f[i-][j][]+1.0*dis[d[i-]][c[i]])*k[i-];
                double _10f=(f[i-][j][]+1.0*dis[c[i-]][c[i]])*(1.0-k[i-]);
                f[i][j][]=min(f[i][j][],_10s+_10f);
                double _01s=(f[i-][j-][]+1.0*dis[c[i-]][d[i]])*k[i];
                double _01f=(f[i-][j-][]+1.0*dis[c[i-]][c[i]])*(1.0-k[i]);
                f[i][j][]=min(f[i][j][],_01s+_01f);
            }
            if(j>=){
                double _11ss,_11ff,_11sf,_11fs;
                _11ss=(f[i-][j-][]+dis[d[i-]][d[i]])*k[i-]*k[i];
                _11ff=(f[i-][j-][]+1.0*dis[c[i-]][c[i]])*(1.0-k[i-])*(1.0-k[i]);
                _11fs=(f[i-][j-][]+1.0*dis[c[i-]][d[i]])*(1.0-k[i-])*k[i];
                _11sf=(f[i-][j-][]+1.0*dis[d[i-]][c[i]])*k[i-]*(1.0-k[i]);
                double tmp=_11ss+_11ff+_11sf+_11fs;
                f[i][j][]=min(f[i][j][],tmp);
            }
        }
    }
    double ans=0x3f3f3f3f;
    for(int i=;i<=m;++i)
        ans=min(ans,min(f[n][i][],f[n][i][]));
    printf("%.2lf\n",ans);
    return ;
}