【BZOJ】2815: [ZJOI2012]灾难

简要题意：

给一个有向无环图，问每个节点删掉之后会导致多少个点不可达。

似乎以前拿来考过....

我们定义一棵树，它满足对应点造成的灭绝值即为当点的子树大小-1

　　按照被捕食者--->捕食者的关系拓扑排序，然后依次建树，建到当前点的时候可以作为当前生物食物的点应当已经在树中了。如果当前点代表的生物要灭亡，很好理解那么可以作为它食物的生物都要灭亡，所以将这个点丢到它的可以作为它食物的生物的所有点的最近$LCA$点之下就可以了，可以用倍增$LCA$来维护动态加点。

　　为了方便，可以令一个点作为所有生产者的根。

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<vector>
 #include<cstdlib>
 #include<cmath>
 #include<cstring>
 using namespace std;
 #define maxn 100100
 #define llg int
 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
 llg n,m,e[maxn],du[maxn],dl[maxn],tpx[maxn],deep[maxn],f[maxn][],size[maxn];
 vector<llg>a[maxn],c[maxn],g[maxn];
 void in(llg x)
 {
     llg d;
     while ()
     {
         scanf("%d",&d);
         if (d==) break;
         du[x]++;
         a[d].push_back(x);
         c[x].push_back(d);
     }
 }

 void TP()
 {
     dl[]=n+;
     llg head=,tail=,x,v,w,cnt=;
     do
     {
         x=dl[++head];
         w=a[x].size();
         tpx[x]=++cnt;
         for (llg i=;i<w;i++)
         {
             v=a[x][i];
             du[v]--;
             if (du[v]==)
             {
                 dl[++tail]=v;
             }
         }
     }while (head!=tail);
 }

 bool cmp(llg a,llg b){return tpx[a]<tpx[b];}

 void make_f(llg x) {for (llg i=;i<=;i++) f[x][i]=f[f[x][i-]][i-];}

 llg find(llg x,llg y)
 {
     if (deep[x]<deep[y]) swap(x,y);
     for (llg i=;i>=;i--)
         if (deep[f[x][i]]>=deep[y])
             x=f[x][i];
     if (x==y) return x;
     for (llg i=;i>=;i--)
         if (f[x][i]!=f[y][i])
             x=f[x][i],y=f[y][i];
     return f[x][];
 }

 void work(llg x)
 {
     llg ro;
     if (c[x].size()==) ro=c[x][];
     else
     {
         llg w=c[x].size();
         ro=find(c[x][],c[x][]);
         for (llg i=;i<w;i++) ro=find(ro,c[x][i]);
     }
     f[x][]=ro; deep[x]=deep[ro]+;
     g[ro].push_back(x);
     make_f(x);
 }

 void dfs(llg x)
 {
     llg w=g[x].size();
     size[x]=;
     for (llg i=;i<w;i++)
     {
         dfs(g[x][i]);
         size[x]+=size[g[x][i]];
     }
 }

 int main()
 {
     yyj("catas");
     cin>>n;
     for (llg i=;i<=n;i++) in(i);
     for (llg i=;i<=n;i++)
         if (du[i]==)
         {
             du[i]++;
             c[i].push_back(n+);
             a[n+].push_back(i);
         }
     TP();
     for (llg i=;i<=n+;i++) e[i]=i;
     sort(e+,e+n+,cmp);
     deep[n+]=;
     for (llg i=;i<=n;i++)
         work(e[i]);
     dfs(n+);
     for (llg i=;i<=n;i++) printf("%d\n",size[i]-);
     return ;
 }

---