bzoj 2815 灭绝树

对于一个食物网（一个DAG），一个物种死亡后，某些物种就必然死亡，求出必然死亡的是那些物种。

灭绝树的另一种含义是：“灭绝树跟节点到节点u的路径上的节点由那些原图中从根节点到节点u的所有路径中都经过了的点“。

 /**************************************************************
     Problem: 2815
     User: idy002
     Language: C++
     Result: Accepted
     Time:424 ms
     Memory:14036 kb
 ****************************************************************/

 #include <cstdio>
 #include <vector>
 #define maxn 65544
 #define maxp 17
 using namespace std;

 int n;
 vector<int> g[maxn], vg[maxn], gg[maxn];
 int indgr[maxn], topo[maxn], topo_cnt;
 int anc[maxn][maxp+], depth[maxn];
 int siz[maxn];

 void make_topo() {
     vector<int> stk;
     for( int i=; i<=n; i++ )
         if( indgr[i]== ) stk.push_back(i);
     while( !stk.empty() ) {
         int u = stk.back();
         topo[++topo_cnt] = u;
         stk.pop_back();
         for( int t=; t<g[u].size(); t++ ) {
             int v = g[u][t];
             indgr[v]--;
             if( indgr[v]== ) stk.push_back(v);
         }
     }
 }

 int lca( int u, int v ) {
     if( depth[u]<depth[v] ) swap(u,v);
     int t = depth[u]-depth[v];
     for( int i=; t; t>>=,i++ )
         if( t& ) u=anc[u][i];
     if( u==v ) return u;
     for( int p=maxp; p>=&&anc[u][]!=anc[v][]; p-- )
         if( anc[u][p]!=anc[v][p] ) u=anc[u][p], v=anc[v][p];
     return anc[u][];
 }

 void dfs( int u ) {
     siz[u] = ;
     for( int t=; t<gg[u].size(); t++ ) {
         int v = gg[u][t];
         dfs(v);
         siz[u] += siz[v];
     }
 }

 void work() {
     make_topo();
     for( int i=; i<=n; i++ ) {
         anc[i][] = i;
         depth[i] = ;
     }
     for( int i=; i<=n; i++ ) {
         int u = topo[i];
         if( vg[u].size()== ) continue;
         int ca = vg[u][];
         for( int t=; t<vg[u].size(); t++ ) {
             int v = vg[u][t];
             ca = lca(ca,v);
         }
         gg[ca].push_back(u);
         anc[u][] = ca;
         depth[u] = depth[ca]+;
         for( int p=; p<=maxp; p++ )
             anc[u][p] = anc[anc[u][p-]][p-];
     }
     dfs(n);
     for( int i=; i<=n-; i++ )
         printf( "%d\n", siz[i]- );
 }

 int main() {
     scanf( "%d", &n );
     for( int v=,u; v<=n; v++ ) {
         while() {
             scanf( "%d", &u );
             if( u== ) break;
             g[u].push_back(v);
             vg[v].push_back(u);
             indgr[v]++;
         }
     }
     for( int u=; u<=n; u++ )
         if( indgr[u]== ) {
             g[n+].push_back(u);
             vg[u].push_back(n+);
             indgr[u]++;
         }
     n++;
     work();
 }