题解：【CF1878E】 Iva & Pav

题目链接

如果一个区间的按位与和在第 \(k\) 位上是 \(1\)，说明这个区间内所有的数在这一位上都是 \(1\)。拆位，对于每一位预处理出每个位置向后连续 \(1\) 段最长到什么位置。对于每次询问将 \(k\) 拆位，如果当前位 \(k\) 为 \(0\)，且更高位不存在 \(1\)，那么如果一个区间按位与和在这一位上为 \(1\) 就可以大于 \(k\)，所以对于这些位置和记录的最长后缀连续长度取最大值即可。找到 \(k\) 二进制位下最高的 \(1\)，从这一位开始，如果 \(k\) 在当前考虑到的位上是 \(1\)，那么区间按位与和也必须是 \(1\)；如果 \(k\) 在这一位是 \(0\)，那么按位与和在这一位选 \(1\) 就赢了。维护一个指针表示为了和 \(k\) 保持相等取到的右端点，如果找到 \(k\) 这一位为 \(0\)，答案就和当前指针和这一位取 \(1\) 的记录的右端点中的较小值再取个最大值即可。

这样可以做到 \(\mathcal O((n + q) \log V)\)，多测清空要多清一点，这种找后继的清少了可能会继承掉上组的答案。

#include<bits/stdc++.h>
#define ld long double
#define ui unsigned int
#define ull unsigned long long
#define int long long
#define eb emplace_back
#define pb pop_back
#define ins insert
#define mp make_pair
#define pii pair<int,int>
#define fi first
#define se second
#define power(x) ((x)*(x))
using namespace std;
 
namespace FastIO
{
    template<typename T=int> inline T read()
    {
        T s=0,w=1; char c=getchar();
        while(!isdigit(c)) {if(c=='-') w=-1; c=getchar();}
        while(isdigit(c)) s=(s*10)+(c^48),c=getchar();
        return s*w;
    }
    template<typename T> inline void read(T &s)
    {
        s=0; int w=1; char c=getchar();
        while(!isdigit(c)) {if(c=='-') w=-1; c=getchar();}
        while(isdigit(c)) s=(s*10)+(c^48),c=getchar();
        s=s*w;
    }
    template<typename T,typename... Args> inline void read(T &x,Args &...args)
    {
        read(x),read(args...);
    }
    template<typename T> inline void write(T x,char ch)
    {
        if(x<0) x=-x,putchar('-');
        static char stk[25]; int top=0;
        do {stk[top++]=x%10+'0',x/=10;} while(x);
        while(top) putchar(stk[--top]);
        if(ch!='~') putchar(ch);
        return;
    }
}
using namespace FastIO;
 
namespace MTool
{   
    #define TA template<typename T,typename... Args>
    #define TT template<typename T>
    static const int Mod=998244353;
    TT inline void Swp(T &a,T &b) {T t=a;a=b;b=t;}
    TT inline void cmax(T &a,T b) {a=max(a,b);}
    TT inline void cmin(T &a,T b) {a=min(a,b);}
    TA inline void cmax(T &a,T b,Args... args) {a=max({a,b,args...});}
    TA inline void cmin(T &a,T b,Args... args) {a=min({a,b,args...});}
    TT inline void Madd(T &a,T b) {a=a+b>=Mod?a+b-Mod:a+b;}
    TT inline void Mdel(T &a,T b) {a=a-b<0?a-b+Mod:a-b;}
    TT inline void Mmul(T &a,T b) {a=a*b%Mod;}
    TT inline void Mmod(T &a) {a=(a%Mod+Mod)%Mod;}
    TT inline T Cadd(T a,T b) {return a+b>=Mod?a+b-Mod:a+b;}
    TT inline T Cdel(T a,T b) {return a-b<0?a-b+Mod:a-b;}
    TT inline T Cmul(T a,T b) {return a*b%Mod;}
    TT inline T Cmod(T a) {return (a%Mod+Mod)%Mod;}
    TA inline void Madd(T &a,T b,Args... args) {Madd(a,Cadd(b,args...));}
    TA inline void Mdel(T &a,T b,Args... args) {Mdel(a,Cadd(b,args...));}
    TA inline void Mmul(T &a,T b,Args... args) {Mmul(a,Cmul(b,args...));}
    TA inline T Cadd(T a,T b,Args... args) {return Cadd(Cadd(a,b),args...);}
    TA inline T Cdel(T a,T b,Args... args) {return Cdel(Cdel(a,b),args...);}
    TA inline T Cmul(T a,T b,Args... args) {return Cmul(Cmul(a,b),args...);}
    TT inline T qpow(T a,T b) {int res=1; while(b) {if(b&1) Mmul(res,a); Mmul(a,a); b>>=1;} return res;}
    TT inline T qmul(T a,T b) {int res=0; while(b) {if(b&1) Madd(res,a); Madd(a,a); b>>=1;} return res;}
    TT inline T spow(T a,T b) {int res=1; while(b) {if(b&1) res=qmul(res,a); a=qmul(a,a); b>>=1;} return res;}
    TT inline void exgcd(T A,T B,T &X,T &Y) {if(!B) return X=1,Y=0,void(); exgcd(B,A%B,Y,X),Y-=X*(A/B);}
    TT inline T Ginv(T x) {T A=0,B=0; exgcd(x,Mod,A,B); return Cmod(A);}
    #undef TT
    #undef TA
}
using namespace MTool;
 
inline void file()
{
    freopen(".in","r",stdin);
    freopen(".out","w",stdout);
    return;
}
 
bool Mbe;
 
namespace LgxTpre
{
    static const int MAX=300010;
    static const int inf=2147483647;
    static const int INF=4557430888798830399;
    
    int T,n,m,a[MAX],l,r,k;
    int mix,suf[35][MAX];
 
    inline void lmy_forever()
    {
    	read(T);
    	while(T--)
    	{
    		read(n),mix=0,a[n+1]=0;
    		for(int i=1;i<=n;++i) read(a[i]),cmax(mix,a[i]);
    		for(int k=0;k<=30;++k) for(int i=1;i<=n+10;++i) suf[k][i]=-1;
    		for(int k=0;k<=__lg(mix);++k) for(int i=n;i;--i) if(a[i]>>k&1) {if(a[i+1]>>k&1) suf[k][i]=suf[k][i+1]; else suf[k][i]=i;} else suf[k][i]=-1;
			read(m);
			for(int i=1;i<=m;++i)
			{
				read(l,k),r=-1;
				for(int j=__lg(mix);j>__lg(k);--j) cmax(r,suf[j][l]);
				int nr=suf[__lg(k)][l];
				for(int j=__lg(k)-1;~j&&~nr;--j) if(k>>j&1) cmin(nr,suf[j][l]); else cmax(r,min(suf[j][l],nr));
				cmax(r,nr),write(r,' ');
			}
			puts("");
		}
    }
}
 
bool Med;
 
signed main()
{
//  file();
    fprintf(stderr,"%.3lf MB\n",abs(&Med-&Mbe)/1048576.0);
    int Tbe=clock();
    LgxTpre::lmy_forever();
    int Ted=clock();
    cerr<<1e3*(Ted-Tbe)/CLOCKS_PER_SEC<<" ms\n";
    return (0-0);
}