ST解决RMQ问题

它的基本方法就是DP。

F[i][j]表示从i往后数2^j个数字中的最大/最小值，即区间[i,i+2^j-1]的最值。

采取二分的思想f[i][j]=max{ f[i][j-1] , f[i+2^(j-1)][j-1] }

然后用DP进行求解即可。

 

查询：[m,n]

先求出一个最大的k满足2^k<=n-m+1

于是我们把[m,n]分成了两个长度均为2^k的区间（部分重叠）[m,m+2^k-1]和[n-2^k+1,n]

接下来返回这两个区间的最值即可

 

算法分析：

由于采用了二分的思想求解，DP的时间复杂度是O(logn)

查询的时间复杂度为O(1)

由于开了一个二维数组f[N][log N]

空间复杂度为O(N*log N)

代码长度短，效率比线段树高，处理静态数据的不二之选。

program pro;
var
        n,m:longint;
        st:array[0..200000,0..20]of longint;
        i,j,l,r:longint;

function max(x,y:longint):longint;
begin
        if x>y then exit(x)
        else exit(y);
end;

function query(l,r:longint):longint;
var
        ii,kk,tt:longint;
begin
        tt:=r-l+1;
        kk:=trunc(ln(tt)/ln(2));
        exit(max(st[l,kk],st[r-(1 shl kk)+1,kk]));
end;

procedure main;
var
        ii,jj,kk:longint;
begin
        kk:=trunc(ln(n)/ln(2))+1;
        for ii:=1 to kk do
        for jj:=1 to n do
        if jj+(1 shl (ii-1))<=n then
        st[jj,ii]:=max(st[jj,ii-1],st[jj+(1 shl (ii-1)),ii-1])
        else st[jj,ii]:=st[jj,ii-1];
end;

begin
        readln(n);
        for i:=1 to n do read(st[i,0]);
        main;
        readln(m);
        for i:=1 to m do
        begin
                readln(l,r);
                writeln(query(l,r));
        end;
end.