HDU6604(Blow up the city DAG支配树模板提)

题目
题意：给你一个DAG，定义中心城市为没有出度的点。一个DAG里会存在多个中心城市。q次询问x,y。问炸掉一个城市使得a,b至少有一个不能到达中心城市得方案数。
T<=10 n<=1e5 m<=2e5 q<=1e5 Time Limit: 5000 MS Memory Limit: 131072 K
思路：将所有的中心城市用一个超级终点root连接起来。就转化为在DAG的支配树上求。
就是a到搭终点的必经点+b到达终点的必经点-a,b共同的必经点。
dep[x]+dep[y]-dep[LCA(x,y)]-1(-1的意思是root这个点不算到答案中去。。这事自己加的)
关键就是求出DAG的支配树。DAG的支配树可以用topo排序和倍增LCA建树。理解一下。

一个点的最近必经点是他所有儿子的LCA.

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
#define en '\n
using namespace std;
typedef long long ll;
const int N=1e5+10,M=2e5+10;
int n,m;
struct Edge{int to,nex;};
struct Graph{
    Edge e[M];int head[N],tot;
    void init(){
        tot=0;
        for(int i=1;i<=n+1;++i) head[i]=0;
    }
    void add(int from,int to){
        e[++tot]=(Edge){to,head[from]};head[from]=tot;
    }
}A,B;
int dep[N],fa[N][25],out[N],a[N],root,K=20;
queue<int>q;
void TOPO(){
    root=n+1;
    for(int i=1;i<=n;++i)
        if(!out[i]) q.push(i),A.add(root,i),B.add(i,root);
    int num=0;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        a[++num]=x;
        for(int i=A.head[x];i;i=A.e[i].nex){
            int y=A.e[i].to;
            if(!(--out[y])) q.push(y);
        }
    }
}
int LCA(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    for(int j=K;j>=0;j--)
        if(dep[x]-(1<<j)>=dep[y]) x=fa[x][j];
    if(x==y) return x;
    for(int j=K;j>=0;j--)
        if(fa[x][j]!=fa[y][j]) x=fa[x][j],y=fa[y][j];
    return fa[x][0];
}
int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        A.init(),B.init();
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;++i) out[i]=0;
        while(m--){
            int u,v;scanf("%d%d",&u,&v);
            A.add(v,u),B.add(u,v);
            ++out[u];
        }
        TOPO();dep[root]=1;
        for(int i=1;i<=n;++i){
            int x=a[i],bf=-1;
            for(int j=B.head[x];j;j=B.e[j].nex){
                int y=B.e[j].to;
                bf=(bf==-1)?y:LCA(bf,y);
            }
            dep[x]=dep[bf]+1,fa[x][0]=bf;
            for(int i=1;i<=K;++i) fa[x][i]=fa[fa[x][i-1]][i-1];
        }
        int q;scanf("%d",&q);
        while(q--){
            int x,y;scanf("%d%d",&x,&y);
            printf("%d\n",dep[x]+dep[y]-dep[LCA(x,y)]-1);
        }
    }
}