POJ 2195 Going Home【最小费用流 二分图最优匹配】

题目大意：一个n*m的地图，上面有一些人man（m）和数量相等的house(H) 图上的距离为曼哈顿距离 问所有人住进一所房子（当然一个人住一间咯）距离之和最短是多少？

思路：一个人一间房，明显是二分图的模型，边权为人和房子的曼哈顿距离，然后算一下最小距离即可 懒得学KM了 最小费用流的经典建图

 

 

#include <stdio.h>

#include <string.h>

#include <algorithm>

#include <queue>

#include <iostream>

#define maxn 40000

#define inf 0x3f3f3f3f

using namespace std;

struct POINT

{

    int x;

    int y;

}a[maxn],b[maxn];

int head[maxn],root[maxn],point[maxn],next[maxn];

int flow[maxn],cost[maxn],pre[maxn],now=0,dist[maxn];

int n,m,h,h2;

char ch[maxn];

int mabs(int x)

{

    return x>0?x:-x;

}

int distanc(int i,int j)

{

    return mabs(a[i].x-b[j].x)+mabs(a[i].y-b[j].y);

}

void add(int x,int y,int v,int c)

{

    next[++now]=head[x];

    head[x]=now;

    point[now]=y;

    flow[now]=v;

    cost[now]=c;

    root[now]=x;

    next[++now]=head[y];

    head[y]=now;

    point[now]=x;

    flow[now]=0;

    cost[now]=-c;

    root[now]=y;

}

int spfa(int s,int t)

{

    memset(pre,0,sizeof(pre));

    for(int i=0;i<=t;i++)dist[i]=inf;

    dist[s]=0;

    int visit[maxn]={0};

    visit[s]=1;

    queue<int>q;

    q.push(s);

    while(!q.empty())

    {

        int u=q.front();

        q.pop();

        visit[u]=0;

        for(int i=head[u];i;i=next[i])

        {

            int k=point[i];

            if(dist[u]+cost[i]<dist[k] && flow[i]!=0)

            {

                pre[k]=i;

                dist[k]=dist[u]+cost[i];

                if(!visit[k])

                {

                    visit[k]=1;

                    q.push(k);

                }

            }

        }

    }

    return dist[t]!=inf;

}

int main()

{

    while(1)

    {

        scanf("%d%d",&n,&m);

        if(n==0&&m==0)break;

        now=0;

        memset(head,0,sizeof(head));

        int ans=0;

        h=h2=0;

        for(int i=1;i<=n;i++)

        {

            scanf("%s",ch+1);

            for(int j=1;j<=m;j++)

            {

                if(ch[j]=='H')a[++h].x=i,a[h].y=j;

                else if(ch[j]=='m')b[++h2].x=i,b[h2].y=j;

            }

        }

        for(int i=1;i<=h;i++)

        {

            for(int j=1;j<=h;j++)

            {

                int u=distanc(i,j);

                add(j,i+h+1,1,u);

            }

        }

        int s=h*2+10,t=h*2+11;

        for(int i=1;i<=h;i++)add(s,i,1,0);

        for(int i=1;i<=h;i++)add(i+h+1,t,1,0);

        while(spfa(s,t))

        {

            int e=pre[t],minx=flow[e];

            while(e)

            {

                if(flow[e]<minx)minx=flow[e];

                e=pre[root[e]];

            }

            e=pre[t];

            while(e)

            {

                flow[e]-=minx;

                flow[((e-1)^1)+1]+=minx;

                e=[root[e]];

            }

            ans+=dist[t]*minx;

        }

        printf("%d\n",ans);

    }

    return 0;

}