POJ1459 Power Network 网络流

　　题中给出了各种点的定义，其实我们只要建立两个虚拟节点就可以了。一个是源点流到所有的p点，一个是汇点，所有的c点流过去，再求最大流就可以了。当然题中的输入数据比较恶心，这里用了字符串的处理方式来解决。

　　代码如下：

#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <queue>
#define S N
#define T (N+1)
#define INF 0x3fffffff
#define MAXN 105
using namespace std;

int N, np, nc, M, flow[MAXN][MAXN], cap[MAXN][MAXN];
int c[MAXN], p[MAXN], Maxflow;

void getint(int &t)
{
    char c;
    while (c = getchar(), c < '0' || c > '9') {};
    t = c - '0';
    while (c = getchar(), c >= '0' && c <= '9') {
        t = t * 10 + c - '0';
    }
}

void bfs()
{
    bool finish = false;
    int pos;
    while (!finish) {
        memset(c, 0, sizeof (c));
        c[S] = INF;  // 虚拟出一个无穷的流从源点流出来
        queue<int>q;
        q.push(S);
        while (!q.empty()) {
            if (c[T]) {
                break;
            }
            pos = q.front();
            q.pop();
            for (int i = 0; i <= N+1; ++i) {
                if (!c[i] && flow[pos][i] < cap[pos][i]) {
                // 前面这个!c[i]很重要，它表示本轮中推送过流量的点不能作为二次推流的节点
                    c[i] = min(cap[pos][i]-flow[pos][i], c[pos]); // 这个值一定是大于零的
                    p[i] = pos;  // 记录路径，以便后面的更新
                    q.push(i);
                }
            }
        }
        if (c[T] == 0) {  // 本次寻找增广路径失败
            finish = true;
        }
        else {
            Maxflow += c[T];  // 流入到汇点的新流量
            pos = T;
            while (pos != S) {  // 没有回溯到源点
                flow[p[pos]][pos] += c[T];
                flow[pos][p[pos]] -= c[T];
                pos = p[pos];
            }
        }
    }
}

int main()
{
    char s[15];
    int x, y, z; // 从0号节点开始，我们可以设置一个N号节点，并建立c到N节点的流量无上限的边
    while (scanf("%d %d %d %d", &N, &np, &nc, &M) == 4) {
        memset(cap, 0, sizeof (cap));
        memset(flow, 0, sizeof (flow));
        Maxflow = 0;
        for (int i = 0; i < M; ++i) {
            getint(x), getint(y), getint(z);
        //    printf("%d %d %d\n", x, y, z);
            cap[x][y] = z;  // 建立边的容量
        }
        for (int i = 0; i < np; ++i) {
            getint(x), getint(y);
            cap[N][x] = y;  // 创建源点
        }
        for (int i = 0; i < nc; ++i) {
            getint(x), getint(y);
            cap[x][N+1] = y;  // 创建汇点
        }
        bfs();
        printf("%d\n", Maxflow);
    }
    return 0;
}