HDU2255奔小康赚大钱-完全带权二分图最佳匹配

这题一打眼是最大费用最大流，但是边太多（完全二分图），会被卡掉。

正解：KM算法

一个不错的博客，但是有些东西还是不是很明白。

http://blog.sina.com.cn/s/blog_691ce2b701016reh.html

关于他写的最后一句，每次要更新stack，实测本题不需要，可能别的地方会用到，以后再讨论。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define inf 0x3f3f3f3f
using namespace std;
int dis[303][303];
int lx[333],ly[333];
int mind[333];
int visx[333];
int visy[333];
int match[333];
int N;
bool dfs(int u)
{
    visx[u] = 1;
    for (int i = 1; i<= N; i++)
    {
        if (visy[i]) continue;
        int t = lx[u] + ly[i] - dis[u][i];
        if (t == 0)
        {
            visy[i] = 1;
            if (match[i]==-1 || dfs(match[i]))
            {
                match[i] = u;
                return true;
            }
        }
        else if (mind[i] > t)
        {
            mind[i] = t;
        }
    }
    return false;
}
int KM()
{
    memset(ly,0,sizeof(lx));
    memset(match,-1,sizeof(match));
    //枚举一方面的点
    for (int i = 1 ;i<=N; i++)
    {   //储存最小的mind
        memset(mind,inf,sizeof(mind));
        memset(visx,0,sizeof(visx));
        memset(visy,0,sizeof(visy));
        while (!dfs(i)) //这里需要多次dfs同一个点，因为更改lx和ly会使得需要再次dfs
        {
            int d = inf;
            for (int j = 1 ; j<=N; j++)
            {
                if (!visy[j]) d = min(d,mind[j]);
            }
            for (int j = 1; j<=N; j++)
            {
                if (visx[j]) lx[j]-=d;
                if (visy[j]) ly[j]+=d;
            }
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
        }
    }
    int ans = 0;
    for (int i= 1;i<= N;i++)
            ans += dis[match[i]][i];
    return ans;
}
int main()
{
    while (scanf ("%d",&N)!=EOF)
    {
        for (int i = 1; i<=N;i++)
        lx[i] = -inf;
        for (int i = 1 ; i<= N; i++)
        {
            for (int j = 1; j <= N; j++)
            {
                scanf("%d",&dis[i][j]);
                lx[i] = max(lx[i],dis[i][j]);
            }
        }
        printf ("%d
",KM());
    }
}