【模板】RMQ算法 && 用RMQ做LCA（倍增+动态规划）

早学了RMQ，今天才写了
其实就是一个比较特别的动规吧(⊙o⊙)…
用来求区间最值的一个算法

嗯…就是一个动态转移方程
dp[i][j]表示以第i个数字为开头，往后2的j次方个数字中最大的那个数字是哪一个
dp[i][j]=min/max (dp[i][j-1],dp[i+ (2的(j-1)次方) ][j-1])
很神奇的一个东西。
也很好实现

查询（x,y）的时候
int k=log(y-x+1) / log(2)
ans=min/max (dp[x][k],dp[y-(2的k次方)+1][k]);

很好理解。
然后我们就在LCA里面证明一下可行性

就是说一个树的DFS序中，入队位置 left_son < father < right_son;
这里写的两个儿子不是直属father的意思，表示的是所有在father以下的节点

然后我们就可以很容易知道left_son 和 right_son 的最近公共祖先应该是
DFS序中，两儿子之间深度最小的那个
很显然，LCA被转化成了一个区间最值的问题。

有些地方可以稍作处理
然后 给自己看的代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define N 500010
using namespace std;
#if 0
Writer: Goes && G.S.M.
Just a game
Enjoy it
#endif

inline int read()
{
    char ch='*';
    while(!isdigit(ch=getchar()));
    int num=ch-'0';
    while(isdigit(ch=getchar()))num=num*10+ch-'0';
    return num;
}//好像是说这个读入优化最快的，确实不用这个在洛谷就会T
struct ss{
	int to,nex;
}edge[N<<1];
int head[N],ecnt;
void add(int x,int y){
	edge[++ecnt]=(ss){y,head[x]};
	edge[++ecnt]=(ss){x,head[y]};
	head[x]=ecnt-1;head[y]=ecnt;
}



int q[N<<1],tot;
int f[N<<1][25];
int deep[N],p[N];
int n,m,root;


void DFS(int pos,int fa){
	q[tot++]=pos;p[pos]=tot-1;
	for(int i=head[pos];i;i=edge[i].nex)
	if(edge[i].to!=fa){
		deep[edge[i].to]=deep[pos]+1;
		DFS(edge[i].to,pos);
		q[tot++]=pos;
	}
	return;
}

void RMQ(){
	for(int i=0;i<tot;i++)	f[i][0]=q[i];
	for(int j=1;(1<<j)<=tot;j++)
		for(int i=0;i+(1<<j)-1<tot;i++)
		if(deep[f[i][j-1]]<deep[f[i+(1<<(j-1))][j-1]])
			f[i][j]=f[i][j-1];
		else f[i][j]=f[i+(1<<(j-1))][j-1];
}
void take(int x,int y)
{
	if(p[y]<p[x]) swap(x,y);
	int k=log(p[y]-p[x]+1)/log(2);

	if(deep[f[p[x]][k]]<deep[f[p[y]-(1<<k)+1][k]])
        printf("%d\n",f[p[x]][k]);
    else printf("%d\n",f[p[y]-(1<<k)+1][k]);
}
int main()
{
	n=read(),m=read(),root=read();
	for(int i=1;i<=n-1;i++)
		add(read(),read());
	
	DFS(root,0);RMQ();
	
	for(int i=1;i<=m;i++)
		take(read(),read());
		
	return 0;
}

差不多就这样吧(⊙o⊙)…