poj 2195 Going Home

http://poj.org/problem?id=2195

任何问题 都是 难了不会 会了不难 难就难在由不会变成会

尤其是刚接触到一个新知识点的时候硬着头皮，耐心地去看去理解，一定能学会，然后你就会发现它

原来并不难

本题是一个最小费用流

详解见代码注释

#include<iostream>

#include<string>

#include<cstring>

#include<algorithm>

#include<queue>

#include<cstdio>



using namespace std;

const int N=1001;//后台数据水了 开1000就过了 其实更小也过不过也不能太小 其实理论上是10000

const int M=100000005;

int flow[N][N];//保存流

int pay[N][N];//保存费用

struct node

{

    int x,y;

    char c;

}mem[10010];//保存输入的m点和H点



void find_flow_pay(int i,int j)

{

   if(mem[i].c==mem[j].c)//两个点必须一个是m一个是H

   return ;

   if(mem[i].c=='m')

   {

       pay[i][j]=abs(mem[j].x-mem[i].x)+abs(mem[j].y-mem[i].y);//求费用

       pay[j][i]=-pay[i][j];//反向费用变负的

       flow[i][j]=1;//正向流为1 反向为0

   }

   else//同上 只不过这个 j点为m而已

   {

       pay[j][i]=abs(mem[j].x-mem[i].x)+abs(mem[j].y-mem[i].y);

       pay[i][j]=-pay[j][i];

       flow[j][i]=1;

   }

}

int Spfa(int n)

{

    bool in[N];//是否在队列中

    int dist[N];//距离

    int f[N];//前驱点

    memset(in,false,sizeof(in));

    for(int i=1;i<=n;++i)

    dist[i]=M;//初始化最大

    dist[0]=0;//超级源点为0

    queue<int>str;

    str.push(0);

    in[0]=true;

    while(!str.empty())//由于没有负环 所以直接找到队列为空就开

    {

        int x=str.front();

        str.pop();

        in[x]=false;

        for(int i=1;i<=n;++i)

        {

            if(flow[x][i]>0&&dist[x ]+pay[x][i]<dist[i])

            {

                dist[i]=dist[x]+pay[x][i];//更新距离 其实就是费用

                f[i]=x;//标记前驱

                if(in[i]==false)//当i不在队列中就进队列 并标记

                {

                    in[i]=true;

                    str.push(i);

                }

            }

        }

    }

    if(dist[n]==M)//到达不了超级终端则返回0

    return 0;

    int k=n;

    while(k!=0)//更新流

    {

        int pre=f[k];

        --flow[pre][k];

        ++flow[k][pre];

        k=pre;

    }

    return dist[n];



}

int main()

{

    int n,m;

    while(cin>>n>>m)

    {

        if(n==0&&m==0)

        break;

        int I=1;

        char ctemp;

        memset(flow,0,sizeof(flow));

        memset(pay,0,sizeof(pay));

        for(int i=1;i<=n;++i)

        {

            getchar();//吃掉回车

            for(int j=1;j<=m;++j)

            {

                scanf("%c",&ctemp);

                if(ctemp!='.')

                {

                    mem[I].c=ctemp;mem[I].x=i;mem[I].y=j;

                    for(int l=1;l<I;++l)

                    {

                        find_flow_pay(l,I);//球两点的流和花费

                    }

                    ++I;

                }

            }

        }

        for(int i=1;i<I;++i)//0 为超级源点 I 为超级终点

        {

            if(mem[i].c=='m')

            {

                flow[0][i]=1;

            }

            else

            {

                flow[i][I]=1;

            }

        }

        int ans=0;

        while(1)

        {

           int k=Spfa(I);

           if(k==0)//无流更新

           break;

           ans=ans+k;

        }

        cout<<ans<<endl;

    }

    return 0;

}