POJ 2195 min_cost_max_flow

给定一个N*M的地图，地图上有若干个man和house，且man与house的数量一致。man每移动一格需花费$1（即单位费用=单位距离），一间house只能入住一个man。现在要求所有的man都入住house，求最小费用。

建图：
s连每个房间 c为1 cost为0
每个人连t c为1 cost为0
每个房间连每个人 c为1 cost为abs（x）+abs（y）

然后跑最小费用流

算法一次敲对了，数组开小了，坑爹啊

#include<cstdio>#include<queue>#include<cstring>#include<algorithm>#include<cmath>using namespace std;const int n_max=205;const int m_max=21505; const int inf=1<<30;struct edge{int u,v,cap,cost,pre;void add(int pu,int pv,int pcap,int pcost,int ppre){u=pu;v=pv;cap=pcap;cost=pcost;pre=ppre;}}e[m_max];int ne;int head[n_max];int min_cost,max_flow;struct point{int x,y;void add(int xx,int yy){x=xx;y=yy;}}hou[n_max],man[n_max];int nhou,nman;void init(){ne=0;memset(head,-1,sizeof(head));}void addedge(int u,int v,int cap,int cost){e[ne].add(u,v,cap,cost,head[u]);head[u]=ne++;e[ne].add(v,u,0,-cost,head[v]);head[v]=ne++;}void min_cost_max_flow(int s,int t){min_cost=max_flow=0;while(true){int p[n_max];memset(p,-1,sizeof(p));int d[n_max];int i;for(i=0;i<n_max;i++){d[i]=inf;}d[s]=0;bool inq[n_max];memset(inq,0,sizeof(inq));queue<int> q;q.push(s);inq[s]=true;int u,v;while(!q.empty()){u=q.front();q.pop();inq[u]=false;for(i=head[u];i!=-1;i=e[i].pre){v=e[i].v;if(e[i].cap&&d[v]>d[u]+e[i].cost){d[v]=d[u]+e[i].cost;p[v]=i;if(!inq[v]){q.push(v);inq[v]=true;}}}}if(d[t]==inf){break;}int nmin=inf;for(i=p[t];i!=-1;i=p[e[i].u]){nmin=min(nmin,e[i].cap);}for(i=p[t];i!=-1;i=p[e[i].u]){e[i].cap -= nmin;e[i^1].cap += nmin;}min_cost += d[t]*nmin;max_flow += nmin;}}int main(){char str[110];int n,m,s,t;while(~scanf("%d%d",&n,&m)){nhou=nman=0;if(n+m==0){break;}int i,j;for(i=0;i<n;i++){scanf("%s",str);for(j=0;j<m;j++){if(str[j]=='H'){hou[++nhou].add(i,j);}else if(str[j]=='m'){man[++nman].add(i,j);}}}s=0;t=nhou+nman+1;init();for(i=1;i<=nhou;i++){for(j=1;j<=nman;j++){int cost=fabs(hou[i].x-man[j].x) + fabs(hou[i].y-man[j].y);addedge(i,nhou+j,1,cost);}}for(i=1;i<=nhou;i++){addedge(s,i,1,0);}for(j=1;j<=nman;j++){addedge(nhou+j,t,1,0);}min_cost_max_flow(s,t);printf("%d\n",min_cost);}return 0;}