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. 在安装mysql数据库的过程中,显示msvcp100.dll丢失?则么办?

    方案一:重装操作系统为windows10专业版 方案二:问题: 解答: 报错原因是VC运行库不全或者没有安装导致,百度搜索VC集合下载安装, 链接:https://pan.baidu.com/s/1U ...

  2. tomcat 编码问题

    默认情况下,tomcat使用的的编码方式:iso8859-1 修改tomcat下的conf/server.xml文件 找到如下代码:    < Connector port="8080 ...

  3. ubuntu 17 编译BTCoin

    一. 安装开发环境 sudo apt-get update sudo apt-get install build-essential libtool autotools-dev autoconf pk ...

  4. 智能家居DIY-空气质量检测篇-获取温度和湿度篇

    目录 智能家居DIY-空气质量检测篇-获取空气污染指数 前言 话说楼主终于升级当爸了,宝宝现在5个月了,宝宝出生的时候是冬天,正是魔都空气污染严重的时候,当时就想搞个自动开启空气净化器,由于种种原因一 ...

  5. mac查看网页时翻页

    1 fn + 上下  翻页 2 command + 上下  一下子到文档头和文档尾部

  6. 项目中nodejs包高效升级插件npm-check-updates

    nodejs包高效升级插件npm-check-updates 最近想升级npm的包 1.//常规的包升级方式/2.npm update (包) 到npm一搜发现了一个很好的升级工具 npm-check ...

  7. C# unicode GBK UTF-8和汉字互转

    界面: 源码: using System; using System.Collections.Generic; using System.ComponentModel; using System.Da ...

  8. 改善程序与设计的55个具体做法 day4

    今天晚上回到小区门口,买了点冬枣,要结账的时候想起来,钥匙没带,落公司了! TNND,没办法再回趟公司,拿了钥匙,来回一个小时,汗~ 条款10:令operator=返回一个reference to * ...

  9. Centos7 配置yum源 安装epel

    一.什么是epel如果既想获得 RHEL 的高质量.高性能.高可靠性,又需要方便易用(关键是免费)的软件包更新功能,那么 Fedora Project 推出的 EPEL(Extra Packages ...

  10. LINQ 学习路程 -- 查询操作 where

    1.where Filtering Operators Description Where Returns values from the collection based on a predicat ...