poj-2195-Going Home最小费用最大流

重新切一遍最小费用最大流~~~

这到题目的数据范围有问题，尽量开大就好了~~

 

#include<stdio.h>

#include<iostream>

#include<algorithm>

#include<string.h>

#include<queue>

using namespace std;

#define INF 999999

struct list

{

	int u;

	int v;

	int next;

}node[1000001];

int num;

int head[1001];

int cost[1001][1001];

void add(int l,int r,int v)

{

	node[num].u=r;

	node[num].v=v;

	node[num].next=head[l];

	head[l]=num++;

}

int nos;

int juli(int x,int y,int m)

{

	return abs(y/m-x/m)+abs(y%m-x%m);

}

int pre[1001];

int bfs()

{

	int visit[1001];

	int dist[1001];

	int i;

	for(i=0;i<=nos;i++)dist[i]=INF;

	memset(visit,0,sizeof(visit));

	memset(pre,-1,sizeof(pre));

	queue<int>q;

	q.push(0);

	visit[0]=1;

	dist[0]=0;

	while(!q.empty())

	{

		int e=q.front();

		q.pop();

		//cout<<"弹"<<" "<<e<<endl;

		visit[e]=0;

		for(i=head[e];i!=-1;i=node[i].next)

		{

			int r=node[i].u;

			if(dist[e]+node[i].v<dist[r]&&cost[e][r])

			{

				pre[r]=e;

				dist[r]=dist[e]+node[i].v;

			//	printf("pre[%d]=%d,dist[%d]=%d\n",r,e,r,dist[r]);

				if(!visit[r])

				{

					q.push(r);

				//	cout<<"入"<<" "<<r<<endl;

					visit[r]=1;

				}

			}

		}

	}

	if(dist[nos]!=INF)return 1;

	return 0;

}

void change()

{

	int minx;

	minx=INF;

	int i;

	for(i=nos;pre[i]!=-1;i=pre[i])

	{

		minx=min(minx,cost[pre[i]][i]);

	}

	for(i=nos;pre[i]!=-1;i=pre[i])

	{

		cost[pre[i]][i]-=minx;

		cost[i][pre[i]]+=minx;

	}

}

int main()

{

	int i,j,m,n;

	char str[100001];

	while(scanf("%d%d%*c",&n,&m)&&(n||m))

	{

		memset(head,-1,sizeof(head));

		memset(cost,0,sizeof(cost));

		int hm,mm;

		hm=mm=0;

		int ms[100001];

		int hs[100001];

		for(i=0;i<n;i++)

		{

			scanf("%s",str);

			for(j=0;j<m;j++)

			{

				if(str[j]=='m')ms[++mm]=i*m+j;

				else if(str[j]=='H')hs[++hm]=i*m+j;

			}

			getchar();

		}

		nos=hm+mm+1;

		for(i=1;i<=mm;i++)

		{

			add(0,i,0);

			add(i,0,0);

			cost[0][i]=1;

			cost[i][0]=0;

		}

		for(j=1;j<=hm;j++)

		{

			add(j+mm,mm+hm+1,0);

			add(mm+hm+1,j+mm,0);

			cost[j+mm][mm+hm+1]=1;

			cost[mm+hm+1][j+mm]=0;

		}

		for(i=1;i<=mm;i++)

		{

			for(j=1;j<=hm;j++)

			{

				int ju=juli(ms[i],hs[j],m);

				add(i,j+mm,ju);

				add(j+mm,i,-ju);

				cost[i][j+mm]=INF;

				cost[j+mm][i]=0;

			}

		}

		while(bfs())

		{

			change();

			//cout<<"------------------------------------"<<endl;

		}

		int sum=0;

		for(i=1;i<=mm;i++)

			for(j=1;j<=hm;j++)

			{

				sum+=(INF-cost[i][j+mm])*juli(ms[i],hs[j],m);

			}

		cout<<sum<<endl;

	}

	return 0;

}