BZOJ1003

[BZOJ1003]ZJOI2006物流运输
神奇的dp。。。

目录

思路

• 
  
• 首先，若我们知道任意一段时间$[i,j]$的最小花费$cost[i][j]$，很容易得出dp方程：
  

$$dp[i]=min_{j=1}^{i-1}\{dp[j]+cost[j+1][i]+k\}$$
表示到第i天时的最小花费(最开始 $dp[i]=cost[1][i]$，就是不该航线的情况)

• 那么我们就暴力跑最短路好了。。。求出$cost[i][j]$就好了，其实挺水的。

代码

#include <cstdio>
#include <queue>
#include <utility>
#include <vector>
#include <functional>
#include <algorithm>

using std::vector;
using std::greater;
using std::priority_queue;
using std::pair;
using std::min;

const int MAXN=105,MAXM=25,INF=1e9;
int n,m,k,e;
int cost[MAXN][MAXN],dp[MAXN];
bool cant[MAXM][MAXN];

priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > > hp;

struct node{
    int dis,he;
    bool use,tag;
}d[MAXM];
struct line{
    int to,nex,w;
}ed[MAXM*MAXM<<1];

inline void addE(int u,int v,int w){
    static int cnt=0;
    ed[++cnt]=(line){v,d[u].he,w};
    d[u].he=cnt;
}

inline int Dijsktra(){
    for(int i=1;i<=m;++i){
        d[i].dis=INF;
        d[i].use=false;
    }
    d[1].dis=0;
    hp.push(pair<int,int>(d[1].dis,1));
    pair<int,int> tmp;
    int &u=tmp.second;
    while(hp.size()){
        tmp=hp.top();
        hp.pop();
        if(d[u].use) continue;
        d[u].use=true;
        for(int i=d[u].he,v,w;i;i=ed[i].nex){
            v=ed[i].to,w=ed[i].w;
            if(d[v].dis>d[u].dis+w && !d[v].tag){
                d[v].dis=d[u].dis+w;
                hp.push(pair<int,int>(d[v].dis,v));
            }
        }
    }
    return d[m].dis;
}

inline void deal(){
    int tmpd;
    scanf("%d",&tmpd);
    for(int i=1,p,a,b;i<=tmpd;++i){
        scanf("%d%d%d",&p,&a,&b);
        for(int j=a;j<=b;++j)
            cant[p][j]=true;
    }
    for(int i=1;i<=n;++i){
        for(int j=i;j<=n;++j){
            for(int k=1;k<=m;++k){
                d[k].tag=false;
                for(int l=i;l<=j;++l)
                    d[k].tag|=cant[k][l];
            }
            cost[i][j]=Dijsktra();
            if(cost[i][j]<INF)
                cost[i][j]*=(j-i+1);
        }
    }
}

inline void DP(){
    for(int i=1;i<=n;++i)
        dp[i]=cost[1][i];
    for(int i=2;i<=n;++i){
        for(int j=1;j<i;++j)
            dp[i]=min(dp[i],dp[j]+cost[j+1][i]+k);
    }
}

int main(){
    scanf("%d%d%d%d",&n,&m,&k,&e);
    for(int i=1,u,v,w;i<=e;++i){
        scanf("%d%d%d",&u,&v,&w);
        addE(u,v,w),addE(v,u,w);
    }
    deal();
    DP();
    printf("%d",dp[n]);
    return 0;
}