RAN-Rendezvous【POI2012】

RAN-Rendezvous【POI2012】

题目描述

给定一个 n个顶点的有向图，每个顶点有且仅有一条出边。

对于顶点 i，记它的出边为 (i,a[i])。

再给出q组询问，每组询问由两个顶点 a、b组成，要求输出满足下面条件的 x、y：

1.从顶点 a沿着出边走x步和从顶点b沿着出边走y步后到达的顶点相同。

2.在满足条件 1的情况下，如果解不唯一，则还需要令 max(x,y)最小。

3.在满足条件 1和 2的情况下，如果解不唯一，则还需要令 min(x,y)最小。

4.在满足条件 1、2和 3的情况下，如果解不唯一，则还需要令 x≥y。

如果不存在满足条件1的 x、y，输出 -1 -1。

输入格式

第一行两个正整数n和q。

第二行n个正整数 a[1],a[2],…,a[n]。

下面 q行，每行两个正整数 a,b，表示一组询问。

输出格式

输出q行，每行两个整数，表示每次询问的x和y

输入输出样例

输入样例1

12 5 4 3 5 5 1 1 12 12 9 9 7 1 7 2 8 11 1 2 9 10 10 5

输出样例1

2 3 1 2 2 2 0 1 -1 -1

说明

【数据范围】

2 <= n, q <= 500000; 1 <= a, b, a[i] <= n。 

【耗时限制】2000ms 【内存限制】256MB

---

以上是题目描述

讲解分析

这道题思路

是真的简单！！！但是代码是真的长......

上交楼主的思路:

第一句话就明显地点明了这是一道基环树的题。并且还是一颗内向基环树。

我们考虑每对x,y的情况：

1.x,y不在一棵基环树上。那么输出-1即可。

2.x,y在同一棵子树上。那么答案就是树上距离，用lca求解即可。

3.x,y在同一棵基环树上，但不在同一棵子树上。由于是有向图，我们可以枚举x走过环到y的距离和y走过环到x的距离，根据题目输出答案即可。

具体流程：

1.首先预处理出环。

2.预处理出树上倍增的fa数组和每个点所在的子树、基环树编号，并且准备好要使用的距离的信息。

对于每个询问，在线回答即可。

时间复杂度为O(N+Q) ~ O((N+Q)logN)

下图是别人帮楼主整理的思路，可能更好懂一些

这是楼主的代码，与思路一同上交

码风清奇，各位见谅

在拿代码的同时，别忘了点三连（关注，点赞，分享）

#include <iostream>
#include <cstdio>
#include <cstring> 
#include <algorithm>
#include <sstream>
#include <cmath>
#include <vector>
#include <queue>
using namespace std;
const int N=5e5+5,M=log2(N)+2;
vector<int> G[N];
int n,m,d[N],p[N][M],cc[N];
int log2n=log2(N)+0.5,root[N];
int num,vis[N],len[N],id[N];
int fa[N],q;
void bfs(int r){
	d[r]=1;
	root[r]=r;
	vis[r]=num;
	queue<int>q;q.push(r);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int i=0;i<G[u].size();i++){
			int v=G[u][i];
			if(!d[v]&&!id[v]){
			    root[v]=r,vis[v]=num;
				d[v]=d[u]+1;
				p[v][0]=u;
				for(int j=1;j<=log2n;j++) p[v][j]=p[p[v][j-1]][j-1];
				q.push(v);
			}
		}
	}
}
int LCA(int x,int y){
	if(d[x]<d[y]) swap(x,y);
	for(int j=log2n;j>=0;j--){
		if(d[p[x][j]]>=d[y]) x=p[x][j];
	}
	if(x==y) return x;
	for(int j=log2n;j>=0;j--){
		if(p[x][j]!=p[y][j]) 
		    x=p[x][j],y=p[y][j];
	}
	return p[x][0];
}
bool better(int x1,int y1,int x2,int y2){
    int max1=max(x1,y1),max2=max(x2,y2),min1=min(x1,y1),min2=min(x2,y2);
    if(max1<max2) return true;
    if(max1==max2&&min1<min2) return true;
    if(max1==max2&&min1==min2&&x1>=y1) return true;
    return false;
}
int main(){
	scanf("%d%d",&n,&q);
	int v;
	for(int i=1;i<=n;i++){
	    scanf("%d",&v);
	    fa[i]=v;
	    G[v].push_back(i);
	}
	for(int i=1;i<=n;i++){
	    if(vis[i]) continue;
	    ++num;
	    int r=i;
	    while(!vis[r]) vis[r]=num,r=fa[r];
	    int cnt=0;
	    while(!id[r]) id[r]=++cnt,cc[cnt]=r,r=fa[r];
	    len[num]=cnt;
	    for(int i=1;i<=cnt;i++) bfs(cc[i]);
	}
	int a,b,x,y;
	for(int i=1;i<=q;i++){
	    scanf("%d%d",&a,&b);
	    if(vis[a]!=vis[b]) printf("-1 -1\n");
	    else if(root[a]==root[b]){
	        int lca=LCA(a,b);
	        printf("%d %d\n",d[a]-d[lca],d[b]-d[lca]);
	    }
	    else{
	        int cnt=len[vis[a]];
	        int ra=root[a],rb=root[b];
	        int x1=d[a]-1,y1=d[b]-1+(id[ra]-id[rb]+cnt)%cnt;
	        int x2=d[a]-1+(id[rb]-id[ra]+cnt)%cnt,y2=d[b]-1;
	        if(better(x1,y1,x2,y2)) printf("%d %d\n",x1,y1);
	        else printf("%d %d\n",x2,y2);
	    }
	}
    return 0;
}