2022-05-01每日刷题打卡

2022-05-01每日刷题打卡

代码源——每日一题

给定一个长度为 n 的数组 a1,a2,…,an。 请你求出下面式子的模1e9+7的值。

$$\sum i=1，n-1\sum j=i+1，n(aiXORaj)$$
输入格式
第一行一个数字 n。

接下来一行 n 个整数 a1,a2,…,an。

输出格式
一行一个整数表示答案。

样例输入

3
1 2 3

样例输出

6

数据规模

所有数据保证 2≤n≤300000,0≤ai<2^60

先解释一下样例：

1^ 2=3,1^3=2;结果就是：3+2=5;

2^3=1；结果就是：5+1=6。

就是当前元素后它后面的所有元素都进行xor运算，并把结果加起来，第一下可能会想到n^2的暴力解法，但这显然是不行的。

让我们把目光放到位运算来看，xor运算是，二进制同一位置上的值不同才为1，不然就为0。然后我们用二进制的形式看看样例，即：001、010、011。

我们先拿1和后面两个进行运算，1的第一个位置上是1，也就说我们想要拿到这个位置的1（xor运算后是1），那么和他进行运算的数，第一个位置应该是0才对，只有数值2满足这个条件，所以我们可以拿下1个2^0；

再看1的第二个位置，是0，也就说我们要找第二个位置是1的数，数值2和3都满足这个条件，所以我们可以拿下2个2^1；

再看第三个位置，是0，我们还是要找为1的数，但是2和3都不满足，所以我们拿下了0个2^2。

（后面还有第四个到第60个位置，为了方便理解我这里不一一距离，反正后面的都是0，都拿不下）。

那么1运算后我们得到的结果就是2*2^ 1+1 * 2^ 0=5；然后我们在用同样方法知道2可以拿下1个2^ 0，此时总结果就是5+1 *2^0=6了。

现在你应该明白了这题的写法，就是算出每一个位置上1和0的个数，然后用我们当前数来去找需要的数，如果我们的第30个位置上是0，那我们就要计算有多少个是第30个位置上的值为1，有多少，说明我们能拿下多少个2^30。最后把这些结果都加起来%1e9+7就是我们要的结果了。

#include<iostream>
using namespace std;
#include<vector>
#include<algorithm>
#include<math.h>
#include<set>
#include<numeric>
#include<string>
#include<string.h>
#include<iterator>
#include<map>
#include<unordered_map>
#include<stack>
#include<list>
#include<queue>
#include<iomanip>

#define endl '\n';
typedef long long ll;
typedef pair<ll, ll>PII;
const int N = 5e4 + 50;
const int MOD = 1e9 + 7;

inline ll read() {
    ll x = 0; char ch = getchar();
    while (ch < '0' || ch > '9') ch = getchar();
    while (ch >= '0' && ch <= '9') x = (x << 3) + (x << 1) + ch - '0', ch = getchar();
    return x;
}

inline void write(ll x) {
    if (x > 9) write(x / 10);
    putchar(x % 10 | '0');
}

int main()
{
    int n;
    n=read();
    vector<ll>v(n);
    ll res = 0, mx = 0, ans = 1, cnt = -1;
    vector<int>mymap1(61), mymap0(61);
    vector<ll>power(62,1);
    for (int i = 0; i < n; i++)v[i] = read();
    for (int i = 1; i <= 61; i++)
    {
        power[i] = power[i - 1] * 2 % MOD;
    }
    for (ll i = n - 1; i >= 0; i--)
    {
        for (ll j = 60; j >= 0; j--)
        {
            ll z = (1LL << j);
            if ((v[i] & z))
            {
                res = (res % MOD + (mymap0[j] * power[j]) % MOD)%MOD;
                mymap1[j]++;
            }
            else
            {
                res = (res % MOD + (mymap1[j] * power[j]) % MOD) % MOD;
                mymap0[j]++;
            }
        }
    }

    write(res);
    return 0;
}