#include<iostream>
using namespace std;
const int maxn=100+5;
int w[maxn][maxn],n,m;
int lx[maxn],ly[maxn];
int Left[maxn];
bool s[maxn],t[maxn];
int a[maxn][maxn];
bool match(int i)
{
s[i]=true;
for(int j=1;j<=n;j++) if(lx[i]+ly[j]==w[i][j]&&!t[j]){
t[j]=true;
if(!Left[j]||match(Left[j])){
Left[j]=i;
return true;
}
}
return false;
}
void update()
{
int a=(1<<30);
for(int i=1;i<=n;i++) if(s[i])
for(int j=1;j<=n;j++) if(!t[j])
a=min(a,lx[i]+ly[j]-w[i][j]);
for(int i=1;i<=n;i++){
if(s[i]) lx[i]-=a;
if(t[i]) ly[i]+=a;
}
}
void KM()
{
for(int i=1;i<=n;i++){
Left[i]=lx[i]=ly[i]=0;
for(int j=1;j<=n;j++)
lx[i]=max(lx[i],w[i][j]);
}
for(int i=1;i<=n;i++){
for(;;)
{
for(int j=1;j<=n;j++) s[j]=t[j]=0;
if(match(i)) break;else update();
}
}
}
int main()
{
int T;
scanf("%d",&T);
for(int i=1;i<=T;i++)
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
scanf("%d",&a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
w[i][j]=0;
for(int k=1;k<=m;k++)
w[i][j]-=abs(a[k][i]-j);
}
KM();
int ans=0;
for(int i=1;i<=n;i++) ans+=w[Left[i]][i];
printf("Case #%d: %d\n",i,-ans);
}
return 0;
}
HDU 3523 Image copy detection(KM)
最新推荐文章于 2023-04-08 17:41:34 发布