【题解】洛谷P1772 物流运输(动态规划 spfa)

对于这道题目我们可以将其分解为两个分问题来解决，首先求在任意一段时间内某些码头被禁用的情况下从1到m的最短路径，然后运用动态规划求出最小值。

对于最短路径，我们可以用spfa来求。因为码头有禁用的情况，所以这里我们开一个ban[i][j]数组，代表第i号码头在第j天被禁用。time[i][j]记录在第i到第j天里从1到m的最短路径。spfa中给出要求的开始时间和末尾时间，初始化can数组代表i号码头是否可以使用为true。接着循环，倘若在第ss到第tt天里，i号码头有禁用的时候，那么can[i]就==false（这里指的是不改变路径的情况，改变路径我们留在动态规划里解决）。接着是spfa的模板，注意在判断时应当加上当前点到目标点有边、当前点can==true，最后得到在天数条件限制下的最短路径。

对于动态规划，我们枚举天数n，初始化f数组代表第i天的最小花费f[i]=time[1][i]*i。接着枚举j从1到i-1,代表在其中的某一个时间里改变了航向，那么我们可以得到状态转移方程f[i]=min(f[i],f[j]+k+(ll)time[j+1][i]*(i-j))。最后输出第n天的花费f[n]即可

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=110; //day
const int maxm=25; //dock
bool can[maxm],ban[maxm][maxn],b[maxm];
int to[maxm][maxm],time[maxn][maxn];
int dis[maxm];
int team[maxm*4];
int s=0,t=0;
int n,m,k,e;
int spfa(int ss,int tt)
{
	s=t=0;
	memset(can,true,sizeof(can));
	memset(b,false,sizeof(b));
	
	for(int i=1;i<=m;i++)
	{
		dis[i]=1e9;
	}dis[1]=0;
	for(int i=1;i<=m;i++)
	{
		for(int j=ss;j<=tt;j++)
		{
			if(ban[i][j]) can[i]=false; 
		}
	}
	team[t]=1;
	t++;
	b[1]=true;
	
	while(s!=t)
	{
		int now=team[s];
		s++;
		s%=m;
		b[now]=false;
		for(int i=1;i<=m;i++)
		{
			if(to[now][i]&&can[i]&&dis[i]>dis[now]+to[now][i])
			{
				dis[i]=dis[now]+to[now][i];
				if(b[i]==false)
				{
					team[t]=i;
					t++;
					t%=m;
					b[i]=true;
				}
			}
		}
	}
	return dis[m];
}
long long f[maxn];
int main()
{
	cin>>n>>m>>k>>e;
	for(int i=1;i<=e;i++)
	{
		int x,y,z;
		cin>>x>>y>>z;
		to[x][y]=z;
		to[y][x]=z; 
	}
	int q; cin>>q;
	
	for(int i=1;i<=q;i++)
	{
		int x,y,z;
		cin>>x>>y>>z;
		for(int j=y;j<=z;j++)
		{
			ban[x][j]=true;
		}
	}
	
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			time[i][j]=spfa(i,j);
			
	for(int i=1;i<=n;i++)
	{
		f[i]=(ll)time[1][i]*i;
		for(int j=1;j<i;j++)
		{
			f[i]=min(f[i],f[j]+k+(ll)time[j+1][i]*(i-j));
		}
	}
	
	cout<<f[n]<<endl;
	return 0;
}