bzoj1003 [ZJOI2006]物流运输

题目

预处理cost[i,j]表示第i天到第j天都跑这条路的最短路。也就是说保证路径上的所有点在i到j天都可以跑的前提下的最短路。
然后dp。f[i]表示前i天的最小费用，那么f[i]=min{f[j]+cost[j+1,i+k]} ，初始值f[i]=cost[1,i]。
先用spfa预处理。然后转移的时候O（1）就好了。

#include<bits/stdc++.h>
#define LL long long
#define inf 10000000
using namespace std;
LL first[21],nex[1001],to[1001],val[1001],size;
LL cost[101][101];
LL dis[21];
LL l[21],r[21];
bool E[21],kx[21][101];
LL f[101];
LL m,n,e,k,d,x,y,z;
void add(LL x,LL y,LL z)
{
    nex[size]=first[x];
    first[x]=size;
    to[size]=y;
    val[size]=z;
    size++;
}
bool pd(int x,int l,int r)
{
    for(int i=l;i<=r;i++)
        if(!kx[x][i])return false;
    return true;
}
LL spfa(LL s,LL t,LL dl,LL dr)
{
    for(int i=1;i<=m;i++)
        dis[i]=inf;
    memset(E,false,sizeof(E));
    dis[s]=0;E[s]=true;
    queue <int> Q;
    Q.push(s);
    while(!Q.empty())
    {
        LL x=Q.front();Q.pop();
        E[x]=false;
        for(LL i=first[x];i!=-1;i=nex[i])
        {
            LL u=to[i];
            if(pd(u,dl,dr))
            {
                if(dis[x]+val[i]<dis[u])
                {
                    dis[u]=dis[x]+val[i];
                    if(!E[u])
                    {
                        E[u]=true;
                        Q.push(u);
                    }
                }
            }
        }
    }
    return dis[t];
}
int main()
{
    scanf("%lld%lld%lld%lld",&n,&m,&k,&e);
    memset(first,-1,sizeof(first));
    for(LL i=1;i<=e;i++)
    {
        scanf("%lld%lld%lld",&x,&y,&z);
        add(x,y,z);
        add(y,x,z);
    }
    memset(kx,true,sizeof(kx));
    scanf("%lld",&d);
    for(LL i=1;i<=d;i++)
    {
        scanf("%lld%lld%lld",&x,&y,&z);
        for(int j=y;j<=z;j++)
            kx[x][j]=false;
    }
    for(LL i=1;i<=n;i++)
    {
        for(LL j=i;j<=n;j++)
        {
                cost[i][j]=spfa(1,m,i,j);
                cost[j][i]=cost[i][j];
        }
    }
//  for(int i=1;i<=n;i++)
//      cost[0][i]=cost[i][0]=cost[1][i];
//  for(int i=0;i<=n;i++)
//  {
//      for(int j=0;j<=n;j++)
//          cerr<<cost[i][j]<<" ";
//      cerr<<endl;
//  }
    memset(f,63,sizeof(f));
    f[0]=0;
    for(int i=1;i<=n;i++)
        for(int j=0;j<i;j++)
            f[i]=min(f[i],f[j]+(i-j)*cost[i][j+1]+k);
    cout<<f[n]-k;
    return 0;
}