题目大意: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;
}