[网络流24题-10]分配问题

最新推荐文章于 2021-05-10 23:32:45 发布

hanyuweining

最新推荐文章于 2021-05-10 23:32:45 发布

阅读量235

点赞数

CC 4.0 BY-SA版权

分类专栏：题解 ————图论———— 网络流

本文链接：https://blog.csdn.net/hanyuweining/article/details/84556429

题解同时被 3 个专栏收录

94 篇文章

订阅专栏

————图论————

31 篇文章

订阅专栏

网络流

25 篇文章

订阅专栏

本文详细介绍了一种解决分配问题的高效算法——最小费用最大流（Minimum Cost Maximum Flow，MFMC），并提供了完整的C++代码实现。通过构建图模型，算法能够求解资源分配问题，以最小的代价达到最大的效益。代码中包含了SPFA算法用于寻找增广路径，并通过调整流量来优化总成本。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

分配问题

裸的MFMC

直接源连人人连工工连汇用1限制配对然后直接MFMC就好啦

终于自己建了个模然后s和t建小了QAQ

附代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define inf 20021225
#define ll long long
using namespace std;

struct edge{int to,lt,f,c,fr;}e[40010];
queue<int> que;bool vis[210];
int in[210],s,t,cnt=1,dis[210],nn,from[210];
void add(int x,int y,int f,int c)
{
	e[++cnt].to=y;e[cnt].fr=x;e[cnt].lt=in[x];e[cnt].f=f;e[cnt].c=c;in[x]=cnt;
	e[++cnt].to=x;e[cnt].fr=y;e[cnt].lt=in[y];e[cnt].f=0;e[cnt].c=-c;in[y]=cnt;
}
bool spfa()
{
	//while(!que.empty())	que.pop();
	for(int i=1;i<=nn;i++)	dis[i]=inf,vis[i]=0;
	dis[s]=0;vis[s]=1;que.push(s);
	while(!que.empty())
	{
		int x=que.front();que.pop();
		//printf("%d ",x);//printf("QAQ");
		for(int i=in[x];i;i=e[i].lt)
		{
			int y=e[i].to;
			if(e[i].f&&dis[y]>dis[x]+e[i].c)
			{
				dis[y]=dis[x]+e[i].c;from[y]=i;
				if(!vis[y])	vis[y]=1,que.push(y);
			}
		}
		vis[x]=0;
	}
	return dis[t]<inf;
}
int flow()
{
	int ff=inf;
	for(int i=t;i!=s;i=e[from[i]].fr)	ff=min(ff,e[from[i]].f);
	for(int i=t;i!=s;i=e[from[i]].fr)	e[from[i]].f-=ff,e[from[i]^1].f+=ff;
	return dis[t]*ff;
}
int mcmf()
{
	int ans=0;//printf("QAQ");
	while(spfa())	ans+=flow();	
	return ans;
}
int c[101][101],n;
int main()
{
	scanf("%d",&n);
	s=2*n+1,t=2*n+2;nn=t;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			scanf("%d",&c[i][j]),add(i,j+n,1,c[i][j]);
	for(int i=1;i<=n;i++)	add(s,i,1,0),add(i+n,t,1,0);
	printf("%d\n",mcmf());
	memset(in,0,sizeof(in));cnt=1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)	add(i,j+n,1,-c[i][j]);
	for(int i=1;i<=n;i++)	add(s,i,1,0),add(i+n,t,1,0);
	printf("%d\n",-mcmf());
	return 0;
}