BZOJ 2815 [ZJOI2012]灾难 （支配树）拓扑+倍增

【问题描述】

阿米巴是小强的好朋友。

阿米巴和小强在草原上捉蚂蚱。小强突然想，如果蚂蚱被他们捉灭绝了，那么吃蚂蚱的小鸟就会饿死，而捕食小鸟的猛禽也会跟着灭绝，从而引发一系列的生态灾难。

学过生物的阿米巴告诉小强，草原是一个极其稳定的生态系统。如果蚂蚱灭绝了，小鸟照样可以吃别的虫子，所以一个物种的灭绝并不一定会引发重大的灾难。

我们现在从专业一点的角度来看这个问题。我们用一种叫做食物网的有向图来描述生物之间的关系：

一个食物网有 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。保证输入的食物网没有环。

传送门

。。查了一下午的错啊woc....

现在才知道这东西叫做支配树。。。

首先理解一下题目之后发现，删除一个点的前提是删除它的所有儿子。

从入度为零的点开始，沿着拓扑序进行，这个似乎很显然。。

问题就是接下来怎么办= =

想了好久，后面发现可以根据原来能够直接删除的点构成的一棵树来做。

就是说，假如一个点v，一开始被u直接到达，且v的入度为1，那么v就是树中u的儿子。

对于在树上的点，答案就是它的子树size-1，

对于不在树上的点，求出它所有的直接来源点（就是反向边的点）LCA总的LCA，然后这个LCA的答案+1.

由于可能有很多入度为零的点，所以建立一个超级源，连接这些点，

然后再按拓扑序搞一下就好了= =

问题是……实现麻烦。。就比如说，求出所有儿子LCA总的LCA，这个过程就有很多特判。。

于是默默调试一个下午。。

就要弃疗的时候，终于看到一线生机。。。。

“支配树：支配点之间连边”

支配点就是说，u~v的所有路径都经过了一个点x，x就是v的支配点，

……所以对于不在树上的点u，求出那个LCA之后，那么LCA便是u的支配点了。。

也就是说u和LCA连边。

然后发现是加点加边、在线LCA，那就倍增，每次加边都倍增一次即可。

建完整棵树后，答案就是每个点的size了（一个点u的支配点是v，u对v有1的贡献）

实现的时候其实还可以归约一下，直接在拓扑排序的过程中实现就好了。

当然为了得出一个点的来源点，还需要记录反向边。

一开始打得很挫（实在太挫了不放上来了噗）= =忍不住了看了看别人的代码。。。

dalao写得好哇

于是用他的方法重写了一遍，时间快了1倍……= =

#include<bits/stdc++.h>
using namespace std;
const int 
	N=65539,
	M=100005;
int n,Eoppcnt,Ecnt[2],sz[N],dep[N];
int into[N],fa[N][17];
struct Edge{
	int next,to;
}E[2][M],Eopp[M];int head[2][N],heado[N];
void add(int o,int u,int v){
	E[o][++Ecnt[o]].next=head[o][u];
	E[o][Ecnt[o]].to=v;
	head[o][u]=Ecnt[o];
	if (o){
		fa[v][0]=u,dep[v]=dep[u]+1;
		for (int i=1;i<=16;i++)
			fa[v][i]=fa[fa[v][i-1]][i-1];
	} else{
		Eopp[++Eoppcnt].next=heado[v];
		Eopp[Eoppcnt].to=u;
		heado[v]=Eoppcnt;
	}
}
int get_LCA(int x,int y){
	if (!x) return y;
	if (dep[x]<dep[y]) swap(x,y);
	int t=dep[x]-dep[y];
	for (int i=0;i<=16;i++)
		if (t&(1<<i)) x=fa[x][i];
	for (int i=16;~i;i--)
		if (fa[x][i]!=fa[y][i])
			x=fa[x][i],y=fa[y][i];
	if (x==y) return x;
		else return fa[x][0];
}
void topo_solve(){
	queue<int>Q;
	while (!Q.empty()) Q.pop();
    dep[++n]=Ecnt[1]=0;
    for (int i=1;i<n;i++)
        if (!into[i]) add(1,n,i),Q.push(i);
    while(!Q.empty()){
        int u=Q.front();Q.pop();
        for (int i=head[0][u];i;i=E[0][i].next){
            int v=E[0][i].to;
            into[v]--;
            if (!into[v]){
                Q.push(v);int x=0;
                for (int j=heado[v];j;j=Eopp[j].next) x=get_LCA(x,Eopp[j].to);
                add(1,x,v);
            }
        }
    }
}
void dfs(int u,int pre){
	sz[u]=1;
	for (int i=head[1][u];i;i=E[1][i].next){
		int v=E[1][i].to;
		if (v==pre) continue;
		dfs(v,u),sz[u]+=sz[v];
	}
}
int main(){
	scanf("%d",&n);
	Ecnt[0]=Eoppcnt=0;int x;
	for (int i=1;i<=n;i++)
		while (scanf("%d",&x),x) add(0,x,i),into[i]++;
	topo_solve(),dfs(n,-1);
	for (int i=1;i<n;i++) printf("%d\n",sz[i]-1);
	return 0;
}