约会 Rendezvous:基环树

提炼：tarjan判环,dfs建树,倍增lca,预处理环两点间距离

我犯的错误：

1.基环树不只有一棵，可以有很多

2.自环不能将其忽略，（对于我的算法）应该将其特殊考虑在算法内

3.代码一定要简洁有力，不能让自己调都恶心

Code

数组含义：

que[N]队列,in_que[N]tarjan判是否入队,dfn[N],low[N],bel[N]点属于的集合

ver[N]根节点的集合编号,to[N],nxt[N],head[N],f[N][20]倍增lca,d[N]深度,bel_rt[N]点属于的树的根

extra[N],Ex[N]特判自环专用，vector<int>vec[N]tarjan点集合

#include<cstdio>
#include<vector>
const int L=<<|;
char buffer[L],*S,*T;
#define getchar() ((S==T&&(T=(S=buffer)+fread(buffer,1,L,stdin),S==T))?EOF:*S++)
using namespace std;
const int N=5e5+;
int n,k,num_tarjan,num_bian,num_huan,top,
que[N],in_que[N],dfn[N],low[N],bel[N],ver[N],to[N],nxt[N],head[N],f[N][],d[N],bel_rt[N],extra[N],Ex[N];
vector<int>vec[N];
int read(){
    int x=;char ch=getchar();
    while(ch<''||ch>'')ch=getchar();
    while(ch>=''&&ch<='')x=x*+ch-'',ch=getchar();
    return x;
}
void add(int x,int y){if(x==y)extra[++extra[]]=x;
    to[++num_bian]=y,nxt[num_bian]=head[x],head[x]=num_bian;
}
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a>b?b:a;}
void tarjan(int x){
    dfn[x]=low[x]=++num_tarjan;
    que[++top]=x;in_que[x]=;
    for(int i=head[x],y;i;i=nxt[i])
        if(!dfn[y=to[i]])tarjan(y),low[x]=min(low[x],low[y]);
        else if(in_que[y])low[x]=min(low[x],dfn[y]);
    if(dfn[x]==low[x]){
        num_huan++;
        int y;
        do{
            y=que[top--];
            in_que[y]=;
            vec[num_huan].push_back(y);
            bel[y]=num_huan;
        }while(y!=x);
        if((int)vec[num_huan].size()>)Ex[++Ex[]]=num_huan;
    }
}
void dfs(int x,int depth,int root){
    bel_rt[x]=root;
    d[x]=depth;
    for(int i=head[x],y;i;i=nxt[i]){
        if(bel[root]!=bel[y=to[i]]&&!d[y=to[i]]){
            f[y][]=x;
            for(int j=;j<=;++j)f[y][j]=f[f[y][j-]][j-];
            dfs(y,depth+,root);
        }
    }
}
int get(int x,int y){
    if(d[x]<d[y])swap(x,y);
    for(int i=;i>=;--i)if(d[f[x][i]]>=d[y])x=f[x][i];
    if(x==y)return y;
    for(int i=;i>=;--i)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
    return f[y][];
}
int main(){
    n=read(),k=read();
    for(int i=;i<=n;++i)add(read(),i);
    for(int i=;i<=n;++i)if(!dfn[i])tarjan(i);
    for(int hua=;hua<=Ex[];++hua){int huan=Ex[hua];
        for(int i=;i<(int)vec[huan].size();++i){
            dfs(vec[huan][i],,vec[huan][i]);
            ver[vec[huan][i]]=i+;
        }
    }
    for(int i=;i<=extra[];++i)
        if(!d[extra[i]]){
            dfs(extra[i],,extra[i]);
            ver[extra[i]]=;
        }
    for(int i=,x,y;i<=k;++i){
        x=read(),y=read();
        if(bel[bel_rt[x]]!=bel[bel_rt[y]])
            printf("-1 -1\n");
        else if(bel_rt[x]==bel_rt[y]){
            int lca=get(x,y);
            printf("%d %d\n",d[x]-d[lca],d[y]-d[lca]);
        }
        else{
            int ai=d[x]-,bi=d[y]-,huan_dis=(int)vec[bel[bel_rt[x]]].size();
            int len_b=ver[bel_rt[x]]-ver[bel_rt[y]],len_a;
            if(len_b<)len_a=-*len_b,len_b=huan_dis-len_a;
            else len_a=huan_dis-len_b;
            int Ma=max(ai+len_a,bi),Mb=max(ai,bi+len_b),ma=min(ai+len_a,bi),mb=min(ai,bi+len_b);
            //条件一
            if(Ma<Mb)
                printf("%d %d\n",ai+len_a,bi);
            else if(Ma>Mb)
                printf("%d %d\n",ai,bi+len_b);
            else{
                //条件二
                if(ma<mb)
                    printf("%d %d\n",ai+len_a,bi);
                else if(ma>mb)
                    printf("%d %d\n",ai,bi+len_b);
                else{
                    //条件三
                    if(ai+len_a>=bi)
                        printf("%d %d\n",ai+len_a,bi);
                    else
                        printf("-1 -1\n");
                }
            }
        }

        }
    return ;
}