hdu 2853Assignment (二分图的最大权匹配)

//hdu2853
//二分图的最大权匹配，图中有一些边，在等权值的条件下有优先选择权。
//因为是选n个进行匹配，我们将每条边的权值扩大（n+1）倍，这样每条边相对大小是一样的
//对于有优先选择权的边，我们将这些边的权值加一，这样权值大的优先，但是加1的结果不会大于n
//最后除以n就可以把这部分忽略掉。答案对n取模得出的就是我们选择有优先选择权边的数量。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
using namespace std;
#define clr(a) memset(a,0,sizeof(a));
const int N=500;
const int INF=1<<28;
int G[500][500];

struct KM{
    int lx[N],ly[N];//顶标
    int visx[N],visy[N];
    int line[N];
    int we[N];//已配对所对应边权值
    int nx,ny;
    int slack;
    bool find(int u){
        visx[u]=true;
        for(int v=1;v<=ny;v++){
            if(!visy[v])//v不在相等子图中
            {
                int t=lx[u]+ly[v]-G[u][v];
                if(t==0){
                    visy[v]=1;
                    if(line[v]==-1||find(line[v])){
                        line[v]=u;
                        we[v]=G[u][v];
                        return 1;
                    }
                }
                else if(t<slack)
                    slack=t;
            }
        }
        return 0;
    }
    bool km(){//返回是否完备匹配
        clr(lx);clr(ly);
        //memset(lx, -INF, sizeof(lx));
        //memset(ly, -INF, sizeof(ly));
        memset(line,-1,sizeof(line));
        for(int i=1;i<=nx;i++){
            for(int j=1;j<=ny;j++){
                if(G[i][j]>lx[i])
                    lx[i]=G[i][j];//初始化为xi所连接的最大权值
            }
        }
        for(int u=1;u<=nx;u++)
        {
            while(true){
                slack=INF;
                clr(visx);clr(visy);
                if(find(u))break;
                if(slack==INF)return 0;
                for(int i=1;i<=nx;i++){
                    if(visx[i])lx[i]-=slack;
                }
                for(int i=1;i<=ny;i++){
                    if(visy[i])ly[i]+=slack;
                }
            }
        }
        return 1;
    }
    int best_match(){
        int res=0;
        km();
        for(int u=1;u<=ny;u++){
            if(line[u]==-1)continue;
            res+=we[u];
        }
        return res;
    }
}kk;

int main(){
    int n,m,u;
    while(scanf("%d%d",&n,&m)!=EOF){
        clr(G);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&u);
                G[i][j]=u*(n+1);
            }
        }
        int sum=0;
        for(int i=1;i<=n;i++){
            scanf("%d",&u);
            sum+=G[i][u]/(n+1);
            G[i][u]++;
        }
        kk.nx=n,kk.ny=m;
        int ans=kk.best_match();
        printf("%d %d\n",n-ans%(n+1),(ans/(n+1))-sum);
      
    }
    return 0;
}