Hdu 2853（KM算法）

Hdu 2853

 

（1）题意：

n个公司匹配m个任务，每个公司接受某一个任务，都会获得一个值价值Eij，

每个公司只能接受一个任务，每个任务也只能被一个公司完成，

最终n个公司要完成n个任务，要求最终的价值之和最大。

现在给出一个当前每个公司完成的任务，询问要怎样修改每个公司的任务，

才能使这个公司的价值之和最大，并且修改的路线数量最小。

 

（2）思路：

很明显是匹配问题，而且要求权值之和最大。

但是限制修改路线数量最小。

解决方法是让每个权值Eij变为Eij*（n+1）倍，如果是原有的路线，就变为Eij*（n+1）+1；

然后求出最佳完备匹配，原来保留的边的数量是ans%(n+1),增加的价值是ans-原有的匹配权值之和。

具体的原理参见这篇文章

下面说一下我自己的理解：

只要考虑两种情况，

一、如果原有匹配与新匹配完全不一样，就将所有原有的边变为最佳匹配的边；

二、如果最佳匹配有很多种，而且原匹配很接近最佳匹配，或者就是最佳匹配，可以找一种策略

在选择最佳匹配时，尽量选取原来的匹配使匹配。可以在不影响结果的前提下尽量小的改变原来的边权值来实现这个策略。

如果将边权值扩大n+1倍，边的权值然后将原有的边权值+1，那1不会对最大权值边的选取产生影响，同时，

还可以尽量走原有的边。

 

（3）代码：
 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int INF = 1e9+10;
const int maxn = 55;
int mp[maxn][maxn],pre[maxn],lx[maxn],ly[maxn],vx[maxn],vy[maxn],cha,n,m;
int bh[maxn];
int MIN(int x,int y){
	return x<y?x:y; 
}
int MAX(int x,int y){
	return x>y?x:y;
}
bool dfs(int u){
	vx[u] = 1;
	for(int v=1;v<=m;v++)
	if(vy[v]==0&&mp[u][v]==lx[u]+ly[v]){
		vy[v] = 1;
		if(pre[v]==-1||dfs(pre[v])){
			pre[v] = u;
			return true;
		}
	}
	else if(vy[v]==0){
		cha = MIN(cha,lx[u]+ly[v]-mp[u][v]);
	}
	return false;
}
void km(int w){
	memset(ly,0,sizeof(ly));
	for(int i=1;i<=n;i++){
		lx[i] = mp[i][1];
		for(int j=2;j<=m;j++) lx[i] = MAX(lx[i],mp[i][j]);
	}
	memset(pre,-1,sizeof(pre));
	for(int i=1;i<=n;i++){
		while(true){
			memset(vx,0,sizeof(vx));
			memset(vy,0,sizeof(vy));
			cha = INF;
			if(dfs(i)) break;
			for(int j=1;j<=n;j++)
			if(vx[j]==1) lx[j] -= cha;
			for(int j=1;j<=m;j++)
			if(vy[j]==1) ly[j] += cha;
		}
	}
	int ans = 0;
	for(int i=1;i<=m;i++)
	if(pre[i]!=-1) ans += mp[pre[i]][i];
	printf("%d %d\n",n-ans%(n+1),ans/(n+1)-w);
}
int main(void){
	while(~scanf("%d%d",&n,&m)){
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++){
				scanf("%d",&mp[i][j]);
				mp[i][j] = mp[i][j]*(n+1);
			}
		int yuan = 0;
		for(int i=1;i<=n;i++){
			int x;scanf("%d",&x);
			yuan += mp[i][x]/(n+1);
			mp[i][x]++;
		}
		km(yuan);
	}
	return 0;
}