【cogs2711】jump，二分答案+倍增套ST表

传送门
思路：
好久没写题解了来爽一发
还是比较有意思的一道题目
首先想到二分答案x，我们枚举每一个点i，点i在x步内所能到达的点显然能形成一个区间$[L_i,R_i]$，那么$[1,L_i-1]$和$[R_i+1,n]$中的点都要走大于x步才能到达，我们只要判断这两个区间内是否存在任意点j，它在x步内所能到达的点的区间不包含点i，即点$i\notin [L_j,R_j]$，这个我们可以搞一个前缀/后缀的min/max来记录所能到达点的区间交的情况，直接$O(1)$判断就可以了
那么现在问题就是如何快速求出点i在x步内所能到达的区间$[L_i,R_i]$
比较简单的想法是预处理每个点在$2^i$步内所能到达点的区间，然后对于每个x，把它二进制分解一下就好了，每次相当于进行区间查询min/max，可以用线段树来做，每个节点开$\log n$个变量记录所辖区间内$2^i$的L/R，空间复杂度$O(n\log n)$，时间复杂度$O(n\log^3 n)$
虽然之前$n=10^5$也有过写$O(n\log^3 n)$AC题目的情况（一道很像“购票”的树上斜率优化），但我还是觉得很虚，所以把线段树改成了ST表，空间复杂度$O(n\log^2 n)$，时间复杂度$O(n\log^2 n)$，内存要开到220M+…
如果采用我这种写法，要注意倍增处理和ST表处理交替进行
好像还有不用二分答案的做法。。。
代码：

#include<cstdio>
#include<iostream>
#define M 100005
using namespace std;
int n;
int a[M],L[17][17][M],R[17][17][M];
int pre[M][2],sub[M][2],b[M][2],bit[M];
int cal(int l,int r,int dis,bool tp)
{
    int t=bit[r-l+1];
    if (tp)
        return max(R[dis][t][l],R[dis][t][r-(1<<t)+1]);
    else
        return min(L[dis][t][l],L[dis][t][r-(1<<t)+1]);
}
bool check(int dis)
{
    --dis;
    for (int t,l,r,tl,tr,i=1;i<=n;++i)
    {
        l=i;r=i;
        t=dis;
        for (;t;t^=t&-t)
        {
            if (l!=1) tl=min(l,cal(l,r,bit[t&-t],0));
            else tl=l;
            if (r!=n) tr=max(r,cal(l,r,bit[t&-t],1));
            else tr=r;
            l=tl;r=tr;
        }
        b[i][0]=l;
        b[i][1]=r;
    }
    pre[1][0]=b[1][0];
    pre[1][1]=b[1][1];
    sub[n][0]=b[n][0];
    sub[n][1]=b[n][1];
    for (int i=2;i<=n;++i)
        pre[i][0]=max(pre[i-1][0],b[i][0]),
        pre[i][1]=min(pre[i-1][1],b[i][1]);
    for (int i=n-1;i;--i)
        sub[i][0]=max(sub[i+1][0],b[i][0]),
        sub[i][1]=min(sub[i+1][1],b[i][1]);
    for (int i=1;i<=n;++i)
    {
        if (b[i][0]>1&&(pre[b[i][0]-1][0]>i||i>pre[b[i][0]-1][1])) return 1;
        if (b[i][1]<n&&(sub[b[i][1]+1][0]>i||i>sub[b[i][1]+1][1])) return 1;
    }
    return 0;

}
main()
{
    freopen("jump2017.in","r",stdin);
    freopen("jump2017.out","w",stdout);
    scanf("%d",&n);
    bit[0]=-1;
    for (int i=1;i<=n;++i)
        scanf("%d",a+i),
        bit[i]=bit[i>>1]+1,
        L[0][0][i]=max(1,i-a[i]),
        R[0][0][i]=min(n,i+a[i]);
    for (int j=1;1<<j<=n;++j)
        for (int k=1;k+(1<<j)-1<=n;++k)
            L[0][j][k]=min(L[0][j-1][k],L[0][j-1][k+(1<<j-1)]),
            R[0][j][k]=max(R[0][j-1][k],R[0][j-1][k+(1<<j-1)]);
    for (int i=1;1<<i<=n;++i)
    {
        for (int p,q,t,k=1;k<=n;++k)
            L[i][0][k]=cal(L[i-1][0][k],R[i-1][0][k],i-1,0),
            R[i][0][k]=cal(L[i-1][0][k],R[i-1][0][k],i-1,1);
        for (int j=1;1<<j<=n;++j)
            for (int k=1;k+(1<<j)-1<=n;++k)
                L[i][j][k]=min(L[i][j-1][k],L[i][j-1][k+(1<<j-1)]),
                R[i][j][k]=max(R[i][j-1][k],R[i][j-1][k+(1<<j-1)]);
    }
    int l=1,r=n,mid=l+r>>1,ans=0;
    for (;l<=r;mid=l+r>>1)
        if (check(mid))
            ans=mid,l=mid+1;
        else
            r=mid-1;
    printf("%d\n",ans);
}