题目大意
给定一个序列,求ai opt aj(
i<j
)的最大值。
opt是and/or/xor
xor
一个数一个数插入进trie中每次查找一发
n log n
and和or
从高位向低位贪心
尽量使高位位运算后结果为1
我们可以枚举一个数作为必须在位运算中出现的数,然后找到一个数与其位运算结果最大。
假如是and的话,对于某一位,该数此位为1,则另一个数如果该位也能为1就强制该位是1,否则不作要求。
假如是or的话,对于某一位,该数此为为0,则另一个数如果该位也能为1就强制该位是1,否则不作要求。
那么如何检验强制某些位为1是否存在这样的数呢?
我们设
f[x]=∑ni=1[x and a[i]=x]
怎么求f[x]?首先f[a[i]]先++啦
然后考虑分治,对[l,r]求f,那么先分治求[l,mid]和[mid+1,r]
然后对于右半部分,刚好和左半部分一个数对应,也就是恰好多1个1,那么就可以直接加上左边对应的f值。具体看代码。
求出f后就知道如何强制了,直接检验强制后是否有(注意and的话强制后其本身也会算1个)
n log n
#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=100000+10,maxd=20,maxx=(1<<20)-1;
int a[maxn],bz[maxn],b[maxn][21],g[maxn*21][2],f[(1<<20)+10];
int i,j,k,l,t,n,m,ca,czy,mx,ans,tot,root;
int read(){
int x=0,f=1;
char ch=getchar();
while (ch<'0'||ch>'9'){
if (ch=='-') f=-1;
ch=getchar();
}
while (ch>='0'&&ch<='9'){
x=x*10+ch-'0';
ch=getchar();
}
return x*f;
}
int calc(int x,int y,int c){
if (c==1) return x&y;
else if (c==2) return x^y;
else return x|y;
}
void insert(int &x,int d,int id){
if (!x){
x=++tot;
g[x][0]=g[x][1]=0;
}
if (d<0) return;
insert(g[x][b[id][d]],d-1,id);
}
int find(int x,int d,int id){
if (d<0) return 0;
if (g[x][1-b[id][d]]) return (1<<d)+find(g[x][1-b[id][d]],d-1,id);
else return find(g[x][b[id][d]],d-1,id);
}
void solve(int l,int r){
if (l==r) return;
int mid=(l+r)/2,i;
solve(l,mid);
solve(mid+1,r);
fo(i,l,mid) f[i]+=f[i+mid+1-l];
}
void work(int x){
int i,j=0,t=0;
fd(i,maxd,0){
if (b[x][i]==1) t+=(1<<i);
else{
if (f[j+(1<<i)]){
j+=(1<<i);
t+=(1<<i);
}
}
}
if (t>ans) ans=t;
}
void doit(int x){
int i,j=0,t=0;
fd(i,maxd,0){
if (b[x][i]==1)
if (f[j+(1<<i)]>1){
j+=(1<<i);
t+=(1<<i);
}
}
if (t>ans) ans=t;
}
int main(){
freopen("maximum.in","r",stdin);freopen("maximum.out","w",stdout);
ca=read();
while (ca--){
n=read();czy=read();
mx=0;
fo(i,1,n){
a[i]=read();
mx=max(mx,a[i]);
}
ans=0;
if (n<=1000){
fo(i,1,n)
fo(j,i+1,n){
t=calc(a[i],a[j],czy);
if (t>ans) ans=t;
}
printf("%d\n",ans);
continue;
}
if (mx<=1024){
fo(i,0,1024) bz[i]=0;
fo(i,1,n) bz[a[i]]++;
fo(i,0,1024)
fo(j,i,1024)
if ((i<j&&bz[i]&&bz[j])||(i==j&&bz[i]>1)){
t=calc(i,j,czy);
if (t>ans) ans=t;
}
printf("%d\n",ans);
continue;
}
fo(i,1,n){
k=a[i];
fo(j,0,maxd){
b[i][j]=k%2;
k/=2;
}
}
if (czy==2){
tot=root=0;
insert(root,maxd,1);
fo(i,2,n){
t=find(root,maxd,i);
if (t>ans) ans=t;
insert(root,maxd,i);
}
printf("%d\n",ans);
continue;
}
fo(i,0,maxx) f[i]=0;
fo(i,1,n) f[a[i]]++;
solve(0,maxx);
if (czy==3){
fo(i,1,n) work(i);
printf("%d\n",ans);
continue;
}
if (czy==1){
fo(i,1,n) doit(i);
printf("%d\n",ans);
continue;
}
}
}