bzoj 1003 物流运输 DP+Spfa

一看到题目，吓了一跳，看到数据范围……
spfa（j，i）表示第j天到第i天可用点的最短路。
spfa（i，j）
if(j!=0)dp[i]=min(dp[i],dp[j]+(i-j)*dis[m]+k);
else dp[i]=min(dp[i],dp[j]+(i-j)*dis[m]);

#include<cstdio>
#include<iostream>
#include<queue>
using namespace std;
int n,m,k,e;
bool use[205][105];
int s[205][105];
struct bian
{   int to,d;
}b[20005];
int fst[105],nxt[20005],tot=1;
void build(int f,int t,int d)
{
    b[++tot].to=t;
    b[tot].d=d;
    nxt[tot]=fst[f];
    fst[f]=tot;
}
int inf=1000000000;
int dis[25];bool vis[25];bool ab[25];
queue<int> q;
void spfa(int x,int t)
{
    for(int i=1;i<=m;i++)
    {
        if(s[i][t]-s[i][x]>0)ab[i]=0;
        else ab[i]=1;
    }
    for(int i=1;i<=m;i++)
        dis[i]=inf;
    dis[1]=0;vis[1]=1;q.push(1);
    int u,v;
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=fst[u];i;i=nxt[i])
        {
            v=b[i].to;
            if(!ab[v])  continue;
            if(dis[v]>dis[u]+b[i].d)
            {
                dis[v]=dis[u]+b[i].d;
                if(!vis[v])
                {
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
}
int dp[1050];
int main()
{
//  freopen("trans.in","r",stdin);
//  freopen("trans.out","w",stdout);
    int p,a,b;
    scanf("%d%d%d%d",&n,&m,&k,&e);
    int u,v,d;
    for(int i=1;i<=e;i++)
    {
        scanf("%d%d%d",&u,&v,&d);
        build(u,v,d);
        build(v,u,d);
    }
    scanf("%d",&d);
    for(int i=1;i<=d;i++)
    {
        scanf("%d%d%d",&p,&a,&b);
        for(int i=a;i<=b;i++)
            use[p][i]=1;
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            s[j][i]=s[j][i-1]+use[j][i];
    for(int i=1;i<=n;i++)
        dp[i]=inf;
    spfa(0,1);
    dp[1]=dis[m];
    for(int i=2;i<=n;i++)
    {
        for(int j=0;j<i;j++)
        {
            spfa(j,i);
            if(dis[m]!=inf)
            {
                if(j!=0)dp[i]=min(dp[i],dp[j]+(i-j)*dis[m]+k);
                else    dp[i]=min(dp[i],dp[j]+(i-j)*dis[m]);
            }
        }
    }
    printf("%d",dp[n]);
//  fclose(stdin);
//  fclose(stdout);
    return 0;
}