题目大意
对于数列a1..n(∀i≠j,ai≠aj),找到一段区间[l,r](l<r),设区间次大值为mx2,最大化max{m2⊕ai|i∈[l,r]}。
1≤n≤50000,0≤ai≤109
题目分析
对于这类题目,一个很通用的思路就是枚举次大值是哪个数,然后求出合法的最大区间,再进行求解。
首先我们讨论如何找出合法区间。
显然要使ai为区间次大值,那区间内只能有一个比它大的。我们求出ai左边比ai大的最右边的数l1,次右边的数l2;同理,在右边找出比ai大的最左边的数r1,次左边的数r2。显然答案区间为(l2,r1)或(l1,r2)。找出这些位置可以用二分加RMQ解决,这里不再累赘。
然后问题转化为给定一个区间,求区间内的数与给定值异或的最大值,显然,如果我们将区间所有数二进制存入Trie中,每次贪心地选择与定值对应位相反的方向走,问题就可以迎刃而解。
但是怎么求给定区间的Trie呢,我们发现这个满足区间减法,直接用可持续化Trie解决即可。
代码实现
#include <iostream>
#include <cstdio>
#include <cctype>
#include <cmath>
using namespace std;
int read()
{
int x=0,f=1;
char ch=getchar();
while (!isdigit(ch))
{
if (ch=='-')
f=-1;
ch=getchar();
}
while (isdigit(ch))
{
x=x*10+ch-'0';
ch=getchar();
}
return x*f;
}
const int D=30;
const int N=50000;
const int S=N*(D+1);
const int LGN=15;
bool buf[D];
struct TRIE
{
int branch[S+1][2];
int root[N+1];
int sum[S+1];
int tot;
int newnode()
{
++tot;
branch[tot][0]=branch[tot][1]=sum[tot]=0;
return tot;
}
void insert(int rt1,int &rt2,int x)
{
rt2=newnode();
branch[rt2][0]=branch[rt1][0],branch[rt2][1]=branch[rt1][1],sum[rt2]=sum[rt1]+1;
if (x==-1)
return;
insert(branch[rt1][buf[x]],branch[rt2][buf[x]],x-1);
}
int query(int rt1,int rt2,int x)
{
if (x==-1)
if (sum[branch[rt2][1^buf[x]]]-sum[branch[rt1][1^buf[x]]])
return 1;
else
return 0;
if (sum[branch[rt2][1^buf[x]]]-sum[branch[rt1][1^buf[x]]])
return (1<<x)+query(branch[rt1][1^buf[x]],branch[rt2][1^buf[x]],x-1);
else
return query(branch[rt1][buf[x]],branch[rt2][buf[x]],x-1);
}
}t;
int a[N+1];
int n,ans;
struct RMQ
{
int num[N+1][LGN+1];
int lgn;
void preparation()
{
lgn=(int)floor(log(n)/log(2));
for (int i=1;i<=n;i++)
num[i][0]=a[i];
for (int j=1;j<=lgn;j++)
for (int i=1;i<=n-(1<<j)+1;i++)
if (num[i][j-1]>num[i+(1<<j-1)][j-1])
num[i][j]=num[i][j-1];
else
num[i][j]=num[i+(1<<j-1)][j-1];
}
int query(int l,int r)
{
int lgrg=(int)floor(log(r-l+1)/log(2));
if (num[l][lgrg]>num[r-(1<<lgrg)+1][lgrg])
return num[l][lgrg];
else
return num[r-(1<<lgrg)+1][lgrg];
}
}rmq;
int search_left(int ed,int x)
{
int ret=-1,l=1,r=ed,mid;
while (l<=r)
{
mid=l+r>>1;
if (rmq.query(mid,ed)>x)
{
ret=mid;
l=mid+1;
}
else
r=mid-1;
}
return ret;
}
int search_right(int ed,int x)
{
int ret=-1,l=ed,r=n,mid;
while (l<=r)
{
mid=l+r>>1;
if (rmq.query(ed,mid)>x)
{
ret=mid;
r=mid-1;
}
else
l=mid+1;
}
return ret;
}
int main()
{
freopen("alo.in","r",stdin);
freopen("alo.out","w",stdout);
n=read();
for (int i=1;i<=n;i++)
a[i]=read();
rmq.preparation();
for (int i=1;i<=n;i++)
{
for (int j=0;j<D;j++)
buf[j]=(a[i]>>j)&1;
t.insert(t.root[i-1],t.root[i],D-1);
}
ans=0;
for (int i=1;i<=n;i++)
{
int l1,l2,r1,r2;
l1=search_left(i-1,a[i]);
if (l1!=-1)
{
l2=search_left(l1-1,a[i]);
if (l2==-1)
l2=0;
}
else
l1=0;
r1=search_right(i+1,a[i]);
if (r1!=-1)
{
r2=search_right(r1+1,a[i]);
if (r2==-1)
r2=n+1;
}
else
r1=n+1;
for (int j=0;j<D;j++)
buf[j]=(a[i]>>j)&1;
if (l1)
ans=max(ans,t.query(t.root[l2],t.root[r1-1],D-1));
if (r1!=n+1)
ans=max(ans,t.query(t.root[l1],t.root[r2-1],D-1));
}
printf("%d\n",ans);
fclose(stdin);
fclose(stdout);
return 0;
}