RMQ问题旨在对于大量的询问,高效的返回一个区间中的最大值或最小值。RMQ为在线算法,对于一个长度为n的序列,预处理需要O(nlogn),查询只需要O(1),但不支持值的修改。
RMQ的实现
RMQ的本质是动态规划,用一个二维数组d[i][j]来存储该序列以ai为首位2^j长度的序列中的最大/最小值,那么d[i][j]就可以由d[i][j-1]和d[i+2^(j-1)][j-1]推出,如图:
算法步骤:
1.通过一个循环把d[i][0]的值全部置为a[i];
2.从小到大枚举j
从小到大枚举i(i+2^j不能大于整个序列长度n)
d[i][j]=max(d[i][j-1],d[i+(1<<(j-1))][j-1]);
核心代码:
int dp() {
for(int i=0; i<n; i++) {
d[i][0]=a[i];
}
for(int j=1; (1<<j)<=n; j++)
for(int i=0; i+(1<<j)<n; i++) {
d[i][j]=max(d[i][j-1],d[i+(1<<(j-1))][j-1]);
}
}
查询:找到最小的k使L+2*k>=R;
输出d[L][K],d[R-k+1][K]的最大值
原则是两个区间可以重复,但不能遗漏
int find_max(int L,int R) {
if (L > R) return 0;
int k=0;
while((1<<(k+1))<=R-L+1)k++;
return max(d[L][k],d[R-(1<<k)+1][k]);
}