HDU 2853 Assignment KM

最新推荐文章于 2021-03-19 16:41:37 发布

ACMmaxx

最新推荐文章于 2021-03-19 16:41:37 发布

阅读量578

点赞数

分类专栏：图论

本文链接：https://blog.csdn.net/acmmaxx/article/details/37316449

版权

图论专栏收录该内容

30 篇文章 0 订阅

订阅专栏

题目大意：n个公司，m个任务（n<=m）。每个公司最多做一个任务，每个任务最多被一个公司做。公司i完成任务j，可以得到价值wij。现在给一种分配方式，要想使得价值和最大，求最优分配方式的价值和与目前的价值的差值和达到最优分配方式，最少改变目前的分配方式的几种关系。

第一个问题是n和m不等，那么加如果虚点，虚点到每个任务的价值都是0。

第二问题是怎么求达到最优分配方式，最少改变目前的分配方式的几种关系？这里有种巧妙的思路，首先，把每个边权乘上k（k>n），最大完备匹配/k就是原图的最大完备匹配。然后，原来的分配方式的边权都+1。这样最大完备匹配/k仍然是原图的最大完备匹配，而且，这样的话，求匹配时候，一定是优先使用原来的匹配边的。

/* hdu 2255 */
//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
//typedef __int64 LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const int INF=100011122;
const double INFF=1e100;
const double eps=1e-8;
const int mod=9999991;
const int NN=305;
const int MM=2000010;
/* ****************** */

int w[NN][NN],lx[NN],ly[NN];
int match[NN],slack[NN],match1[NN];
bool visx[NN],visy[NN];

bool find(int x,int n)
{
    visx[x]=true;
    int i,t;
    for(i=1;i<=n;i++)
    {
        if(visy[i])
            continue;

        t=lx[x]+ly[i]-w[x][i];
        if(t==0)
        {
            visy[i]=true;
            if(match[i]==-1 || find(match[i],n))
            {
                match[i]=x;
                return true;
            }
        }
        else
            slack[i]=min(slack[i],t);
    }
    return false;
}

int KM(int n)
{
    int i,j,t;
    for(i=1;i<=n;i++)
    {
        match[i]=-1;
        ly[i]=0;
        lx[i]=-INF;
        for(j=1;j<=n;j++)
            lx[i]=max(lx[i],w[i][j]);
    }

    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
            slack[j]=INF;

        for(;;)
        {
            for(j=1;j<=n;j++)
                visx[j]=visy[j]=false;

            if( find(i,n) ) break;

            t=INF;
            for(j=1;j<=n;j++)
                if(!visy[j])
                    t=min(t,slack[j]);
            for(j=1;j<=n;j++)
            {
                if(visx[j])
                    lx[j]-=t;
                if(visy[j])
                    ly[j]+=t;
            }
        }
    }

    t=0;
    for(i=1;i<=n;i++)
    {
        t+=w[ match[i] ][i];
    }
    return t;
}

int main()
{
    int i,j,n,m,ans,ans1,t;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
            {
                scanf("%d",&w[i][j]);
                w[i][j]*=55;
            }
        for(i=n+1;i<=m;i++)
            for(j=1;j<=m;j++)
                w[i][j]=0;

        ans1=0;
        memset(match1,-1,sizeof(match1));
        for(i=1;i<=n;i++)
        {
            scanf("%d",&t);
            match1[t]=i;
            ans1+=w[i][t];
            w[i][t]++;
        }

        ans=KM(m);

        t=0;
        for(i=1;i<=m;i++)
        {
            if(match1[i]!=-1 && match1[i]!=match[i])
                t++;
        }

        printf("%d %d\n",t,ans/55-ans1/55);
    }
    return 0;
}