hdu - 2853 二分图最优匹配 优先使用原匹配边（☆）

题意：给定一个二分图，N个点对应M个点，两两之间存在一组关系，每组关系一个权值。题目中了给定了一个匹配方案，现在要求满足这组关系中的最大的匹配权值在原方案上增长了多少？并且还要求出在原匹配方案上改变(最少)多少条边才能够得到这个最大匹配？

        左右不对称最优匹配问题 N <= M ，直接用KM算法可求得最优匹配，问题就在于有原始匹配边，尽量使用原有匹配

方法：由于最优匹配最终必有N条边， 所以将原图所有的边的权值乘以（N+1）扩大（N + 1）倍，并且如果边是原匹配中的一条边就再 加 1     所以边权就变为了 （N+1）的倍数或（N+1）的倍数 + 1      

所以最终得到的最优权值     除以（N+1）就是最优权值解   、 %（N+1）就是复用旧边的条数

老边加1保证了在最优解唯一时不丢失最优解，在最优解不唯一时保证了优先用原匹配边

这样就实现尽量使用原匹配边的情况下求最优匹配

链接：hdu 2853

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <set>
#include <map>
#include <stack>
#include <queue>
#define INF 0x3f3f3f3f

using namespace std;

const int inf = 0x3f3f3f3f;
const int maxn = 1005;

int n, nx, ny, w[maxn][maxn];
int link[maxn], lx[maxn], ly[maxn], slack[maxn];
bool visx[maxn], visy[maxn];

bool dfs(int x){
    visx[x] = true;
    for(int y = 1; y <= ny; y++){
        if(visy[y]) continue;
        int t = lx[x] + ly[y] - w[x][y];
        if(t == 0){
            visy[y] = true;
            if(link[y] == -1 || dfs(link[y])){
                link[y] = x;
                return true;
            }
        }
        else
            if(slack[y] > t){ //不在相等子图中slack取最小的
                slack[y] = t;
            }
    }
    return false;
}

int KM(){
    //nx = ny = n;
    memset(link, -1, sizeof(link));
    memset(ly, 0, sizeof(ly));
    for(int i = 1; i <= nx; i++){ //lx 初始化为与它关联边中最大的
        lx[i] = -INF;
        for(int j = 1; j <= ny; j++){
            if(w[i][j] > lx[i]) lx[i] = w[i][j];
        }
    }

    for(int x = 1; x <= nx; x++){
        for(int i = 1; i <= ny; i++){
            slack[i] = INF;
        }
        while(true){
            memset(visx, false, sizeof(visx));
            memset(visy, false, sizeof(visy));

            //若成功(找到了增广路)，则该点增广完毕，下一个点进入增广
            if(dfs(x)) break;

            //若失败，则需要改变一些点的顶标，使得图中可行边的数量增加
            //(1)将所有在增广轨中的 X 方点的标号 全部减去一个常数 d ；
            //(2)将所有在增广轨中的 Y 方点的标号 全部加上一个常数 d ；
            int d = INF;
            for(int i = 1; i <= ny; i++){
                if(!visy[i] && d > slack[i]) d = slack[i];
            }
            for(int i = 1; i <= nx; i++){
                if(visx[i]) lx[i] -= d;
            }
            for(int i = 1; i <= ny; i++){
                if(visy[i]) ly[i] += d;
                else slack[i] -= d;
            }
        }
    }
    int res = 0;
    for(int i = 1; i <= ny; i++){
        if(link[i] > -1) res += w[link[i]][i];
    }
    return res;
}

int main ()
{
    int t, kcase = 0;
    int N, M, E;
    while(cin >> N >> M) {
        for(int i = 0; i < maxn; i++) {
            for(int j = 0; j < maxn; j++) {
                w[i][j] = -INF;
            }
        }
        for(int i = 1; i <= N; i++) {
            for(int j = 1; j <= M; j++) {
                scanf("%d", &E);
                w[i][j] = E * (N + 2);
            }
        }
        int sum = 0;
        for(int i = 1; i <= N; i++) {
            scanf("%d", &E);
            sum += w[i][E];
            w[i][E]++;
        }
        sum /= (N + 2);
        n = N;
        nx = N;
        ny = M;
        int mm = KM();
        printf("%d %d\n", n - mm % (N + 2), mm / (N + 2) - sum);
        /*int ans = 0, flag = 0;
        for(int i = 1; i <= n; i++) {
            if(w[link[i]][i] == -INF) {
                flag = 1;
                break;
            }
            ans += w[link[i]][i];
        }
        printf("Case %d: ", ++kcase);
        if(flag){
            puts("-1");
        }
        else {
            printf("%d\n", mm);
        }*/
    }
    return 0;
}