P1772 [ZJOI2006]物流运输

传送门（洛谷）

---

算法：最短路与动态规划结合

---

$c1[id][j]$表示$id$号港口在$j$天的时候是否是正常运行

$co[i][j]$表示从$i$到$j$天走同一条路线的最小花费

$f[i]$表示前$i$天的最小花费


转移方程：

$f[i]=min(f[i],f[j]+co[j+1][i]\ast (i-j)+delta)$

也就是在第$j+1$天进行换航线的转移

接下来主要是求$co[i][j]$
这个我们只需要在算第$i$天到第$j$天的时候把所有在这段时间类的所有不合法的港口标记一下就$OK$，然后跑最短路就行了

Code

#include<bits/stdc++.h>
#include<queue>
#define ll long long
#define rep(i,a,b) for(int (i)=(a);(i)<=(b);(i)++)
using namespace std;
const int maxn=1e4+10;
const int maxm=1e3+10;
int m,n,delta,k,cnt,d;
int head[maxn],c1[maxm][maxm],co[maxm][maxm],vis[maxn],close[maxn],dis[maxn];
ll f[maxn];

template <class t> inline void read(t &x) {
	int f=1;x=0;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
	while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
	x*=f;
}

struct node{
	int v,w,nex;
}e[maxn<<1];

void add(int u,int v,int w) {
	e[++cnt].v=v;
	e[cnt].w=w;
	e[cnt].nex=head[u];
	head[u]=cnt;
}

void readdata() {
	read(n),read(m),read(delta),read(k);
	rep(i,1,k) {
		int u,v,w;
		read(u),read(v),read(w);
		add(u,v,w);
		add(v,u,w);
	}
	read(d);
	rep(i,1,d) {
		int id,x,y;
		read(id),read(x),read(y);
		rep(j,x,y) c1[id][j]=true;
	}
}

void spfa() {
	memset(dis,0x3f,sizeof(dis));
	memset(vis,false,sizeof(vis));
	queue<int > q;
	dis[1]=0;
	q.push(1);
	while(q.size()) {
		int u=q.front();q.pop();
		vis[u]=false;
		for(int i=head[u];i;i=e[i].nex) {
			int v=e[i].v,w=e[i].w;
			if(close[v]) continue;
			if(dis[v]>dis[u]+w) {
				dis[v]=dis[u]+w;
				if(!vis[v]) {
					q.push(v);
					vis[v]=true;
				}
			}
		}
	}
}

void work() {
	rep(i,1,n)
	 rep(j,1,n) {
	 	memset(close,false,sizeof(close));
	 	rep(r,i,j)
	 	 rep(l,1,m) 
	 		if(c1[l][r]) close[l]=true;
	 		spfa();
	 		co[i][j]=dis[m];
	 }
	memset(f,0x3f,sizeof(f));
	rep(i,1,n) {
		f[i]=(ll)co[1][i]*i;
		for(int j=i-1;j>=0;j--) {
			f[i]=min(f[i],f[j]+(ll)co[j+1][i]*(i-j)+(ll)delta);
		}
	}
	printf("%lld\n",f[n]);
}

int main() {
	freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}