序列异或和

版权声明:QQ:1716638489 https://blog.csdn.net/hypHuangYanPing/article/details/80612781
【题意】给定n个数的数组A和数组B,求所有A[i]+B[j]的异或和(1<=i,j<=n)。n<=200000。
【算法】二分+模拟
【题解】将答案分成(A[i]+B[j]-A[i]^B[j])的异或和 以及 A[i]^B[j]的异或和,即单独考虑进位(后面部分很好算)。
二进制题目必须拆位,通过进位使第k位+1的数对必须满足 ( A[i] & ((1<<k)-1) ) + ( B[i] & ((1<<k)-1) ) >= (1<<k)。
首先预处理cx[k][i]=B[i] & ((1<<k)-1),然后对所有cx[k]排序。(如果是累加预处理的话先全部处理出来再排序)
枚举数位k,然后枚举A[i] & ((1<<k)-1),在cx[k-1]中二分出满足要求的数对个数%2后累加进答案。
最后异或和部分就很好算啦。
还是需要一些线性基基础的吧
复杂度O(n*28*log n)。
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+7;
int n,A[maxn],B[maxn],cx[40][maxn],dx[maxn],ans=0;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&A[i]);
    for(int i=1;i<=n;i++) scanf("%d",&B[i]);
    for(int k=0;k<28;k++){
        if(!k) for(int i=1;i<=n;i++)cx[k][i]=((1<<k)&B[i]);
        else  for(int i=1;i<=n;i++) cx[k][i]=cx[k-1][i]^((1<<k)&B[i]);
    }
    for(int k=0;k<28;k++) sort(cx[k]+1,cx[k]+n+1);
    for(int k=1;k<=28;k++){
        for(int i=1;i<=n;i++){
            dx[i]^=(1<<(k-1))&A[i];
            int x=lower_bound(cx[k-1]+1,cx[k-1]+n+1,((1<<k)-dx[i]))-cx[k-1];
            ans^=(((n-x+1)%2)<<k);
        }
    }
    int a,b,c,d;
    for(int k=0;k<=28;k++){
        a=b=c=d=0;
        for(int i=1;i<=n;i++)
            if((1<<k)&A[i])b++;
            else a++;
        for(int i=1;i<=n;i++)
            if((1<<k)&B[i])d++;
            else c++;
        ans^=(((1ll*a*d%2+1ll*b*c%2)%2)<<k);
    }
    printf("%d\n",ans);
    return 0;
}

猜你喜欢

转载自blog.csdn.net/hypHuangYanPing/article/details/80612781