POJ3436 ACM Computer Factory(最大流)

题目链接。

分析：

题意很难懂。

大体是这样的：给每个点的具体情况，1.容量 2。进入状态 3.出去状态。求最大流。

因为有很多点，所以如果一个点的出去状态满足另一个点的进入状态，则这两个点可以连一条边。容量为两者容量的较小值。

再建立一个超源、一个超汇。让超源与所有进入状态全为0或者不全为0但只包含0和2的点连边，同时让所有出去状态全部为1的与超汇连边。

然后求最大流.

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

const int maxn = ;
const int INF = (<<);

int in[maxn][], cap[maxn][maxn], n, flow[maxn][maxn];

int EK(int s, int t) {
    queue<int> q;
    int p[maxn], a[maxn];

    int f = ;

    memset(flow, , sizeof(flow));

    while(true) {
        memset(a, , sizeof(a));
        a[s] = INF;
        q.push(s);
        while(!q.empty()) {
            int u = q.front(); q.pop();
            for(int v = ; v < n; v++) if(!a[v] && cap[u][v] > flow[u][v]) {
                p[v] = u; q.push(v);
                a[v] = min(a[u], cap[u][v]-flow[u][v]);
            }
        }

        if(a[t] == ) break;
        for(int u=t; u != s; u = p[u]) {
            flow[p[u]][u] += a[t];
            flow[u][p[u]] -= a[t];
        }
        f += a[t];
    }
    return f;
}

int main() {
    int p;

    while(scanf("%d%d", &p, &n) == ) {

        memset(cap, , sizeof(cap));

        for(int i=; i<*p+; i++) {        //初始化超源，超汇
            in[][i] = ;
            in[n+][i] =;
        }

        for(int i=; i<=n; i++) {           //输入数据
            for(int j=; j<*p+; j++) {
                scanf("%d", &in[i][j]);
            }
        }

        n+=;                               //点数+2

        for(int i=; i<n; i++) {            //将可以连通的点，记录,算出每条边的容量
            for(int j=; j<n; j++) {
                if(i == j) continue;

                bool flag = true;
                for(int k=; k<=p; k++) {
                    if( !((in[j][k] == ) || (in[i][k+p] == in[j][k])) ) {
                        flag = false;
                    }
                }
                if(flag && i == ) cap[i][j] = in[j][];
                else if(flag && j == n-) cap[i][j] = in[i][];
                else if(flag) cap[i][j] += min(in[i][], in[j][]);
            }
        }

        printf("%d ", EK(, n-));          //增广路算法

        int cnt = ;                        //计数
        for(int i=; i<n-; i++) {
            for(int j=; j<n-; j++) {
                if(flow[i][j] > ) cnt++;
            }
        }
        printf("%d\n", cnt);

        for(int i=; i<n-; i++) {          //输出
            for(int j=; j<n-; j++) {
                if(flow[i][j] > ) printf("%d %d %d\n", i, j, flow[i][j]);
            }
        }
    }

    return ;
}