POJ1273：Drainage Ditches——题解

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

题目大意：

n点m边网络流，求1-n最大流。

——————————————

网络流板子，切了。

#include <cstdio>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=;
const int INF=;
inline int read(){
    int X=,w=;char ch=;
    while(ch<''||ch>''){w|=ch=='-';ch=getchar();}
    while(ch>=''&&ch<='')X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
struct node{
    int next;
    int to;
    int w;
}edge[*maxn];
int head[*maxn],cnt=-;
void add(int u,int v,int w){//u起点v终点w容量
    cnt++;
    edge[cnt].to=v;
    edge[cnt].w=w;
    edge[cnt].next=head[u];
    head[u]=cnt;
    return;
}
int lev[maxn],cur[maxn];//lev层数，cur[i]为以i为起点的边的编号
bool bfs(int m){//强制1为源点，m为汇点
    int dui[m],r=;//队列和右指针
    for(int i=;i<=m;i++){//初始化
        lev[i]=-;
        cur[i]=head[i];
    }
    dui[]=,lev[]=;
    int u,v;//u起点v终点
    for(int l=;l<=r;l++){//左指针
        u=dui[l];
        for(int e=head[u];e!=-;e=edge[e].next){
            v=edge[e].to;
            if(edge[e].w>&&lev[v]==-){//1.能走 2.未分层
                lev[v]=lev[u]+;
                r++;
                dui[r]=v;//v入队
                if(v==m)return ;//分层完毕
            }
        }
    }
    return ;//无法分层
}
int dinic(int u,int flow,int m){//u当前点，flow为下面的点能够分配多大的流量,m终点
    if(u==m)return flow;//终点直接全流入
    int res=,delta;//res实际流量
    for(int &e=cur[u];e!=-;e=edge[e].next){//'&'相当于cur[u]=e;即流满的边不会再被扫一次
        int v=edge[e].to;
        if(edge[e].w>&&lev[u]<lev[v]){//只能从低层往高层流
            delta=dinic(v,min(edge[e].w,flow-res),m);
            if(delta>){//如果增广
                edge[e].w-=delta;//正向边容量减少
                edge[e^].w+=delta;//反向边仍量增加（暗示退流）
                res+=delta;//扩张流量增加
                if(res==flow)break;//可流的都流完了，及时跳出
            }
        }
    }
    if(res!=flow)lev[u]=-;//没流完，说明以后不能从这个点流出任何流量，那就不需要这个点了
    return res;
}
int main(){
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF){
    cnt=-;
    memset(head,-,sizeof(head));
    for(int i=;i<=n;i++){
        int s=read();
        int e=read();
        int c=read();
        add(s,e,c);
        add(e,s,);
    }
    int ans=;
    while(bfs(m)==)ans+=dinic(,INF,m);
    printf("%d\n",ans);
    }
    return ;
}