D11 树链剖分 P3379【模板】最近公共祖先(LCA)
视频链接:D11 树链剖分 P3379【模板】最近公共祖先(LCA)_哔哩哔哩_bilibili
// 树链剖分 O(mlogn) #include <iostream> #include <cstring> #include <algorithm> #include <vector> using namespace std; const int N=500010; int n,m,s,a,b; vector<int> e[N]; int fa[N],son[N],dep[N],siz[N],top[N]; void dfs1(int u,int f){ //搞fa,son,dep fa[u]=f;siz[u]=1;dep[u]=dep[f]+1; for(int v:e[u]){ if(v==f) continue; dfs1(v,u); siz[u]+=siz[v]; if(siz[son[u]]<siz[v])son[u]=v; } } void dfs2(int u,int t){ //搞top top[u]=t; //记录链头 if(!son[u]) return; //无重儿子 dfs2(son[u],t); //搜重儿子 for(int v:e[u]){ if(v==fa[u]||v==son[u])continue; dfs2(v,v); //搜轻儿子 } } int lca(int u,int v){ while(top[u]!=top[v]){ if(dep[top[u]]<dep[top[v]])swap(u,v); u=fa[top[u]]; } return dep[u]<dep[v]?u:v; } int main(){ scanf("%d%d%d",&n,&m,&s); for(int i=1; i<n; i++){ scanf("%d%d",&a,&b); e[a].push_back(b); e[b].push_back(a); } dfs1(s,0); dfs2(s,s); while(m--){ scanf("%d%d",&a,&b); printf("%d\n",lca(a,b)); } return 0; }
// 树链剖分 O(mlogn) #include <iostream> #include <cstring> #include <algorithm> #include <vector> using namespace std; const int N=500010; int n,m,s,a,b; int h[N],to[N<<1],ne[N<<1],idx; int fa[N],son[N],dep[N],siz[N],top[N]; void add(int u,int v){ to[++idx]=v,ne[idx]=h[u],h[u]=idx; } void dfs1(int u,int f){ //搞fa,son,dep fa[u]=f;siz[u]=1;dep[u]=dep[f]+1; for(int i=h[u];i;i=ne[i]){ int v=to[i]; if(v==f) continue; dfs1(v,u); siz[u]+=siz[v]; if(siz[son[u]]<siz[v])son[u]=v; } } void dfs2(int u,int t){ //搞top top[u]=t; //记录链头 if(!son[u]) return; //无重儿子 dfs2(son[u],t); //搜重儿子 for(int i=h[u];i;i=ne[i]){ int v=to[i]; if(v==fa[u]||v==son[u])continue; dfs2(v,v); //搜轻儿子 } } int lca(int u,int v){ while(top[u]!=top[v]){ if(dep[top[u]]<dep[top[v]])swap(u,v); u=fa[top[u]]; } return dep[u]<dep[v]?u:v; } int main(){ scanf("%d%d%d",&n,&m,&s); for(int i=1; i<n; i++){ scanf("%d%d",&a,&b); add(a,b); add(b,a); } dfs1(s,0); dfs2(s,s); while(m--){ scanf("%d%d",&a,&b); printf("%d\n",lca(a,b)); } return 0; }
分类:
D 图论
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!