【bzoj2741】[FOTILE模拟赛]L 可持久化Trie树+分块

题目描述

FOTILE得到了一个长为N的序列A,为了拯救地球,他希望知道某些区间内的最大的连续XOR和。
即对于一个询问,你需要求出max(Ai xor Ai+1 xor Ai+2 ... xor Aj),其中l<=i<=j<=r。
为了体现在线操作,对于一个询问(x,y):
l = min ( ((x+lastans) mod N)+1 , ((y+lastans) mod N)+1 ).
r = max ( ((x+lastans) mod N)+1 , ((y+lastans) mod N)+1 ).
其中lastans是上次询问的答案,一开始为0。

输入

第一行两个整数N和M。
第二行有N个正整数,其中第i个数为Ai,有多余空格。
后M行每行两个数x,y表示一对询问。

输出

共M行,第i行一个正整数表示第i个询问的结果。

样例输入

3 3
1 4 3
0 1
0 1
4 3

样例输出

5
7
7


题解

可持久化Trie树+分块

把区间异或和转化为前缀相异或的形式,问题就转化为:在 $[l-1,r]$ 中选出 $i$ 和 $j$ ,使得 $sum_i\text{^}sum_j$ 最大。

如果确定一个端点,另一个端点在某区间内的话,可以使用可持久化Trie树来解决。

本题两个端点都是只有区间范围限制。考虑分块,预处理出 $f[x][y]$ 表示从第 $x$ 块到第 $y$ 块的答案(这里我傻逼了,可以处理出块到端点的答案)。然后对于询问,整块直接取出答案,零碎的部分再使用可持久化Trie树暴力即可。

时间复杂度 $O(n\sqrt n\log n)$

注意本题有一个坑点:输入的x和y加上lastans后可能会爆int,因此需要预先取模再相加。

#include <cmath>
#include <cstdio>
#include <algorithm>
#define N 12010
using namespace std;
int a[N] , sum[N] , root[N] , c[N * 32][2] , si[N * 32] , tot , f[115][115];
inline void insert(int x , int &y , int v)
{
int i , u = y = ++tot;
bool t;
si[y] = si[x] + 1;
for(i = 1 << 30 ; i ; i >>= 1)
t = v & i , c[u][t] = ++tot , c[u][t ^ 1] = c[x][t ^ 1] , x = c[x][t] , u = c[u][t] , si[u] = si[x] + 1;
}
inline int query(int x , int y , int v)
{
int i , ans = 0;
bool t;
for(i = 1 << 30 ; i ; i >>= 1)
{
t = v & i;
if(si[c[x][t ^ 1]] == si[c[y][t ^ 1]]) x = c[x][t] , y = c[y][t];
else ans += i , x = c[x][t ^ 1] , y = c[y][t ^ 1];
}
return ans;
}
int main()
{
int n , m , i , j , k , b , l , r , x , y , ans = 0;
scanf("%d%d" , &n , &m);
for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &a[i]) , sum[i + 1] = sum[i] ^ a[i];
b = (int)sqrt(n + 1);
for(i = 1 ; i <= n + 1 ; i ++ ) insert(root[i - 1] , root[i] , sum[i]);
for(i = 1 ; i <= n / b + 1 ; i ++ )
{
for(j = i ; j <= n / b + 1 ; j ++ )
{
f[i][j] = f[i][j - 1];
for(k = (j - 1) * b + 1 ; k <= j * b && k <= n + 1 ; k ++ )
f[i][j] = max(f[i][j] , query(root[(i - 1) * b] , root[k] , sum[k]));
}
}
while(m -- )
{
scanf("%d%d" , &l , &r) , l = (l % n + ans % n) % n + 1 , r = (r % n + ans % n) % n + 1;
if(l > r) swap(l , r);
x = (l - 1) / b + 1 , y = (r ++ ) / b + 1 , ans = f[x + 1][y - 1];
if(x == y)
for(i = l ; i <= r ; i ++ )
ans = max(ans , query(root[l - 1] , root[r] , sum[i]));
else
{
for(i = l ; i <= x * b ; i ++ ) ans = max(ans , query(root[l - 1] , root[r] , sum[i]));
for(i = r ; i > (y - 1) * b ; i -- ) ans = max(ans , query(root[l - 1] , root[r] , sum[i]));
}
printf("%d\n" , ans);
}
return 0;
}
上一篇:HTML5标签:使用canvas元素在网页上绘制线条和圆(1)


下一篇:-_-#【Canvas】导出在元素上绘制的图像