hdu 2853

虚拟赛一开始lyf就对我说这是一道匹配的题目，我一看明显裸的最优匹配，敲完提交wrong，

题目要求改变尽量少的公司，就是如果遇到相等的权值，优先选择跟他原来匹配的，KM匹配是按序号大小来的，如果一个公司原来匹配的序号较大，前面有权值相等的点时，KM就会选择前面的点参加匹配。想了好长时间不知道怎么去优先选择原来匹配的边，

最后想着如果把原来匹配的边变得大一些的话，就可以，但是变大的话就会影响最优匹配的总值，而且变大的话还会影响原来比他大的权值，所以就是所有的权值都得扩大，我想到的是都*100，原来匹配的边再加1，因为最多选50条边，也就是最多有50个01相加，不会超过一百，得到的答案除以一百，就把加的1都去掉了，只要扩大的倍数大于n就可以，，

 

#include<stdio.h>

#include<string.h>

#define N 55

#define inf 0x3fffffff

int map[N][N],lx[N],ly[N],sx[N],sy[N],d[N],n,m,match[N],Kmatch[N],tch[N];

int find(int x)

{

    sx[x]=1;

    for(int i=1;i<=m;i++)

    {

        if(sy[i]==1)continue;

        int temp=lx[x]+ly[i]-map[x][i];

        if(temp==0)

        {

            sy[i]=1;

            if(match[i]==-1||find(match[i])==1)

            {

                match[i]=x;

                Kmatch[x]=i;

                return 1;

            }

        }

        else d[i]=d[i]>temp?temp:d[i];

    }

    return 0;

}

int KM()

{

    int i,j,k,min,sum;

    memset(ly,0,sizeof(ly));

    memset(match,-1,sizeof(match));

    memset(Kmatch,-1,sizeof(Kmatch));

    for(i=1;i<=n;i++)

    {

        lx[i]=map[i][1];

        for(j=2;j<=m;j++)

            if(lx[i]<map[i][j])

                lx[i]=map[i][j];

    }

    for(i=1;i<=n;i++)

    {

        for(j=1;j<=m;j++)

            d[j]=inf;

        while(1)

        {

            memset(sx,0,sizeof(sx));

            memset(sy,0,sizeof(sy));

            if(find(i)==1)break;

            min=inf;

            for(k=1;k<=m;k++)

                if(sy[k]==0&&min>d[k])

                    min=d[k];

            for(j=1;j<=n;j++)

                if(sx[j]==1)lx[j]-=min;

            for(j=1;j<=m;j++)

                if(sy[j]==1)ly[j]+=min;

        }

    }

    sum=0;

    for(i=1;i<=n;i++)

        sum+=map[i][Kmatch[i]];

    return sum;

}

int main()

{

    int i,j,k,ans,sum;

    while(scanf("%d%d",&n,&m)!=-1)

    {

        for(i=1;i<=n;i++)

            for(j=1;j<=m;j++)

                map[i][j]=-inf;

        for(i=1;i<=n;i++)

            for(j=1;j<=m;j++)

            {

                scanf("%d",&map[i][j]);

                map[i][j]*=100;

            }

        ans=0;

        for(i=1;i<=n;i++)

        {

            scanf("%d",&k);

            tch[i]=k;

            ans+=map[i][k];

            map[i][k]++;

        }

        ans/=100;

        sum=KM()/100;

        k=0;

        for(i=1;i<=n;i++)

            if(tch[i]!=Kmatch[i])

                k++;

        printf("%d %d\n",k,sum-ans);

    }

    return 0;

}