转移方程很好想 f[i]=f[j]+cost[j+1][i]*(i-j)+k,其中f[i]表示前i天的最小花费,cost[i][j]表示从第i天到第j天可行的最短路,可预处理出来。
注意我们在选择第一条路线时,本不用加k,但我们的方程是看做从0转化了一次路线,已经加上了k,所以最后答案应为f[n]-k。
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<queue>
#define inf 10000000
using namespace std;
int n,m,x[100000],y[100000],z[100000],k,e,d,jishu=0;
int xuan[500],ying[500],num[500],f[200],cost[200][200],head[50],dis[50];
int xx,yy,zz;
bool flag[500];
queue<int> q;
struct bian
{
int qi,zhong,w,next;
};
bian c[200000];
bool check(int i,int k,int j)
{
if(ying[j]<i) return 1;
if(xuan[j]>k) return 1;
return 0;
}
void add(int x,int y,int z)
{
c[++jishu].qi=x;
c[jishu].zhong=y;
c[jishu].w=z;
c[jishu].next=head[x];
head[x]=jishu;
}
void spfa(int x)
{
dis[x]=0;
flag[x]=1;
q.push(x);
while(!q.empty())
{
int t=q.front();
q.pop();
for(int i=head[t];i;i=c[i].next)
if(dis[c[i].zhong]>dis[t]+c[i].w)
{
dis[c[i].zhong]=dis[t]+c[i].w;
if(!flag[c[i].zhong])
{
flag[c[i].zhong]=1;
q.push(c[i].zhong);
}
}
flag[t]=0;
}
}
int main()
{
cin>>n>>m>>k>>e;
for(int i=0;i<=n;++i)
f[i]=inf;
for(int i=1;i<=e;++i)
cin>>x[i]>>y[i]>>z[i];
cin>>d;
for(int j=1;j<=d;++j)
cin>>num[j]>>xuan[j]>>ying[j];
for(int i=1;i<=n;++i)
{
for(int k=i;k<=n;++k)
{
bool vis[50];
for(int j=1;j<=m;++j)
vis[j]=1;
vis[1]=vis[m]=1;
jishu=0;
for(int j=1;j<=m;++j)
dis[j]=inf;
memset(flag,0,sizeof(flag));
memset(head,0,sizeof(head));
for(int j=1;j<=d;++j)
if(!check(i,k,j))
vis[num[j]]=0;
for(int j=1;j<=e;++j)
{
if(vis[x[j]]&&vis[y[j]])
{
add(x[j],y[j],z[j]);
add(y[j],x[j],z[j]);
}
}
spfa(1);
cost[i][k]=dis[m];
}
}
f[0]=0;
for(int i=1;i<=n;++i)
for(int j=0;j<=i-1;++j)
{
if(cost[j+1][i]!=inf&&f[j]!=inf)
f[i]=min(f[i],f[j]+cost[j+1][i]*(i-j)+k);
}
cout<<f[n]-k;
return 0;
}