预处理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;
}