BZOJ4316: 小C的独立集 圆方树+树形DP

版权声明：xgc原创文章，未经允许不得转载。 https://blog.csdn.net/xgc_woker/article/details/82501449

Description
求一个仙人掌无向图最大独立集。

---

Sample Input
5 6
1 2
2 3
3 1
3 4
4 5
3 5

---

Sample Output
2

---

先建圆方树。
建好后做树形DP。
对于每一个环，将它拆环分两种情况讨论。
其他的点直接树形DP。
圆方树，唔

---

#include <cstdio>
#include <cstring>

using namespace std;
int _max(int x, int y) {return x > y ? x : y;}

struct edge {
    int x, y, next;
} e[130000], c[130000]; int len, cnt, last[51000], lt[110000];
int u, n, id, f[51000][2], fa[51000], dfn[51000];
bool cir[51000];
int sta[51000][2];

void ins(int x, int y) {
    e[++len].x = x; e[len].y = y;
    e[len].next = last[x]; last[x] = len;
}

void add(int x, int y) {
    c[++cnt].x = x; c[cnt].y = y;
    c[cnt].next = lt[x]; lt[x] = cnt;
}

void dfs(int x) {
    dfn[x] = ++id;
    for(int k = last[x]; k; k = e[k].next) {
        int y = e[k].y;
        if(y != fa[x]) {
            if(!dfn[y]) {
                cir[x] = 0;
                fa[y] = x; dfs(y);
                if(!cir[x]) add(x, y), add(y, x);
            } else if(dfn[y] < dfn[x]) {
                int i = x; n++;
                add(i, n), add(n, i);
                do {
                    i = fa[i];
                    add(n, i), add(i, n);
                    cir[i] = 1;
                } while(i != y);
            }
        }
    }
}

void treedp(int x, int fa) {
    if(x <= u) {
        f[x][0] = 0, f[x][1] = 1;
        for(int k = lt[x]; k; k = c[k].next) {
            int y = c[k].y;
            if(y != fa) {
                treedp(y, x);
                if(y <= u) {
                    f[x][0] += _max(f[y][0], f[y][1]);
                    f[x][1] += f[y][0];
                }
            }
        }
    } else {
        int tp = 0;
        for(int k = lt[x]; k; k = c[k].next) {
            int y = c[k].y;
            if(y != fa) treedp(y, x);
        } for(int k = lt[x]; k; k = c[k].next) sta[++tp][0] = f[c[k].y][0], sta[tp][1] = f[c[k].y][1];
        for(int i = tp - 1; i; i--) sta[i][0] += _max(sta[i + 1][1], sta[i + 1][0]), sta[i][1] += sta[i + 1][0];
        f[fa][0] = sta[1][0];
        for(int i = 1; i < tp; i++) sta[i][0] -= _max(sta[i + 1][1], sta[i + 1][0]), sta[i][1] -= sta[i + 1][0];
        sta[tp][1] = -999999999;
        for(int i = tp - 1; i; i--) sta[i][0] += _max(sta[i + 1][1], sta[i + 1][0]), sta[i][1] += sta[i + 1][0];
        f[fa][1] = sta[1][1];
    }
}

int main() {
    int m; scanf("%d%d", &n, &m); u = n;
    for(int i = 1; i <= m; i++) {
        int x, y; scanf("%d%d", &x, &y);
        ins(x, y), ins(y, x);
    } dfs(1);
    treedp(1, 0);
    printf("%d\n", _max(f[1][0], f[1][1]));
    return 0;
}