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
#include<iostream>#include<stdio.h>#include<string.h>#include<algorithm>using namespace std;#define INF 99999999#define maxn 305int lx[maxn],ly[maxn];//顶标int Match[maxn];//记录匹配值int visx[maxn],visy[maxn];int w[maxn][maxn];//权值int slack[maxn];//slack为修改量int ans,n;bool findPath(int x)//寻找最优解{    int temp;    visx[x]=1;    for(int y=1; y<=n; y++)    {        if(visy[y])continue;        //temp=w[x][y]-lx[x]-ly[y];//用此TLE        //if(temp==0)        if(w[x][y]==lx[x]+ly[y])//说明是相等子图        {            visy[y]=1;            if(!Match[y]||findPath(Match[y]))            {                Match[y]=x;                return true;            }        }        else        slack[y]=min(slack[y],lx[x]+ly[y]-w[x][y]);    }    return false;}void km(){    memset(Match,0,sizeof(Match));    //初始化顶标    memset(lx,0,sizeof(lx));    memset(ly,0,sizeof(ly));    /*for(int i=1; i<=n; i++)    {        lx[i]=-1;//若求最小权匹配了lx[i]=INF        ly[i]=0;    }*/    for(int i=1; i<=n; i++)        for(int j=1; j<=n; j++)            //if(lx[i]<w[i][j])//若求最小权匹配lx[i]>w[i][j]                lx[i]=max(lx[i],w[i][j]);//    for(int x=1; x<=n; x++)    {        for(int i=1;i<=n;i++)        slack[i]=INF;        while(1)        {            memset(visx,0,sizeof(visx));            memset(visy,0,sizeof(visy));            if(findPath(x))break;            int tmp=INF;            for(int i=1;i<=n;i++)            {                if(!visy[i])                {                    if(tmp>slack[i])                    tmp=slack[i];                }            }            if(tmp==INF)return ;            for(int i=1; i<=n; i++)            {                if(visx[i]) lx[i]-=tmp;//若求最小权匹配改为lx[i]+=slack;                if(visy[i]) ly[i]+=tmp;//ly[i]-=slack;                else                slack[i]-=tmp;            }        }    }}int main(){    while(scanf("%d",&n)!=EOF)    {        ans=0;        for(int i=1; i<=n; i++)            for(int j=1; j<=n; j++)                scanf("%d",&w[i][j]);        km();        for(int i=1; i<=n; i++)            ans+=w[Match[i]][i];        printf("%d\n",ans);    }    return 0;}