hdu 2586 How far away ? lca求最短路
来源:互联网 发布:剑三丐姐脸型数据 编辑:程序博客网 时间:2024/05/18 22:15
给定一棵n个结点的树,有m个询问,问a,b之间的最短路是多少?。
设dis[x]为树根到x节点的路程,所以只要求出a,b的lca,答案就是dis[a]+dis[b]-2*dis[lca(a,b)]
dis的记录在遍历树的时候可以求得。
#include <iostream>#include <cstdio>#include <cstdlib>#include <cstring>#include <algorithm>#include <vector>#include <map>#include <cmath>#include <set>#include <queue>using namespace std;const int INF=1e9+10;const double EPS = 1e-10; typedef long long ll;const int mod=1e9+7;const int N = 40010; const int M = 25;int n,m;int dir[N];int dp[2*N][M]; //这个数组记得开到2*N,因为遍历后序列长度为2*n-1 bool vis[N]; struct edge { int u,v,w,next; }e[2*N]; int tot,head[N]; inline void add(int u ,int v ,int w ,int &k) { e[k].u = u; e[k].v = v; e[k].w = w; e[k].next = head[u]; head[u] = k++; u = u^v; v = u^v; u = u^v; e[k].u = u; e[k].v = v; e[k].w = w; e[k].next = head[u]; head[u] = k++; } int ver[2*N],R[2*N],first[N]; //ver:节点编号 R:深度 first:点编号位置 void dfs(int u ,int dep) { vis[u] = true; ver[++tot] = u; first[u] = tot; R[tot] = dep; for(int k=head[u]; k!=-1; k=e[k].next) if( !vis[e[k].v] ) { int v = e[k].v,w = e[k].w; dir[v] = dir[u] + w; dfs(v,dep+1); ver[++tot] = u; R[tot] = dep; } }//返回的不是最小值,而是最小值时对应的下标void ST(int n) { for(int i=1;i<=n;i++) dp[i][0] = i; for(int j=1;(1<<j)<=n;j++) { for(int i=1;i+(1<<j)-1<=n;i++) { int a = dp[i][j-1] , b = dp[i+(1<<(j-1))][j-1]; dp[i][j] = R[a]<R[b]?a:b; } } }//中间部分是交叉的。 int RMQ(int l,int r) { int k=0; while((1<<(k+1))<=r-l+1) k++; int a = dp[l][k], b = dp[r-(1<<k)+1][k]; //保存的是编号 return R[a]<R[b]?a:b; } int LCA(int u ,int v) { int x = first[u] , y = first[v]; if(x > y) swap(x,y); int res = RMQ(x,y); return ver[res]; } void init1(){ memset(head,-1,sizeof(head)); memset(vis,false,sizeof(vis));}void init2(){tot=0;dfs(1,1);ST(2*N-1);}int main(){//freopen("out.txt","w",stdout);int t;scanf("%d",&t);while(t--){int num=0;init1();scanf("%d %d",&n,&m);for(int i=0;i<n-1;i++){int u,v,w; scanf("%d%d%d",&u,&v,&w); add(u,v,w,num);}dir[1]=0;init2();while(m--){int u,v;scanf("%d %d",&u,&v);int k=LCA(u,v);printf("%d\n",dir[u]+dir[v]-2*dir[k]);}} return 0;}
const int N = 40010; const int M = 25;int dp[2*N][M]; //这个数组记得开到2*N,因为遍历后序列长度为2*n-1 bool vis[N]; struct edge { int u,v,w,next; }e[2*N]; int tot,head[N]; inline void add(int u ,int v ,int w ,int &k) { e[k].u = u; e[k].v = v; e[k].w = w; e[k].next = head[u]; head[u] = k++; u = u^v; v = u^v; u = u^v; e[k].u = u; e[k].v = v; e[k].w = w; e[k].next = head[u]; head[u] = k++; } int ver[2*N],R[2*N],first[N]; //ver:节点编号 R:深度 first:点编号位置 void dfs(int u ,int dep) { vis[u] = true; ver[++tot] = u; first[u] = tot; R[tot] = dep; for(int k=head[u]; k!=-1; k=e[k].next) if( !vis[e[k].v] ) { int v = e[k].v,w = e[k].w; dfs(v,dep+1); ver[++tot] = u; R[tot] = dep; } }//返回的不是最小值,而是最小值时对应的下标,n=2*N-1void ST(int n) { for(int i=1;i<=n;i++) dp[i][0] = i; for(int j=1;(1<<j)<=n;j++) { for(int i=1;i+(1<<j)-1<=n;i++) { int a = dp[i][j-1] , b = dp[i+(1<<(j-1))][j-1]; dp[i][j] = R[a]<R[b]?a:b; } } }//中间部分是交叉的。 int RMQ(int l,int r) { int k=0; while((1<<(k+1))<=r-l+1) k++; int a = dp[l][k], b = dp[r-(1<<k)+1][k]; //保存的是编号 return R[a]<R[b]?a:b; } int LCA(int u ,int v) { int x = first[u] , y = first[v]; if(x > y) swap(x,y); int res = RMQ(x,y); return ver[res]; } void init1(){ memset(head,-1,sizeof(head)); memset(vis,false,sizeof(vis));}void init2(){tot=0;dfs(1,1);ST(2*N-1);}
0 0
- hdu 2586 How far away ? lca求最短路
- 【LCA求最短距离】hdu 2586 How far away ?
- hdu 2586 How far away ?(LCA 求两点距离)
- hdu 2586 How far away ?(离线tarjan求LCA)
- HDU 2586How far away ? tarjan算法求LCA
- HDU 2586 How far away? LCA模板
- HDU 2586 - How far away ? (LCA)
- hdu 2586 How far away ?(LCA)
- HDU--2586--How far away ?【LCA】
- HDU 2586 How far away ?(LCA)
- hdu 2586 How far away ?(LCA)
- hdu 2586 How far away ?(LCA)
- HDU - 2586 How far away ?(LCA)
- How far away ? (hdu 2586 LCA)
- hdu 2586 How far away ?(LCA)
- HDU---2586-How far away(LCA)
- HDU 2586 How far away ?(LCA)
- HDU 2586How far away ?(LCA)
- jquery checked 操作多选
- Android 网络框架Volley源码讲解
- xUtils框架
- Sping基础
- java 常用方法记录
- hdu 2586 How far away ? lca求最短路
- SASS学习
- redis部署和启动
- 配置文件ehcache.xml详解(3)—<ehcache>
- jetBrains 系列产品激活序列号
- webpack入门学习2
- Activity和多个fragments的切换及数据的传递
- 一天搞定CSS: 清除浮动(float)--13
- 时间对话框DatePickerDialog