HDU-2255-KM算法

                           奔小康赚大钱

Problem Description

传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革：重新分配房子。
这可是一件大事,关系到人民的住房问题啊。村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住（如果有老百姓没房子住的话，容易引起不安定因素），每家必须分配到一间房子且只能得到一间房子。
另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.(当然是在他们的经济范围内).现在这个问题就是村领导怎样分配房子才能使收入最大.(村民即使有钱购买一间房子但不一定能买到,要看村领导分配的).

 

Input

输入数据包含多组测试用例，每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量)，接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。

 

Output

请对每组数据输出最大的收入值，每组的输出占一行。

 

Sample Input

2100 1015 23

 

Sample Output

123

一些概念：

设G=[V,E]是一个无向图，M属于E是G的若干条边的集合，如果M中的任意两条边都没有公共端点，就称M是一个匹配。

设P是图G的一条路径，如果该路径的任意两条相邻边一定是一条属于M而另一条不属于M，就称P是一条交错轨。

交错轨中有两个端点都是未盖点（不在M中构成边的点）的交错轨叫可增广轨。

若两个未盖点之间只有一条边则该边为可增广轨。因为一条单边为交错轨，且两端点均为未盖点。

一条关于M的增广路径的长度必为奇数，且路上的第一条边和最后一条边都不属于M。

若匹配M中不存在可增广轨，M一定是最大匹配。

最大匹配方案可能有很多，但最大匹配数是一样的。

如果一个二分图，X部和Y部的顶点数相等，若存在一个匹配包含X部与Y部的所有顶点，则称为完美匹配。

先服用

https://www.cnblogs.com/wenruo/p/5264235.html#3479945

再服用

http://blog.csdn.net/c20180630/article/details/71080521

https://www.cnblogs.com/chenhuan001/archive/2013/05/05/3060449.html

效果更好

开始的时候用C++一直超时，浪费那么多时间，，，，，

#include<stdio.h>
#include<string.h>
#define INF 0x3f3f3f3f
int G[305][305],match[305],n;
int visx[305],visy[305];
int dbx[305],dby[305],slack[305];
int dfs(int x)
{
	visx[x]=1;
	for(int y=0;y<n;y++)
	{
		if(visy[y]) continue;
		int c=dbx[x]+dby[y]-G[x][y];
		if(c==0)
		{
			visy[y]=1;
			if(match[y]==-1||dfs(match[y]))
			{
				match[y]=x;
				return 1;
			}
		}
		else if(c<slack[y])
			slack[y]=c;
	}
	return 0;
}
int KM()
{
	memset(match,-1,sizeof(match));
	memset(dbx,0,sizeof(dbx));
	memset(dby,0,sizeof(dby));
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			if(G[i][j]>dbx[i])
				dbx[i]=G[i][j];
	for(int x=0;x<n;x++)
	{
		memset(slack,INF,sizeof(slack));
		while(1)
		{
			memset(visx,0,sizeof(visx));
			memset(visy,0,sizeof(visy));
			if(dfs(x))break;
			int d=INF;
			for(int i=0;i<n;i++)
				if(!visy[i]&&slack[i]<d)
					d=slack[i];
			for(int i=0;i<n;i++)
			{
				if(visx[i]) dbx[i]-=d;
				if(visy[i]) dby[i]+=d;
				else slack[i]-=d;
			}
		}
	}
	int ans=0;
	for(int i=0;i<n;i++)
		if(match[i]>-1)
		ans+=G[match[i]][i];
	return ans;
}
int main()
{
	while(scanf("%d",&n)!=EOF)
	{
		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		scanf("%d",&G[i][j]);
		printf("%d\n",KM());
	}
}