【景区导游——LCA】

题目

代码

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 1e5 + 10;
const int M = 2 * N;
int p[N][18], d[N], a[N];
ll dis[N][18]; //注意这里要开long long
int h[N], e[M], ne[M], idx, w[M];
int n, k;
void add(int a, int b, int c)
{
    e[idx] = b, ne[idx] = h[a], h[a] = idx, w[idx++] = c;
}
void dfs(int u)
{
    for (int i = h[u]; ~i; i = ne[i])
    {
        int j = e[i];
        if (d[j])
            continue;
        d[j] = d[u] + 1;
        p[j][0] = u;
        dis[j][0] = w[i];
        for (int k = 1; k <= 17; k++)
        {
            p[j][k] = p[p[j][k - 1]][k - 1];
            dis[j][k] = dis[j][k - 1] + dis[p[j][k - 1]][k - 1];
        }
        dfs(j);
    }
}
ll lca(int a, int b)
{
    ll retv = 0;

    if (d[a] < d[b])
        swap(a, b);
    for (int i = 17; i >= 0; i--)
    {
        if (d[p[a][i]] >= d[b])
        {
            retv += dis[a][i];
            a = p[a][i];
        }
    }

    if (a == b)
        return retv;

    for (int i = 17; i >= 0; i--)
    {
        if (p[a][i] != p[b][i])
        {
            retv += dis[a][i];
            retv += dis[b][i];
            a = p[a][i];
            b = p[b][i];
        }
    }

    retv += dis[a][0];
    retv += dis[b][0];

    return retv;
}
int main()
{
    memset(h, -1, sizeof h);

    cin >> n >> k;
    for (int i = 1; i < n; i++)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a, b, c);
        add(b, a, c);
    }
    
    d[1] = 1;
    dfs(1);
    
    ll tmp = 0;
    cin >> a[1];
    for (int i = 2; i <= k; i++)
    {
        cin >> a[i];
        tmp += lca(a[i - 1], a[i]);
    }

    for (int i = 1; i <= k; i++)
    {
        ll ans = tmp;
        if (i == 1)
            ans -= lca(a[1], a[1 + 1]);
        else if (i == k)
            ans -= lca(a[k - 1], a[k]);
        else
        {
            ans -= lca(a[i - 1], a[i]);
            ans -= lca(a[i], a[i + 1]);
            ans += lca(a[i - 1], a[i + 1]);
        }

        cout << ans << ' ';
    }

    return 0;
}
相关推荐
那个村的李富贵16 小时前
CANN加速下的AIGC“即时翻译”:AI语音克隆与实时变声实战
人工智能·算法·aigc·cann
power 雀儿16 小时前
Scaled Dot-Product Attention 分数计算 C++
算法
琹箐16 小时前
最大堆和最小堆 实现思路
java·开发语言·算法
renhongxia117 小时前
如何基于知识图谱进行故障原因、事故原因推理,需要用到哪些算法
人工智能·深度学习·算法·机器学习·自然语言处理·transformer·知识图谱
坚持就完事了17 小时前
数据结构之树(Java实现)
java·算法
算法备案代理17 小时前
大模型备案与算法备案,企业该如何选择?
人工智能·算法·大模型·算法备案
赛姐在努力.17 小时前
【拓扑排序】-- 算法原理讲解,及实现拓扑排序,附赠热门例题
java·算法·图论
野犬寒鸦19 小时前
从零起步学习并发编程 || 第六章:ReentrantLock与synchronized 的辨析及运用
java·服务器·数据库·后端·学习·算法
霖霖总总19 小时前
[小技巧66]当自增主键耗尽:MySQL 主键溢出问题深度解析与雪花算法替代方案
mysql·算法
rainbow688919 小时前
深入解析C++STL:map与set底层奥秘
java·数据结构·算法