bzoj 2815

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

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

【输入格式】 
输入文件 catas.in 的第一行是一个正整数 N，表示生物的种数。生物从 1 标
号到 N。 
接下来 N 行，每行描述了一个生物可以吃的其他生物的列表，格式为用空
格隔开的若干个数字，每个数字表示一种生物的标号，最后一个数字是 0 表示列
表的结束。 
【输出格式】 
输出文件 catas.out 包含N行，每行一个整数，表示每个生物的灾难值。 
【样例输入】 
5 
0 
1 0 
1 0 
2 3 0 
2 0 
【样例输出】 
4 
1 
0 
0 
0 
【样例说明】 
样例输入描述了题目描述中举的例子。 
【数据规模】 
对 50%的数据，N ≤ 10000。 
对 100%的数据，1 ≤ N ≤ 65534。 
输入文件的大小不超过 1M。保证输入的食物网没有环。 

 这是题目（破bzoj居然没有），白爷爷讲过，然而写起来gg；

 作者wxh

#include<bits/stdc++.h>

using namespace std;

const int maxn=66666;

int n,cnt,dep[maxn],q[maxn],h,tail,f[maxn][20],ans[maxn],head[maxn],in[maxn];

struct edge
{
	int to,nxt;
}e[maxn<<1];

vector <int> vec[maxn];

void addedge(int x,int y)
{
	e[++cnt].to=y;
	e[cnt].nxt=head[x];
	head[x]=cnt;
}

void topu()
{
	for(int i=1;i<=n;i++)
		if(!in[i])
			q[++tail]=i;
	while(h^tail)
	{
		int now=q[++h];
		for(int i=0;i<vec[now].size();i++)
		{
			in[vec[now][i]]--;
			if(!in[vec[now][i]])
				q[++tail]=vec[now][i];
		}
	}
}

void pushup(int &x,int d)
{
	for(int i=16;i>=0;i--)
		if(dep[x]-d>=(1<<i))
			x=f[x][i];
}

int lca(int x,int y)
{
	if(x==-1) return y;
	if(dep[x]<dep[y]) swap(x,y);
	if(dep[x]!=dep[y])
		pushup(x,dep[y]);
	if(x==y) return x;
	for(int i=16;i>=0;i--)
	{
		if(f[x][i]!=f[y][i])
		{
			x=f[x][i];
			y=f[y][i];
		}
	}
	return f[x][0];
}

void push(int x)
{
	for(int i=0;i<16;i++)
	{
		if(!f[x][i]) break;
		f[x][i+1]=f[f[x][i]][i];
	}
}

void build()
{
	for(int i=n;i;i--)
	{
		int now=q[i],fa=-1;
		for(int j=0;j<vec[now].size();j++)
		{
			fa=lca(fa,vec[now][j]);
		}
		if(fa==-1) fa=0;
		addedge(fa,now);
		dep[now]=dep[fa]+1;
		f[now][0]=fa;
		push(now);
	}
}

void dfs(int x)
{
	int tmp=1;
	for(int i=head[x];i;i=e[i].nxt)
	{
		dfs(e[i].to);
		tmp+=ans[e[i].to];
	}
	ans[x]=tmp;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		int x;
		while(scanf("%d",&x)&&x)
		{
			in[x]++;
			vec[i].push_back(x);
		}
	}
	topu();
	build();
	dfs(0);
	for(int i=1;i<=n;i++)
		printf("%d\n",ans[i]-1);
	return 0;
}