POJ 1273 & & HDU 1532 Drainage Ditches (최대 흐름 - DINIC)

설명 FJ 는 물 을 연못 에서 개울 로 끌 어 올 려 야 한다. 그 사이 에 M 개의 점 과 N 개의 수로 가 있 고 수로 의 연결 점 이 있 으 며 수 로 는 유량 제한 이 있다. 연못 에서 얼마나 많은 물 을 개울 로 흘 릴 수 있 느 냐 고 물 었 다.Input 여러 조 의 입력, 각 조 의 입력 첫 번 째 행위 두 개의 정수 M 과 N 은 각각 수로 수 와 점 수 를 표시 합 니 다. 그 다음 에 M 줄 은 줄 마다 세 개의 정수 a, b 와 c 는 a 점 에서 b 점 까지 유량 이 c 인 수로 가 있 음 을 표시 합 니 다. 파일 끝 에 출력 은 각 조 의 사례 에 대해 연못 (번호 1) 에서 개울 (번호 N) 까지 출력 할 수 있 습 니 다.최대 유량 Sample Input 5, 4, 1, 2, 40, 1, 4, 20, 2, 3, 30, 4, 10 Sample Output 50 Solution 최대 스 트림 누 드 문제

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
#define maxn 222 
#define maxm 555
#define INF 0x3f3f3f3f 
int head[maxn],cur[maxn],d[maxn],st[maxm],s,e,no,n;
struct point
{
    int u,v,flow,next;
    point(){};
    point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)
{
    p[no]=point(x,y,head[x],z); 
    head[x]=no++;
    p[no]=point(y,x,head[y],0); 
    head[y]=no++;
}
void init()
{
    memset(head,-1,sizeof(head));
    no=0;
}
bool bfs()
{
    int i,x,y;
    queue<int>q;
    memset(d,-1,sizeof(d));
    d[s]=0; 
    q.push(s);
    while(!q.empty())
    {
        x=q.front();    
        q.pop();
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&& d[y = p[i].v]<0)
            {
                d[y]=d[x]+1;
                if(y==e)    
                    return true;
                q.push(y);
            }
        }
    }
    return false;
}
int dinic()
{
    int i,loc,top,x=s,nowflow,maxflow=0;
    while(bfs()){
        for(i=s;i<=e;i++)   
            cur[i]=head[i];
        top=0;
        while(true)
        {
            if(x==e)
            {
                nowflow=INF;
                for(i=0;i<top;i++)
                {
                    if(nowflow>p[st[i]].flow)
                    {
                        nowflow=p[st[i]].flow;
                        loc=i;
                    }
                }
                for(i=0;i<top;i++)
                {
                    p[st[i]].flow-=nowflow;
                    p[st[i]^1].flow+=nowflow;
                }
                maxflow+=nowflow;
                top=loc;    
                x=p[st[top]].u;
            }
            for(i=cur[x];i!=-1;i=p[i].next)
                if(p[i].flow&&d[p[i].v]==d[x]+1) 
                    break;
            cur[x]=i;
            if(i!=-1)
            {
                st[top++]=i;
                x=p[i].v;
            }
            else 
            {
                if(!top)    
                    break;
                d[x]=-1;
                x=p[st[--top]].u;
            }
        }
    }
    return maxflow;
}
int main()
{
    int m;
    while(~scanf("%d%d",&m,&n))
    {
        init();//    
        s=1;//      
        e=n;//      
        int u,v,c;
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&u,&v,&c);
            add(u,v,c);
        }
        printf("%d
",dinic());
    }
    return 0;
}