链接:https://ac.nowcoder.com/acm/contest/370/I
来源:牛客网
时间限制:C/C++ 2秒,其他语言4秒
空间限制:C/C++ 262144K,其他语言524288K
64bit IO Format: %lld
题目描述
Rinne 最近学习了位运算相关的知识,她想运用自己学习的知识发明一个加密算法。
首先她有一个源数组 A,还有一个密钥数组 B,现在她想生成加密后的数组 C。
她发明的方法是:当计算CiCi的时候,首先将 CiCi 赋值为Ci−1Ci−1,然后加上 AiAi 分别与每一个满足 j < i 的 BjBj 异或后的和,然后加上 BiBi 分别与每一个满足 j < i 的 AjAj 异或后的和,最后加上 AiAi 与 BiBi 的异或和。
形式化的讲,关于 CiCi 的递推式为以下式子:
C0=0C0=0
Ci=Ci−1+AixorBi+(∑i−1j=1(AixorBj+AjxorBi))Ci=Ci−1+AixorBi+(∑j=1i−1(AixorBj+AjxorBi))
现在她想用程序来实现这个过程,你能帮帮她吗?由于输出可能太大,你只需要输出每个 CiCi 模 109+7109+7的结果即可。
输入描述:
第一行一个整数 N,表示数组 A 和 B 的长度。 第二行 N 个整数表示数组 A。 第三行 N 个整数表示数组 B。
输出描述:
输出一行 N 个整数,表示加密后的数组 C。
示例1
输入
复制
10 65605 70259 77306 43823 61443 98602 9261 7662 46394 83019 81393 5966 61479 24259 92528 96132 35859 47981 11702 71736
输出
复制
15796 166270 623824 1132402 1650729 2445262 3256941 4150718 5106184 6353038
备注:
N≤10^5,ai≤10^9,N≤10^5,ai≤10^9
题意:见题面;
解法:异或操作两个操作数当且仅当该位不同,异或后才为1,否则为零,所以考虑按位处理,即分别记录a,b数组前 i 个数中第 j 位 为1的个数bit_a[i][j],bit_b[i][j];具体看代码吧
附代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const double eps=1.0e-5;
const int maxn=200000+10;
const ll mod=1e9+7;
int n,bit_a[100005][40],bit_b[100005][40],a[100005],b[100005];
ll c[100005];
int main()
{
//cout <<((1<<16)+(1<<6)+5);
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 i=1;i<=n;i++){
for(int j=0;j<=31;j++){
bit_a[i][j]=bit_a[i-1][j]+((a[i]&(1<<j))?1:0);
bit_b[i][j]=bit_b[i-1][j]+((b[i]&(1<<j))?1:0);
}
}
for(int i=1;i<=n;i++){
c[i]=c[i-1];
for(int j=0;j<=31;j++){
if(a[i]&(1<<j)) c[i]=(c[i]+(ll)(i-bit_b[i][j])*(1<<j))%mod;
else c[i]=(c[i]+(ll)bit_b[i][j]*(1<<j))%mod;
if(b[i]&(1<<j)) c[i]=(c[i]+(ll)(i-bit_a[i][j])*(1<<j))%mod;
else c[i]=(c[i]+(ll)bit_a[i][j]*(1<<j))%mod;
}
c[i]=((c[i]-(ll)(a[i]^b[i]))%mod+mod)%mod;
}
for(int i=1;i<=n;i++) printf("%lld%c",c[i],i==n?'\n':' ');
}