poj2195——Going Home

每次选择最小费用可增广路进行增广改进，直到不存在增广路为掌握spfa算法、最小费用最大流、用最大流求二分匹配！理解最小费用最大流原理： 
每次选择最小费用可增广路进行增广改进，直到不存在增广路为止！然后调整流网络！
#include<stdio.h>#include<string.h>#include<math.h>#define max 205#define maxcost 200005int house[max][2];int man[max][2];int cost[max][max];int g[max][max];int flow[max][max];int pre[max],dist[max],max_f[max];//pre[] 记录结点的父节点 ; max_f[] 记录从源点到当前点最多可以增加的容量，即限制值  int d[max];int abs(int x){    if(x>=0)        return x;    else        return -x;}bool spfa(int num,int source,int sink){    int queue[max*max],front=0,rear=0,i,j;    memset(pre,-1,sizeof(pre));    memset(d,0,sizeof(d));    d[source]=1;    for(i=0;i<num;i++)    {        dist[i]=maxcost;    }    dist[source]=0;    queue[rear]=source;rear++;    max_f[source]=maxcost;    while(front!=rear)    {        int temp=queue[front];front++;            for(i=0;i<num;i++)            {                if(((g[temp][i]-flow[temp][i])>0)&&(dist[i]>(dist[temp]+cost[temp][i])))                {                    dist[i]=dist[temp]+cost[temp][i];                    max_f[i]=(max_f[temp]<(g[temp][i]-flow[temp][i])?max_f[temp]:(g[temp][i]-flow[temp][i]));                    pre[i]=temp;                    if(d[i]==0)                    {                        queue[rear]=i;rear++;                        d[i]=1;                    }                }                                }        d[temp]=0;        }    if(pre[sink]!=-1)        return true;    return false;}int sum_cost,sum_flow;void min_max(int num,int source,int sink){    sum_cost=0;sum_flow=0;    memset(flow,0,sizeof(flow));    while(spfa(num,source,sink))    {        int k=sink;        while(pre[k]>=0)        {            flow[pre[k]][k]+=max_f[sink];            flow[k][pre[k]]=-flow[pre[k]][k];//反对称性            k=pre[k];         }        sum_cost+=dist[sink];        sum_flow+=max_f[sink];    }}main(){    int n,m,i,j;    char temp[max];    int house_n,man_n,source,sink;    while(scanf("%d%d",&n,&m)!=EOF)    {        house_n=0;        man_n=0;        if(n==0&&m==0)            break;        for(i=0;i<n;i++)        {            scanf("%s",temp);            for(j=0;j<m;j++)            {                if(temp[j]=='H')                {                    house[house_n][0]=i;                    house[house_n][1]=j;                    house_n++;                }                else if(temp[j]=='m')                {                    man[man_n][0]=i;                    man[man_n][1]=j;                    man_n++;                }            }        }        memset(g,0,sizeof(g));        memset(cost,0,sizeof(cost));        for(i=0;i<man_n;i++)        {            for(j=0;j<house_n;j++)            {                cost[i][j+man_n]=abs(man[i][0]-house[j][0])+abs(man[i][1]-house[j][1]);                cost[j+man_n][i]=-cost[i][j+man_n];//负费用，表示回流会减小费用                g[i][j+man_n]=1;            }            g[man_n+house_n][i]=1;            g[man_n+i][man_n+house_n+1]=1;        }        min_max(man_n+house_n+2,man_n+house_n,man_n+house_n+1);        printf("%d/n",sum_cost);    }}
         g[man_n+house_n][i]=1;//构造超级源点
         g[man_n+i][man_n+house_n]=1;//构造超级汇点