额,大一暑假写的RMQ,还是比较干净的,翻出来存起来以后备用。
#include "bits/stdc++.h"
using namespace std;
#define inf 10009
#define INF 999999999
//Sparse-Table 学习手敲:求区间最小值,nlogn预处理,1查询,不可动态维护
//ST是一个打表的递推式 ,用2的k次来拓展
//d(i,j)表示i开始的长度为2^j的最小值,j->0,i随意,此处也为0
//d(i,j+1)=min( d(i,j),d(i+2^j,j) ) 此列由前一列某两行推出
//方便程序的书写,用 d(i,j)=min( d(i,j-1),d(i+2^(j-1),j-1) )
int a[inf],d[inf][100];
int n;
void init_d() //预处理nlogn,计算d数组
{
int i,j;
for(i=0;i<n;i++)d[i][0]=a[i]; //2的0次=1
for(j=1;(1<<j)<=n;j++) //j<=n,因为表示长度,下标为n-1,长度为n
for(i=0;i+(1<<j)-1<n;i++) //要减一,此处是长度,要转为下标
d[i][j]=min( d[i][j-1],d[i+(1<<(j-1))][j-1] );
}
int cal_mini(int l,int r)
{
int k=0;
while((1<<(k+1))<=r-l+1)k++;
return min( d[l][k],d[r-(1<<k)+1][k] );
}
int main()
{
int i,j;
int m,l,r;
scanf("%d",&n);
for(i=0;i<n;i++)scanf("%d",&a[i]);
init_d();
int k=0;
while((1<<(k+1))<=n)k++;
for(i=0;i<n;i++)
{
for(j=0;j<=k;j++)
printf("%d ",d[i][j]);
printf("\n");
}
scanf("%d",&m);
while(m--)
{
scanf("%d%d",&l,&r);
printf("%d\n",cal_mini(l,r));
}
return 0;
}