[BZOJ2815][ZJOI2012]灾难（倍增lca+top）

题目：我是超链接

题解：

啊其实可以1A的......但因为我写top的时候手滑.......

依次考虑每个生物i：假设我们已经构建好了排序在i之前的生物组成的“灭绝树”了。 

假设i的食物有x[0]~x[k]（x[0]~x[k]在拓扑排序中比i靠前）， 

很显然，只有x[0]~x[k]在树上的公共祖先的灭绝会导致i灭绝，否则i一定可以找到能让它活下来的食物。 

于是，我们可以把i挂在x[0]~x[k]的最近公共祖先下面。 

处理完所有的生物，我们得到的树就是整个图的灭绝树了。之后dfs计算一下每个节点对应的子树大小，子树大小-1即危险程度

建一棵灭绝树是为了更快查询lca+统计答案

代码：

#include <cstdio>
#include <queue>
#include <cstring>
#include <iostream>
#define sz 31
#define N 65540
using namespace std;
int tot,nxt[N*4],point[N*4],v[N*4],tot1,nxt1[N*4],point1[N*4],v1[N*4],tot2,nxt2[N*4],point2[N*4],v2[N*4],n,in
,ran
,h
;
int f
[sz],ans
;
void addline(int x,int y){++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;}
void addline1(int x,int y){++tot1; nxt1[tot1]=point1[x]; point1[x]=tot1; v1[tot1]=y;}
void addline2(int x,int y){++tot2; nxt2[tot2]=point2[x]; point2[x]=tot2; v2[tot2]=y;}
void topsort()
{
queue<int>q;
for (int i=1;i<=n;i++)
if (!in[i]) q.push(i);
int num=0;
while (!q.empty())
{
int now=q.front(); q.pop();
ran[++num]=now;
for (int i=point[now];i;i=nxt[i])
{
in[v[i]]--;
if (!in[v[i]]) q.push(v[i]);
}
}
}
int lca(int x,int y)
{
if (h[x]<h[y]) swap(x,y);
int k=h[x]-h[y];
for (int i=0;i<sz;i++)
if ((k>>i)&1) x=f[x][i];
if (x==y) return x;
for (int i=sz-1;i>=0;i--)
if (f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
return f[x][0];
}
void dfs(int now)
{
for (int i=point2[now];i;i=nxt2[i])
{
dfs(v2[i]);
ans[now]+=ans[v2[i]];
}
ans[now]++;
}
int main()
{
int i,j;
scanf("%d",&n);
for (i=1;i<=n;i++)
{
int x;
scanf("%d",&x);
while (x) {in[i]++;addline1(i,x);addline(x,i);scanf("%d",&x);}
}
topsort();
for (i=1;i<=n;i++)
{
int x=v1[point1[ran[i]]];
for (j=point1[ran[i]];j;j=nxt1[j])
x=lca(x,v1[j]);
addline2(x,ran[i]);
h[ran[i]]=h[x]+1;
f[ran[i]][0]=x;
for (j=1;j<sz;j++) f[ran[i]][j]=f[f[ran[i]][j-1]][j-1];
}
dfs(0);
for (i=1;i<=n;i++) printf("%d\n",ans[i]-1);
}