bzoj 3166: [Heoi2013]Alo （可持久化trie+线段树）

3166: [Heoi2013]Alo

Time Limit: 20 Sec  Memory Limit: 256 MB
Submit: 873  Solved: 415
[Submit][Status][Discuss]

Description

Welcome to ALO ( Arithmetic and Logistic Online)。这是一个VR MMORPG ，
如名字所见，到处充满了数学的谜题。
现在你拥有n颗宝石，每颗宝石有一个能量密度，记为ai，这些宝石的能量
密度两两不同。现在你可以选取连续的一些宝石（必须多于一个）进行融合，设为  ai, ai+1, …, a j，则融合而成的宝石的能量密度为这些宝石中能量密度的次大值
与其他任意一颗宝石的能量密度按位异或的值，即，设该段宝石能量密度次大值
为k，则生成的宝石的能量密度为max{k xor ap | ap ≠ k , i ≤ p ≤ j}。 
现在你需要知道你怎么选取需要融合的宝石，才能使生成的宝石能量密度最大。 

Input

第一行，一个整数 n，表示宝石个数。 
第二行， n个整数，分别表示a1至an，表示每颗宝石的能量密度，保证对于i ≠ j有 ai ≠ aj。 
 

Output

输出一行一个整数，表示最大能生成的宝石能量密度。 

Sample Input

5
9 2 1 4 7

Sample Output

14

HINT

【样例解释】 

选择区间[1,5]，最大值为 7 xor 9。 

 

 

对于 100%的数据有 1 ≤ n ≤ 50000, 0 ≤ ai ≤ 10^9

Source

加强型数据By Hta

[Submit][Status][Discuss]

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#define N 150000#define bit 30using namespace std;int l[N],r[N],a[N],b[N],ch[50003*32][3],size[50003*32];int n,m,tr[N],root[N],sz;int cmp(int x,int y){    return a[x]>a[y];}void update(int now){    tr[now]=tr[now<<1]+tr[now<<1|1];}void pointchange(int now,int l,int r,int pos){    if (l==r) {        tr[now]=1;        return;    }    int mid=(l+r)/2;    if (pos<=mid) pointchange(now<<1,l,mid,pos);    else pointchange(now<<1|1,mid+1,r,pos);    update(now);}int query(int now,int l,int r,int ll,int rr){    if (ll>rr) return 0;    if (ll<=l&&r<=rr) return tr[now];    int mid=(l+r)/2; int ans=0;    if (ll<=mid) ans+=query(now<<1,l,mid,ll,rr);    if (rr>mid) ans+=query(now<<1|1,mid+1,r,ll,rr);    return ans;}int findr(int now,int l,int r,int ll,int rr,int size){    if (ll>rr) return n+1;    if (l==r) return l;    int mid=(l+r)/2;    int  t=query(now,l,r,mid+1,rr);    if (t>=size) return findr(now<<1|1,mid+1,r,max(mid+1,ll),rr,size);    else return findr(now<<1,l,mid,ll,min(mid,rr),size-t);}int findl(int now,int l,int r,int ll,int rr,int size){    if (ll>rr) return 1;    if (l==r) return l;    int mid=(l+r)/2;    int t=query(now,l,r,ll,mid);    if (t>=size) return findl(now<<1,l,mid,ll,min(rr,mid),size);    else return findl(now<<1|1,mid+1,r,max(mid+1,ll),rr,size-t); }void insert(int i,int x){    int pre=root[i-1];    root[i]=++sz;    int now=root[i];    for (int i=bit-1;i>=0;i--) {        int t=(x>>i)&1;        size[now]=size[pre]+1;        ch[now][t]=++sz;        ch[now][t^1]=ch[pre][t^1];        now=ch[now][t];        pre=ch[pre][t];    }    size[now]=size[pre]+1;}int find(int i,int j,int x){    if (j-i+1<2) return 0;    int l=root[i]; int r=root[j+1];    int ans=0;    for (int i=bit-1;i>=0;i--) {        int t=((x>>i)&1)^1;        if (size[ch[r][t]]-size[ch[l][t]]==0) t^=1;        else ans+=1<<i;        l=ch[l][t];        r=ch[r][t];    }    return ans;}int main(){    scanf("%d",&n);    for (int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=i,l[i]=1,r[i]=n;    sort(b+1,b+n+1,cmp);    for (int i=1;i<=n;i++) {        pointchange(1,1,n,b[i]);        int t=query(1,1,n,1,b[i]-1);        //if (t==0) l[b[i]]=b[i];        //else         if (t>=2) l[b[i]]=findr(1,1,n,1,b[i]-1,2)+1;        t=query(1,1,n,b[i]+1,n);        //if (t==0) r[b[i]]=b[i];        //else         if (t>=2) r[b[i]]=findl(1,1,n,b[i]+1,n,2)-1;        //cout<<b[i]<<" "<<l[b[i]]<<" "<<r[b[i]]<<endl;    }    insert(1,0);    for (int i=1;i<=n;i++){        insert(i+1,a[i]);    }    int ans=0;    for (int i=1;i<=n;i++) {        if (i==b[1]) continue;    //  cout<<a[i]<<" "<<l[i]<<" "<<r[i]<<endl;        ans=max(ans,find(l[i],r[i],a[i]));    }    printf("%d\n",ans);}