BZOJ4316: 小C的独立集

建出圆方树在上面做树形dp
对于一个圆点，
f[i][0/1]表示点i是否选，i的子树内的最大独立集
对于一个环，我们用他的方点代表除了父亲剩余的链
g[i][0/1][0/1]表示i这个环（不包括环的父亲），两端是否选的最大独立集

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1e9
using namespace std;

inline void up(int &x,const int &y){if(x<y)x=y;}
inline void down(int &x,const int &y){if(x>y)x=y;}
const int maxn = 110000;

int n,m;
int e[maxn][2];
struct edge{int y,nex;}a[maxn<<1],b[maxn<<1];
int len,fir[maxn],len2,fir2[maxn];
inline void ins(const int x,const int y){a[++len]=(edge){y,fir[x]};fir[x]=len;}
inline void ins2(const int x,const int y){b[++len2]=(edge){y,fir2[x]};fir2[x]=len2;}

vector<int>V[maxn]; int cnt;
int T[maxn],tp;
int did,dfn[maxn],low[maxn];
void tarjan(const int x,const int pre)
{
    dfn[x]=low[x]=++did; T[++tp]=x;
    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(k!=(pre^1))
    {
        if(!dfn[y])
        {
            tarjan(y,k); down(low[x],low[y]);
            if(low[y]==dfn[x])
            {
                V[++cnt].push_back(x);
                int la=0,ti=1;
                while(la!=y) V[cnt].push_back(la=T[tp--]);
            }
            else if(low[y]>dfn[x]) tp--;
        }
        else down(low[x],dfn[y]);
    }
}
void build()
{
    for(int i=1;i<=m;i++)
    {
        int x=e[i][0],y=e[i][1];
        if(dfn[x]>dfn[y]) swap(x,y);
        if(low[y]>dfn[x]) ins2(x,y),ins2(y,x);
    }
    for(int i=1;i<=cnt;i++)
    {
        for(int j=0;j<V[i].size();j++) ins2(V[i][j],n+i),ins2(n+i,V[i][j]);
    }
}
int f[maxn][2],g[maxn][2][2],tg[2][2];
void dp(const int x,const int fa)
{
    if(x<=n)
    {
        f[x][0]=0,f[x][1]=1;
        for(int k=fir2[x],y=b[k].y;k;k=b[k].nex,y=b[k].y) if(y!=fa)
        {
            dp(y,x);
            if(y>n) f[x][0]+=f[y][0],f[x][1]+=f[y][1];
            else f[x][0]+=max(f[y][0],f[y][1]),f[x][1]+=f[y][0];
        }
    }
    else
    {
        int i=x-n;
        for(int j=1;j<V[i].size();j++) dp(V[i][j],x);
        g[x][0][0]=f[V[i][1]][0],g[x][1][1]=f[V[i][1]][1];
        for(int j=2;j<V[i].size();j++)
        {
            const int y=V[i][j];
            tg[0][0]=max(g[x][0][1],g[x][0][0])+f[y][0];
            tg[0][1]=g[x][0][0]+f[y][1];
            tg[1][0]=max(g[x][1][0],g[x][1][1])+f[y][0];
            tg[1][1]=g[x][1][0]+f[y][1];
            for(int ti=0;ti<2;ti++) for(int tj=0;tj<2;tj++) g[x][ti][tj]=tg[ti][tj];
        }
        for(int i=0;i<2;i++) for(int j=0;j<2;j++) up(f[x][0],g[x][i][j]);
        f[x][1]=g[x][0][0];
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    len=1;
    for(int i=1;i<=m;i++)
    {
        int x,y; scanf("%d%d",&x,&y);
        e[i][0]=x,e[i][1]=y;
        ins(x,y); ins(y,x);
    }
    did=cnt=tp=0; tarjan(1,0);
    len2=0; build();

    dp(1,0);
    printf("%d\n",max(f[1][0],f[1][1]));

    return 0;
}