【ZOJ 3642】 Just Another Information Sharing Problem

一道网络流的模板题。

结果比赛的时候没A，感觉应该没问题的，当时就怀疑自己是不是模板敲错了，debug了好久，虽然是发现了几个地方的疏漏，但还是没A，然后。。。。然后我就去吃饭了- - 

后来听说其他队有用网络流A掉了，而且是套模板的，我就感觉应该是哪个小地方写错，就没冲动去查代码。。。

然后，回校第一天下午，没事干，又去查了下错误，把代码完完整整看了一遍，终于发现自己写错了一个地方- -|| 我就不明白了，当时为什么没有RE？！

<code>

#include<stdio.h>
#include<iostream>
#include<map>
#include<algorithm>
#include<string.h>
#include<assert.h>
#include<math.h>
using namespace std;
#define M 500
#define oo 300

struct node {
    int ax[15];
    int a, x;
} pnt[M];
int mark[2 * M], next[2 * M], head[2 * M], layer[2 * M], Q[M * 2];
int tot;

struct ss {
    int f, to, next, cap;
} edge[M*M];

void _add(int f, int to, int cap) {
    edge[tot].f = f;
    edge[tot].to = to;
    edge[tot].cap = cap;
    edge[tot].next = head[f];
    head[f] = tot++;
}

void add(int f, int to, int cap) {
    _add(f, to, cap);
    _add(to, f, 0);
}

bool BFS(int begin, int end) {
    int i, j, l, h, k, y;
    for (i = begin; i <= end; i++) layer[i] = -1;
    layer[begin] = 0; //深度
    l = h = 0;
    Q[l++] = begin;
    while (h < l) {
        k = Q[h++];
        for (i = head[k]; i != -1; i = edge[i].next) {
            y = edge[i].to;
            if (edge[i].cap > 0 && layer[y] == -1) {
                layer[y] = layer[k] + 1;
                if (y == end) return true;
                Q[l++] = y;
            }
        }
    }
    return false;
}

int DFS(int x, int exp, int end) { //exp 最小流量
    mark[x] = 1;
    if (x == end)return exp;
    int y, temp, i;
    for (i = next[x]; i != -1; i = edge[i].next) {
        y = edge[i].to;
        if (edge[i].cap > 0 && layer[y] == layer[x] + 1 && !mark[y]) {
            if ((temp = (DFS(y, min(exp, edge[i].cap), end))) > 0) {
                edge[i].cap -= temp;
                edge[i^1].cap += temp;
                return temp;
            }
        }
    }
    return 0;
}

int Dinic_flow(int begin, int end) {
    int i, ans = 0, flow;
    while (BFS(begin, end)) {
        for (i = begin; i <= end; i++)next[i] = head[i];
        while (1) {
            for (i = begin; i <= end; i++) mark[i] = 0; //mark用来标记是否访问过
            flow = DFS(begin, oo, end);
            if (flow == 0)break; //flow==0 表示已经找不到新的扩展路
            ans += flow;
        }
    }
    return ans;
}
map<int, int>mat;
int aa[M * 12], len;
int sot;

void deal(int n, int m) {
    mat.clear();
    sot = 1;
    int pre = -1;
    for (int i = 0; i < len; i++) {
        if (pre != aa[i]) {
            mat[aa[i]] = sot++;
            pre = aa[i];
        }
    }
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < pnt[i].a; j++) {
            assert(mat[pnt[i].ax[j]]);
            pnt[i].ax[j] = mat[pnt[i].ax[j]];
        }
    }
    memset(head, -1, sizeof (head));
    tot = 0;
    for (int i = 1; i < sot; i++) {
        int flag = 0;
        for (int j = 0; j < pnt[m].a; j++) {
            if (i == pnt[m].ax[j]) {
                flag = 1;
                break;
            }
        }
        if (flag == 0) {
            add(0, i, 1);
        }
    }
    for (int i = 0; i < n; i++) {
        if (i == m) continue;
        for (int j = 0; j < pnt[i].a; j++) {
            int kk = pnt[i].ax[j];
            add(kk, sot + i, 1);
        }
    }
    for (int i = 0; i < n; i++) {
        if (i == m)continue;
        add(sot + i, n + sot, pnt[i].x);
    }
}

int main() {
    int n, m, i, j;
    while (scanf("%d", &n) != EOF) {
        len = 0;
        for (i = 0; i < n; i++) {
            scanf("%d%d%d", &pnt[i].a, &pnt[i].x, &pnt[i].x);
            for (j = 0; j < pnt[i].a; j++) {
                scanf("%d", &pnt[i].ax[j]);
                aa[len++] = pnt[i].ax[j];
            }
        }
        sort(aa, aa + len);
        scanf("%d", &m);
        m--;
        deal(n, m);
        int res = Dinic_flow(0, sot + n) + pnt[m].a;
        printf("%d\n", res);
    }
    return 0;
}

</code>

不知为何，这几次zju月赛的题目偏向简单，以前都只能做1、2题的，现在一般都是4、5题。

上次跟聪哥吃饭，有听他说他们月赛好像是内部练习赛的题目放上去的。。。

不过总得来说，数据还是很强的~这点非常有感触- -