ZOJ 1734 Power Network （EK算法 模板题）

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=1734

题目大意：

电力发电与传输问题。若干个发电站与若干个消费点。

题目分析： 

初学EK算法，模板提。
不断的按bfs层次去收集流量。
寻找增广路径的方法不确定导致时间复杂度不确定。

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=120;
const int ub=1e6;
const double inf=1e-4;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}

/*
题目大意：
电力发电与传输问题。

题目分析：
初学EK算法，模板提。
不断的按bfs层次去收集流量。
寻找增广路径的方法不确定导致时间复杂度不确定。
*/

int n,np,nc,m,cap[maxn][maxn];
int from,to,value;

int EK(int s,int t){
    queue<int> Q;
    int flow[maxn][maxn],ret=0;
    int minv[maxn],pre[maxn];///最小流量,前继链接
    mst(flow,0);///局部变量要初始化
    while(1){
        Q.push(n);
        mst(minv,0),minv[n]=mod;
        while(!Q.empty()){
            int tp=Q.front();Q.pop();
            for(int i=0;i<=n+1;i++) if(cap[tp][i]>flow[tp][i]&&!minv[i]){
                pre[i]=tp;
                Q.push(i);
                minv[i]=min(minv[tp],cap[tp][i]-flow[tp][i]);
            }
        }
        if(minv[t]==0) break;///无法再约束
        for(int i=t;i!=s;i=pre[i]){
            flow[i][pre[i]]-=minv[t];
            flow[pre[i]][i]+=minv[t];
        }
        ret+=minv[t];
    }
    return ret;
}

int main(){
    while(scanf("%d %d %d %d",&n,&np,&nc,&m)!=EOF){
        mst(cap,0);
        while(m--){
            scanf(" (%d,%d)%d",&from,&to,&value);
            cap[from][to]=value;
        }
        while(np--){
            scanf(" (%d)%d",&from,&value);
            cap[n][from]=value;
        }
        while(nc--){
            scanf(" (%d)%d",&to,&value);
            cap[to][n+1]=value;
        }
        ///初始网络构建完毕
        printf("%d\n",EK(n,n+1));
    }
    return 0;
}