POJ1273 Drainage Ditches（网络流&增广路）

/*PKU 1273 Drainage Ditches
  Edmonds-Karp Maximum Flow  O(VE^2)*/
#include <cstdio>
#include <queue>

using namespace std;

const int N = 210;

int ad[N][N],flow[N][N];
int pre[N],inc[N];
int n,m;

bool BFS(int ad[N][N],int src,int des)
{
    queue<int>q;
    int i,j;
   
    memset(pre,-1,sizeof(pre));
    q.push(src); inc[src] = 10000000; pre[src] = 1;
    while(!q.empty())
    {
        i = q.front();
        q.pop();
        for(j = 1; j <= n; j++)  //遍历所有与i相邻的节点
        {
            if(pre[j] == -1 && ad[i][j] > flow[i][j]) //ad是这条边上的总容量(即这条边上还可以让流增大)
            {
                inc[j]=inc[i]<ad[i][j] - flow[i][j]?inc[i]:ad[i][j] - flow[i][j]; //inc是到i节点为止增加的流量的最小值
       //一条路径上能流过的流量取决于这条路径上能流的最小的边上的流量
                pre[j] = i;  //标记j的前一个节点为i，记录路径信息
                if(j == des)  //如果到达汇点则返回
       return true;
                q.push(j);  //否则继续搜索下面的点（边）
            }
        }
    }
    return false;
}

int Edmonds_Karp(int ad[N][N],int src,int des)
{
    int i,ans = 0;
    while(BFS(ad,src,des))  //如果这条路径能够达到汇点（即上面的函数是遇到汇点返回的）
    {
        for(i = des; i != src; i = pre[i]) //从汇点到源点
        {
            flow[pre[i]][i] += inc[des]; 
    //找到一条增广路，然后更改这条路上的所有边的剩余流量，inc是这条路径上面本次增广的流量
            flow[i][pre[i]]  = - flow[pre[i]][i];  

        }
        ans += inc[des];  //最后要得到的是所有路径的最大流量的总和（即这个网络的最大传输量）
    }
    return ans;
}

int main()
{
    int i,j,k;
    while(scanf("%d%d",&m,&n)==2)
    {
        memset(ad,0,sizeof(ad));
        memset(flow,0,sizeof(flow));
        while(m--)
        {
            scanf("%d%d%d",&i,&j,&k);  //建图并且赋最大容量
            ad[i][j] += k;
        }
        printf("%d/n",Edmonds_Karp(ad,1,n));
    }
    return 0;
}