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