【题解】NOIP-2016 换教室

题目描述

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

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

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

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

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

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

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

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

输入输出格式

输入格式：

第一行四个整数 n,m,v,e n , m , v , e 。 n n 表示这个学期内的时间段的数量； m m 表示牛牛最多可以申请更换多少节课程的教室； v v 表示牛牛学校里教室的数量； e e 表示牛牛的学校里道路的数量。

第二行 n n 个正整数，第 i i （1≤i≤n) （ 1 ≤ i ≤ n ) 个正整数表示 ci c i ​，即第 i i 个时间段牛牛被安排上课的教室；保证 1≤ci≤v 1 ≤ c i ≤ v 。

第三行 n n 个正整数，第 i i （1≤i≤n1） （ 1 ≤ i ≤ n 1 ） 个正整数表示 di d i ​，即第 i i 个时间段另一间上同样课程的教室；保证 1≤di≤v1 1 ≤ d i ≤ v 1 。

第四行 n n 个实数，第 i i （1≤i≤n） （ 1 ≤ i ≤ n ） 个实数表示 ki k i ​，即牛牛申请在第 i i 个时间段更换教室获得通过的概率。保证 0≤ki≤1 0 ≤ k i ​ ≤ 1 。

接下来 e e 行，每行三个正整数 aj a j , bj b j , wj w j ​，表示有一条双向道路连接教室 aj a j ​ , bj b j ​，通过这条道路需要耗费的体力值是 wj w j ​；保证 1≤aj,bj≤v,1≤wj≤100 1 ≤ a j , b j ≤ v , 1 ≤ w j ≤ 100

保证 1≤n≤2000 1 ≤ n ≤ 2000 ， 0≤m≤20000 0 ≤ m ≤ 20000 ， 1≤v≤300 1 ≤ v ≤ 300 ， 0≤e≤900000 0 ≤ e ≤ 900000

保证通过学校里的道路，从任何一间教室出发，都能到达其他所有的教室。

保证输入的实数最多包含 3 3 位小数。

输出格式：

输出一行，包含一个实数，四舍五入精确到小数点后恰好 2 2 位，表示答案。你的输出必须和标准输出完全一样才算正确。

测试数据保证四舍五入后的答案和准确答案的差的绝对值不大于 4×10−3 4 × 10 − 3 （如果你不知道什么是浮点误差，这段话可以理解为：对于大多数的算法，你可以正常地使用浮点数类型而不用对它进行特殊的处理）

输入输出样例

输入样例#1：

3 2 3 3
2 1 2
1 2 1
0.8 0.2 0.5
1 2 5
1 3 3
2 3 1

输出样例#1：

2.80

说明

【提示】

道路中可能会有多条双向道路连接相同的两间教室。 也有可能有道路两端连接
的是同一间教室。
2.请注意区分n,m,v,e的意义, n不是教室的数量, m不是道路的数量。
特殊性质1:图上任意两点 ai a i ​, bi b i ​, ai≠bi a i ​ ≠ b i ​间,存在一条耗费体力最少的路径只包含一条道路。
特殊性质2:对于所有的 1≤i≤n 1 ≤ i ≤ n ， ki=1 k i = 1 。

题目概要

给定无向连通图，换教室的概率，求最小期望位移

简要思路

Floyd求图论部分，dp部分概率&期望

Floyd就不讲了

dp方程贼恶心

f[i][j][0] f [ i ] [ j ] [ 0 ] 表示前 i i 天换了 j j 次教室，第 i i 天未换教室的期望

f[i][j][1] f [ i ] [ j ] [ 1 ] 表示前 i i 天换了 j j 次教室，第 i i 天换教室的期望

方程 1 1 号：

f[i][j][0]=min(p[i−1]dis[d[i−1]][c[i]]+(1−p[i−1])dis[c[i−1]][c[i]]+f[i−1][j][1] f [ i ] [ j ] [ 0 ] = m i n ( p [ i − 1 ] d i s [ d [ i − 1 ] ] [ c [ i ] ] + ( 1 − p [ i − 1 ] ) d i s [ c [ i − 1 ] ] [ c [ i ] ] + f [ i − 1 ] [ j ] [ 1 ]

,dis[c[i−1]][c[i]]+f[i−1][j][0]); , d i s [ c [ i − 1 ] ] [ c [ i ] ] + f [ i − 1 ] [ j ] [ 0 ] ) ;

方程 2 2 号：(方程太长写一行显示不出来)

f[i][j][1]=min(f[i−1][j−1][0]+p[i]dis[c[i−1]][d[i]]+(1−p[i])dis[c[i−1]][c[i]], f [ i ] [ j ] [ 1 ] = m i n ( f [ i − 1 ] [ j − 1 ] [ 0 ] + p [ i ] d i s [ c [ i − 1 ] ] [ d [ i ] ] + ( 1 − p [ i ] ) d i s [ c [ i − 1 ] ] [ c [ i ] ] ,

f[i−1][j−1][1]+p[i−1]p[i]dis[d[i−1]][d[i]]+p[i−1](1−p[i])dis[d[i−1]][c[i]]+ f [ i − 1 ] [ j − 1 ] [ 1 ] + p [ i − 1 ] p [ i ] d i s [ d [ i − 1 ] ] [ d [ i ] ] + p [ i − 1 ] ( 1 − p [ i ] ) d i s [ d [ i − 1 ] ] [ c [ i ] ] +

(1−p[i−1])p[i]dis[c[i−1]][d[i]]+(1−p[i−1])(1−p[i])dis[c[i−1]][c[i]]); ( 1 − p [ i − 1 ] ) p [ i ] d i s [ c [ i − 1 ] ] [ d [ i ] ] + ( 1 − p [ i − 1 ] ) ( 1 − p [ i ] ) d i s [ c [ i − 1 ] ] [ c [ i ] ] ) ;

代码

#include
using namespace std;
#define cl1(x) memset(x,-1,sizeof(x))
#define rg register
#define INF 2000000001

template  inline void read(_Tp &x){char c11=getchar();x=0;
    while(c11<'0'||c11>'9')c11=getchar();while(c11>='0'&&c11<='9'){x=x*10+c11-'0';c11=getchar();}
}

const int maxn=2020,maxv=305,maxe=90050;
int n,m,e,v;
int c[maxn],d[maxn];
double p[maxn];
int dis[maxn][maxn];
double f[maxn][maxn][2];

void init();

void floyd(){
    for(rg int k=1;k<=v;++k)
    for(rg int i=1;i<=v;++i)if(i!=k)
    for(rg int j=1;j<=v;++j)if(i!=j and j!=k)
    if(dis[i][j]>dis[i][k]+dis[k][j])
    dis[i][j]=dis[i][k]+dis[k][j];
}

void dp(){
    for(rg int i=1;i<=n;++i)
    for(rg int j=0;j<=m;++j)
    f[i][j][0]=f[i][j][1]=INF;
    f[1][0][0]=f[1][1][1]=0;
    for(rg int i=2;i<=n;++i)
    for(rg int j=0;j<=i and j<=m;++j){
        int cc=dis[c[i-1]][c[i]],cd=dis[c[i-1]][d[i]],dc=dis[d[i-1]][c[i]],dd=dis[d[i-1]][d[i]];
        if(j<=i-1)f[i][j][0] = min( p[i-1]*dc + (1-p[i-1])*cc + f[i-1][j][1] , cc+f[i-1][j][0]);
        if(j)f[i][j][1] = min( f[i-1][j-1][0] + p[i]*cd + (1-p[i])*cc , f[i-1][j-1][1] + p[i-1]*p[i]*dd + p[i-1]*(1-p[i])*dc + (1-p[i-1])*p[i]*cd + (1-p[i-1])*(1-p[i])*cc);
    }
    double ans=INF;
    for(rg int j=0;j<=m&&j<=n;++j) ans=min(ans,min(f[n][j][0],f[n][j][1]));
    printf("%.2lf",ans);
}

int main(){
    init();
    floyd();
    dp();
    return 0;
}

void init(){
    read(n);read(m);read(v);read(e);
    for(rg int i=1;i<=n;++i)read(c[i]);
    for(rg int i=1;i<=n;++i)read(d[i]);
    for(rg int i=1;i<=n;++i)scanf("%lf",&p[i]);
    int A,B,C;
    cl1(dis);
    memset(dis,0x3f3f3f3f,sizeof(dis));
    for(rg int i=1;i<=v;++i)dis[i][i]=0;
    for(rg int i=1;i<=e;++i){
        read(A);read(B);read(C);
        if(A==B)continue;
        else dis[A][B]=dis[B][A]=min(dis[A][B],C);
    }
}