HDU 1533 Going Home
2015-05-18 11:09
323 查看
HDU 1533 Going Home
题目:HDU 1533
题意:
m代表人,H代表家,将全部人安排最短的路线回家。
思路:
将“人”与每个“家”连一条线,并且以两者的横纵坐标的差值之和的负值为边权。
(即为-(abs(x1-x2)+abs(y1-y2))),因为KM算法是求最大值,所以取负值,输出-ans就行
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=1005;
int g[maxn][maxn];
bool visitx[maxn], visity[maxn];
int lx[maxn], ly[maxn];
int slack[maxn];
int match[maxn];
int n,m;
char s[maxn][maxn];
int nn,mm;
struct jud
{
int x,y;
}H[maxn],M[maxn];
bool Hungary(int u) //匈牙利算法
{
visitx[u] = true;
for(int i = 1; i <= n; ++i)
{
if(visity[i])
continue;
if(lx[u] + ly[i] == g[u][i])
{
visity[i] = true;
if(match[i] == -1 || Hungary(match[i]))
{
match[i] = u;
return true;
}
}
else //不在相等子图
slack[i] = min(slack[i], lx[u] + ly[i] - g[u][i]);
}
return false;
}
void KM()
{
int temp;
memset(match, -1, sizeof(match));
memset(lx, 0, sizeof(lx)); //初始化顶标
memset(ly, 0, sizeof(ly)); //ly[i]为0
for(int i = 1; i <= n; ++i) //lx[i]为权值最大的边
for(int j = 0; j <= m; ++j)
lx[i] = max(lx[i], g[i][j]);
for(int i = 1; i <= n; i++) //对n个点匹配
{
for(int j = 1; j <= n; j++)
slack[j] = INT_MAX;
while(1)
{
memset(visitx, false, sizeof(visitx));
memset(visity, false, sizeof(visity));
if(Hungary(i)) //匹配成功
break;
else //匹配失败,找最小值
{
temp = INT_MAX;
for(int j = 1; j <= n; j++)
if(!visity[j])
if(temp > slack[j])
temp = slack[j];
for(int j = 1; j <= n; j++) //更新顶标
{
if(visitx[j])
lx[j] -= temp;
}
for(int j=1;j<=m;++j)
{
if(visity[j])
ly[j] += temp;
else
slack[j] -= temp;
}
}
}
}
int ans = 0;
for(int i = 1; i <= m; i++) //权值相加
ans += g[match[i]][i];
printf("%d\n", -ans);
}
int main()
{
// freopen("in.txt","r",stdin);
while(scanf("%d %d", &nn,&mm) != EOF&&(nn||mm))
{
n=m=1;
for(int i=0;i<nn;i++)
{
scanf("%s",s[i]);
}
for(int i=0;i<nn;i++)
{
for(int j=0;j<mm;j++)
{
if(s[i][j]=='H')
H
.x=i,H[n++].y=j;
if(s[i][j]=='m')
M[m].x=i,M[m++].y=j;
}
}
n--;m--;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
g[i][j]=-(abs(H[i].x-M[j].x)+abs(H[i].y-M[j].y));
}
}
KM();
}
return 0;
}
题目:HDU 1533
题意:
m代表人,H代表家,将全部人安排最短的路线回家。
思路:
将“人”与每个“家”连一条线,并且以两者的横纵坐标的差值之和的负值为边权。
(即为-(abs(x1-x2)+abs(y1-y2))),因为KM算法是求最大值,所以取负值,输出-ans就行
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=1005;
int g[maxn][maxn];
bool visitx[maxn], visity[maxn];
int lx[maxn], ly[maxn];
int slack[maxn];
int match[maxn];
int n,m;
char s[maxn][maxn];
int nn,mm;
struct jud
{
int x,y;
}H[maxn],M[maxn];
bool Hungary(int u) //匈牙利算法
{
visitx[u] = true;
for(int i = 1; i <= n; ++i)
{
if(visity[i])
continue;
if(lx[u] + ly[i] == g[u][i])
{
visity[i] = true;
if(match[i] == -1 || Hungary(match[i]))
{
match[i] = u;
return true;
}
}
else //不在相等子图
slack[i] = min(slack[i], lx[u] + ly[i] - g[u][i]);
}
return false;
}
void KM()
{
int temp;
memset(match, -1, sizeof(match));
memset(lx, 0, sizeof(lx)); //初始化顶标
memset(ly, 0, sizeof(ly)); //ly[i]为0
for(int i = 1; i <= n; ++i) //lx[i]为权值最大的边
for(int j = 0; j <= m; ++j)
lx[i] = max(lx[i], g[i][j]);
for(int i = 1; i <= n; i++) //对n个点匹配
{
for(int j = 1; j <= n; j++)
slack[j] = INT_MAX;
while(1)
{
memset(visitx, false, sizeof(visitx));
memset(visity, false, sizeof(visity));
if(Hungary(i)) //匹配成功
break;
else //匹配失败,找最小值
{
temp = INT_MAX;
for(int j = 1; j <= n; j++)
if(!visity[j])
if(temp > slack[j])
temp = slack[j];
for(int j = 1; j <= n; j++) //更新顶标
{
if(visitx[j])
lx[j] -= temp;
}
for(int j=1;j<=m;++j)
{
if(visity[j])
ly[j] += temp;
else
slack[j] -= temp;
}
}
}
}
int ans = 0;
for(int i = 1; i <= m; i++) //权值相加
ans += g[match[i]][i];
printf("%d\n", -ans);
}
int main()
{
// freopen("in.txt","r",stdin);
while(scanf("%d %d", &nn,&mm) != EOF&&(nn||mm))
{
n=m=1;
for(int i=0;i<nn;i++)
{
scanf("%s",s[i]);
}
for(int i=0;i<nn;i++)
{
for(int j=0;j<mm;j++)
{
if(s[i][j]=='H')
H
.x=i,H[n++].y=j;
if(s[i][j]=='m')
M[m].x=i,M[m++].y=j;
}
}
n--;m--;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
g[i][j]=-(abs(H[i].x-M[j].x)+abs(H[i].y-M[j].y));
}
}
KM();
}
return 0;
}
相关文章推荐
- hdu 1533 Going Home(最小费用最大流)
- HDU1533-Going Home 简单KM
- hdu 1533 Going Home (最小费用最大流)
- HDU 1533 Going Home 最小费用流
- HDU--1533--Going Home--KM算法
- HDU 1533 Going Home(KM算法模板样例)
- hdu 1533 Going Home 最小费用最大流
- hdu 1533 || poj 2195 Going Home (最小费用最大流)
- HDU-1533 Going Home
- HDU 1533 Going Home【km应用】
- HDU 1533 Going Home【最小费用流|KM算法】
- HDU 1533 Going Home【KM算法】
- POJ 2195 && HDU 1533 Going Home(最小费用最大流-mcmf)
- hdu 1533 Going Home(zkw费用流)
- hdu 1533 Going Home (最大权匹配,KM算法)
- hdu 1533 Going Home(二分KM算法)
- HDU 1533 && POJ 2195 Going Home 最小费用流模板题
- HDU-1533 Going Home(二分图匹配)
- hdu 1533 Going Home(最小权值之和or最小费用流)
- hdu(1533)Going Home(poj2195)