BZOJ2815: [ZJOI2012]灾难（支配树）

传送门

题解：
支配树裸题了。。

写这篇博客只是为了提供一份支配树的代码。。

#include<bits/stdc++.h>using namespace std;inline int rd(){    char ch=getchar();int i=0,f=1;    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}    while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=getchar();}    return i*f;}inline void W(int x){    static int buf[50];    if(!x){putchar('0');return;}    if(x<0){putchar('-');x=-x;}    while(x){buf[++buf[0]]=x%10;x/=10;}    while(buf[0])putchar(buf[buf[0]--]+'0');}inline void gmin(int &x,int y){(x>y)?(x=y):0;}const int N=1e5+50;int rt,n,sdom[N],idom[N],id[N],dfn[N],mnsd[N],tot,anc[N],fa[N];vector<int>pre[N],buc[N],edge[N];struct DominatorTree{    int sze[N],v[N],nt[N],g[N],ec;    inline void add(int x,int y){nt[++ec]=g[x];g[x]=ec;v[ec]=y;}    inline void dfs(int x){        sze[x]=1;        for(int j=g[x];j;j=nt[j])dfs(v[j]),sze[x]+=sze[v[j]];    }}g;inline void dfs(int x,int f=0){    dfn[x]=++tot;id[tot]=x;fa[x]=f;    sdom[x]=mnsd[x]=anc[x]=x;    for(int e=edge[x].size()-1;e>=0;e--){        int v=edge[x][e];        if(dfn[v])continue;        dfs(v,x);    }}inline void getanc(int p){    if(anc[p]!=p)getanc(anc[p]);    if(dfn[sdom[mnsd[anc[p]]]]<dfn[sdom[mnsd[p]]])mnsd[p]=mnsd[anc[p]];    anc[p]=anc[anc[p]];}inline int eval(int p){  //获取到根节点路径上sdom最小的点    getanc(p);    return mnsd[p];}int main(){    n=rd();rt=n+1;    for(int i=1;i<=n;i++){        int cnt=0,x=rd();        while(x){            pre[i].push_back(x);            edge[x].push_back(i);            ++cnt;x=rd();        }        if(!cnt){            pre[i].push_back(rt);            edge[rt].push_back(i);        }    }    dfs(rt);fa[rt]=rt;    for(int i=tot;i>=1;i--){        int u=id[i];        for(int j=pre[u].size()-1;j>=0;j--){            int t=eval(pre[u][j]);            if(dfn[sdom[t]]<dfn[sdom[u]])sdom[u]=sdom[t];        }        anc[u]=fa[u];int f=fa[u];        buc[sdom[u]].push_back(u);        for(int j=buc[f].size()-1;j>=0;j--){ //此时该节点所有子树以父节点为半支配点的点路径全部已经处理出            int v=buc[f][j];int t=eval(v);            if(dfn[sdom[t]]<dfn[sdom[v]])idom[v]=t;            else idom[v]=sdom[v];        }        buc[f].clear();    }    for(int i=2;i<=tot;i++){        int u=id[i];        if(idom[u]!=sdom[u])idom[u]=idom[idom[u]];        g.add(idom[u],u);    }    g.dfs(id[1]);    for(int i=1;i<=n;i++){        W(g.sze[i]-1);putchar('\n');    }}