Description
给出 n n 个非负整数和 m m 组查询,每次查询给出一区间,查询满足 l≤i≤j≤r l ≤ i ≤ j ≤ r 且 ai a i ^ ai+1 a i + 1 ^ ... . . . ^ aj=k a j = k 的 (i,j) ( i , j ) 个数
Input
第一行三个整数 n,m,k n , m , k ,之后输入 n n 个整数,最后 m m 行每行一组查询
Output
对于每组查询,输出查询结果
Sample Input
6 2 3
1 2 1 1 0 3
1 6
3 5
Sample Output
7
0
Solution
记 si=a1 s i = a 1 ^ a2 a 2 ^ ... . . . ^ ai a i ,则查询转化为求 l≤i≤j≤r,si−1 l ≤ i ≤ j ≤ r , s i − 1 ^ sj=k s j = k 的个数,把所有查询左端点减一,则查询转化为求 l≤i<j≤r,si l ≤ i < j ≤ r , s i ^ sj=k s j = k 的个数,固定右端点 j j ,则满足条件的左端点为且 si s i = sj s j ^ k k ,固定左端点同理,那么只要记录中满足 l≤i≤r,si=x l ≤ i ≤ r , s i = x 的 i i 个数,那么当右端点从 r r 变成时,对答案的贡献为 num[sr+1 n u m [ s r + 1 ^ k] k ] ,右端点减一或者左端点加一减一同理,故每次区间端点移动一个位置只需要 O(1) O ( 1 ) 的修改和更新贡献,用莫队算法 O(nn−−√) O ( n n ) 即可求解
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 1111111
int n,m,k,a[maxn],sum[maxn],num[maxn],pos[maxn];
ll ans;
struct node
{
int l,r,id;
ll ans;
bool operator<(const node &b)const
{
if(pos[l]!=pos[b.l])return l<b.l;
return r<b.r;
}
}q[maxn];
int cmp(node x,node y)
{
return x.id<y.id;
}
void And(int x)
{
x=sum[x];
ans+=num[x^k];
num[x]++;
}
void Plus(int x)
{
x=sum[x];
num[x]--;
ans-=num[x^k];
}
int main()
{
while(~scanf("%d%d%d",&n,&m,&k))
{
memset(num,0,sizeof(num));
int mm=(int)sqrt(n);
sum[0]=0;
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
sum[i]=sum[i-1]^a[i];
pos[i]=(i-1)/mm+1;
}
for(int i=0;i<m;i++)
{
scanf("%d%d",&q[i].l,&q[i].r);
q[i].id=i,q[i].l--;
}
sort(q,q+m);
ans=0;
num[0]=1;
int l=0,r=0;
for(int i=0;i<m;i++)
{
while(r<q[i].r)And(++r);
while(r>q[i].r)Plus(r--);
while(l<q[i].l)Plus(l++);
while(l>q[i].l)And(--l);
q[i].ans=ans;
}
sort(q,q+m,cmp);
for(int i=0;i<m;i++)printf("%I64d\n",q[i].ans);
}
return 0;
}