设表示前
天的最优解,容易得到状态转移方程
其中表示从第
天到第
天走同一条路的最优解
时间复杂度
注意到第天与第
天到第
天的路径可能是相同的,这时
就比最优解多了一个
但是实际上总能找到一个,使第
天到第
天走同一条路且与第
天不同,所以不需要特别判断
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=105,maxm=25,maxe=405,INF=0x3f3f3f3f;
int n,k,f[maxn];bool tp[maxm],lt[maxm],lmt[maxn][maxm];
int m,e,tot,son[maxe],w[maxe],nxt[maxe],lnk[maxm],que[maxm],dst[maxm];bool vis[maxm];
inline int read(){
int ret=0,fh=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')fh=-fh;ch=getchar();}
while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();
return ret*fh;
}
void add_e(int x,int y,int z){son[++tot]=y,w[tot]=z,nxt[tot]=lnk[x],lnk[x]=tot;}
int SPFA(bool *V){
memset(dst,0x3f,sizeof(dst));
memset(vis,0,sizeof(vis));
que[1]=1,dst[1]=0,vis[1]=1;
int hd=0,tl=1;
while(hd!=tl){
int x=que[hd=(hd+1)%maxm];
vis[x]=0;
for(int j=lnk[x];j;j=nxt[j]){
int v=son[j];
if(V[v]) continue;
if(dst[x]+w[j]<dst[v]){
dst[v]=dst[x]+w[j];
if(!vis[v]){
vis[v]=1;
que[tl=(tl+1)%maxm]=v;
}
}
}
}
return dst[m];
}
int main(){
freopen("P1003.in","r",stdin);
freopen("P1003.out","w",stdout);
n=read(),m=read(),k=read(),e=read();
for(int i=1;i<=e;i++){
int x=read(),y=read(),z=read();
add_e(x,y,z),add_e(y,x,z);
}
for(int d=read();d;d--){
int p=read(),l=read(),r=read();
for(int i=l;i<=r;i++) lmt[i][p]=1;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++) tp[j]|=lmt[i][j],lt[j]=0;
int x=SPFA(tp);
if(x!=INF) f[i]=x*i;
else f[i]=INF;
for(int j=i;j>1;j--){
for(int t=1;t<=m;t++) lt[t]|=lmt[j][t];
int x=SPFA(lt);
if(x==INF) break;
f[i]=min(f[i],f[j-1]+x*(i-j+1)+k);
}
}
printf("%d\n",f[n]);
return 0;
}