【拓扑排序+LCA构造Dominator Tree】BZOJ2815 [ZJOI2012]灾难

题面在这里

对每一个生产者建立一个共同的食物（太阳233）

可以对整个图拓扑排序

然后按倒拓扑序构造Dominator Tree

显然，对于一个点，它的所有食物的LCA就是其在Dominator Tree上的father

然后就好了

示例程序：

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int red(){
int res=0,f=1;char ch=nc();
while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=nc();}
while ('0'<=ch&&ch<='9') res=res*10+ch-48,ch=nc();
return res*f;
}

const int maxn=65500,maxe=1000005;
int n,g[maxn],f[maxn][17],dep[maxn],w[maxn];
int tot,lnk[maxn],nxt[maxe],son[maxe];
inline void add(int x,int y){ g[y]++;
son[++tot]=y;nxt[tot]=lnk[x];lnk[x]=tot;
}
int que[maxn];
void topo(){
int hed=0,til=0;
for (int i=1;i<=n;i++)
if (!g[i]) que[++til]=i;
while (hed!=til)
for (int j=lnk[que[++hed]];j;j=nxt[j]){
g[son[j]]--;
if (!g[son[j]]) que[++til]=son[j];
}
}
int LCA(int a,int b){
if (dep[a]<dep[b]) swap(a,b);
for (int j=log2(n+1);j>=0;j--)
if (dep[f[a][j]]>=dep[b]) a=f[a][j];
if (a==b) return a;
for (int j=log2(n+1);j>=0;j--)
if (f[a][j]!=f[b][j]) a=f[a][j],b=f[b][j];
return f[a][0];
}
int main(){
n=red();
for (int i=1;i<=n;i++)
for (int x=red();x;x=red()) add(i,x);
topo();
for (int i=0;(1<<i)<=n+1;i++) f[n+1][i]=n+1;dep[n+1]=1;
for (int i=n;i>=1;i--){
int x=-1;
if (!lnk[que[i]]) x=n+1;
for (int j=lnk[que[i]];j;j=nxt[j])
if (x<0) x=son[j];else x=LCA(x,son[j]);
f[que[i]][0]=x;dep[que[i]]=dep[x]+1;
for (int j=1;(1<<j)<=n+1;j++)
f[que[i]][j]=f[f[que[i]][j-1]][j-1];
}
for (int i=1;i<=n;i++) w[f[que[i]][0]]+=w[que[i]]+1;
for (int i=1;i<=n;i++) printf("%d\n",w[i]);
return 0;
}