uoj#158. 【清华集训2015】静态仙人掌

http://uoj.ac/problem/158

预处理dfs序，询问转为区间1的个数，用可持久化bitset预处理出所有可能的修改对应哪些位置，然后用一个bitset维护当前每个点的状态，修改时可以用xor实现

#include<bits/stdc++.h>
int _(){
    int x=,c=getchar();
    while(c<)c=getchar();
    while(c>)x=x*+c-,c=getchar();
    return x;
}
const int N=;
typedef unsigned long long u64;
typedef u64 block[];
typedef block *bits[];
block pool[N*],*pp=pool;
int n,m,q,fa[N],son[N],dep[N],dl[N],dr[N],cb;
std::vector<int>e[N];
int dfn[N],mxd[N],tk=;
bits f[N],fl[N],fr[N];
void init(){
    std::queue<int>q;
    q.push();dep[]=;
    while(q.size()){
        int w=q.front();q.pop();
        for(unsigned i=;i<e[w].size();++i){
            int u=e[w][i];
            if(!dep[u])q.push(u),fa[u]=w,dep[u]=dep[w]+;
            else if(dep[u]==dep[w]&&w<u){
                for(int x=w,y=u,x0=u,y0=w;x!=y;x=fa[x],y=fa[y]){
                    dl[x]=dr[y]=;
                    son[x]=x0,son[y]=y0;
                    x0=x,y0=y;
                }
            }
        }
    }
}
void dfs(int w){
    dfn[w]=++tk;
    for(unsigned i=;i<e[w].size();++i){
        int u=e[w][i];
        if(dep[u]>dep[w]&&u!=son[w])dfs(u);
    }
    mxd[w]=tk;
    if(dep[son[w]]>dep[w]&&son[w])dfs(son[w]);
}
void cpy(bits a,bits b,int x){
    for(int i=;i<cb;++i)a[i]=b[i];
    memcpy(a[x>>]=++pp,b[x>>],sizeof(block));
    (*pp)[x>>&]|=1llu<<(x&);
}
void get(int w){
    if(f[w][])return;
    get(fa[w]);
    cpy(f[w],f[fa[w]],dfn[w]);
}
void getr(int);
void getl(int w){
    if(fl[w][])return;
    int u=dr[w]?son[w]:fa[w];
    if(u==fa[w]&&w!=son[u]){
        dl[u]?getr(u):getl(u);
        cpy(fl[w],dl[u]?fr[u]:fl[u],dfn[w]);
    }else{
        getl(u);
        cpy(fl[w],fl[u],dfn[w]);
    }
}
void getr(int w){
    if(fr[w][])return;
    int u=dl[w]?son[w]:fa[w];
    if(u==fa[w]&&w!=son[u]){
        dl[u]?getr(u):getl(u);
        cpy(fr[w],dl[u]?fr[u]:fl[u],dfn[w]);
    }else{
        getr(u);
        cpy(fr[w],fr[u],dfn[w]);
    }
}
u64 ans[];
int test(int x){
    return ans[x>>]>>(x&)&;
}
int q0(int l,int r){
    int l1=l>>,r1=r>>,s=;
    if(l1==r1)for(int i=l;i<=r;++i)s+=test(i);
    else{
        for(int i=l;(i>>)==l1;++i)s+=test(i);
        for(int i=r;(i>>)==r1;--i)s+=test(i);
        for(int i=l1+;i<r1;++i){
            u64 x=ans[i];
            s+=__builtin_popcount(x&0xffffffff);
            s+=__builtin_popcount(x>>);
        }
    }
    return r-l+-s;
}
void rev(bits w){
    u64*A=ans,*B;
    for(int i=;i<cb;++i){
        B=*(w[i]);
        for(int j=;j<;j+=){
            A[]^=B[];A[]^=B[];A[]^=B[];A[]^=B[];
            A[]^=B[];A[]^=B[];A[]^=B[];A[]^=B[];
            A+=,B+=;
        }
    }
}
int main(){
    n=_();m=_();q=_();
    cb=(n+>>)+;
    for(int i=;i<cb;++i)f[][i]=fl[][i]=fr[][i]=pp;
    for(int i=,a,b;i<=m;++i){
        a=_();b=_();
        e[a].push_back(b);
        e[b].push_back(a);
    }
    init();
    dfs();
    for(int i=;i<=n;++i){
        get(i);
        getl(i);
        getr(i);
    }
    for(int i=;i<=q;++i){
        int op=_(),x=_();
        if(op==)rev(f[x]);
        else if(op==)rev(dl[x]?fr[x]:fl[x]);
        else printf("%d\n",q0(dfn[x],mxd[x]));
    }
    return ;
}