hdu6031 Innumerable Ancestors

题目

答案明显满足二分性质，所以我们可以先二分出答案的深度，然后判断的话，把A集合用倍增跳到那个深度，再对B集合做同样的操作，判断两个祖先集合的交是否为空集，只要用map就可以呢。

还有另外一个解法，就是深度最大的LCA在DFS序中与它们两点距离最近，然后贪一贪就好了。

#include
#define N 500000
using namespace std;
int n,q,m,x,y,la,lb,a[N+5],b[N+5];
int dep[N+5],f[N+5][25];
int first[N+5],nxt[N+5],to[N+5],siz;
int l,r,mid,ans;
inline void link(int x,int y)
{
    nxt[siz]=first[x];
    first[x]=siz;
    to[siz++]=y;
}
inline void reset()
{
    memset(first,-1,sizeof(first));
    memset(f,0,sizeof(f));
    siz=0;
}
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 read()
{
    int x=0,b=1;
    char c=nc();
    for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
    for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
    return x*b;
}
inline void write(int x)
{
    if(x==0)putchar('0');
    else
    {
        char buf[15];
        int len=0;
        if(x<0)putchar('-'),x=-x;
        while(x)buf[++len]=x%10+'0',x/=10;
        for(int i=len;i>=1;i--)putchar(buf[i]);
    }
    putchar(' ');
}
inline void pre(int x,int fa)
{
    dep[x]=dep[fa]+1;
    for(int i=1;i<=20;i++)if(f[x][i-1])f[x][i]=f[f[x][i-1]][i-1];else break;
    for(int i=first[x];i!=-1;i=nxt[i])
    {
        int u=to[i];
        if(u==fa)continue;
        f[u][0]=x;
        pre(u,x);
    }
}
inline int getf(int x,int d)
{
    if(d<0)return -1;
    if(d==0)return x;
    for(int i=20;i>=0;i--)
        if(d>=(1<1<return x;
}
inline bool check(int d)
{
    set <int> E;
    for(int i=1;i<=la;i++)
    {
        int dis=dep[a[i]]-d;
        int ret=getf(a[i],dis);
        if(ret==-1)continue;
        E.insert(ret);
    }
    for(int i=1;i<=lb;i++)
    {
        int dis=dep[b[i]]-d;
        int ret=getf(b[i],dis);
        if(ret==-1)continue;
        if(E.count(ret))return true;
    }
    return false;
}
int main()
{
    freopen("in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        reset();
        for(int i=1;iscanf("%d%d",&x,&y),link(x,y),link(y,x);
        pre(1,0);
        for(int i=1;i<=m;i++)
        {
            scanf("%d",&la);for(int j=1;j<=la;j++)scanf("%d",&a[j]);
            scanf("%d",&lb);for(int j=1;j<=lb;j++)scanf("%d",&b[j]);
            l=0,r=n;
            while(l<=r)
            {
                mid=(l+r)/2;
                if(check(mid))ans=mid,l=mid+1;
                else r=mid-1;
            }
            printf("%d\n",ans);
        }
    }
    return 0;
}