【BZOJ3166】[Heoi2013]Alo
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
9 2 1 4 7
Sample Output
14
HINT
【样例解释】
选择区间[1,5],最大值为 7 xor 9。
对于 100%的数据有 1 ≤ n ≤ 50000, 0 ≤ ai ≤ 10^9
题解:如果我们知道了对于一个数i,它是那些区间的次大值,然后直接用可持久化Trie树+贪心就OK了,现在问题是i到底是那些区间的次大值。
一个不难想的思路就是将所有数排序,然后一个一个扔到原数组里(我一开始从小到大排的序,结果发现没法做)
所以我们将所有数从大到小排序,然后一个一个扔回到数组里,这是他的前驱和后继都是刚好比它大的数,那么它(前驱的前驱,后继的后继)这段区间里的所有包含它的子区间的次大值一定就是这个数。具体实现方法:用set维护一下每个数的位置就好了
别忘了用set之前,先无脑的把0和n+1都扔到set里去
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #include <set> using namespace std; const int maxn=50010; int n,tot,ans; int ch[maxn*100][2],s[maxn*100],v[maxn],ls[maxn],rs[maxn],rt[maxn]; set<int> ms; struct node { int val,org; }p[maxn]; int rd() { int ret=0,f=1; char gc=getchar(); while(gc<'0'||gc>'9') {if(gc=='-')f=-f; gc=getchar();} while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar(); return ret*f; } void ins(int x,int &y,int num) { y=++tot; int i,u=y,a; for(i=1<<30;i;i>>=1) { a=(num&i)>0; ch[u][a]=++tot,ch[u][a^1]=ch[x][a^1],u=ch[u][a],x=ch[x][a],s[u]=s[x]+1; } } int query(int a,int b,int num) { int i,sum=0,d; for(i=1<<30;i;i>>=1) { d=!(num&i); if(s[ch[b][d]]-s[ch[a][d]]) a=ch[a][d],b=ch[b][d],sum|=i; else a=ch[a][d^1],b=ch[b][d^1]; } return sum; } bool cmp(node a,node b) { return a.val>b.val; } int main() { n=rd(); int i,j; ins(0,rt[0],0); for(i=1;i<=n;i++) v[i]=rd(),p[i].val=v[i],p[i].org=i,ins(rt[i-1],rt[i],v[i]); sort(p+1,p+n+1,cmp); set<int>::iterator it; ms.insert(0),ms.insert(n+1); for(i=1;i<=n;i++) { if(i>1) { it=ms.upper_bound(p[i].org),it++; if(it==ms.end()) rs[p[i].org]=n; else rs[p[i].org]=(*it)-1; it--,it--; if(it==ms.begin()) ls[p[i].org]=0; else it--,ls[p[i].org]=(*it); } ms.insert(p[i].org); } for(i=1;i<=n;i++) { if(i==p[1].org) continue; ans=max(ans,query(rt[ls[i]],rt[rs[i]],v[i])); } printf("%d",ans); return 0; }