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

;
int lx
, ly
, visx
, visy
,linkx
, linky
;
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
;
int ffk
;
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 ;
}