BZOJ1003: [ZJOI2006]物流运输（洛谷P1772）

最短路 DP

BZOJ题目传送门
洛谷题目传送门

首先spfa预处理从$i$到$j$天用同一路线的最短路，记为$w[i][j]$。
$dp[i]$表示前 i i <script id="MathJax-Element-27" type="math/tex">i</script>天的最小代价。
那么有转移方程dp[i]=min(dp[i],dp[j]+w[j+1][i]+k)
代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 25
#define T 105
using namespace std;
struct edge{
    int next,to,dis;
}ed[N*N<<1];
int h[N],w[T][T],dp[T],n,m,c,k,t;
bool lmt[N][T],f[N];
void addedge(int x,int y,int z){
    ed[++k].next=h[x],ed[k].to=y,ed[k].dis=z,h[x]=k;
}
int spfa(){
    int d[N],que[N*N*N],r=0,w=1;
    memset(d,0x3f,sizeof(d));
    memset(que,0,sizeof(que));
    que[1]=1,d[1]=0;
    while (r<w){
        int x=que[++r]; f[x]=false;
        for (int i=h[x],v;i;i=ed[i].next)
            if (d[v=ed[i].to]>d[x]+ed[i].dis){
                d[v]=d[x]+ed[i].dis;
                if (!f[v]) que[++w]=v,f[v]=true;
            }
    }
    return d[n];
}
int main(){
    scanf("%d%d%d%d",&t,&n,&c,&m);
    memset(w,0x3f,sizeof(w));
    for (int i=1;i<=m;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedge(x,y,z),addedge(y,x,z);
    }
    int d; scanf("%d",&d);
    for (int i=1;i<=d;i++){
        int x,l,r;
        scanf("%d%d%d",&x,&l,&r);
        for (int j=l;j<=r;j++)
            lmt[x][j]=true;
    }
    for (int i=1;i<=t;i++)
        for (int j=i;j<=t;j++){
            memset(f,false,sizeof(f));
            for (int k=1;k<=n;k++)
                for (int p=i;p<=j;p++)
                    f[k]|=lmt[k][p];
            w[i][j]=spfa();
        }
    for (int i=1;i<=t;i++)
        for (int j=i;j<=t;j++)
            if (w[i][j]<1e9) w[i][j]*=(j-i)+1;
    for (int i=1;i<=t;i++) dp[i]=w[1][i];
    for (int i=2;i<=t;i++)
        for (int j=1;j<i;j++)
            dp[i]=min(dp[i],dp[j]+w[j+1][i]+c);
    return printf("%d\n",dp[t]),0;
}