Golden Tiger Claw UVA - 11383（km原理）

这题使我对km多了一些看法

写给自己看。。

km结束后bx[i] + by[j] == w[i][j], 所以所有bx与by的和即为w的和

而且记住bx[i] + by[j] >= w[i][j] 这个式子 那么bx与by的最小值 即为km结束后的值

#include <iostream>

#include <cstdio>

#include <cstring>

#include <iostream>

#include <queue>

#include <algorithm>

#include <vector>

#define mem(a, b) memset(a, b, sizeof(a))

using namespace std;

const int maxn = , INF = 0x7fffffff;

int usedx[maxn], usedy[maxn], w[maxn][maxn], bx[maxn], by[maxn], cx[maxn], cy[maxn];

int nx, ny, n, minn, min_value;

bool dfs(int u)

{

    usedx[u] = ;

    for(int i=; i<=ny; i++)

    {

        if(usedy[i] == -)

        {

            int t = bx[u] + by[i] - w[u][i];

            if(t == )

            {

                usedy[i] = ;

                if(cy[i] == - || dfs(cy[i]))

                {

                    cy[i] = u;

                    cx[u] = i;

                    return ;

                }

            }

            else if(t > )

                minn = min(minn, t);

        }

    }

    return ;

}

int km()

{

    mem(cx, -);

    mem(cy, -);

    mem(bx, -);

    mem(by, );

    for(int i=; i<=nx; i++)

        for(int j=; j<=ny; j++)

            bx[i] = max(bx[i], w[i][j]);

    for(int i=; i<=nx; i++)

    {

        while()

        {

            minn = INF;

            mem(usedx, -);

            mem(usedy, -);

            if(dfs(i)) break;

            for(int j=; j<=nx; j++)

                if(usedx[j] != -) bx[j] -= minn;

            for(int j=; j<=ny; j++)

                if(usedy[j] != -) by[j] += minn;

        }

    }

    min_value = ;

    for(int i=; i<=n; i++)

        if(cx[i] != -)

            min_value += w[i][cx[i]];

    return min_value;

}

int main()

{

    while(~scanf("%d",&n))

    {

        for(int i=; i<=n; i++)

            for(int j=; j<=n; j++)

                scanf("%d",&w[i][j]);

        nx = ny = n;

        int ans = km();

    for(int i=; i<=nx; i++)

    {

        if(i != ) printf(" ");

        printf("%d", bx[i]);

    }

    printf("\n");

    for(int i=; i<=ny; i++)

    {

        if(i != ) printf(" ");

        printf("%d", by[i]);

    }

    printf("\n");

    cout<< ans <<endl;

    }

    return ;

}