HDU-2586 How far away ?(Tarjian_LCA)

博客围绕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 ] dis[u] + dis[v] - 2 * dis[rt] dis[u]+dis[v]2dis[rt]
我用的离线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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值