ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Codeforces1466 E. Apollo versus Pan(位运算)

2021-01-01 12:02:03  阅读:210  来源: 互联网

标签:versus xj Codeforces1466 xi xk int sum cin Apollo


题意:

给定长度为n的非负序列x,要求计算:
在这里插入图片描述
答案对1e9+7取模。

数据范围:n<=5e5,0<=x[i]<=260

解法:

∑ i = 1 n ∑ j = 1 n ∑ k = 1 n ( x i & x j ) ∗ ( x j ∣ x k ) = ∑ i = 1 n ∑ j = 1 n ( x i & x j ) ∗ ∑ k = 1 n ( x j ∣ x k ) = ∑ j = 1 n ( ∑ i = 1 n ( x i & x j ) ∗ ∑ k = 1 n ( x j ∣ x k ) ) = ∑ j = 1 n f ( j ) ∗ g ( j ) 按 位 计 算 贡 献 即 可 \sum_{i=1}^n\sum_{j=1}^n\sum_{k=1}^n(x_{i}\&x_{j})*(x_{j}|x_{k})\\ =\sum_{i=1}^n\sum_{j=1}^n(x_{i}\&x_{j})*\sum_{k=1}^n(x_{j}|x_{k})\\ =\sum_{j=1}^n(\sum_{i=1}^n(x_{i}\&x_{j})*\sum_{k=1}^n(x_{j}|x_{k}))\\ =\sum_{j=1}^nf(j)*g(j)\\ 按位计算贡献即可 i=1∑n​j=1∑n​k=1∑n​(xi​&xj​)∗(xj​∣xk​)=i=1∑n​j=1∑n​(xi​&xj​)∗k=1∑n​(xj​∣xk​)=j=1∑n​(i=1∑n​(xi​&xj​)∗k=1∑n​(xj​∣xk​))=j=1∑n​f(j)∗g(j)按位计算贡献即可

code:

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxm=2e6+5;
const int mod=1e9+7;
int a[maxm];
int cnt[66];
int n;
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);
    int T;cin>>T;
    while(T--){
        cin>>n;
        for(int i=1;i<=n;i++)cin>>a[i];
        int ans=0;
        for(int j=0;j<60;j++)cnt[j]=0;
        for(int i=1;i<=n;i++)for(int j=0;j<60;j++)cnt[j]+=(a[i]>>j&1);
        for(int j=1;j<=n;j++){
            int f=0,g=0;
            for(int p=0;p<60;p++){
                if(a[j]>>p&1){
                    f+=(1ll<<p)%mod*cnt[p]%mod;
                    g+=(1ll<<p)%mod*n%mod;
                    f%=mod,g%=mod;
                }else{
                    g+=(1ll<<p)%mod*cnt[p]%mod;
                    g%=mod;
                }
            }
            ans+=f*g%mod;
            ans%=mod;
        }
        cout<<ans<<endl;
    }
    return 0;
}

标签:versus,xj,Codeforces1466,xi,xk,int,sum,cin,Apollo
来源: https://blog.csdn.net/weixin_44178736/article/details/112058018

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有