Good Bye 2020 E（位运算+多重和

E. Apollo versus Pan
题意：
给定一个长度为 $n$ 的序列，然后对下列式子求和，答案对 $10^9+7$ 取模
$$\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n(x_i\&x_j)\ast (x_j|x_k)$$
思路：
题解
先交换求和号顺序
$$\sum _{j=1}^n\sum _{i=1}^n\sum _{k=1}^n(x_i\&x_j)\ast (x_j|x_k)$$
提出公因式，对于最后一个项来说，$i$ 和 $j$ 都是固定的，可以看成常数提出来
$$\sum _{j=1}^n\sum _{i=1}^n((x_i\&x_j)\ast \sum _{k=1}^n(x_j|x_k))$$
对于第二项的每一项，$j$ 是固定的，$k$ 每次都是从 $1$ 到 $k$，也是固定的，所以可以提出来第三项，然后变成两项的乘积
$$\sum _{j=1}^n(\sum _{i=1}^n(x_i\&x_j))\ast (\sum _{k=1}^n(x_j|x_k))$$
然后拆位预处理，每次计算每一位的贡献就好了
code：

#include<bits/stdc++.h>
#define endl '\n'
#define ll long long
#define ull unsigned long long
#define ld long double
#define all(x) x.begin(), x.end()
#define eps 1e-6
using namespace std;
const int maxn = 2e6 + 9;
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
ll n, m;
ll d[maxn];

void work()
{
	ll cnt[70] = {0};
	cin >> n;
	for(int i = 1; i <= n; ++i){
		cin >> d[i];
		for(int j = 0; j < 60; ++j) if(d[i] & (1ll << j))
			++cnt[j];
	}
	ll ans = 0;
	for(int i = 1; i <= n; ++i){
		ll a = 0, b = 0;
		for(int j = 0; j < 60; ++j){
			ll x = (1ll << j) % mod;
			if(d[i] & (1ll << j))
				(a += cnt[j] * x % mod) %= mod, (b += n * x % mod) %= mod;
			else 
				(b += cnt[j] * x % mod) %= mod;
		}
		(ans += a * b % mod) %= mod;
	}
	cout << ans << endl;
}

int main()
{
	ios::sync_with_stdio(0);
	int TT;cin>>TT;while(TT--)
	work();
	return 0;
}