【NOIP2016】换教室-最短路+DP

最新推荐文章于 2025-07-11 17:05:45 发布

原创最新推荐文章于 2025-07-11 17:05:45 发布 · 262 阅读

1 ·

CC 4.0 BY-SA版权

-------DP------- 同时被 3 个专栏收录

114 篇文章

订阅专栏

-------图论-------

99 篇文章

订阅专栏

概率与期望

17 篇文章

订阅专栏

这篇博客介绍了NOIP2016比赛中的一道题目，涉及大学课程申请和教室更换的问题。学生牛牛需要根据概率申请更换教室，以期望最小化在教室间移动耗费的体力值。文章详细解析了题意，输入输出格式，并提供了题解思路及代码实现，主要运用了最短路和动态规划的方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门:洛谷-换教室

题意

    对于刚上大学的牛牛来说，他面临的第一个问题是如何根据实际情况申请合适的课程。
    在可以选择的课程中，有 $2n$ 节课程安排在 $n$ 个时间段上。在第 $i（1 \leq i \leq n )$ 个时间段上，两节内容相同的课程同时在不同的地点进行，其中，牛牛预先被安排在教室 $c_i$ 上课，而另一节课程在教室 $d_i$ 进行。
在不提交任何申请的情况下，学生们需要按时间段的顺序依次完成所有的 $n$ 节安排好的课程。如果学生想更换第 $i$ 节课程的教室，则需要提出申请。若申请通过，学生就可以在第 $i$ 个时间段去教室 $d_i$ 上课，否则仍然在教室 $C_i$ 上课。
    由于更换教室的需求太多，申请不一定能获得通过。通过计算，牛牛发现申请更换第 $i$ 节课程的教室时，申请被通过的概率是一个已知的实数 $k_i$ ，并且对于不同课程的申请，被通过的概率是互相独立的。
    学校规定，所有的申请只能在学期开始前一次性提交，并且每个人只能选择至多 $m$ 节课程进行申请。这意味着牛牛必须一次性决定是否申请更换每节课的教室，而不能根据某些课程的申请结果来决定其他课程是否申请；牛牛可以申请自己最希望更换教室的 $m$ 门课程，也可以不用完这 $m$ 个申请的机会，甚至可以一门课程都不申请。
    因为不同的课程可能会被安排在不同的教室进行，所以牛牛需要利用课间时间从一间教室赶到另一间教室。
    牛牛所在的大学有 $v$ 个教室，有 $e$ 条道路。每条道路连接两间教室，并且是可以双向通行的。由于道路的长度和拥堵程度不同，通过不同的道路耗费的体力可能会有所不同。当第 $i（1 \leq i \leq n-1）$ 节课结束后，牛牛就会从这节课的教室出发，选择一条耗费体力最少的路径前往下一节课的教室。
    现在牛牛想知道，申请哪几门课程可以使他因在教室间移动耗费的体力值的总和的期望值最小，请你帮他求出这个最小值。

输入

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

第二行 $n$ 个正整数，第 $i （1 \leq i \leq n ）$ 个正整数表示 $c_i$ ,即第 $i$ 个时间段牛牛被安排上课的教室；保证 $1 \le c_i \le v$ 。
第三行 $n$ 个正整数，第 $i（1 \leq i \leq n）$ 个正整数表示 $d_i$ ，即第 $i$ 个时间段另一间上同样课程的教室；保证 $1 \le d_i \le v$ 。
第四行 $n$ 个实数，第 $i （1 \leq i \leq n）$ 个实数表示 $k_i$ ，即牛牛申请在第 $i$ 个时间段更换教室获得通过的概率。保证 $0 \le k_i \le 1$ 。
接下来 $e$ 行，每行三个正整数 $a_j, b_j, w_j$ ，表示有一条双向道路连接教室 $a_j, b_j$ ，通过这条道路需要耗费的体力值是 $w_j$ ；保证 $1 \le a_j, b_j \le v$ ， $1 \le w_j \leq100$ 。
保证 $1 \leq n \leq 2000$ ， $0 \leq m \leq 2000$ ， $1 \leq v \leq 300$ ， $0 \leq e \leq 90000$ 。

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

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

输出

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

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

题解

先跑一遍floyd，再dp，f[i][j][0]表示到第i阶段且当前教室不变的最小值。
水题。就是调试了很久。

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define db double
#define INF 0x7fffffff
using namespace std;
typedef long long ll;
const db eps=1e-5;
const int MX=2e3+10;
const int N=305;
const int M=2e5+10;
int n,m,V,E,tot;
int lc[MX<<2];ll dis[MX][MX];
db pro[MX],f[MX][MX][2],anss;
struct P{
    int tag,dd;
    bool operator <(const P& u)const{
       return dd<u.dd;
    }
}op;
priority_queue<P>Q;

inline int read()
{
    char ch=getchar();int x=0,tw=1;
    while(ch<'0' || ch>'9') {if(ch=='-') tw=-1;ch=getchar();}
    while(ch>='0' && ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*tw;
}

inline db min(db x,db y)
{
    return (x-y)>eps? y:x;
}

inline db max(db x,db y)
{
    return (x-y)>eps? x:y;
}

inline void dj()
{
    for(int kk=1;kk<=V;kk++)
     for(int i=1;i<=V;i++)
      for(int j=1;j<i;j++){
        if(dis[i][kk]+dis[kk][j]<dis[i][j]){
            dis[i][j]=dis[j][i]=dis[i][kk]+dis[kk][j];
        }
     }
}

int main(){
    //freopen("xx.txt","r",stdin);
    //freopen("wht.txt","w",stdout);
    n=read();m=read();V=read();E=read();
    for(int i=1;i<=V;i++)
      for(int j=1;j<=V;j++)
        dis[i][j]=INF;
    for(int i=1;i<=V;i++) dis[i][i]=0;
    for(int i=0;i<n;i++)
      for(int j=0;j<=m;j++)
        f[i][j][0]=f[i][j][1]=INF;
    f[0][0][0]=0;f[0][1][1]=0;//this mistake ori->f[0][0][0]=f[0][0][1]=0;
    for(int i=0;i<n;i++) lc[i<<1]=read();
    for(int i=0;i<n;i++) lc[i<<1|1]=read();
    for(int i=0;i<n;i++) scanf("%lf",&pro[i]);
    for(int i=1;i<=E;i++){
        int x=read(),y=read(),z=read();
        if(x==y) continue;
        if(dis[x][y]>z) dis[x][y]=dis[y][x]=z;
    }
    dj();
    for(int i=0;i<(n-1);i++){
        int qw= m>(i+1)? (i+1):m;
        for(int j=0;j<=qw;j++){
            if(j<m){
              db possibility=pro[i]*pro[i+1]*dis[lc[(i<<1)|1]][lc[((i+1)<<1)|1]]+pro[i]*(1-pro[i+1])*dis[lc[i<<1|1]][lc[(i+1)<<1]]+(1-pro[i])*pro[i+1]*dis[lc[i<<1]][lc[((i+1)<<1)|1]]+(1-pro[i])*(1-pro[i+1])*dis[lc[i<<1]][lc[(i+1)<<1]];
             if(j) f[i+1][j+1][1]=min(f[i+1][j+1][1],f[i][j][1]+possibility);
             if(j<=i) f[i+1][j+1][1]=min(f[i+1][j+1][1],f[i][j][0]+pro[i+1]*dis[lc[i<<1]][lc[((i+1)<<1)|1]]+(1-pro[i+1])*dis[lc[i<<1]][lc[(i+1)<<1]]);  
            }
            if(j<=i) f[i+1][j][0]=min(f[i+1][j][0],f[i][j][0]+dis[lc[i<<1]][lc[(i+1)<<1]]);
            if(j) f[i+1][j][0]=min(f[i+1][j][0],f[i][j][1]+pro[i]*dis[lc[i<<1|1]][lc[(i+1)<<1]]+(1-pro[i])*dis[lc[i<<1]][lc[(i+1)<<1]]);
            }
    }
    int ss= m>(n-1)? (n-1):m;
    anss=min(f[n-1][0][0],f[n-1][0][1]);
    for(int i=1;i<=ss;i++)
      anss=min(anss,min(f[n-1][i][0],f[n-1][i][1]));
    printf("%.2f\n",anss);
    return 0;
}