【ZJOI2012】【BZOJ2815】灾难 (catas) {拓扑+倍增lca}

• 【题目描述】
  阿米巴是小强的好朋友。
  阿米巴和小强在草原上捉蚂蚱。小强突然想，如果蚂蚱被他们捉灭绝了，那 么吃蚂蚱的小鸟就会饿死，而捕食小鸟的猛禽也会跟着灭绝，从而引发一系列的生态灾难。
  学过生物的阿米巴告诉小强，草原是一个极其稳定的生态系统。如果蚂蚱灭绝了，小鸟照样可以吃别的虫子，所以一个物种的灭绝并不一定会引发重大的灾难。
  我们现在从专业一点的角度来看这个问题。我们用一种叫做食物网的有向图 来描述生物之间的关系：一个食物网有 N个点，代表 N 种生物，如果生物 x 可以吃生物 y，那么从 y 向 x 连一个有向边。这个图没有环。
  图中有一些点没有连出边，这些点代表的生物都是生产者，可以通过光合作用来生存；而有连出边的点代表的都是消费者，它们必须通过吃其他生物来生存。
  如果某个消费者的所有食物都灭绝了，它会跟着灭绝。
  我们定义一个生物在食物网中的“灾难值”为，如果它突然灭绝，那么会跟 着一起灭绝的生物的种数。

  

  如果小强和阿米巴把草原上所有的羊都给吓死了，那么狼会因为没有食物而灭绝，而小强和阿米巴可以通过吃牛、牛可以通过吃草来生存下去。所以，羊的灾难值是1。但是，如果草突然灭绝，那么整个草原上的5种生物都无法幸免，所以，草的灾难值是4。
  给定一个食物网，你要求出每个生物的灾难值。

• 【Sample Input】
  （输入文件 catas.in 的第一行是一个正整数 N，表示生物的种数。生物从1标号到N。
  接下来 N 行，每行描述了一个生物可以吃的其他生物的列表，格式为用空格隔开的若干个数字，每个数字表示一种生物的标号，最后一个数字是 0 表示列 表的结束。）
  5
  0
  1 0
  1 0
  2 3 0
  2 0

• 【Sample Output】
  （输出文件 catas.out 包含N行，每行一个整数，表示每个生物的灾难值。）
  4
  1
  0
  0
  0

• 【数据规模】
  对 50%的数据，N ≤ 10000。
  对 100%的数据，1 ≤ N ≤ 65534。
  输入文件的大小不超过 1M。保证输入的食物网没有环。

---

【题解】拓扑+倍增
{50%}反向建边拓扑排序后暴力做就可以了；
{100%}首先依然是拓扑排序。接下来思考的问题就是：什么条件下一种生物会灭绝？
<1>生物本身灭绝了；
<2>此生物的所有食物都灭绝了。
考虑情况<2>，要使得此生物的所有食物灭绝，这些食物的公共祖先必然灭绝。换句话说，如果某种生物的所有食物的公共祖先灭绝，此生物就会灭绝。根据这个性质，找最近公共祖先，建立新图。这里用到树上倍增的算法。
（整理）对于生物 x，找其所有食物的最近公共组先lca，从lca向x连一条有向边。建图完毕后dfs遍历新图即可。

---

#include <cstdio>
#include <iostream>
#define N 65540
struct edge{ int to,nxt;}e[10*N],e1[10*N],e2[10*N]; 
int n,cnt,cnt1,cnt2,d[N],de[N],f[N][20],
    ans[N],fi[N],fi1[N],fi2[N],q[N]; 
    void add(int u,int v)
    {
        e[++cnt].to=v;e[cnt].nxt=fi[u];fi[u]=cnt; 
    }
    void add1(int u,int v)
    {
        e1[++cnt1].to=v;e1[cnt1].nxt=fi1[u];fi1[u]=cnt1; 
    }
    void add2(int u,int v)
    {
        e2[++cnt2].to=v;e2[cnt2].nxt=fi2[u];fi2[u]=cnt2; 
    }
    void topology()
    {
        int h=1,t=0;
        for (int i=1;i<=n;++i)
            if (!d[i]) q[++t]=i;
        for (;h<=t;++h)
            for (int i=fi1[q[h]];i;i=e1[i].nxt)
                if (!--d[e1[i].to]) q[++t]=e1[i].to;
    }
    int findlca(int u,int v)
    {
        if (de[u]<de[v]) std::swap(u,v);
        for (int i;de[u]>de[v];u=f[u][i-1]) 
            for (i=1;de[f[u][i]]>de[v];++i);
        for (int i;u!=v;u=f[u][i-1],v=f[v][i-1])
            for (i=1;f[u][i]!=f[v][i];++i);
        return u;
    }
    void dfs(int x)
    {
        ans[x]=1;
        for (int i=fi2[x];i;i=e2[i].nxt)
        {
            dfs(e2[i].to);
            ans[x]+=ans[e2[i].to];
        }
    } 
int main()
{
    scanf("%d\n",&n);
    for (int i=1;i<=n;++i)
    {
        int x;
        for (scanf("%d",&x);x;scanf("%d",&x))
        {
            add(i,x);add1(x,i);++d[i];
        }
        scanf("\n");
    }
    topology();
    for (int i=1;i<=n;++i)
    {
        int x=e[fi[q[i]]].to;
        for (int j=fi[q[i]];j;j=e[j].nxt)
            x=findlca(x,e[j].to);
        add2(x,q[i]);de[q[i]]=de[x]+1;
        f[q[i]][0]=x;
        for (int k=x,l=0;f[k][l];k=f[k][l++])
            f[q[i]][l+1]=f[k][l];
    }
    dfs(0);
    for (int i=1;i<=n;++i) printf("%d\n",--ans[i]);
    return 0; 
}

【题外话】O(∩_∩)O这题竟然能1A~~~惊呆了开心嘿嘿