#include<iostream>
#include<vector>
using namespace std;
const int maxn=100+5;
struct point{int x,y;point(int x=0,int y=0):x(x),y(y) {}};
vector<point> st,en;
char map[maxn][maxn];
int w[maxn][maxn],n;
int lx[maxn],ly[maxn];
int Left[maxn];
bool s[maxn],t[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 r,c;
while(scanf("%d%d",&r,&c)&&(r+c))
{
st.clear();en.clear();
for(int i=1;i<=r;i++) scanf("%s",map[i]+1);
for(int i=1;i<=r;i++){
for(int j=1;j<=c;j++){
if(map[i][j]=='m') st.push_back(point(i,j));
if(map[i][j]=='H') en.push_back(point(i,j));
}
}
for(int i=0;i<st.size();i++)
{
for(int j=0;j<en.size();j++)
{
w[j+1][i+1]=-(abs(en[j].x-st[i].x)+abs(en[j].y-st[i].y));
}
}
n=st.size();
KM();
int ans=0;
for(int i=1;i<=n;i++) ans+=-w[Left[i]][i];
printf("%d\n",ans);
}
return 0;
}
HDU 1533 Going Home(最佳二分匹配)
最新推荐文章于 2017-01-09 20:09:38 发布