Codeforces 1045C DFS 树 + 贪心 + 倍增

最新推荐文章于 2022-04-17 20:00:17 发布

SHOHOKUKU

最新推荐文章于 2022-04-17 20:00:17 发布

阅读量250

点赞数

分类专栏：图论基本算法文章标签：算法

本文链接：https://blog.csdn.net/neweryyy/article/details/123015757

版权

基本算法同时被 2 个专栏收录

252 篇文章 1 订阅

订阅专栏

图论

205 篇文章 3 订阅

订阅专栏

题意

传送门 Codeforces 1045C Hyperspace Highways

题解

图中任一简单环节点构成的子图一定是完全图。构造 $D F S$ 树，由于非树边仅可能为返祖边，所以，若任一节点位于任一简单环中，那么环上的节点是这个节点的祖先节点或孙子节点。

$L C A (u, v)$ 一定是路径 $(u, v)$ 的必经点（否则会出现非返祖边的非树边），故分别考虑 $u - L C A (u, v)$ 与 $v - L C A (u, v)$ 的路径即可；位于同一个环中（考虑最大环）的节点一步能到达的深度最小的祖先节点相同，那么贪心地走这个祖先节点即可，若一步恰好超出 $L C A$ ，根据完全图性质，则一定也可以一步走到 $L C A$ 。

总时间复杂度 $O\big((N+Q)\log N)\big)$ 。

#include <bits/stdc++.h>
using namespace std;
#define rep(i, l, r) for (int i = l, _ = r; i < _; ++i)
#define _rep(i, r, l) for (int i = r - 1, _ = l; i >= _; --i)
#define pb push_back
typedef long long ll;
const int MAXN = 1E5 + 5, MAXLG = 18;
int N, M, Q;
int dep[MAXN], lg[MAXN], par[MAXLG][MAXN], f[MAXLG][MAXN];
vector<int> G[MAXN];

void dfs(int u, int p, int d)
{
    dep[u] = d, par[0][u] = f[0][u] = p;
    for (auto &v : G[u])
        if (v != p && dep[v] != -1)
        {
            int &t = f[0][u];
            if (t != -1 && dep[v] < dep[t])
                t = v;
        }
    rep(i, 0, lg[N])
    {
        par[i + 1][u] = par[i][u] == -1 ? -1 : par[i][par[i][u]];
        f[i + 1][u] = f[i][u] == -1 ? -1 : f[i][f[i][u]];
    }
    for (auto &v : G[u])
        if (dep[v] == -1)
            dfs(v, u, d + 1);
}

int lca(int u, int v)
{
    if (dep[u] < dep[v])
        swap(u, v);
    int d;
    while ((d = dep[u] - dep[v]) > 0)
        u = par[lg[d]][u];
    if (u == v)
        return u;
    _rep(i, lg[N] + 1, 0)
    {
        if (par[i][u] != par[i][v])
            u = par[i][u], v = par[i][v];
    }
    return par[0][u];
}

int get(int u, int p)
{
    if (dep[u] == dep[p])
        return 0;
    int res = 0;
    _rep(i, lg[N] + 1, 0)
    {
        if (f[i][u] != -1 && dep[f[i][u]] > dep[p])
            res += 1 << i, u = f[i][u];
    }
    return res + 1;
}

int main()
{
    ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    cin >> N >> M >> Q;
    rep(i, 0, M)
    {
        int u, v;
        cin >> u >> v;
        --u, --v;
        G[u].pb(v), G[v].pb(u);
    }
    lg[0] = -1;
    rep(i, 1, N + 1)
        lg[i] = lg[i >> 1] + 1;
    memset(dep, -1, sizeof(dep));
    dfs(0, -1, 0);
    while (Q--)
    {
        int u, v;
        cin >> u >> v;
        --u, --v;
        int p = lca(u, v);
        int res = get(u, p) + get(v, p);
        cout << res << '\n';
    }
    return 0;
}