题面
首先有个性质,就是以某个点为开始\结束的所有前缀\后缀or和只有
log
log
种,虽然很显然,但要是没想到这题就别做了。
考虑分块,设块大小为
B
B
,每个块记录,
su
s
u
,
pr
p
r
三个数组,分别表示块内长度为
i
i
的or和最大值,所有不同的后缀or和,所有不同的前缀or和。三个都可以在求出。
然后修改的时候直接重构上述东西。
询问的时候先考虑答案在块内的情况,从左往右扫块,如果当前块有比答案更优的就更新,复杂度
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<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#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<int,int>
#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;i<len;i++)
chkmax(mx[i+1],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;
}