POJ 1459 构图+最大流（Edmond_karp模版）

题目链接： http://poj.org/problem?id=1459

题目大意：

　　给定一个有向图，n个点，m条边( 1<=n<=100, m<=n^2 )，点有三类： np个 power station，nc个consumer，其余的是普通的点，其中np个power station可以产生power供给网络传递，nc个consumer消耗能量，每个power station有一个最大的供应值value，每个consumer有一个最大的消耗值value，问最终可以消耗能量的最大值;

分析：

　　我的构图方式： 虚拟一个源点S（S=0）,其余输入的点都+1，虚拟一个汇点T（T=n+1）,S向所有的power station引一条边，边权（容量）为相应的power station的value值，所有的consumer向T引一天边，边权（容量）为相应consumer的value值;

　　之后用Edmond_karp模版求最大流即可;

注意：

　　我开始想用%*来处理一些多余的空格，回车的输入，后来失败了，后来读入m个串，np个串，nc个串来解决问题（下面的代码就是这样），然后又看了大众的代码，发现都是这样输入的：

scanf(" (%d,%d)%d", ...); // 即括号前面加了空格;

scanf(" (%d)%d", ...);

代码：

/*1459    Accepted    212K    375MS    C++    2908B    2012-04-23 19:59:30*/
#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;

#define mpair make_pair
#define pii pair<int,int>
#define MM(a,b) memset(a,b,sizeof(a));
typedef long long lld;
typedef unsigned long long u64;
template<class T> bool up_max(T& a,const T& b){return b>a? a=b,1 : 0;}
template<class T> bool up_min(T& a,const T& b){return b<a? a=b,1 : 0;}
#define maxn 110
#define inf 2100000000
int n, np, nc, m, S, T;
int mark[maxn]; // 0: station;  1: consumer; 2: anthter;
// all station should be linked to the sourse; and all consumer should be linked to the sink;
int g[maxn][maxn];
void Read(){
    int i,j,u,v,w;
    char ch[100];
    for(i=0;i<=n+1;++i)
        for(j=0;j<=n+1;++j)
            g[i][j]= 0;

    for(i=1;i<=m;++i){
        //scanf( "%*[ ]%*['\n'](%d,%d)%d", &u, &v, &w );
        //scanf( "%*[ ](%d,%d)%d", &u, &v, &w );
        scanf("%s", ch);
        u= 0;
        for(j=1;ch[j]!=',';++j)
            u= u*10 + ch[j]-'0';
        v= 0;
        for(++j;ch[j]!=')';++j)
            v= v*10 + ch[j]-'0';
        w= 0;
        for(++j;ch[j];++j)
            w= w*10 + ch[j]-'0';
        ++u, ++v;
        g[u][v]= w;
    }


    for(i=1;i<=np;++i){
        //scanf("%*[ ](%d)%d", &u, &w);
        //scanf("%*[ ]%*['\n'](%d)%d", &u, &w);
        scanf("%s", ch);
        u= 0;
        for(j=1;ch[j]!=')';++j)
            u= u*10 + ch[j]-'0';
        w= 0;
        for(++j;ch[j];++j)
            w= w*10 + ch[j]-'0';
        ++u;
        g[0][u]= w;
    }
    for(i=1;i<=nc;++i){
        //scanf("%*[ ](%d)%d", &u, &w);
        //scanf("%*[ ]%*['\n'](%d)%d", &u, &w);
        scanf("%s", ch);
        u= 0;
        for(j=1;ch[j]!=')';++j)
            u= u*10 + ch[j]-'0';
        w= 0;
        for(++j;ch[j];++j)
            w= w*10 + ch[j]-'0';
        ++u;
        g[u][T]= w;
    }
}

int pre[maxn];
int que[maxn];
bool vis[maxn];
bool bfs(){
    MM( vis, 0 );
    int head= 0, tail= 0;
    que[tail++]= S;
    vis[S]= 1;
    while(head<tail){
        int u= que[head++];
        for(int i=1;i<=T;++i){
            if( g[u][i] && !vis[i] ){
                vis[i]= 1;
                pre[i]= u;
                if(i==T) return 1;
                que[tail++]= i;
            }
        }
    }
    return 0;
}
int Edmond_karp(){
    int i,ans= 0;
    while( bfs() ){
        int t= inf;
        for(i= T; i!= S; i= pre[i])
            up_min( t, g[pre[i]][i] );
        ans+= t;
        for(i=T; i!=S; i= pre[i]){
            g[pre[i]][i]-= t;
            g[i][pre[i]]+= t;
        }
    }
    return ans;
}

int main()
{
    while( scanf("%d%d%d%d", &n, &np, &nc, &m) != EOF ){
        S= 0, T= n+1;
        Read();
        printf("%d\n", Edmond_karp() );
    }
}