hdu 2586 How far away ? lca求最短路
2017-05-10 17:59
567 查看
给定一棵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
;
int dp[2*N][M]; //这个数组记得开到2*N,因为遍历后序列长度为2*n-1
bool vis
;
struct edge
{
int u,v,w,next;
}e[2*N];
int tot,head
;
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
;
//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;
}
设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
;
int dp[2*N][M]; //这个数组记得开到2*N,因为遍历后序列长度为2*n-1
bool vis
;
struct edge
{
int u,v,w,next;
}e[2*N];
int tot,head
;
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
;
//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 ; struct edge { int u,v,w,next; }e[2*N]; int tot,head ; 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 ; //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-1 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); }
相关文章推荐
- HDU2586 How far away ?(最近公共祖先lca,离线Tarjan,最短路)
- 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 在线算法(倍增法)/离线算法(Tarjan算法)
- How far away ? - HDU 2586 LCA 或水题
- 【LCA最近公共祖先】HDU 2586 How far away ?
- hdu 2874 Connections between cities hdu 2586 How far away ? LCA
- LCA和RMQ练习,水题(LCA模板题):【hdu 2586】 How far away ?
- 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 ?(DFS/LCA/RMQ)
- hdu 2586 How far away ? (LCA)
- hdu 2586 How far away ? Lca的模板了、、
- hdu 2586 How far away ?(LCA离线)