[ZJOI2006] 物流运输（最短路+dp）

[ZJOI2006] 物流运输

题目描述

物流公司要把一批货物从码头 A 运到码头 B。由于货物量比较大，需要 $n$ 天才能运完。货物运输过程中一般要转停好几个码头。

物流公司通常会设计一条固定的运输路线，以便对整个运输过程实施严格的管理和跟踪。由于各种因素的存在，有的时候某个码头会无法装卸货物。这时候就必须修改运输路线，让货物能够按时到达目的地。

但是修改路线是—件十分麻烦的事情，会带来额外的成本。因此物流公司希望能够订一个 $n$ 天的运输计划，使得总成本尽可能地小。

输入格式

第一行是四个整数 $n,m,k,e$。$n$ 表示货物运输所需天数，$m$ 表示码头总数，$k$ 表示每次修改运输路线所需成本，$e$ 表示航线条数。

接下来 $e$ 行每行是一条航线描述，包括了三个整数，依次表示航线连接的两个码头编号以及航线长度。其中码头 A 编号为 $1$，码头 B 编号为 $m$。单位长度的运输费用为 $1$。航线是双向的。

再接下来一行是一个整数 $d$，后面的 $d$ 行每行是三个整数 $p,a,b$。表示编号为 $p$ 的码头在 $[a,b]$ 天之内无法装卸货物。同一个码头有可能在多个时间段内不可用。但任何时间都存在至少一条从码头 A 到码头 B 的运输路线。

输出格式

包括了一个整数表示最小的总成本。
总成本为 $n$ 天运输路线长度之和 $+k\times$ 改变运输路线的次数。

样例 #1

样例输入 #1

5 5 10 8
  1 2 1
  1 3 3
  1 4 2
  2 3 2
  2 4 4
  3 4 1
  3 5 2
  4 5 2
  4
  2 2 3
  3 1 1
  3 3 3
  4 4 5

样例输出 #1

32

提示

【数据范围】
对于 $100\%$ 的数据，$1\le n\le 100$，$1\le m\le 20$, $1\le k\le 500$, $1\le e\le 200$。

【样例输入说明】

上图依次表示第 $1$ 至第 $5$ 天的情况，阴影表示不可用的码头。

【样例输出说明】

前三天走 $1\to 4\to 5$，后两天走 $1\to 3\to 5$，这样总成本为 $(2+2)\times 3+(3+2)\times 2+10=32$。

_NOI导刊2010提高（01）

思路

这道题的最短路相信大家都能看的出来。

• 但这道题最难的就是它有的港口时间性关闭，那我也不知道我什么时候转向，而且还要求最优的方案，那dp合适不过了。
  总结：当题目要求转折，且不确定什么时候转且要求最优方案，dp合适不过了。
• 因此，我们用f[i]表示在前i天的最小花费。
• 那么对于时间范围内关闭的我们该怎么维护呢？我们可以用个二维数组进行维护，比如：close[i][j]表示i点j时间关闭。然后当close是有值的状态，那么我们的occpy就得设为有值状态。（这样就把二维的转化成一维状态）。
• 那我们的状态转移是什么：f[i]=min(f[i],f[j]+k+(long long)cost[j+1][i]*(i-j));//换边和不换边，其中：j表示断点 意为从第j天更换航路 分为1-j + j-i + 更换费用 。
• 细节：dist不能加long long，要不然会出问题：就是memset 0x3f，这样的dist很大，如果dist本身每被更新，然后再乘以一个数，那就溢出了，就变成负数。

代码

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N = 1010,M = 2*N;

int e[M],ne[M],w[M],h[N],idx;
bool st[N];
long long cost[N][N];
long long f[N],occpy[N];
int close[N][N],dist[N];//港口关闭时间，close[i][j]表示i点j时间关闭
long long n,m,k,p,d;

void add(int a,int b,int c){
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

long long spfa(int x,int y){
    queue<int>q;
    memset(st,0,sizeof st);
    memset(dist,0x3f,sizeof dist);
    memset(occpy,0,sizeof occpy);
    
    for(int i=1;i<=m;i++){
        for(int j=x;j<=y;j++){
            if(close[i][j])occpy[i]=1;//从x到y有一天关闭也算关闭 保证共同性 
        }
    }
    
    q.push(1);
    st[1]=true;
    dist[1]=0;
    
    while(q.size()){
        int t=q.front();
        q.pop();
        st[t]=false;
        
        for(int i=h[t];~i;i=ne[i]){
            int j=e[i];
            
            if(occpy[j])continue;
            
            if(dist[j]>dist[t]+w[i]){
                dist[j]=dist[t]+w[i];
                if(!st[j]){
                    st[j]=true;
                    q.push(j);
                }
            }
        }
    }
    
    return dist[m];
}

signed main(){
    cin>>n>>m>>k>>p;
    
    memset(h,-1,sizeof h);
    
    for(int i=1;i<=p;i++){
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
        add(b,a,c);
    }
    
    cin>>d;
    
    for(int i=1;i<=d;i++){
        int p,a,b;
        cin>>p>>a>>b;
        for(int j=a;j<=b;j++){
            close[p][j]=1;
        }
    }
    
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            //因为题目有时间范围，所以这边要i-j（预处理）
            cost[i][j]=spfa(i,j);//用cost[i][j] 处理出从第i天到第j天共同的最短路的长度 
        }
    }
    
    
    for(int i=1;i<=n;i++){
        f[i]=(long long)cost[1][i]*i;
        for(int j=1;j<i;j++){
            f[i]=min(f[i],f[j]+k+(long long)cost[j+1][i]*(i-j));//换边和不换边
            //j表示断点  意为从第j天更换航路 分为1-j + j-i + 更换费用 
            // cout<<(i-j)<<endl;
        }
    }
    
    cout<<f[n];
    
    return 0;
    
    
}