HDU 2255 奔小康赚大钱 (KM算法求带权匹配)

xiaoxiao2021-02-28  64

奔小康赚大钱 HDU - 2255 

传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子。  这可是一件大事,关系到人民的住房问题啊。村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住(如果有老百姓没房子住的话,容易引起不安定因素),每家必须分配到一间房子且只能得到一间房子。 

另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.(当然是在他们的经济范围内).现在这个问题就是村领导怎样分配房子才能使收入最大.(村民即使有钱购买一间房子但不一定能买到,要看村领导分配的). 

Input 输入数据包含多组测试用例,每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量),接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。  Output 请对每组数据输出最大的收入值,每组的输出占一行。  Sample Input 2 100 10 15 23 Sample Output 123 中文题面,很容易知道意思。 题目就是让我们求一个二分图的最大匹配。二分图一边是村名,一边是房屋。

这道题目给的完全图,跑一遍KM算法就可以了

#include<iostream> #include<cstdio> #include<cstring> #include<functional> #include<algorithm> #define maxn 305 #define inf 0x3f3f3f3f using namespace std; int n; int link[maxn][maxn]; bool visx[maxn],visy[maxn];//标记访问点集X和点集Y int ex[maxn],ey[maxn];//点集X和点集Y的期望 int match[maxn];//匹配情况 int slack[maxn];//本轮未参与匹配的点集Y需要提高的期望 int dfs(int u)//找增广路过程 { visx[u]=1; for(int v=1;v<=n;v++) { if(link[u][v]&&!visy[v]) { int gap=ex[u]+ey[v]-link[u][v]; if(gap==0)//双方期望和等于边权的时候才匹配 { visy[v]=1; if(match[v]==-1||dfs(match[v])) { match[v]=u; return 1; } } else slack[v]=min(slack[v],gap); } } return 0; } int KM() { int i,j; memset(match,-1,sizeof(match)); memset(ex,0,sizeof(ex)); memset(ey,0,sizeof(ey)); for(i=1;i<=n;i++)//点集X初始化的期望为相关的边中边权最大的一个 { ex[i]=link[i][1]; for(j=2;j<=n;j++) ex[i]=max(ex[i],link[i][j]); } for(i=1;i<=n;i++) { fill(slack,slack+n+1,inf); while(1) { memset(visx,0,sizeof(visx)); memset(visy,0,sizeof(visy)); if(dfs(i))//找到匹配才退出 break; int d=inf;//找不到匹配就改变匹配的点的期望 for(j=1;j<=n;j++) if(!visy[j]) d=min(d,slack[j]); for(j=1;j<=n;j++) { if(visx[j])//本轮参与匹配的点集X降低期望 ex[j]-=d; if(visy[j])//本轮参与匹配的点集Y提高期望 ey[j]+=d; else slack[j]-=d; } } } int sum=0; for(i=1;i<=n;i++)//最大匹配数目肯定为N sum+=link[ match[i] ][i]; return sum; } int main() { int i,j,ans; while(scanf("%d",&n)!=EOF) { for(i=1;i<=n;i++) for(j=1;j<=n;j++) scanf("%d",&link[i][j]); ans=KM(); printf("%d\n",ans); } }

转载请注明原文地址: https://www.6miu.com/read-62831.html

最新回复(0)