HDU-2586 How far away ？（Tarjian_LCA)

最新推荐文章于 2022-12-04 12:21:01 发布

stormjing7

最新推荐文章于 2022-12-04 12:21:01 发布

阅读量280

点赞数

CC 4.0 BY-SA版权

分类专栏：解题报告 HDU LCA

本文链接：https://blog.youkuaiyun.com/GodJing007/article/details/90414867

解题报告同时被 3 个专栏收录

178 篇文章

订阅专栏

HDU

21 篇文章

订阅专栏

LCA

5 篇文章

订阅专栏

博客围绕HDU-2586题目展开，题目是给定一棵树及边权值，求多组两点最短路径。分析指出可先求出每个节点到根节点的权值，设u、v的LCA为rt，其最短路径为dis[u]+dis[v]−2∗dis[rt]，采用离线Tarjian算法，用链式前向星存树和查询。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

HDU-2586 How far away ？

题目

给一颗树，及其每条边权值，求多组两点最短路径，

分析

树上两点最短路径可以先把每个节点到跟节点的权值dis求出，设rt为u，v的LCA。
u，v最短路径： $d i s [u] + d i s [v] - 2 * d i s [r t]$
我用的离线Tarjian算法，所以要先把查询存起来，这里用的链式前向星存树，和查询。
（用在线应该方便许多）

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
#define d(x) cout << (x) << endl
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#pragma GCC diagnostic error "-std=c++11"
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
const int N = 4e4 + 10;
const int M = 2e2 + 10;

int t, n, m;
int cnt, num;
struct edge{
    int next, to, w;
} e[N];
int head[N];

struct node{
    int next, to;
    int index;  //查询序号
} ask[N];
int h[N];

int fa[N];      //并查集
int rt[N];      //求根节点
int vis[N];     //标记数组
int dis[N];     //节点到根节点权值
int ans[N];     //最后答案

void add(int x, int y, int w)
{
    e[cnt].to = y;
    e[cnt].w = w;
    e[cnt].next = head[x];
    head[x] = cnt++;
}
void addask(int u, int v, int index)
{
    ask[num].to = v;
    ask[num].next = h[u];
    ask[num].index = index;
    h[u] = num++;
    ask[num].to = u;
    ask[num].next = h[v];
    ask[num].index = index;
    h[v] = num++;
}
void init()
{
    cnt = num = 0;
    memset(head, -1, sizeof(head));
    memset(h, -1, sizeof(h));
    memset(vis, 0, sizeof(vis));
    memset(rt, 0, sizeof(rt));
    memset(dis, 0, sizeof(dis));
    for (int i = 1; i <= n; i++){
        fa[i] = i;
    }
}
int find(int x)
{
    return fa[x] = fa[x] == x ? x : find(fa[x]);
}

int tarjian(int rt)
{
    for (int i = head[rt]; i != -1; i = e[i].next){
        int x = e[i].to;
        dis[x] = dis[rt] + e[i].w;
        tarjian(x);
        fa[x] = rt;
    }
    vis[rt] = 1;
    for (int i = h[rt]; i != -1; i = ask[i].next){
        int x = ask[i].to;
        if(vis[x]){
            ans[ask[i].index] = dis[rt] + dis[x] - 2 * dis[find(x)];
        }
    }
}
void solve()
{
    int x, y, w, r;
    scanf("%d%d", &n, &m);
    init();
    for (int i = 1; i < n; i++){
        scanf("%d%d%d", &x, &y, &w);
        add(x, y, w);
        rt[y]++;
    }
    for(int i = 1; i <= n; i++){
        if(!rt[i]){
            r = i;  break;
        }
    }
    for (int i = 1; i <= m; i++){
        scanf("%d%d", &x, &y);
        addask(x, y, i);
    }
    tarjian(r);
    for(int i = 1; i <= m; i++){
        printf("%d\n", ans[i]);
    }
}
int main()
{
    scanf("%d", &t);
    while(t--){
        solve();
    }
    return 0;
}