链接:https://ac.nowcoder.com/acm/contest/5600/J
来源:牛客网
异或和之和
时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 262144K,其他语言524288K
64bit IO Format: %lld
题目描述
给一个数组,数组内有 个正整数。
求这些数任取3个数异或运算后求和的值。
有几个个三元组,计算这些三元组内部异或,之后求和。(具体操作可以见样例描述)
由于该值可能过大,输出其对 取模的值。
输入描述:
输出描述:
任取三个数、三元组内部位异或后求和对取模的值。
示例1
输入
复制
4
3 4 5 6
输出
复制
10
说明
共有4个三元组:{3,4,5}、{3,4,6}、{3,5,6}、{4,5,6}
相加为。
**这题异或和,遇到这种题我们暴力肯定不行的,我们想像像异或,按位与,按位或等等都是按照每一位进行处理的。这一题我们也可以对每一位进行处理,将每一位上的1统计起来,2^18最多64位。
只有1才是有值的,为1的情况就是
1xor1xor1
1xor0xor0
这两种情况
那么我们就是要求组合数了,比如说总共6个数,第3位1的个数有4个。
那么所有的组合就是c(4,3)+c(6-4,2)*4;
那么组合数怎么求呐。
我们知道按照公式就是n!/m!(n-m)!
因为数据过大,我们要对其进行取余操作,众所周知,取余加法,乘法,减法都适合,唯独减法不可以。
那么我们怎么办呐,有大佬就想出来逆元这个概念,将除法变成乘法就可以进行取余了。
逆元怎么求的,读者可以自行查阅一下下。
具体解释看代码
**
#include<bits/stdc++.h>
using namespace std;
#define ld long double
typedef long long ll;
const ll p=1e9+7;
int vs[65];
ll quick_pow(ll a,ll b)//快速幂
{
ll res=1;
while(b)
{
if(b&1)
{
res=res*a%p;
}
a=a*a%p;
b>>=1;
}
return res;
}
ll inv(ll a)//a的逆元,在p位质数的情况下才有,详细情况应该是费马小定理推出的,不过扩展欧几里得好像也可以求
{
return quick_pow(a,p-2);
}
ll zuheshu(ll a,ll b)
{
if(a<b)
return 0;
ll res=1;
for(ll i=0;i<b;i++)
{
res=res*(a-i)%p;
res=res*inv(i+1)%p;
}
return res;
}
int main(){
int n;
scanf("%d",&n);
for(int i=0;i<n;i++)
{
ll s;
int k=0;
scanf("%lld",&s);
while(s)//将这个数的每一位1统计起来
{
if(s&1)
{
vs[k]++;
}
k++;s>>=1;
}
}
ll res=0;
for(int i=0;i<=64;i++)
{//第i位,就需要乘上2^(i-1);
//zuheshu(n-vs[i],2)*vs[i]%p)这里就是上面说的1xor0xor0的情况,就是在所有的0的情况下选择2个,然后在选择一个1就是其种数
res+=(1LL<<i)%p*(zuheshu(vs[i],3)%p+zuheshu(n-vs[i],2)*vs[i]%p)%p;
res%=p;
}
cout<<res<<endl;
}
知识点:位运算
思路:做这一类题有个技巧,就是每一位分别去处理,1e18对应二进制的64位。
可以先统计出每一位的1的个数和0的个数,那么异或和就可以用组合数学的方式求出来了,因为异或为1一定是:
11 XORXOR 11 XORXOR 11
或者
11 XORXOR 00 XORXOR 00
这两种情况中的一种,分别处理即可。(别忘了对应位最后乘以 2^i )