//给一个n*m的图,
//m表示人,h表示房子
//问全部人走回家的最小步数
//每一个人仅仅能进一间房
//非常明显的最大带权匹配
//每一个人到每每间房的距离即为权值
//因为是求最小,仅仅要改为负权即可
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<cmath>
using namespace std ;
const int maxn = 110 ;
const int inf = 0x3f3f3f3f ;
int map[maxn][maxn] ;
int lx[maxn] ,ly[maxn] ;
int slack[maxn] , match[maxn] ;
int visx[maxn] , visy[maxn] ;
char str[maxn][maxn] ;
int mx[maxn] , my[maxn] ;
int hx[maxn] , hy[maxn] ;
int n , m ;
bool find(int x)
{
visx[x] = 1 ;
for(int i = 1;i <= n;i++)
{
if(visy[i])continue ;
int tmp = lx[x] + ly[i] - map[x][i] ;
if(tmp == 0)
{
visy[i] = 1 ;
if(match[i] == -1 || find(match[i]))
{
match[i] = x ;
return true ;
}
}
else slack[i] = min(slack[i] , tmp) ;
}
return false ;
}
int KM()
{
memset(match , -1 , sizeof(match)) ;
memset(ly , 0 , sizeof(ly)) ;
for(int i = 1;i <= n;i++)
{
lx[i] = -inf ;
for(int j = 1;j <= n;j++)
lx[i] = max(lx[i] , map[i][j]) ;
}
for(int i = 1;i <= n;i++)
{
for(int j = 1;j <= n;j++)
slack[j] = inf ;
while(1)
{
memset(visx , 0 ,sizeof(visx)) ;
memset(visy , 0 , sizeof(visy)) ;
if(find(i))break ;
int d = inf ;
for(int j = 1;j <= n;j++)
if(!visy[j])
d = min(slack[j] , d) ;
for(int j = 1;j <= n;j++)
if(visx[j])
lx[j] -= d ;
for(int j = 1;j <= n;j++)
if(visy[j])
ly[j] += d ;
else
slack[j] -= d ;
}
}
int ans = 0 ;
for(int j = 1;j <= n;j++)
ans += map[match[j]][j];
return ans ;
}
int main()
{
// freopen("in.txt" , "r" , stdin) ;
while(scanf("%d%d" , &n , &m) && (n+m))
{
int len_m = 0 , len_h = 0 ;
for(int i = 1;i <= n;i++)
{
scanf("%s" , &str[i][1]) ;
for(int j = 1;j <= m;j++)
if(str[i][j] == 'm')
{
mx[++len_m] = i ;
my[len_m] = j ;
}
else if(str[i][j] == 'H')
{
hx[++len_h] = i ;
hy[len_h] = j ;
}
}
memset(map , 0 ,sizeof(map)) ;
for(int i = 1;i <= len_m;i++)
for(int j = 1;j <= len_h;j++)
map[i][j] = (int)(-abs((double)(mx[i]-hx[j])) - abs((double)(my[i]-hy[j]))) ;
n = len_m ;
int ans = KM() ;
cout<<-ans<<endl;
}
}

hdu1533Going Home KM算法的更多相关文章

  1. 匈牙利算法与KM算法

    匈牙利算法 var i,j,k,l,n,m,v,mm,ans:longint; a:..,..]of longint; p,f:..]of longint; function xyl(x,y:long ...

  2. 【HDU2255】奔小康赚大钱-KM算法

    Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Description ...

  3. HDU2255-奔小康赚大钱-二分图最大权值匹配-KM算法

    二分图最大权值匹配问题.用KM算法. 最小权值的时候把权值设置成相反数 /*-------------------------------------------------------------- ...

  4. KM算法及其优化的学习笔记&&bzoj2539: [Ctsc2000]丘比特的烦恼

    感谢  http://www.cnblogs.com/vongang/archive/2012/04/28/2475731.html 这篇blog里提供了3个链接……基本上很明白地把KM算法是啥讲清楚 ...

  5. poj 2195 KM算法

    题目链接:http://poj.org/problem?id=2195 KM算法模板~ 代码如下: #include "stdio.h" #include "string ...

  6. hdu 2255 奔小康赚大钱--KM算法模板

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2255 题意:有N个人跟N个房子,每个人跟房子都有一定的距离,现在要让这N个人全部回到N个房子里面去,要 ...

  7. HDU(2255),KM算法,最大权匹配

    题目链接 奔小康赚大钱 Time Limit: 1000/1000MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Su ...

  8. 二分图 最大权匹配 km算法

    这个算法的本质还是不断的找增广路: KM算法的正确性基于以下定理:若由二分图中所有满足A[i]+B[j]=w[i,j]的边(i,j)构成的子图(称做相等子图)有完备匹配,那么这个完备匹配就是二分图的最 ...

  9. hdu 2255 奔小康赚大钱 KM算法

    看到这么奇葩的题目名我笑了,后来这么一个裸的KM调了2小时我哭了…… 这是个裸的KM算法,也没什么多说的,主要是注意多组数据时,每次都要把各种数组清空啊,赋值啊什么的,反正比较麻烦.至于为什么调了2小 ...

随机推荐

  1. BZOJ 1391 网络流

    vis[0]没有清零查一年- //By SiriusRen #include <cstdio> #include <cstring> #include <algorith ...

  2. Metasploit的攻击实例讲解----辅助扫描工具

    不多说,直接上干货! 怎么弹出来这个呢,连续按两次tab. msf > use auxiliary/scanner/ Display all possibilities? (y or n) us ...

  3. zookeeper伪分布安装配置

    1.下载路径为:http://mirrors.cnnic.cn/apache/zookeeper/stable/ 2.安装: 第一步 解压zookeeper压缩包: 进入 zookeeper安装目录 ...

  4. PHP+MySQL登录注册,完整版,详细注释

    纯手写打造. 下载地址:http://download.csdn.net/detail/qq_33599520/9779970 项目结构: 下面是代码: <!DOCTYPE html> & ...

  5. 《剑指offer》栈的压入、弹出序列

    一.题目描述 输入两个整数序列,第一个序列表示栈的压入顺序,请判断第二个序列是否为该栈的弹出顺序.假设压入栈的所有数字均不相等.例如序列1,2,3,4,5是某栈的压入顺序,序列4,5,3,2,1是该压 ...

  6. Java中将String转json对象

    import org.json.simple.JSONObject; import org.json.simple.parser.JSONParser; import org.json.simple. ...

  7. ubuntu 同时安装anaconda2和anaconda3

    说明:先根据Ubuntu预装的python2.7来安装Anaconda2,然后将Anaconda3作为其环境安装在envs文件夹下. 重要提示:有一些软件需要py2.7的环境,比如XX-Net, 最好 ...

  8. HNU 12933 Random Walks Catalan数 阶乘求逆元新技能

    一个Catalan数的题,打表对每个数都求一次逆元会T,于是问到了一种求阶乘逆元的打表新方法. 比如打一个1~n的阶乘的逆元的表,假如叫inv[n],可以先用费马小定理什么的求出inv[n],再用递推 ...

  9. 思考一下activity的启动模式

    在android里,有4种activity的启动模式.分别为:"standard" (默认) "singleTop" "singleTask" ...

  10. UML之序列图

    一 序列图概述: 序列图主要用于展示对象之间交互的顺序. 序列图将交互关系表示为一个二维图.纵向是时间轴,时间沿竖线向下延伸. 横向轴代表了在协作中各独立对象的类元角色.类元角色用生命线表示.当对象存 ...