BZOJ1003: [ZJOI2006]物流运输（SPFA+DP）

设表示前天的最优解，容易得到状态转移方程

其中表示从第天到第天走同一条路的最优解

时间复杂度

注意到第天与第天到第天的路径可能是相同的，这时就比最优解多了一个

但是实际上总能找到一个，使第天到第天走同一条路且与第天不同，所以不需要特别判断

#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;
}