【bzoj 1003】物流运输（最短路+DP）

传送门biu~
用s[i][j]表示从i天到j天不改变运输路线的最小成本，f[i]代表到第i天所需最小的成本，显然f[i]可以由方程f[i]=min(s[1][i],min{f[j]+s[j][i]+k})推出。

#include<bits/stdc++.h>
using namespace std;
int n,m,k,e;
int head[25],nex[405],to[405],len[405],tp;
int spfa[105][105],mark[25][105],dis[25],in[25];
int s[105][105],f[105];
bool bj[25];
inline void add(int x,int y,int l){
    nex[++tp]=head[x];
    head[x]=tp;
    to[tp]=y;
    len[tp]=l;
}
inline void init(){
    scanf("%d%d%d%d",&n,&m,&k,&e);
    for(int i=1;i<=e;++i){
        int x,y,l;
        scanf("%d%d%d",&x,&y,&l);
        add(x,y,l);
        add(y,x,l);
    }
    int d;
    scanf("%d",&d);
    for(int i=1;i<=d;++i){
        int p,a,b;
        scanf("%d%d%d",&p,&a,&b);
        for(int j=a;j<=b;++j){
            mark[p][j]=1;
        }
    }
}
queue<int>q;
inline int SPFA(int a,int b){
    memset(bj,0,sizeof(bj));
    for(int i=a;i<=b;++i){
        for(int j=1;j<=m;++j){
            if(mark[j][i])      bj[j]=1;
        }
    }
    for(int i=2;i<=m;++i)       dis[i]=10000000;
    dis[1]=0;q.push(1);in[1]=1;
    while(!q.empty()){
        int x=q.front();q.pop();in[x]=0;
        for(int i=head[x];i;i=nex[i]){
            if(bj[to[i]])   continue;
            if(dis[x]+len[i]<dis[to[i]]){
                dis[to[i]]=dis[x]+len[i];
                if(!in[to[i]])  q.push(to[i]),in[to[i]]=1;
            }
        }
    }
    return dis[m]*(b-a+1);
}
inline void solve(){
    for(int i=1;i<=n;++i){
        for(int j=i;j<=n;++j){
            s[i][j]=SPFA(i,j);
        }
    }
    for(int i=1;i<=n;++i){
        f[i]=s[1][i];
        for(int j=1;j<=i-1;++j)
            f[i]=min(f[i],f[j]+s[j+1][i]+k);
    }
    printf("%d",f[n]);
}
int main(){
    init();
    solve();
    return 0;
}