POJ1459 - Power Network

最新推荐文章于 2018-02-28 11:28:33 发布

VisJiao

最新推荐文章于 2018-02-28 11:28:33 发布

阅读量217

点赞数

分类专栏： ----OnlineJudge---- POJ ----图论---- 最大流

本文链接：https://blog.csdn.net/VicJiao/article/details/78813599

版权

----OnlineJudge---- 同时被 3 个专栏收录

37 篇文章 0 订阅

订阅专栏

----图论----

13 篇文章 0 订阅

订阅专栏

POJ

6 篇文章 0 订阅

订阅专栏

原题链接

题意简述

原题看了好几遍才看懂…
给出一个 $n(n\leq100)$ 个点， $m(m\leq n^2)$ 条边的有向图。 $n$ 个点中有 $np$ 个源点， $nc$ 个汇点，每个源点和汇点都有流出上限和流入上限。求最大流。

题解

建一个真 · 源点和一个真 · 汇点。真 · 源点向所有原源点连容量等于其上限的边，所有原汇点向真 · 汇点连容量等于其上限的边。跑一遍最大流即可。

Code

//Power Network
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
inline char gc()
{
    static char now[1<<16],*S,*T;
    if(S==T) {T=(S=now)+fread(now,1,1<<16,stdin); if(S==T) return EOF;}
    return *S++;
}
bool dgt(char ch) {return '0'<=ch&&ch<='9';}
inline int read()
{
    int x=0; char ch=gc();
    while(ch!=EOF&&!dgt(ch)) ch=gc();
    while(dgt(ch)) x=x*10+ch-'0',ch=gc();
    return x;
}
int const N=1e2+10;
int n,np,nc,m; int s,t;
int cnt,h[N];
struct edge{int v,c,nxt;} ed[N*N<<1];
void edAdd(int u,int v,int c)
{
    cnt++; ed[cnt].v=v,ed[cnt].c=c; ed[cnt].nxt=h[u],h[u]=cnt;
    cnt++; ed[cnt].v=u,ed[cnt].c=0; ed[cnt].nxt=h[v],h[v]=cnt;
}
int dpt[N]; int q[N],op,cl;
bool bfs()
{
    op=cl=0; memset(dpt,0,sizeof dpt);
    dpt[q[++cl]=s]=1;
    while(op<cl)
    {
        int u=q[++op]; if(u==t) break;
        for(int i=h[u];i;i=ed[i].nxt)
        {
            int v=ed[i].v,c=ed[i].c;
            if(dpt[v]==0 && c) dpt[q[++cl]=v]=dpt[u]+1;
        }
    }
    return dpt[t]?true:false;
}
int fill(int u,int in)
{
    if(u==t || in==0) return in;
    int out=0;
    for(int i=h[u];i;i=ed[i].nxt)
    {
        int v=ed[i].v,c=ed[i].c;
        if(dpt[v]!=dpt[u]+1 || !c) continue;
        int flow=fill(v,min(in-out,c));
        if(!flow) dpt[v]==0;
        else out+=flow,ed[i].c-=flow,ed[i^1].c+=flow;
    }
    return out;
}
int main()
{
    while(true)
    {
        n=read(),np=read(),nc=read(); m=read();
        if(gc()==EOF) break;
        s=0,t=n+1; cnt=1; memset(h,0,sizeof h);
        for(int i=1;i<=m;i++) {int u=read()+1,v=read()+1,c=read(); edAdd(u,v,c);}
        for(int i=1;i<=np;i++) {int u=read()+1,p=read(); edAdd(0,u,p);}
        for(int i=1;i<=nc;i++) {int u=read()+1,c=read(); edAdd(u,t,c);}
        int ans=0;
        while(bfs()) ans+=fill(s,0x7FFFFFFF);
        printf("%d\n",ans);
    }
    return 0;
}