2791: [Poi2012]Rendezvous

Time Limit: 25 Sec  Memory Limit: 128 MB
Submit: 259  Solved: 160
[Submit][Status][Discuss]

Description

给定一个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。

Input

第一行两个正整数n和q (n,q<=500,000)。

第二行n个正整数a[1],a[2],...,a[n] (a[i]<=n)。

下面q行,每行两个正整数a,b (a,b<=n),表示一组询问。

Output

输出q行,每行两个整数。

Sample Input

12 5

4 3 5 5 1 1 12 12 9 9 7 1

7 2

8 11

1 2

9 10

10 5


Sample Output

2 3

1 2

2 2

0 1

-1 -1


HINT

Source

[Submit][Status][Discuss]

HOME
Back


题解

我是很久没有好好写一篇题解了。

首先分析题意,是要沿着出边跳到相同的节点,说明要求的是lca。

然后说一下基环树的套路:

  1. 树上倍增可以照常写,注意使用的条件就行了。
  2. 一遍拓扑排序把环找出来,这时候连通且未标记的肯定都在同一个环内。
  3. 以环上任意节点为起点开始做一个环标号,环长度前缀和,还有环总长。具体可以直接用出边跳。
  4. 把不在环上的点的出边反向用来建树,以环为根处理深度,环根标号

这样询问时的处理套路:

  1. 不在同一连通块内,即环根标号对应节点的环标号不同,那么lca不存在。
  2. 在同一棵树内,即环根标号相同,直接树上lca解决。
  3. 这题我们主要解决的情况——在同一连通块内,但是不在同一棵树内,即环根标号对应节点的环标号相同,但环根标号不同。

    x,y肯定要调到环上,然后分为两种情况,x,y分别跳。用预处理的3解决。

这样询问时间复杂度可以做到\(O(1)\),不过这题我的简便写法是\(O(\log n)\)的。总时间复杂度\(((n+m)\log n)\)

#include<bits/stdc++.h>
#define co const
#define il inline
#define rg register
template<class T>T read(){
rg T data=0,w=1;rg char ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-') w=-w;
for(;isdigit(ch);ch=getchar()) data=data*10+ch-'0';
return data*w;
}
template<class T>T read(T&x){
return x=read<T>();
}
typedef long long ll;
using namespace std;
typedef pair<int,int> pii;
#define x first
#define y second co int N=5e5+1;
int n,m,t,f[N][20],deg[N];
int pos[N],cnt,len[N],s[N];
int d[N],id[N];
vector<int> e[N];
queue<int> q; int lca(int x,int y){
if(d[x]>d[y]) swap(x,y);
for(int i=19;i>=0;--i)
if(d[f[y][i]]>=d[x]) y=f[y][i];
if(x==y) return x;
for(int i=19;i>=0;--i)
if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
return f[x][0];
}
pii cmp(co pii&a,co pii&b){
if(max(a.x,a.y)<max(b.x,b.y)) return a;
if(max(a.x,a.y)>max(b.x,b.y)) return b;
if(min(a.x,a.y)<min(b.x,b.y)) return a;
if(min(a.x,a.y)>min(b.x,b.y)) return b;
return a.x>=a.y?a:b;
}
int main(){
read(n),read(m);
for(int i=1;i<=n;++i) ++deg[read(f[i][0])];
// prework
t=log(n)/log(2);
for(int i=1;i<=t;++i)
for(int x=1;x<=n;++x) f[x][i]=f[f[x][i-1]][i-1];
// topsort
for(int i=1;i<=n;++i)if(!deg[i]) q.push(i);
while(q.size()){
int x=q.front();q.pop();
if(!--deg[f[x][0]]) q.push(f[x][0]);
}
// circle
for(int i=1;i<=n;++i)if(deg[i]&&!pos[i]){
++cnt;
for(int j=i;!pos[j];j=f[j][0])
pos[j]=cnt,s[j]=++len[cnt];
}
// bfs
for(int i=1;i<=n;++i){
if(pos[i]) id[i]=i,q.push(i);
else e[f[i][0]].push_back(i);
}
while(q.size()){
int x=q.front();q.pop();
for(unsigned i=0;i<e[x].size();++i){
int y=e[x][i];
d[y]=d[x]+1,id[y]=id[x];
q.push(y);
}
}
while(m--){
int x=read<int>(),y=read<int>();
if(pos[id[x]]!=pos[id[y]]) puts("-1 -1");
else if(id[x]==id[y]){
int p=lca(x,y);
printf("%d %d\n",d[x]-d[p],d[y]-d[p]);
}
else{
pii a,b;
int sx=s[id[x]],sy=s[id[y]],now=len[pos[id[x]]];
a.x=d[x]+(sy-sx+now)%now,a.y=d[y];
b.x=d[x],b.y=d[y]+(sx-sy+now)%now;
pii ans=cmp(a,b);
printf("%d %d\n",ans.x,ans.y);
}
}
return 0;
}

LG3533 [POI2012]RAN-Rendezvous的更多相关文章

  1. Solution -「基环树」做题记录

    写的大多只是思路,比较简单的细节和证明过程就不放了,有需者自取. 基环树简介 简单说一说基环树吧.由名字扩展可得这是一类以环为基础的树(当然显然它不是树. 通常的表现形式是一棵树再加一条非树边,把图画 ...

  2. [BZOJ2791][Poi2012]Rendezvous

    2791: [Poi2012]Rendezvous Time Limit: 25 Sec  Memory Limit: 128 MBSubmit: 95  Solved: 71[Submit][Sta ...

  3. 【BZOJ2791】[Poi2012]Rendezvous 倍增

    [BZOJ2791][Poi2012]Rendezvous Description 给定一个n个顶点的有向图,每个顶点有且仅有一条出边.对于顶点i,记它的出边为(i, a[i]).再给出q组询问,每组 ...

  4. 【BZOJ 2791】 2791: [Poi2012]Rendezvous (环套树、树链剖分LCA)

    2791: [Poi2012]Rendezvous Description 给定一个n个顶点的有向图,每个顶点有且仅有一条出边.对于顶点i,记它的出边为(i, a[i]).再给出q组询问,每组询问由两 ...

  5. 「POI2012」约会 Rendezvous

    #2691. 「POI2012」约会 Rendezvous 这题我简直不想说什么了,什么素质,卡常数…… “每个顶点有且仅有一条出边”,所以是一道基环树的题,首先tarjan缩点,在缩完点后的图上求a ...

  6. bzoj 2791 [Poi2012]Rendezvous 基环森林

    题目大意 给定一个n个顶点的有向图,每个顶点有且仅有一条出边. 对于顶点i,记它的出边为(i, a[i]). 再给出q组询问,每组询问由两个顶点a.b组成,要求输出满足下面条件的x.y: 从顶点a沿着 ...

  7. [BZOJ2791]:[Poi2012]Rendezvous(塔尖+倍增LCA)

    题目传送门 题目描述 给定一个有n个顶点的有向图,每个顶点有且仅有一条出边.每次询问给出两个顶点${a}_{i}$和${b}_{i}$​​,求满足以下条件的${x}_{i}$和${y}_{i}$:   ...

  8. [Poi2012]Rendezvous

    题目描述 给定一个n个顶点的有向图,每个顶点有且仅有一条出边. 对于顶点i,记它的出边为(i, a[i]). 再给出q组询问,每组询问由两个顶点a.b组成,要求输出满足下面条件的x.y: 从顶点a沿着 ...

  9. POI2012题解

    POI2012题解 这次的完整的\(17\)道题哟. [BZOJ2788][Poi2012]Festival 很显然可以差分约束建图.这里问的是变量最多有多少种不同的取值. 我们知道,在同一个强连通分 ...

随机推荐

  1. 使用3DES+Base64来加密传输iOS应用数据

    本文转载至 http://www.erblah.com/post/objective-c/shi-yong-3des-base64lai-jia-mi-chuan-shu-iosying-yong-s ...

  2. 【BZOJ5018】[Snoi2017]英雄联盟 背包

    [BZOJ5018][Snoi2017]英雄联盟 Description 正在上大学的小皮球热爱英雄联盟这款游戏,而且打的很菜,被网友们戏称为「小学生」.现在,小皮球终于受不了网友们的嘲讽,决定变强了 ...

  3. iOS 运行时详解

    注:本篇文章转自:http://www.jianshu.com/p/adf0d566c887 一.运行时简介 Objective-C语言是一门动态语言,它将很多静态语言在编译和链接时期做的事放到了运行 ...

  4. zabbix_get 命令介绍

    zabbix_get 是 zabbix 服务端的一个命令,用于检测 agent 端的配置是否正确,可以很方便地知道 key 是否能正常获取到数据,在测试自定义监控的时候特别有用 [root@crazy ...

  5. Bootstrap导航栏头部错位问题

    代码: <section class="header"> <div class="container"> <div class=& ...

  6. 2017-2018-1 20179209《Linux内核原理与分析》第十周作业

    设备与模块 设备分类 块设备 块设备可以以块为单位寻址,块大小随设备不同而不同:设备通常支持重定位操作,也就是对数据的随机访问.块设备的例子有外存,光盘等. 字符设备 字符设备不可寻址,仅供数据的流式 ...

  7. JDK动态代理连接池

    JDK动态代理   1 什么是JDK动态代理 刚刚写ItcastConnection时爽么?因为Connection中的方法太多了,每个都要写,所以很累吧.累点到是没什么,可以完成功能就是好的.但是不 ...

  8. pip-grep

    Pip-pop pip-grep主要是用于方便查看Requirements.txt中那些模块是安装了的.也就是通过输入的然后模块名称然后在Requirements.txt中进行查询.里面比较难的就是d ...

  9. HackerRank - beautiful-binary-string 【字符串】

    题意 给出一个 N 位的 01 串 然后 每次 改动 可以将其中的 (0 -> 1) 或者 (1 -> 0) 然后 求 最少几次 改动 使得 这个 01 串 当中 不存在 连续的 010 ...

  10. Data Structure Array: Maximum of all subarrays of size k

    http://www.geeksforgeeks.org/maximum-of-all-subarrays-of-size-k/ #include <iostream> #include ...