其他分享
首页 > 其他分享> > Rendezvous

Rendezvous

作者:互联网

题面

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

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

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

从顶点a沿着出边走x步和从顶点b沿着出边走y步后到达的顶点相同。
在满足条件1的情况下max(x,y)最小。
在满足条件1和2的情况下min(x,y)最小。
在满足条件1、2和3的情况下x>=y。
如果不存在满足条件1的x、y,输出-1 -1。

输入格式
第一行两个正整数n和q。

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

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

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

数据范围

n,q≤500000,
a[i]≤n,
a,b≤n

输入样例:

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

输出样例:

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

题解

这道题和那道基环树森林是类似的, 从在子树上 和 在环上 考虑

先用dfs, 跑同一个基环树上的所有点, 并找出环

我们使用倍增的方法, 处理环上的每棵子树

对于两个节点, 且二者不在同一颗基环树上, 明显 -1, -1

然后就是倍增找父节点,

直接通过 较大深度的节点 向上走 走到x 的, x y同一颗子树, 答案明显 二者深度分别减去 lca深度

然后是 两者在同意深度, 向上同时走

1.最后按通常情况 f[x][0], 应该是lca, 但是我们由于x, y可能在不同的子树上,

所以当 f[x][0] != 0, 是 lca = f[x][0], 答案和上一样, x y同一颗子树

2.f[x][0] == 0, 明显 x, y都走到了自己所在子树的根节点, 也就是环上

那么按照题意, 选择是让 x -> y, 还是 y -> x, 计算答案即可

#include <bits/stdc++.h>
#define all(n) (n).begin(), (n).end()
#define se second
#define fi first
#define pb push_back
#define mp make_pair
#define sqr(n) (n)*(n)
#define rep(i,a,b) for(int i=a;i<=(b);++i)
#define per(i,a,b) for(int i=a;i>=(b);--i)
#define IO ios::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
typedef pair<ll, ll> PLL;
typedef vector<int> VI;
typedef double db;

const int N = 5e5 + 5;

int n, m, _, k, t;
int h[N], to[N << 1], ne[N << 1], tot;
int f[N][20], dep[N];
int a[N], idx[N], idcnt, idc[N], sizc[N];
int cir[N], cid;
bool v[N];

void add(int u, int v) {
    ne[++tot] = h[u]; h[u] = tot; to[tot] = v;
}

void bfs(int s) {
    queue<int> q; q.push(s); dep[s] = 1;
    while (!q.empty()) {
        int x = q.front(); q.pop();
        for (int i = h[x]; i; i = ne[i]) {
            if (!(i & 1)) continue;
            int y = to[i];
            if (v[y] || dep[y]) continue;
            dep[y] = dep[x] + 1; f[y][0] = x;

            for (int j = 1; j <= t; ++j) f[y][j] = f[f[y][j - 1]][j - 1];
            q.push(y);
        }
    }
}

PII lca(int x, int y) {
    if (idx[x] != idx[y]) return { -1, -1 };
    int a = x, b = y; bool flag = 0;

    if (dep[x] > dep[y]) swap(x, y), flag = 1;
    per(i, t, 0) if (dep[f[y][i]] >= dep[x]) y = f[y][i];
    if (x == y) return { dep[a] - dep[x], dep[b] - dep[x] };

    per(i, t, 0) if (f[x][i] != f[y][i]) x = f[x][i], y = f[y][i];
    if (f[x][0]) return { dep[a] - dep[f[x][0]], dep[b] - dep[f[x][0]] };

    x = idc[x], y = idc[y];
    if (flag) swap(y, x);
    PII p = { dep[a] - 1, dep[b] - 1 };
    int bd = y > x ? sizc[idx[a]] - y + x : x - y;
    int ad = y > x ? y - x : sizc[idx[a]] - x + y;

    if (max(p.fi + ad, p.se) < max(p.se + bd, p.fi)) p.fi += ad;
    else if (max(p.se + bd, p.fi) < max(p.fi + ad, p.se)) p.se += bd;
    else if (min(p.se, p.fi + ad) > min(p.fi, p.se + bd)) p.se += bd;
    else p.fi += ad;

    return p;
}

int dfsc(int u, int bian) {
    idx[u] = idcnt; int s = 0;
    for (int i = h[u]; i; i = ne[i]) {
        int y = to[i];
        if (bian == i + 1 >> 1 || v[y]) continue;
        if (!idx[y]) s = max(s, dfsc(y, i + 1 >> 1));
        else cir[idc[s = u] = ++cid] = u, v[u] = 1;
    }
    return s;
}

void work(int u) {
    cid = 0; ++idcnt;
    for (u = a[dfsc(u, 0)];u != cir[1]; u = a[u]) 
        cir[idc[u] = ++cid] = u, v[u] = 1;
    sizc[idcnt] = cid;
    rep(i, 1, cid) bfs(cir[i]);
}

int main() {
    IO; cin >> n >> m;
    t = log2(n - 1) + 1;
    rep(i, 1, n) cin >> a[i], add(a[i], i), add(i, a[i]);
    rep(i, 1, n) if (!idx[i]) work(i);

    rep(i, 1, m) {
        int x, y; cin >> x >> y;
        PII p = lca(x, y);
        cout << p.fi << ' ' << p.se << '\n';
    }
    return 0;
}

标签:bd,int,dep,Rendezvous,fi,se,define
来源: https://www.cnblogs.com/2aptx4869/p/13584425.html