【题目描述】
阿米巴是小强的好朋友。
阿米巴和小强在草原上捉蚂蚱。小强突然想,如果蚂蚱被他们捉灭绝了,那 么吃蚂蚱的小鸟就会饿死,而捕食小鸟的猛禽也会跟着灭绝,从而引发一系列的生态灾难。
学过生物的阿米巴告诉小强,草原是一个极其稳定的生态系统。如果蚂蚱灭绝了,小鸟照样可以吃别的虫子,所以一个物种的灭绝并不一定会引发重大的灾难。
我们现在从专业一点的角度来看这个问题。我们用一种叫做食物网的有向图 来描述生物之间的关系:一个食物网有 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~~~惊呆了开心嘿嘿