[2018雅礼集训1-23]ete 分块

题面
首先有个性质，就是以某个点为开始\结束的所有前缀\后缀or和只有 log log 种，虽然很显然，但要是没想到这题就别做了。
考虑分块，设块大小为 B B ，每个块记录 mxi m x i ， su s u ， pr p r 三个数组，分别表示块内长度为 i i 的or和最大值，所有不同的后缀or和，所有不同的前缀or和。三个都可以在 O(32∗B) O ( 32 ∗ B ) 求出。
然后修改的时候直接重构上述东西。
询问的时候先考虑答案在块内的情况，从左往右扫块，如果当前块有比答案更优的就更新，复杂度 O(nB) O ( n B ) 。
再考虑跨过块的情况，同样从左往右扫块，考虑右端点再当前块中的答案。首先需要维护好当前块左端点之前的后缀or和，就是再考虑完一块之后把当前数组中所有值or上当前整块or和，然后再加入当前块的后缀or和，去重即可。然后用这些or和去与当前块的前缀拼成答案，考虑后缀从大到小扫，前缀从小到大扫，这样这部分就能做到 O(32) O ( 32 ) ，然后总的就是 O(32∗nB) O ( 32 ∗ n B )
取 B=N−−√ B = N ，那么复杂度为 O(NN−−√∗32) O ( N N ∗ 32 ) 。
代码：

#include
#include
#include
#include
#define fs first
#define sc second
#define chkmax(x,y) x=max(x,y)
#define chkmin(x,y) x=min(x,y)
#define pii pair
#define N 50010
#define ID(x) ((x-1)/B+1)
using namespace std;
const int B=240;
int n,q,a[N];
pii d[70];
int read()
{
    int x=0;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x;
}
void uniq(pii *b,int &size)
{
    int re=0;
    for(int i=1;i<=size;i++)
        if(i==size||b[i].fs>b[i+1].fs) b[++re]=b[i];
    size=re;    
}
struct block
{
    int mx[250],np,ns,qwq,l,r,len;
    pii su[35],pr[35];
    void init(int lx,int rx){l=lx;r=rx;len=r-l+1;}
    void rebuild()
    {
        memset(mx,0,sizeof(mx));
        ns=np=0;
        for(int i=l;i<=r;i++)
        {
            for(int j=1;j<=ns;j++)
                su[j].fs|=a[i];
            su[++ns]=make_pair(a[i],i);
            uniq(su,ns);
            for(int j=1;j<=ns;j++)
                chkmax(mx[i-su[j].sc+1],su[j].fs);  
        }
        for(int i=1;i1],mx[i]);
        for(int i=l;i<=r;i++)
            if(!np||(pr[np].fs|a[i])>pr[np].fs) pr[np+1]=make_pair(pr[np].fs|a[i],i-l+1),np++;  
        qwq=pr[np].fs;  
    }
}blk[250];
int main()
{
    n=read();q=read();
    for(int i=1;i<=n;i++)
        a[i]=read();
    for(int i=1;i<=ID(n);i++)
    {
        blk[i].init(i*B-B+1,min(i*B,n)),blk[i].rebuild();   
    }
    while(q--)
    {
        int opt=read(),x=read(),y;
        if(opt==1) 
        {
            y=read(),a[x]=y,blk[ID(x)].rebuild();
        }
        else
        {
            int ans=n+1,num=0;

            for(int i=1;i<=ID(n);i++)
                while(ans&&blk[i].mx[min(ans-1,blk[i].len)]>=x) ans=min(ans-1,blk[i].len);  
            for(int i=1;i<=ID(n);i++)
            {
                int top=1;
                for(int j=1;j<=blk[i].np;j++)
                    for(;top<=num&&(blk[i].pr[j].fs|d[top].fs)>=x;top++) 
                        chkmin(ans,i*B-B-d[top].sc+1+blk[i].pr[j].sc);
                for(int j=1;j<=num;j++)
                    d[j].fs|=blk[i].qwq;
                for(int j=1;j<=blk[i].ns;j++)
                    d[++num]=blk[i].su[j];
                uniq(d,num);        
            }   
            printf("%d\n",ans>n?-1:ans);
        }
    }
    return 0;
}