KM(n^4)模板

#include<stdio.h>
const int maxn=301;
int lx[maxn],ly[maxn],a[maxn][maxn],vx[maxn],vy[maxn],match[maxn],n;
bool dfs(int x){
	vx[x]=1;
	for(int i=1;i<=n;i++)
		if(!vy[i] && lx[x]+ly[i]==a[x][i]){
			vy[i]=1;
			if(!match[i] || dfs(match[i])){
				match[i]=x;
				return 1;
			}
		}
	return 0;
}
int main(){
	int i,j,k,m;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		int max=0;
		for(j=1;j<=n;j++){
			scanf("%d",&a[i][j]);
			max=max>a[i][j]?max:a[i][j];
		}
		lx[i]=max;		
	}
	for(i=1;i<=n;i++){			
		while(1){		
			for(j=1;j<=n;j++){vx[j]=0;vy[j]=0;}		
			if(dfs(i))break;
			int min=100000000;
			for(j=1;j<=n;j++)
			    if(vx[j])
			    	for(k=1;k<=n;k++)
			    		if(!vy[j])
			    			min=min>(lx[j]+ly[k]-a[j][k])?min:lx[j]+ly[k]-a[j][k];
			if(min==100000000){
				printf("-1\n");
				return 0;
			}
			for(j=1;j<=n;j++)
				if(vx[j])lx[j]-=min;
			for(j=1;j<=n;j++)
				if(vy[j])ly[j]+=min;			
		}	
	}
	int ans=0;
	for(i=1;i<=n;i++)ans+=a[match[i]][i];
	printf("%d\n",ans);	
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值