2742 聚会 (最近公共祖先LCA)

发布于 2021-09-02 12:56

题目描述

Y岛风景美丽宜人,气候温和,物产丰富。

Y岛上有N个城市(编号1,2,…,N),有N-1条城市间的道路连接着它们。

每一条道路都连接某两个城市。

幸运的是,小可可通过这些道路可以走遍Y岛的所有城市。

神奇的是,乘车经过每条道路所需要的费用都是一样的。

小可可,小卡卡和小YY经常想聚会,每次聚会,他们都会选择一个城市,使得3个人到达这个城市的总费用最小。

由于他们计划中还会有很多次聚会,每次都选择一个地点是很烦人的事情,所以他们决定把这件事情交给你来完成。

他们会提供给你地图以及若干次聚会前他们所处的位置,希望你为他们的每一次聚会选择一个合适的地点。

输入

第一行两个正整数,N和M,分别表示城市个数和聚会次数。N≤500000,M≤500000

后面有N-1行,每行用两个正整数A和B表示编号为A和编号为B的城市之间有一条路。

再后面有M行,每行用三个正整数表示一次聚会的情况:小可可所在的城市编号,小卡卡所在的城市编号以及小YY所在的城市编号。

输出

一共有M行,每行两个数Pos和Cost,用一个空格隔开,表示第i次聚会的地点选择在编号为Pos的城市,总共的费用是经过Cost条道路所花费的费用。

样例输入

6 4

1 2

2 3

2 4

4 5

5 6

4 5 6

6 3 1

2 4 4

6 6 6

样例输出

5 2

2 5

4 1

6 0

 思路

先做树上的倍增,求出f[x][i]表示x的第2^i祖先。

对于询问a,b,c,聚会地点肯定是其中一个lca,将三种情况都尝试,距离最小的就是答案。

例如以pos=lca(a,b),求a,b,c到pos的距离为dep[a]+dep[b]+dep[c]-dep[pos]-2*dep[lca(pos,c)];

需要吐槽的是,f[M][20]这样写不超时,f[20][M]这样写就超时...

#include<bits/stdc++.h>using namespace std;const int M=500005;int n,m,pos,cost,ans1,ans2;int tot,head[M],Next[M*2],vet[M*2];int L,dep[M],f[M][20]; //写成f[20][M]就超时...void add(int a,int b){    Next[++tot]=head[a],vet[tot]=b;    head[a]=tot;}void dfs(int x,int pre){    for(int i=head[x]; i; i=Next[i]){        int y=vet[i]; if(y==pre) continue;        f[y][0]=x; dep[y]=dep[x]+1;        for(int j=1; j<L; j++) f[y][j]=f[f[y][j-1]][j-1];        dfs(y,x);    }}int lca(int x,int y){    if(dep[x]<dep[y]) swap(x,y);    for(int i=L-1; i>=0; i--)        if(dep[f[x][i]]>=dep[y]) x=f[x][i];    if(x==y) return x;    for(int i=L-1; i>=0; i--)        if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];    return f[x][0];}int main(){    scanf("%d%d",&n,&m);    while((1<<L)<n) L++;    for(int i=2; i<=n; i++){        int a,b; scanf("%d%d",&a,&b);        add(a,b); add(b,a);    }    dep[1]=1; dfs(1,0);    for(int i=1; i<=m; i++){        int a,b,c,t; scanf("%d%d%d",&a,&b,&c);        pos=lca(a,b);  t=lca(pos,c);        cost=dep[a]+dep[b]+dep[c]-dep[pos]-2*dep[t];        ans1=pos; ans2=cost;        pos=lca(a,c); t=lca(pos,b);        cost=dep[a]+dep[b]+dep[c]-dep[pos]-2*dep[t];        if(cost<ans2){ans1=pos; ans2=cost;}        cost=0;        pos=lca(b,c); t=lca(pos,a);        cost=dep[a]+dep[b]+dep[c]-dep[pos]-2*dep[t];        if(cost<ans2){ans1=pos; ans2=cost;}        printf("%d %d\n",ans1,ans2);    }    return 0;}

本文来自网络或网友投稿,如有侵犯您的权益,请发邮件至:aisoutu@outlook.com 我们将第一时间删除。

相关素材