luogu P4014 分配问题

关于构图:

首先,这是一道最小费用最大流和最大费用最大流(就是最小费用最大流,但建边时的费用取负,就可以保证得到一个最小的费用,再取反,就得到最大费用)的题目,比较裸。

构图:

1.首先有一个超级源点和汇点:st=0,ed=2*n+1。

2.源点向每一个人建一条流量为1,费用为0的边,表示每一个人只能选一份工作。

3.每一份工作向汇点建一条流量为1,费用为0边,表示每一份工作只能有一个人来完成。

4.每一个人向每一个工作建一条流量为1,费用为x的边,表示当前这个人完成这份工作的费用x。

跑两边费用流(最小和最大)即可。


关于代码:

#include<cstdio>
#include<queue>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std;
queue<int> f;
	int n,m,len,st,ed;
	struct node{int x,y,c,d,next;} a[10000];
	int b[200][200],last[10000],dis[10000],pre[10000],pos[10000],p[10000];
	bool bz[10000];
void ins(int x,int y,int c,int d)
{
	a[++len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;a[len].next=last[x];last[x]=len;
}
void init(int x)
{
	len=-1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
	{
		ins(st,i,1,0),ins(i,st,0,0);
		ins(i+n,ed,1,0),ins(ed,i+n,0,0);
		for(int j=1;j<=n;j++)
			ins(i,j+n,1,b[i][j]*x),ins(j+n,i,0,-b[i][j]*x);
	}
	
}
bool spfa()
{
	memset(bz,true,sizeof(bz));
	bz[st]=false;
	memset(dis,63,sizeof(dis));
	dis[st]=0;
	p[st]=INF;
	f.push(st);
	while(!f.empty())
	{
		int x=f.front();
		bz[x]=true;
		for(int i=last[x];i>=0;i=a[i].next)
		{
			int y=a[i].y;
			if(a[i].c>0&&dis[y]>dis[x]+a[i].d)
			{
				dis[y]=dis[x]+a[i].d;
				pos[y]=x;
				pre[y]=i;
				p[y]=min(p[x],a[i].c);
				if(bz[y])
				{
					f.push(y);
					bz[y]=false;
				}
			}
		}
		f.pop();
	}
	return dis[ed]<1061109567;
}
int flow()
{
	int ans=0;
	while(spfa())
	{
		ans+=p[ed]*dis[ed];
		for(int i=ed;i!=st;i=pos[i])
		{
			a[pre[i]].c-=p[ed];
			a[pre[i]^1].c+=p[ed];
		}
	}
	return ans;
}
int main()
{
	scanf("%d",&n);
	st=0,ed=2*n+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			scanf("%d",&b[i][j]);
	init(1);
	printf("%d\n",flow());
	init(-1);
	printf("%d",-flow());
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值