【最短路+dp】P1772 [ZJOI2006]物流运输

这道题目由于观察到边数和节点数都较少,所以求最短路的方式就有很多(选择了spfa)

可以先用N^{2}次spfa求出path[i][j] 即i-j天的最短路

然后就利用dp dp[i]表示前i天的最短路和,递推公式为dp[i]=min(dp[i],dp[j]+path[j+1][i]*(i-j)+k) j为0-i-1

 

代码

#include<bits/stdc++.h>
using namespace std;
const int inf=0x3f3f3f3f;
int n,m,k,e,d;
struct edge
{
	int to,nxt,v;
}G[10005];
int head[30],vis[30],dis[30],valid[30][105],now[30];
int cnt;
long long path[105][105],dp[105];
void add(int x,int y,int z)
{
	G[++cnt].nxt=head[x];G[cnt].v=z; G[cnt].to=y; head[x]=cnt;
}
void spfa()
{
	for(int i=1;i<=m;i++)
		dis[i]=inf,vis[i]=0;
	queue<int> q;
	q.push(1); dis[1]=0;
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		vis[x]=0;
		for(int i=head[x];i;i=G[i].nxt)
		{
			int to=G[i].to;
			if(now[to]) continue;
			if(dis[to]>dis[x]+G[i].v)
			{
				dis[to]=dis[x]+G[i].v;
				if(!vis[to])
				{
					vis[to]=1;
					q.push(to);
				}
			}
		}
	}
}
int main()
{
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
	scanf("%d%d%d%d",&n,&m,&k,&e);
	int x,y,z;
	for(int i=1;i<=e;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z); add(y,x,z);
	}
	scanf("%d",&d);
	for(int i=1;i<=d;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		for(int j=y;j<=z;j++) valid[x][j]=1;
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			memset(now,0,sizeof(now));
			for(int k=i;k<=j;k++)
				for(int a=1;a<=m;a++)
					if(valid[a][k]) now[a]=1;
			spfa();
			path[i][j]=dis[m];
		}
	memset(dp,inf,sizeof(dp));
	for(int i=1;i<=n;i++)
	{
		dp[i]=path[1][i]*i;
		for(int j=i-1;j>=0;j--)
			dp[i]=min(dp[i],dp[j]+path[j+1][i]*(i-j)+k);
	}
	printf("%lld\n",dp[n]);
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值