Good Bye 2020 E. Apollo versus Pan

最新推荐文章于 2022-09-20 23:33:48 发布

edini je sam

最新推荐文章于 2022-09-20 23:33:48 发布

阅读量127

点赞数 3

分类专栏： Codeforces 文章标签： acm竞赛

本文链接：https://blog.csdn.net/yfy20020106/article/details/116353025

版权

Codeforces 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

原题传送门
相关题意请了解后再来看此篇题解 (~~绝不是偷懒~~

正篇

题目让我们求的是 $\displaystyle \sum_{i=1}^n \sum_{j=1}^n \sum_{k=1}^n (x_i \, \& \, x_j) \cdot (x_j \, | \, x_k)$ 。这个式子看起来比较复杂，有三个变量，但是不难发现，这里变量 $\displaystyle i$ 和变量 $\displaystyle j$ 绑定，而变量 $\displaystyle j$ 又与变量 $\displaystyle k$ 绑定，且题目数据范围非常大，我们不太可能枚举两个这样的大范围变量，故考虑对和式化简，以便我们处理。注意到对每一个给定的 $\displaystyle j$ 有 $\displaystyle n$ 个值要跟它取或运算并求和，于是我们将和式进行第一步转化得到如下式子：
$\displaystyle \sum_{i=1}^n \sum_{j=1}^n(x_i \, \& \, x_j) \sum_{k=1}^n (x_j \, | \, x_k)$
但是光化简这一步，我们的和式依旧很难处理，这里的求和顺序是先给定一个 $\displaystyle i$ ，然后枚举 $\displaystyle j$ 在对每一个给定的 $\displaystyle j$ 对 $\displaystyle k$ 进行求和，形成了一个类似链一样的传递，我们不妨在思考一下，前面的两个和式意义是什么？其实就是枚举了 $\displaystyle (i,j)$ 对，那么我们是否可以考虑将 $\displaystyle j$ 转化为"最高优先级"，对于每一个 $\displaystyle i$ 有 $\displaystyle n$ 个 $\displaystyle j$ 可以与其配对，那么反过来，对于每一个 $\displaystyle j$ 同样有 $\displaystyle n$ 个 $\displaystyle i$ 与其配对。于是和式又可以进一步转化：
$\displaystyle\sum_{j=1}^n (\sum_{i=1}^n (x_i \, \& \, x_j) \sum_{k=1}^n (x_j \, | \, x_k))$
对于这个式子，我们就相对好处理很多了，我们只需要能够快速求得 $\displaystyle x_j$ 与整个序列的"&“运算和，以及 $\displaystyle x_j$ 与整个序列的”|“运算和。由于都是位运算，整体考虑，肯定是很困难的，我们不妨把每一位拿出来考虑，毕竟最多也只有60位，枚举一下还是可以接受的。
我们单独思考每一位，一个数的二进制表示的某一位无非两种情况，要么是0，要么就是1。这样我们进行分类：
$\displaystyle (1)$ 如果 $\displaystyle x_j$ 的第 $\displaystyle p$ 位是0，那么这一位的”&“运算一定是0，对答案无贡献，再考虑”|“运算，由于这一位是1，那么它或上任意值都是那个数本身，那么这一位的贡献便是所有数在这一位的贡献。即：
$\displaystyle\sum_{i=1}^n (x_{ip}\, \& \, x_{jp})=0$
$\displaystyle\sum_{k=1}^n (x_{jp} \, | \, x_{kp})=\sum_{i=1}^n x_{ip}$
$\displaystyle (2)$ 如果 $\displaystyle x_j$ 的第 $\displaystyle p$ 位是1，那么这一位的”&“运算和便是所有数在这一位的贡献，在考虑”|"运算，由于1或上任意数都是1，所以这一位将贡献n个1。即：
$\displaystyle\sum_{i=1}^n (x_{ip}\, \& \, x_{jp})=\sum_{i=1}^n x_{ip}$
$\displaystyle\sum_{k=1}^n (x_{jp} \, | \, x_{kp})=n$
注意上述求和的结果是1的个数。

那么这篇题解就到此为止了，以下是我丑陋的代码：

//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<stack>
#include<cmath>
using namespace std;
#define ACCPECT return 0;
const int p = 1e9+7;
const int N = 5e5+10;
typedef long long ll;
typedef pair<ll,ll> PII;
int t;
ll n;
ll x[N];
ll cnt[61];
void solve(){
   scanf("%lld",&n);
   for(int i = 1; i <= 60; ++i) cnt[i] = 0;
   for(int i = 1; i <= n ; ++i)scanf("%lld",&x[i]);
   for(int j = 1; j <= n; ++j){//预处理每一位1的个数
        for(int i = 0; i < 60; ++i){
            if((x[j] >> (i*1ll)) & 1)++cnt[i+1];
        }
   }
   ll ans = 0;
   ll sum_and = 0, sum_or = 0;
   for(int i = 1; i <= n; ++i){
       sum_and = 0, sum_or = 0;
       for(int j = 0; j < 60; ++j){
           if((x[i] >> (j*1ll)) & 1){
               sum_and = sum_and + ((1ll << j)%p) * cnt[j+1] % p;
               sum_and %=p;
               //cout<<j<<":::"<<sum_and<<'\n';
               sum_or = sum_or + ((1ll << j)%p) * n % p;
               sum_or %= p;
               //cout<<j<<":::"<<sum_or<<'\n';
           } else {
                sum_or += ((1ll << j)%p) * cnt[j+1] % p;
                sum_or %= p;
           }
       }
       ans = (ans + (sum_and * sum_or) % p) % p;
   }
   cout<<ans<<'\n';
}
signed main(){
#ifdef ACM_LOCAL
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
#endif
    t = 1;
    scanf("%d",&t);
    while(t--){
        solve();
    }


    fclose(stdin);
    fclose(stdout);
    //--------
    ACCPECT;
    //--------
}