[BZOJ2815][ZJOI2012]灾难 拓扑排序+lca

若A能吃掉B，那么B向A连边。建一个虚根连向所有入度为0的点，拓扑排序一发。
考虑这么一棵树，假如某个节点灭绝，那么它的子树也全部灭绝。构造出这样一棵树就好了。
假设拓扑中前k-1项的树已经构造好了，考虑第k项该接在哪。显然，应该接在原图中所有向它有连边的点在树中的lca下。
构造出来统计一发子树大小就行了。
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=100010;
int n,rd[maxn],a[maxn],st[maxn],fa[maxn][20],dep[maxn],size[maxn];
struct edge
{
    int t;
    bool rev;
    edge *next; 
}*con[maxn];
int read()
{
    int x=0;
    char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x;
}
void ins(int x,int y)
{
    rd[y]++;
    edge *p=new edge;
    p->t=y;
    p->rev=0;
    p->next=con[x];
    con[x]=p;
    p=new edge;
    p->t=x;
    p->rev=1;
    p->next=con[y];
    con[y]=p;
}
int getlca(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    for(int k=19;k>=0;k--)
        if(dep[fa[x][k]]>=dep[y]) x=fa[x][k];

    if(x==y) return x;
    for(int k=19;k>=0;k--)
        if(fa[x][k]!=fa[y][k]) {x=fa[x][k];y=fa[y][k];}
    return fa[x][0];        
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        int x=read();
        while(x!=0) {ins(x,i);x=read();}
    }
    for(int i=1;i<=n;i++)
        if(rd[i]==0) ins(0,i);
    int top=1;st[1]=a[0]=0;
    for(int i=0;i<=n;i++)
    {
        a[i]=st[top--];
        for(edge *p=con[a[i]];p!=NULL;p=p->next)
        if(!p->rev)
        {
            rd[p->t]--;
            if(rd[p->t]==0) st[++top]=p->t;       
        }
    }
    memset(dep,0,sizeof(dep));
    memset(fa,0,sizeof(fa));
    for(int i=1;i<=n;i++)
    {
        int lca=-1;
        for(edge *p=con[a[i]];p!=NULL;p=p->next)
            if(p->rev)
            {
                if(lca==-1) lca=p->t;
                else lca=getlca(p->t,lca);
            }   
        fa[a[i]][0]=lca;
        dep[a[i]]=dep[lca]+1;
        for(int k=1;(1<<k)<=n;k++)
            fa[a[i]][k]=fa[fa[a[i]][k-1]][k-1];

    }
    memset(size,0,sizeof(size));
    for(int i=n;i>=1;i--)
    {
        size[a[i]]++;
        size[fa[a[i]][0]]+=size[a[i]];
    }
    for(int i=1;i<=n;i++)
        printf("%d\n",size[i]-1);
    return 0;
}