POJ 2400 最小權匹配

吐槽：首先，這道題的輸入居然是錯的。要將上下兩個矩陣的位置換一下才可以出樣例，也就是上面那個矩陣是employee對Supervisor的打分，下面那個矩陣才是Supervisor對employee的打分。 題意：給出兩個矩陣，分別是employee對supervision的打分和supervision對employee的打分。當然矩陣中給出的不是分數，而是進來的先後順序，第一個進來的分數就是1，第二個。。。類推，然後分數越低對這個部門越喜歡，同理下一個矩陣。 然後叫你求出，使得他們都最滿意的方案，並且輸出平均不滿意度，這個平均不滿意度就是，假設a這個人到b這個部分，他的分數是1，但是完備匹配後，他被分配到了c，他對c的分數是2，那麼他的不滿意度就是1，然後平均不滿意度就是不滿意度的總和/(2 * n) 。 思路：建圖，對於每個employe和supervision，連一條邊，權值是他們互相的分數的總和。 然後就是求一次最小權匹配，因為權值越小他們越滿意。最後輸出平均不滿意度的時候，為什麼要 - 2 * n 呢，假設所有的employee和supervision都匹配到了他們權值最小的邊，那麼這條邊的值就是2，所以如果當前是的匹配結果的權值是2 * n 的話，那麼所有人都是在他最滿意的位置上，所以不滿意度就是0. CODE：  

#include <set>  
#include <map>  
#include <stack>  
#include <cmath>  
#include <queue>  
#include <cstdio>  
#include <string>  
#include <vector>  
#include <iomanip>  
#include <cstring>  
#include <iostream>  
#include <algorithm>  
#define Max 2505  
#define FI first  
#define SE second  
#define ll long long  
#define PI acos(-1.0)  
#define inf 0x3fffffff  
#define LL(x) ( x << 1 )  
#define bug puts("here")  
#define PII pair<int,int>  
#define RR(x) ( x << 1 | 1 )  
#define mp(a,b) make_pair(a,b)  
#define mem(a,b) memset(a,b,sizeof(a))  
#define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )  
  
using namespace std;  
  
inline void RD(int &ret) {  
    char c;  
    int flag = 1 ;  
    do {  
        c = getchar();  
        if(c == '-')flag = -1 ;  
    } while(c < '0' || c > '9') ;  
    ret = c - '0';  
    while((c=getchar()) >= '0' && c <= '9')  
        ret = ret * 10 + ( c - '0' );  
    ret *= flag ;  
}  
  
#define N 22  
int n ;  
int Map[N][N] ;  
int lx[N] , ly[N] , visx[N] , visy[N] ,linkx[N] , linky[N] ;  
int find(int now){  
    visx[now] = 1 ;  
    for (int i = 1 ; i <= n ; i ++ ){  
        if(!visy[i]){  
            int fk = Map[now][i] - lx[now] - ly[i] ;  
            if(!fk){  
                visy[i] = 1 ;  
                if(linky[i] == -1 || find(linky[i])){  
                    linky[i] = now ;  
                    linkx[now] = i ;  
                    return 1 ;  
                }  
            }  
        }  
    }  
    return 0 ;  
}  
int KM(){  
    mem(linky ,-1) ;  
    mem(ly ,0) ;  
    for (int i = 1 ; i <= n ; i ++ ){  
        lx[i] = inf ;  
        for (int j = 1 ; j <= n ; j ++ )lx[i] = min(lx[i] , Map[i][j]) ;  
    }  
    for (int i = 1 ; i <= n ; i ++ ){  
        while(1){  
            mem(visx, 0) ;mem(visy, 0) ;  
            if(find(i))break ;  
            int d = inf ;  
            for (int j = 1 ; j <= n ; j ++ )  
                if(visx[j])  
                    for (int k = 1 ; k <= n ; k ++ )  
                        if(!visy[k])  
                            d = min(d , Map[j][k] - lx[j] - ly[k]) ;  
            for (int j = 1 ; j <= n ; j ++ ){  
                if(visx[j])lx[j] += d ;  
                if(visy[j])ly[j] -= d ;  
            }  
        }  
    }  
    int ans = 0 ;  
    for (int i = 1 ; i <= n ; i ++ ){  
        if(linky[i] != -1)ans += Map[linky[i]][i] ;  
    }  
    return ans ;  
}  
int cnt , ans ;  
bool vis[N] ;  
int ffk[N] ;  
void dfs(int dp , int sum){  
    if(sum > ans)return ;  
    if(dp > n){  
        if(sum == ans){  
            printf("Best Pairing %d\n",++ cnt) ;  
            for (int i = 1 ; i <= n ; i ++ ){  
                printf("Supervisor %d with Employee %d\n",i , ffk[i]) ;  
            }  
        }  
        return ;  
    }  
    else  
    for (int i = 1 ; i <= n ; i ++ ){  
        if(!vis[i]){  
            vis[i] = 1 ;  
            ffk[dp] = i ;  
            dfs(dp + 1 , sum + Map[dp][i]) ;  
            vis[i] = 0 ;  
        }  
    }  
}  
int main() {  
    int t ;  
    cin >> t ;  
    int ca = 0 ;  
    while(t -- ){  
        cin >> n ;  
        cnt = 0 ;  
        mem(Map ,0) ;  
        for (int i = 1 ; i <= n ; i ++ ){  
            for (int j = 1 ; j <= n ; j ++ ){  
                int fk ;  
                RD(fk) ;  
                Map[fk][i] += j ;  
            }  
        }  
        for (int i = 1 ; i <= n ; i ++ ){  
            for (int j = 1 ; j <= n ; j ++ ){  
                int fk ;  
                RD(fk) ;  
                Map[i][fk] += j ;  
            }  
        }  
        mem(vis ,0) ;  
        ans = KM() ;  
        printf("Data Set %d, Best average difference: %f\n",++ ca ,( ans - n * 2 ) * 0.5 / n) ;  
        dfs(1 , 0) ;  
        puts("") ;  
    }  
    return 0 ;  
}