最大值

18 篇文章 1 订阅
8 篇文章 0 订阅

题目大意

给定一个序列,求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;
        }
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值