一、最近公共祖先(Least Common Ancestors)
对于有根树T的两个结点u、v,最近公共祖先LCA(T,u,v)表示一个结点x,满足x是u、v的祖先且x的深度尽可能大。另一种理解方式是把T理解为一个无向无环图,而LCA(T,u,v)即u到v的最短路上深度最小的点。
这里给出一个LCA的例子:
例一
对于T=<V,E>
V={1,2,3,4,5}
E={(1,2),(1,3),(3,4),(3,5)}
则有:
LCA(T,5,2)=1
LCA(T,3,4)=3
LCA(T,4,5)=3
二、RMQ问题(Range Minimum Query)
RMQ问题是指:对于长度为n的数列A,回答若干询问RMQ(A,i,j)(i,j<=n),返回数列A中下标在[i,j]里的最小值下标。这时一个RMQ问题的例子:
例二
对数列:5,8,1,3,6,4,9,5,7 有:
RMQ(2,4)=3
RMQ(6,9)=6
RMQ问题与LCA问题的关系紧密,可以相互转换,相应的求解算法也有异曲同工之妙。下面给出LCA问题向RMQ问题的转化方法。
对树进行深度优先遍历,每当“进入”或回溯到某个结点时,将这个结点的深度存入数组E最后一位。同时记录结点i在数组中第一次出现的位置(事实上就是进入结点i时记录的位置),记做R[i]。如果结点E[i]的深度记做D[i],易见,这时求LCA(T,u,v),就等价于求E[RMQ(D,R[u],R [v])],(R[u]<R[v])。例如,对于第一节的例一,求解步骤如下:
数列E[i]为:1,2,1,3,4,3,5,3,1
R[i]为:1,2,4,5,7
D[i]为:0,1,0,1,2,1,2,1,0
于是有:
LCA(T,5,2) = E[RMQ(D,R[2],R[5])] = E[RMQ(D,2,7)] = E[3] = 1
LCA(T,3,4) = E[RMQ(D,R[3],R[4])] = E[RMQ(D,4,5)] = E[4] = 3
LCA(T,4,5) = E[RMQ(D,R[4],R[5])] = E[RMQ(D,5,7)] = E[6] = 3
易知,转化后得到的数列长度为树的结点数的两倍加一,所以转化后的RMQ问题与LCA问题的规模同次。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define maxn 100010
int cnt,id;
int head[maxn];
bool vis[maxn];
int dep[maxn*2+1], E[maxn*2+1], R[maxn];
int f[maxn*2+1][20],d[50]; //f[] is RMQ, d[i] is 2^i
struct Edge
{
int to,next,weight;
}edges[maxn]; //邻接表
void init()
{
cnt=id=0;
memset(vis,false,sizeof(vis));
memset(head,-1,sizeof(head));
}
void insert(int a, int b, int weight)
{
edges[cnt].to=b;
edges[cnt].next=head[a];
edges[cnt].weight=weight;
head[a]=cnt++;
}
void DFS(int u, int d)
{
vis[u]=1;
R[u]=id;E[id]=u;dep[id++]=d;
for(int i = head[u]; i != -1; i=edges[i].next)
{
int v=edges[i].to;
if(!vis[v])
{
DFS(v,d+1);
E[id]=u;dep[id++]=d;
}
}
}
void InitRMQ(const int &id,int n)
{
d[0]=1;
for(int i = 1; i < n; i++)d[i]=2*d[i-1];
for(int i = 0; i < id; i++)f[i][0]=i;
int k=int(log(double(n))/log(2.0))+1;
for(int j = 1; j < k; j++)
for(int i = 0; i < id; i++)
{
if(i+d[j-1]-1<id)
f[i][j]=dep[f[i][j-1]]>dep[f[i+d[j-1]][j-1]]?f[i+d[j-1]][j-1]:f[i][j-1];
else break;
}
}
int Query(int x, int y)
{
int k;
k=int(log(double(y-x+1))/log(2.0));
return dep[f[x][k]]>dep[f[y-d[k]+1][k]]?f[y-d[k]+1][k]:f[x][k];
}
void Answer()
{
int Q;scanf("%d",&Q);
for(int i = 0; i < Q; i++)
{
int x,y;
scanf("%d%d",&x,&y); //查询x,y的LCA
x=R[x];y=R[y];
if(x>y)swap(x,y);
printf("%d\n",E[Query(x,y)]);
}
}
int main()
{
int t;scanf("%d",&t);
while(t--)
{
init();
int n,m;scanf("%d%d",&n,&m);
for(int i = 0; i < m; i++)
{
int a,b,c;scanf("%d%d%d",&a,&b,&c);
insert(a,b,c);
}
DFS(1,0);
InitRMQ(id,n);
Answer();
}
return 0;
}