【codechef】Nikitosh and xor（线段树）

画家 Nikitosh 有一个含 N 个元素的数组 A，下标从 1 开始。他想要找出下面式子的最大值：

(A[l1] ⊕ A[l1 + 1] ⊕ · · · ⊕ A[r1]) + (A[l2] ⊕ A[l2 + 1] ⊕ · · · ⊕ A[r2])

其中，1 ≤ l1 ≤ r1 < l2 ≤ r2 ≤ N。式中 x ⊕ y 表示 x 和 y 的 按位异或

因为 Nikitosh 是一个画家而不是一个数学家，他需要你帮忙解决这个问题。

输入格式

输入数据第一行包含一个整数 N，表示数组中的元素个数。

第二行包含 N 个整数 A1, A2, . . . , AN。

输出格式

输出一行包含给定表达式可能的最大值

样例输入

51 2 3 1 2

样例输出

6

样例解释：满足条件的 (l1, r1, l2, r2) 有：(1, 2, 3, 3),(1, 2, 4, 5),(3, 3, 4, 5)。

先贴过来吧。。。。。。

#include<bits/stdc++.h>
#define MAX(x,y) (((x) < (y)) ? (y) : (x))
using namespace std;

int v1[400005],v2[400005];
int a[400005];
inline void Read(int *a)
{
    register char c=0;
    while (c<33) c=getchar_unlocked();
    *a=0;
    while (c>33)
    {
        *a=*a*10+c-'0';
        c=getchar_unlocked();
    }
}
struct Xor 
{
	int data;
  	int nd[20000000][2];
};

Xor rt;
void init()
{
	rt.data = 1;
    rt.nd[0][0] = rt.nd[0][1] = -1;
}
void insert(int number) 
{
   	int nd = 0;
   	for (int i = 31; i >= 0; --i) 
	{
    	int right = (number & (1 << i)) != 0;
      	if (rt.nd[nd][right] == -1) 
		{
	        rt.nd[nd][right] = rt.data;
	        rt.nd[rt.data][0] = -1;
			rt.nd[rt.data][1] = -1;
	        rt.data++;
    	}
    	nd = rt.nd[nd][right];
   	}
}
int query(int number) 
{
    int nd = 0;
    int MAXXOR = 0;
    for (int i = 31; i >= 0; --i) 
	{
    	int right = (number & (1 << i)) != 0;
      	int flag = 1;
      	if (rt.nd[nd][right ^ 1] == -1) 
		{
	        right ^= 1;
	        flag = 0;
      	}
      	if (flag)
        MAXXOR += (1 << i);
      	if (rt.nd[nd][right ^ 1] == -1) while (true) {}
      	nd = rt.nd[nd][right ^ 1];
    }
    return	MAXXOR;
}
int main() 
{
    int n;
	//scanf("%d",&n);
	Read(&n);
    init();
    int i;

    insert(0);
    int MMAX = 0;
    int xr1 = 0,xr2;
    
    for (i=0;i<n;i++) 
	{
    	//scanf("%d",&a[i]);
    	Read(&a[i]);
      	xr1 = xr1^a[i];
      	MMAX = max(MMAX, query(xr1));
      	v1[i]=MMAX;
      	insert(xr1);
    }
    int ans = 0;
    init();
    insert(0);
    MMAX = 0;
    xr1 = 0;
    for (int i = n-1; i >=0; --i) 
	{
  		xr1 = xr1 ^ a[i];
      	MMAX = max(MMAX, query(xr1));
      	v2[n-i-1]=MMAX;
      	insert(xr1);
      	if(i!=0)
        ans =max(v2[n-i-1]+v1[i-1],ans);
    }
    //reverse(v2.begin(),v2.end());
    
    /*for(i=0; i<n-1; i++) 
	{
       ans = max(v1[i]+v2[n-2-i],ans);
    }*/
	printf("%d\n",ans);
  	return 0;
}