A - 春の湊に舟の影 HDU - 2586 LCA/RMQ模板

//#include<bits/stdc++.h>
//#pragma comment(linker, "/STACK:1024000000,1024000000") 
#include<stdio.h>
#include<algorithm>
#include<queue>
#include<string.h>
#include<iostream>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<iomanip>
using namespace std;
#define ll long long
#define pb push_back

const int maxn=4e4+9;
const int maxq=210;
int n,m,root;

struct EDGE{int v;int d;};
vector<EDGE>G[maxn];

int vs[maxn*2-1];
int dep[maxn*2-1];
int depth[maxn];
int id[maxn];

void addedge(int u,int v,int d){
	G[u].pb((EDGE){v,d});
	G[v].pb((EDGE){u,d});
}

void dfs(int u,int fa,int d,int &k){
	id[u]=k;
	vs[k]=u;dep[k++]=d;
	for(int i=0;i<G[u].size();i++){
		EDGE &e=G[u][i];
		if(e.v==fa)continue;
		dfs(e.v,u,d+e.d,k);	
		vs[k]=u;dep[k++]=d;
	}
}
void init_dfs(int n){
	int k=0;
	dfs(root,-1,0,k);
}

int dp[maxn<<1][20];
int mm[maxn<<1];
void init_rmq(int n,int dep[]){
	for(int i=0;i<n;i++){
		dp[i][0]=i;	
	}
	for(int j=1;(1<<j)<=n;j++)
		for(int i=0;i+(1<<j)-1<n;i++){
			int a=dp[i][j-1],b=dp[i+(1<<(j-1))][j-1];
			if(dep[a]<=dep[b]){
				dp[i][j]=a;
			}else dp[i][j]=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];
	if(dep[a]<=dep[b])return a;
	else return b;	
}

int lca(int u,int v){
	return vs[rmq(min(id[u],id[v]),max(id[u],id[v])+1)];
}

int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++)G[i].clear();
		for(int i=1;i<n;i++){
			int u,v,d;
			scanf("%d%d%d",&u,&v,&d);
			addedge(u,v,d);
		}
		root=1;
		init_dfs(n);

		init_rmq(n*2-1,dep);

		for(int i=1;i<=m;i++){
			int u,v;
			scanf("%d%d",&u,&v);
			printf("%d\n",dep[id[u]]+dep[id[v]]-2*dep[id[lca(u,v)]]);
		}
		if(T)puts("");
	}
}

与传统RMQ不同的是DP数组的意义,

反正可以用dep数组查找到对应深度,RMQ就没必要保存最值了,直接保存dfs序列号

posted @ 2017-07-17 11:44  Drenight  阅读(122)  评论(0编辑  收藏  举报