【网络流24题】运输问题 题解

题目传送门

题目大意: n n n 个仓库 m m m 个商店,第 i i i 个仓库有 a i a_i ai 个货物,第 j j j 个商店需要 b j b_j bj 个货物,第 i i i 个仓库向第 j j j 个商店运输 1 1 1 个货物的代价为 c i , j c_{i,j} ci,j a , b a,b a,b 满足 ∑ a i = ∑ b j \sum a_i=\sum b_j ai=bj,求出满足所有商店需求的最小费用和最大费用。

题解

裸费用流,最大费用流可以通过费用取相反数变成最小费用流。

具体建图就是源点向仓库建,流量为 a i a_i ai 费用为 0 0 0,仓库向商店建,流量无限费用为 c i , j c_{i,j} ci,j,最后商店向汇点建,流量 b i b_i bi 费用为 0 0 0

代码如下:

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 210
#define inf 999999999

int n,m,S,T,a[maxn],b[maxn],c[maxn][maxn];
struct edge{int x,y,z,cost,next;};
edge e[maxn*maxn<<1];
int first[maxn],len;
void buildroad(int x,int y,int z,int cost)
{
	e[++len]=(edge){x,y,z,cost,first[x]};
	first[x]=len;
}
int q[maxn],st,ed,f[maxn],fa[maxn];
bool v[maxn];
bool SPFA()
{
	for(int i=1;i<=T;i++)f[i]=-inf;
	memset(fa,0,sizeof(fa));
	st=1;ed=2;q[st]=S;v[S]=true;f[S]=0;
	while(st!=ed)
	{
		int x=q[st++];st=st>T?1:st;v[x]=false;
		for(int i=first[x];i;i=e[i].next)
		{
			int y=e[i].y;
			if(e[i].z&&f[y]<f[x]+e[i].cost)
			{
				f[y]=f[x]+e[i].cost;fa[y]=i;
				if(!v[y])v[q[ed++]=y]=true,ed=ed>T?1:ed;
			}
		}
	}
	return fa[T];
}
void work(int type)
{
	memset(first,0,sizeof(first));len=1;
	for(int i=1;i<=n;i++)buildroad(S,i,a[i],0),buildroad(i,S,0,0);
	for(int i=1;i<=m;i++)buildroad(i+n,T,b[i],0),buildroad(T,i+n,0,0);
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)
	buildroad(i,j+n,inf,c[i][j]*type),buildroad(j+n,i,0,-c[i][j]*type);
	int ans=0;
	while(SPFA())
	{
		int min_flow=inf;
		for(int i=fa[T];i;i=fa[e[i].x])min_flow=min(min_flow,e[i].z);
		for(int i=fa[T];i;i=fa[e[i].x])
		e[i].z-=min_flow,e[i^1].z+=min_flow,ans+=min_flow*e[i].cost;
	}
	printf("%d\n",ans*type);
}

int main()
{
	scanf("%d %d",&n,&m);S=n+m+1,T=S+1;
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	for(int i=1;i<=m;i++)scanf("%d",&b[i]);
	for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)scanf("%d",&c[i][j]);
	work(-1);work(1);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值