ST表——————一失足成千古恨系列2

在此先祝自己这个系列写的越少越好qwq（保证不超过4篇（flag已立））

考试原题：（绝壁是看完复联出的）

第一反应：线段树？？不对，是st表。嗯，没错。哎，st表咋写来着？？完了凉了。

结果：写暴搜的都有60分，结果我爆了0 qwq

           80-->0，与键盘无缘嘤嘤嘤

好了开始说正事

ST表用来干什么的？

给定一个区间，求最值。上面那道题是典型的模板题（虽然有毒瘤数据会卡掉st表，但我们这里不讨论）

复杂度：预处理：O(nlogn),询问：O(1)

先说预处理。

st表示一个二维数组，其中st[i][j]表示区间[i,i+2^j-1]中的最值，这里我们就拿最大值举例。显然st[i][0]=a[i],因为[i,i+1-1]中只有i一个点。那么我们怎么推其他的st[i][j]呢？

我们先来看一张图

这样我们就可以推出来：st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]。为什么呢？因为i+2^(j-1)+2^(j-1)-1=i+2*(2^(j-1))-1=i+2^j-1

这样，我们就可以预处理出来所有的st[i][j]了

int st[maxn][23]//一般第二维不会超过20，这里是应对毒瘤数据范围
void init()
{
    for(int i=1;i<=n;i++)
        st[i][0]=a[i];
    for(int  i=1;(1<<i)<=n;i++)//这里实际上枚举上面说的j
    {for(int j=0;j+(1<<i)-1<n;j++)
       {st[j][i]=max(st[j][i-1],st[j+(1<<(i-1))][i-1]);//这里为什么j是第一维呢？这里我们max()里面变的只有第一维，所以我们要先处理第一维，才能保证后面正常进行（就是个顺序问题）
        }
    } 
}

查询给的区间可不是卡好2^k的，所以我们应该怎么查询呢？我们第二维表示的是i+2^j-1,我们看到这里有个次方，我们可以玄学的考虑一下log。（以2为底）

   这里区间的长度显然是2^j,我们已经考虑到log，那不妨对区间长度len=2^j来个log.（向下取整），设t=log(len)/log(2)向下取整，那么2^t一定<len/2。我们发现，l+2^t一定是在[l,r]的中间往右，而r-2^t一定在中间往左。[l,l+2^t-1]与[r-2^t+1,r]这两个区间是可以完全覆盖[l,r]，而且还不会有超出[l,r]的部分，所以我们就可以从st[l][t]和st[r-2^t+1][t]中选一个最大值。

Update

---

 

 由于窝之前的代码是手糊的，所以出了锅。误导了您真的十分十分抱歉（我发4再也不手糊代码了qwq）

真·能AC板子的代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#define pa pair<int,int>
#include<ctime>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=2147483647;
inline ll read()
{
    char ch=getchar();
    ll x=0; 
    bool f=0;
    while(ch<'0'||ch>'9')
    {
        if(ch=='-') f=1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<3)+(x<<1)+(ch^48);
        ch=getchar();
    }
    return f?-x:x;
}
int n,m,a[1000009],st[1000009][23];
int lg[100009];
int qry(int l,int r)
{
    int qwq=log((double)(r-l+1))/log(2.0);
    return max(st[l][qwq],st[r-(1<<qwq)+1][qwq]);
}
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
    {
        a[i]=read();
        st[i][0]=a[i];
    }
    for(int i=1;(1<<i)<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            st[j][i]=max(st[j][i-1],st[j+(1<<(i-1))][i-1]);
        }
    }
    while(m--)
    {
        int l=read(),r=read();
        printf("%d\n",qry(l,r));
    }
}