天天看点

ST表(RMQ) 总结

f[i][j] 表示 以i号点为起点 的长度为 2^j 次方 终点为i+2^j-1 的最大或者最小值

打表时间复杂度 O(Nlog2(N)) 查询O(1)

打表代码:

for(int i=1;i<=n;i++) 
    {
        scanf("%d",&f[i][0]);
    }
    for(int j=1;j<=log2(n);j++)
    {
        for(int i=1;i+(1<<j)-1<=n;i++)
        {
            f[i][j]=max(f[i][(j-1)],f[i+(1<<j-1)][j-1]);
        }
    }      

将区间分为两段 [i ,i+2^(j-1)-1 ] 和[2^(j-1)+i,i+2^j-1]

查询代码:

scanf("%d%d",&l,&r);
        int k=log2(r-l+1);
        printf("%d\n",max(f[l][k],f[r-(1<<k)+1][k]));      

k 为区间长度的log2 值  查询 l到l+2^k-1 和r-(2^k)+1 到 r

关于二维 RMQ 

查询一个矩阵里面的最大最小值

 例题 POJ https://vjudge.net/problem/POJ-2019

code:

//
#include<iostream>
#include<cstdio>
#include<cmath> 
using namespace std;
int f1[260][260][40],f2[260][260][40];
int n,b,kk;
int RMQ(int x,int y,int c)
{
    int k=log2(c);
    int xx=x+c-1;
    int yy=y+c-1;
    int i=x,j=y;
    int t1=min(f2[i][j][k],f2[xx-(1<<k)+1][j][k]);
    t1=min(t1,min(f2[xx-(1<<k)+1][yy-(1<<k)+1][k],f2[i][yy-(1<<k)+1][k]));
    int t2=max(f1[xx-(1<<k)+1][yy-(1<<k)+1][k],f1[i][yy-(1<<k)+1][k]);
    t2=max(max(f1[i][j][k],f1[xx-(1<<k)+1][j][k]),t2);
    return t2-t1;
    
}
int main()
{
    scanf("%d%d%d",&n,&b,&kk);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            scanf("%d",&f1[i][j][0]);
            f2[i][j][0]=f1[i][j][0];
        }
    for(int k=1;k<=log2(n);k++)
    {
        for(int i=1;i+(1<<k)-1<=n;i++)
            for(int j=1;j+(1<<k)-1<=n;j++)
            {
                int t1=min(f2[i][j][k-1],f2[i+(1<<k-1)][j][k-1]);
                int t2=min(f2[i][j+(1<<k-1)][k-1],f2[i+(1<<k-1)][j+(1<<k-1)][k-1]);
                f2[i][j][k]=min(t1,t2);
                 t1=max(f1[i][j][k-1],f1[i+(1<<k-1)][j][k-1]);
                 t2=max(f1[i][j+(1<<k-1)][k-1],f1[i+(1<<k-1)][j+(1<<k-1)][k-1]);
                f1[i][j][k]=max(t1,t2);
            }
    }
    int x,y;
    while(kk--)
    {
        scanf("%d%d",&x,&y);
        printf("%d\n",RMQ(x,y,b));
    }
}      

转载于:https://www.cnblogs.com/OIEREDSION/p/11330747.html