标签：Golden 10 int py UVA11383 KM link ans include

题目

给定 \(n*n\) 的矩阵，现在给每行安排一个权值 \(x_i\)，给每列安排一个权值 \(y_j\)，

使得 \(x_i+y_j\geq a_{i,j}\)，并且使 \(\sum_{i=1}^nx_i+y_i\) 最小。

分析

学过KM算法的话，就应该知道可以将 \(x_i\) 和 \(y_i\) 当成顶标，并且当 \(x_i+y_j=a_{i,j}\) 时取得最小值，

那么就转换成二分图最大权匹配，直接跑KM算法将匹配的边权和加起来就是答案。

代码

#include <cstdio>
#include <cctype>
#include <cmath>
#include <queue>
using namespace std;
const int N=511;
bool vx[N],vy[N]; queue<int>q;
int slack[N],lx[N],ly[N],G[N][N];
int px[N],py[N],link[N],n,ans;
int iut(){
	int ans=0,f=1; char c=getchar();
	while (!isdigit(c)) f=(c=='-')?-f:f,c=getchar();
	while (isdigit(c)) ans=ans*10+c-48,c=getchar();
	return ans*f;
}
void print(int ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
int min(int a,int b){return a<b?a:b;}
int max(int a,int b){return a>b?a:b;}
void adjust(int y){
	for (int _y;y;y=_y){
		_y=px[link[y]];
		px[link[y]]=y;
		py[y]=link[y]; 
	}
}
void bfs(int st){
	for (int i=1;i<=n;++i) slack[i]=0x3f3f3f3f,vx[i]=vy[i]=0;
	while (!q.empty()) q.pop();
	q.push(st);
	while (1){
		while (!q.empty()){
			int x=q.front();
			vx[x]=1,q.pop();
			for (int y=1;y<=n;++y)
			if (!vy[y]&&slack[y]>lx[x]+ly[y]-G[x][y]){
				slack[y]=lx[x]+ly[y]-G[x][y],link[y]=x;
				if (!slack[y]){
					vy[y]=1;
					if (!py[y]) {adjust(y); return;}
					    else q.push(py[y]);
				}
			}
		}
		int mn=0x3f3f3f3f;
		for (int i=1;i<=n;++i)
		    if (!vy[i]) mn=min(mn,slack[i]);
		for (int i=1;i<=n;++i){
			if (vx[i]) lx[i]-=mn;
		    if (vy[i]) ly[i]+=mn;
		        else slack[i]-=mn;
		}
		for (int i=1;i<=n;++i)
		if (!vy[i]&&!slack[i]){
			vy[i]=1;
			if (!py[i]) {adjust(i); return;}
			    else q.push(py[i]);
		}
	}
}
void KM(){
	for (int i=1;i<=n;++i){
		link[i]=ly[i]=px[i]=py[i]=0,lx[i]=-0x3f3f3f3f;
		for (int j=1;j<=n;++j)
		    lx[i]=max(lx[i],G[i][j]);
	}
	for (int i=1;i<=n;++i) bfs(i);
}
int main(){
	while (scanf("%d",&n)==1){
		for (int i=1;i<=n;++i)
		for (int j=1;j<=n;++j) G[i][j]=iut();
		KM(),ans=0;
		for (int i=1;i<=n;++i) ans+=lx[i],ans+=ly[i];
		for (int i=1;i<=n;++i) print(lx[i]),putchar(i==n?10:32);
		for (int i=1;i<=n;++i) print(ly[i]),putchar(i==n?10:32);
		print(ans),putchar(10);
	}
	return 0;
}

标签：Golden,10,int,py,UVA11383,KM,link,ans,include
来源： https://www.cnblogs.com/Spare-No-Effort/p/15957069.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。