poj 2195 最小费用最大流 EK+SPFA

传送门

题意：给你一个矩阵，里面有数量相等的H和m，代表房子和人，每个人要找个房子，一个人到一个房子的费用是他走的步数，一个房子只能容纳一个人，求所有人都到房子的最小费用。

个人心得：第一次做最小费用最大流，原本不会费用流的，借着这题就是想学习一下，网上各种博客讲的都好神啊。由于本人最大流只会EK算法（最慢的），所以对于网上大牛的博客各种看不懂，直到看到了这个，里面说就是把EK算法的BFS过程改成SPFA，于是顿悟，好吧，就是这点差别。我理解了下她的代码，然后完成了人生第一道费用流。不过目前掌握的算法效率不高，有待学习新的算法。

思路：感觉没啥好说的，主要就是建图，添加超级源点和汇点，然后就是按照EK+SPFA的思路做了。

#include<iostream>#include<cstdio>#include<cstring>#include<queue>#include<cmath>#define maxn 2000using namespace std;int m,n;int hnum,mnum,hx[200],hy[200],mx[200],my[200];int c[210][210],cost[210][210],f[210][210],d[210],pre[210];bool v[210];char str[200];int abs(int x){    return x>0?x:-x;}int EK(int s,int t){    int ans=0;    memset(f,0,sizeof(f));    memset(v,0,sizeof(v));    queue<int>q;    while(1)    {        for(int i=0;i<=t;i++)d[i]=maxn;        d[s]=0;        q.push(s);        v[s]=1;        while(!q.empty())        {            int a=q.front();            q.pop();            v[a]=0;            for(int i=s;i<=t;i++)            {                if(c[a][i]-f[a][i]>0&&d[i]>d[a]+cost[a][i])                {                    pre[i]=a;                    d[i]=d[a]+cost[a][i];                    if(!v[i])                    {                        q.push(i);                        v[i]=1;                    }                }            }        }        if(d[t]==maxn)break;        int mi=maxn;        for(int i=t;i!=s;i=pre[i])        {            if(c[pre[i]][i]-f[pre[i]][i]<mi)mi=c[pre[i]][i]-f[pre[i]][i];        }        for(int i=t;i!=s;i=pre[i])        {            f[pre[i]][i]+=mi;            f[i][pre[i]]-=mi;            ans+=cost[pre[i]][i];        }    }    return ans;}int main(){    while(scanf("%d%d",&m,&n)&&n&&m)    {        hnum=mnum=0;        memset(c,0,sizeof(c));        memset(cost,0,sizeof(cost));        for(int i=0;i<m;i++)        {            scanf("%s",str);            for(int j=0;j<n;j++)            {                if(str[j]=='H')                {                    hx[++hnum]=i;                    hy[hnum]=j;                }                else if(str[j]=='m')                {                    mx[++mnum]=i;                    my[mnum]=j;                }            }        }        for(int i=1;i<=mnum;i++)        {            c[0][i]=1;            for(int j=1;j<=hnum;j++)            {                int hh=mnum+j;                c[i][hh]=1;                cost[i][hh]=abs(mx[i]-hx[j])+abs(my[i]-hy[j]);                cost[hh][i]=-cost[i][hh];            }        }        for(int i=1;i<=hnum;i++)        {            c[i+mnum][hnum+mnum+1]=1;        }        printf("%d\n",EK(0,hnum+mnum+1));    }    return 0;}