//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;
}
hdu 2853Assignment (二分图的最大权匹配)
最新推荐文章于 2020-09-21 22:47:56 发布