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 30
using 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);
}