本文搬运自本人初中博客园博客,若图片加载不出来,可到原文查看:https://www.cnblogs.com/zhangtingxi/p/15791222.html

题目链接

题目

给定一个含 NN 个元素的数组 AA,下标从 11 开始。请找出下面式子的最大值:

(A[l1]A[l1+1]A[r1])+(A[l2]A[l2+1]A[r2])(A[l_1]⨁A[l_1+1]⨁…⨁A[r_1])+(A[l_2]⨁A[l_2+1]…⨁A[r_2]),其中 1l1r1<l2r2N1≤l_1≤r_1<l_2≤r_2≤Nxyx⨁y 表示 xxyy 的按位异或。

思路

枚举分割点 ii,则此时最大值为左边的最大值 lmxlmx 和右边的最大值 rmxrmx 的和。

至于 lmxlmxrmxrmx 怎么求,以 lmxlmx 为例,lmxi=max(lmxi1,find(ai))lmx_i=\max(lmx_{i-1}, find(a_i))

find(ai)find(a_i) 的过程本质上就是在trie树上找离 aia_i 最远的点。

aia_i 表示异或前缀和。

总结

A玩感觉自己的思路好像错了,应该是把前缀和丢到trie树上,但是却对了。

后来改成前缀和之后,也AC了,说明前缀和是真的可以。

也侧面说明了数据的弱。

把前缀和丢trie树的时候,记得一定要丢0。

另一方面,对于这种分两段的题,我想不到用前缀和,只有分一段才想到。

其实对于分两段,可以前后各做一次前缀和并统计,然后枚举中间点合并,理论时间复杂度与分一段的复杂度是一样的。

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
// Problem: 1473:【例题3】Codechef REBXOR
// Contest: SSOIER
// URL: http://ybt.ssoier.cn:8088/problem_show.php?pid=1473
// Memory Limit: 131 MB
// Time Limit: 1000 ms
//
// Powered by CP Editor (https://cpeditor.org)

#include<bits/stdc++.h>
using namespace std;
//#define int long long
inline 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<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
//#define M
//#define mo
#define N 400010
struct node
{
int s[2];
}d[N*30];
int n, m, i, j, k;
int a[N], b[N], lmx[N], rmx[N];
int cnt, ans;

void trie(int x)
{
int i, j, p=1;
for(i=30; i>=0; --i)
{
j= (x& (1<<i)) > 0;
if(!d[p].s[j]) d[p].s[j]=++cnt;
p=d[p].s[j];
}
}

int dfs(int dep, int x, int k, int now)
{
if(dep<0) return now;
int i=((x&(1<<dep))>0), j=i^1;
if(d[k].s[j]) return dfs(dep-1, x, d[k].s[j], now|(1<<dep));
else return dfs(dep-1, x, d[k].s[i], now);
}

signed main()
{
// freopen("tiaoshi.in", "r", stdin);
// freopen("tiaoshi.out", "w", stdout);
n=read();
for(i=1; i<=n; ++i) b[i]=read();
memset(d, 0, sizeof(d));
cnt=1; trie(0);
for(i=1; i<=n; ++i)
{
a[i]=(a[i-1]^b[i]);
// printf("%d ", a[i]);
trie(a[i]);
lmx[i]=max(lmx[i-1], dfs(30, a[i], 1, 0));
}
// printf("\n");
memset(d, 0, sizeof(d));
cnt=1; trie(0);
for(i=n; i>=1; --i)
{
a[i]=(a[i+1]^b[i]);
// printf("%d ", a[i]);
trie(a[i]);
rmx[i]=max(rmx[i+1], dfs(30, a[i], 1, 0));
}
// for(i=1; i<=n; ++i) printf("%d ", lmx[i]); printf("\n");
// for(i=1; i<=n; ++i) printf("%d ", rmx[i]); printf("\n");
for(i=1; i<n; ++i) ans=max(ans, lmx[i]+rmx[i+1]);
// printf("\n");
printf("%d", ans);
return 0;
}