HDU 2853 && HDU 3315

http://acm.hdu.edu.cn/showproblem.php?pid=2853

题意：给一个n-m二分图，边权用一个n*m的矩阵表示，给出初始匹配，求二分图完美匹配相比初始匹配改变了几条边以及改变的数值

这类题的主要思想是增加原配边的权值，但又不影响最后结果。

步骤1：观察顶点数，每条边乘一个大于顶点数的数v

步骤2：对于原配边，每边加1（注意步骤2可以保证km()/v的结果与原结果相同）

步骤3：求完美匹配，答案为res，改变的边数=n-res%v（res%v表示完美匹配中有多少边属于原匹配），要求的完美匹配答案=res/v

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N=;
const int INF=0x3f3f3f3f;
int nx,ny;
int linker[N],lx[N],ly[N],slack[N];
int visx[N],visy[N],w[N][N];
int DFS(int x)
{
    visx[x]=;
    for(int y=;y<=ny;y++){
        if(visy[y])
            continue;
        int tmp=lx[x]+ly[y]-w[x][y];
        if(tmp==){
            visy[y]=;
            if(linker[y]==- || DFS(linker[y])){
                linker[y]=x;
                return ;
            }
        }else if(slack[y]>tmp){
            slack[y]=tmp;
        }
    }
    return ;
}
int change ;
int KM()
{
    int i,j;
    memset(linker,-,sizeof(linker));
    memset(ly,,sizeof(ly));
    for(i=;i<=nx;i++)
        for(j=,lx[i]=-INF;j<=ny;j++)
            if(w[i][j]>lx[i])
                lx[i]=w[i][j];
    for(int x=;x<=nx;x++){
        for(i=;i<=ny;i++)
            slack[i]=INF;
        while(){
            memset(visx,,sizeof(visx));
            memset(visy,,sizeof(visy));
            if(DFS(x))
                break;
            int d=INF;
            for(i=;i<=ny;i++)
                if(!visy[i] && d>slack[i])
                    d=slack[i];
            for(i=;i<=nx;i++)
                if(visx[i])
                    lx[i]-=d;
            for(i=;i<=ny;i++)
                if(visy[i])
                    ly[i]+=d;
                else
                    slack[i]-=d;
        }
    }
    int res=;
    for(i=;i<=ny;i++)
        if(linker[i]!=-)
            res+=w[linker[i]][i];
    change=nx-res% ;
    return res/;
}
int main()
{
    while(~scanf("%d%d",&nx,&ny))
    {
        for(int i=;i<=nx;i++)
            for(int j=;j<=ny;j++)
            {
                scanf("%d",&w[i][j]);
                w[i][j]*= ;
            }
        int res= ;
        for(int i= ;i<=nx ;i++)
        {
            int x ;
            scanf("%d",&x) ;
            res+=w[i][x] ;
            w[i][x]++ ;
        }
        int ans=KM();
        printf("%d %d\n",change,ans-res/);
    }
    return ;
}

http://acm.hdu.edu.cn/showproblem.php?pid=3315

与2853一样思路的题目，区别是这题图没有直接给，要先算一下

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N=;
const int INF=0x3f3f3f3f;
int n,nx,ny;
int linker[N],lx[N],ly[N],slack[N];
int visx[N],visy[N],w[N][N];
int DFS(int x)
{
    visx[x]=;
    for(int y=;y<=ny;y++){
        if(visy[y])
            continue;
        int tmp=lx[x]+ly[y]-w[x][y];
        if(tmp==){
            visy[y]=;
            if(linker[y]==- || DFS(linker[y])){
                linker[y]=x;
                return ;
            }
        }else if(slack[y]>tmp){
            slack[y]=tmp;
        }
    }
    return ;
}
int change ;
int KM()
{
    int i,j;
    memset(linker,-,sizeof(linker));
    memset(ly,,sizeof(ly));
    for(i=;i<=nx;i++)
        for(j=,lx[i]=-INF;j<=ny;j++)
            if(w[i][j]>lx[i])
                lx[i]=w[i][j];
    for(int x=;x<=nx;x++){
        for(i=;i<=ny;i++)
            slack[i]=INF;
        while(){
            memset(visx,,sizeof(visx));
            memset(visy,,sizeof(visy));
            if(DFS(x))
                break;
            int d=INF;
            for(i=;i<=ny;i++)
                if(!visy[i] && d>slack[i])
                    d=slack[i];
            for(i=;i<=nx;i++)
                if(visx[i])
                    lx[i]-=d;
            for(i=;i<=ny;i++)
                if(visy[i])
                    ly[i]+=d;
                else
                    slack[i]-=d;
        }
    }
    int res=;
    for(i=;i<=ny;i++)
        if(linker[i]!=-)
            res+=w[linker[i]][i];
    if(res<=)return - ;
    change=nx-res% ;
    return res/;
}
int V[N],H1[N],A1[N],H2[N],A2[N] ;
int main()
{
    while(~scanf("%d",&n),n)
    {
        nx=ny=n ;
        for(int i= ;i<=n ;i++)
            scanf("%d",&V[i]) ;
           for(int i= ;i<=n ;i++)
               scanf("%d",&H1[i]) ;
        for(int i= ;i<=n ;i++)
               scanf("%d",&H2[i]) ;
        for(int i= ;i<=n ;i++)
               scanf("%d",&A1[i]) ;
        for(int i= ;i<=n ;i++)
               scanf("%d",&A2[i]) ;
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
            {
                int x ;
                x=H1[i]/A2[j] ;
                if(H1[i]%A2[j])x++ ;
                if(x*A1[i]>=H2[j])w[i][j]=V[i]* ;
                else w[i][j]=(-V[i]*) ;
            }
        for(int i= ;i<=n ;i++)
            w[i][i]++ ;
        int ans=KM();
        if(ans==-)puts("Oh, I lose my dear seaco!") ;
        else printf("%d %.3lf%%\n",ans,(nx-change)*100.0/n);
    }
    return ;
}