poj 2195 （最小费用最大流）

恩，不明白C++为何编译不过。。只能交G++

初步理解最小费用最大流。就是利用费用求出最短路，然后沿着最短路，找出最大流，再更改网络流，算出最小费用。

#include<stdio.h>
#include<string.h>
#define MX 100000
int pos_m,pos_h,k,n,m,vs,ans,vt,dis[10010],que[10000],vis[10000],list[10000],pre[1000];
struct Egde{
	int u,v,pre,next,c,f,w;
}eg[100000];
struct Man{
	int x,y;
}man[1000];
struct House{
	int x,y;
}house[1000];
int min(int a,int b){
	return a<b?a:b;
}
int abs(int a){
	if(a>0)
		return a;
	return -a;
}
void add(int u,int v,int c,int w){  //建图
	eg[k].u=u;eg[k].v=v;eg[k].pre=k+1;eg[k].next=list[u];eg[k].c=c;eg[k].f=0;eg[k].w=w;list[u]=k++;
	eg[k].u=v;eg[k].v=u;eg[k].pre=k-1;eg[k].next=list[v];eg[k].c=0;eg[k].f=0;eg[k].w=-w;list[v]=k++;
}
bool spfa(){    //最小费用的 最短路
	int qs=0,qe=1;
	for(int i=0;i<=2*pos_m+1;i++){
		dis[i]=MX;
		vis[i]=false;
	}
	dis[0]=0;
	que[0]=0;
	vis[0]=true;
	while(qs<qe){
		int u=que[qs++];
		for(int i=list[u];i!=-1;i=eg[i].next){
			int v=eg[i].v;
			if(eg[i].c-eg[i].f&&dis[v]>dis[u]+eg[i].w){
				dis[v]=dis[u]+eg[i].w;
				pre[v]=i;
				if(!vis[v]){
					vis[v]=true;
					que[qe++]=v;
				}
			}
		}
		vis[u]=false;
	}
	if(dis[2*pos_m+1]==MX){
		return false;
	}
	return true;
}
void find(){  //在最短路上找最大流
	int sum=MX,p;
	for(int u=2*pos_m+1;u!=0;u=eg[eg[p].pre].v){  //找寻最大流
		p=pre[u];
		sum=min(sum,eg[p].c-eg[p].f);
	}
	for(int u=2*pos_m+1;u!=0;u=eg[eg[p].pre].v){  //更改网络流量，并计算最小费用于 ans
		p=pre[u];
		eg[p].f+=sum;
		eg[eg[p].pre].f-=sum;
		ans+=sum*eg[p].w;
	}
}
int main(){
	while(scanf("%d%d",&n,&m)&&n&&m){
		pos_h=0,pos_m=0;
		k=0;
		for(int i=1;i<=n;i++){
			char s[110];
			scanf("%s",s+1);
			for(int j=1;j<=m;j++){
				if(s[j]=='H'){
					house[pos_h].x=i;
					house[pos_h++].y=j;
				}
				if(s[j]=='m'){
					man[pos_m].x=i;
					man[pos_m++].y=j;
				}
			}
		}
		memset(list,-1,sizeof(list));
		for(int i=0;i<pos_m;i++){
			for(int j=0;j<pos_h;j++){
				add(i+1,pos_m+j+1,1,abs(man[i].x-house[j].x)+abs(man[i].y-house[j].y));
			}
		}
		for(int i=1;i<=pos_m;i++){
			add(0,i,1,0);
			add(i+pos_m,2*pos_m+1,1,0);
		}
		ans=0;
		while(spfa()){
			find();
		}
		printf("%d\n",ans);
	}
	return 0;
}