【模板】LCA

#include<cstdio> #include<cstring> #include<algorithm> #include<set> #define int long long using namespace std; const int N=2e5+5; inline int read() { int X=0; bool flag=1; char ch=getchar(); while(ch<'0'||ch>'9') {if(ch=='-') flag=0; ch=getchar();} while(ch>='0'&&ch<='9') {X=(X<<1)+(X<<3)+ch-'0'; ch=getchar();} if(flag) return X; return ~(X-1); } int n,m,tot,ans; int head[N*2],nxt[N*2],to[N*2],w[N*2],cnt; int f[N][20],d[N],dis[N],dfn[N]; void add(int x,int y,int z) { to[++cnt]=y,w[cnt]=z,nxt[cnt]=head[x],head[x]=cnt; } void dfs(int x,int fath) { dfn[x]=++tot; f[x][0]=fath; for(int i=1;i<20;i++) { f[x][i]=f[f[x][i-1]][i-1]; if(!f[x][i]) break; } for(int i=head[x];i;i=nxt[i]) { int v=to[i]; if(v==fath) continue; dis[v]=dis[x]+w[i]; d[v]=d[x]+1; dfs(v,x); } } int Lca(int x,int y) { if(d[x]<d[y]) swap(x,y); for(int i=19;i>=0;i--) { if(f[x][i]&&d[f[x][i]]>=d[y]) x=f[x][i]; } if(x==y) return x; for(int i=19;i>=0;i--) { if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i]; } return f[x][0]; } int dist(int x,int y) { int lc=Lca(x,y); return dis[x]+dis[y]-2*dis[lc]; } signed main() { n=read(); for(int i=1;i<n;i++) { int x=read(),y=read(),z=read(); add(x,y,z); add(y,x,z); } }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530378.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(3)  评论(0编辑  收藏  举报  
编辑推荐:
· 记一次.NET内存居高不下排查解决与启示
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
· 没有源码,如何修改代码逻辑?
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示