牛客国庆集训派对Day4.B.异或求和(按位统计)
刷牛客一战到底做到的,感觉还挺有趣...
\(Description\)
求给定\(n\)及序列\(A_i\),求$$\sum_{i\lt j\lt k}(A_i\oplus A_j)(A_j\oplus A_k)(A_i\oplus A_k)$$
\(n\leq10^5,\ A_i\leq10^9\)。
\(Solution\)
首先有一个\(O(n\log^2n)\)的做法,代码这里有,看来是跑不过去..(不知道提交记录里有没有过的)
来自国庆正睿dls课件。
能跑过的做法:
注意,口胡的,不保证正确性= =(但是算法是对的)(没办法啊找不到题解,还比较赶时间...)
当然还是想拆开乘法按位统计。
考虑枚举每一个异或结果\(2^x\)的贡献:$$Ans=\sum_{i=0}{29}\sum_{j=0}{29}\sum_{k=0}{29}2{i+j+k}\times ?$$
乘上多少呢?现在我们需要统计\(A_i\oplus A_j\)在第\(i\)为\(1\),且\(A_j\oplus A_k\)在第\(j\)位为\(1\),且\(A_k\oplus A_i\)在第\(k\)位为\(1\)的方案数(不要弄混...\(A_i\)中的\(i\)就是下标,外面的\(i\)是枚举的\(2^i\))。
考虑枚举\(A_i\)第\(i\)位是\(0\)还是\(1\),设为\(a\),那么\(A_j\)的第\(i\)位是\(a\oplus1\);同理枚举\(A_j\)的第\(j\)位是\(b\),那么\(A_k\)的第\(j\)位是\(b\oplus1\);同理枚举\(A_k\)的第\(k\)位的\(c\),那么\(A_i\)的第\(k\)位是\(c\oplus1\)。(这么打累死我了...)
那么合法的\(A_i\)就是,第\(i\)位为\(a\)且第\(k\)位为\(c\oplus1\)的数字,\(A_j,A_k\)同理...
所以记\(cnt[i][j][a][b]\)表示第\(i\)位为\(a\),第\(j\)位为\(b\)的\(A_x\)有多少个,乘起来就OK了。
预处理\(cnt\)的复杂度是\(O(n\log^2n)\),常数很小。最后求和的复杂度是\(O(2^3\log^3n)\)。
答案最后除个\(6\)。
//97ms 1892KB
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define mod 998244353
#define inv6 166374059
#define MAXIN 500000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=31;
int cnt[N][N][2][2],pw[N<<2];
char IN[MAXIN],*SS=IN,*TT=IN;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-48,c=gc());
return now;
}
inline int Calc(int i,int j,int x,int y)
{
if(i>j) std::swap(i,j), std::swap(x,y);
return cnt[i][j][x][y];
}
int main()
{
int n=read();
for(int i=1; i<=n; ++i)
{
int x=read();
for(int a=0; a<N; ++a)
for(int b=a; b<N; ++b) ++cnt[a][b][x>>a&1][x>>b&1];
}
pw[0]=1;
for(int i=1; i<90; ++i) pw[i]=pw[i-1]<<1, pw[i]>=mod&&(pw[i]-=mod);
LL ans=0;
for(int i=0; i<N; ++i)
for(int j=0; j<N; ++j)
for(int k=0; k<N; ++k)
for(int a=0; a<2; ++a)
for(int b=0; b<2; ++b)
for(int c=0; c<2; ++c)
ans+=1ll*pw[i+j+k]*Calc(i,k,a,c^1)%mod*Calc(i,j,a^1,b)%mod*Calc(j,k,b^1,c)%mod;
printf("%lld\n",ans%mod*inv6%mod);
return 0;
}
牛客国庆集训派对Day4.B.异或求和(按位统计)的更多相关文章
- 牛客国庆集训派对Day4 J-寻找复读机
链接:https://www.nowcoder.com/acm/contest/204/J 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 1048576K,其他语言20 ...
- 牛客国庆集训派对Day4 I-连通块计数(思维,组合数学)
链接:https://www.nowcoder.com/acm/contest/204/I 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 1048576K,其他语言20 ...
- 2018 牛客国庆集训派对Day4 - H 树链博弈
链接:https://ac.nowcoder.com/acm/contest/204/H来源:牛客网 题目描述 给定一棵 n 个点的树,其中 1 号结点是根,每个结点要么是黑色要么是白色 现在小 Bo ...
- 牛客国庆集训派对Day4 Solution
A 深度学习 puts(n) #include <bits/stdc++.h> using namespace std; int main() { double n; while ( ...
- 线性基求交(2019牛客国庆集训派对day4)
题意:https://ac.nowcoder.com/acm/contest/1109/C 问你有几个x满足A,B集合都能XOR出x. 思路: 就是线性基求交后,有几个基就是2^几次方. #defin ...
- 牛客国庆集训派对Day6 A Birthday 费用流
牛客国庆集训派对Day6 A Birthday:https://www.nowcoder.com/acm/contest/206/A 题意: 恬恬的生日临近了.宇扬给她准备了一个蛋糕. 正如往常一样, ...
- 2019牛客国庆集训派对day5
2019牛客国庆集训派对day5 I.Strange Prime 题意 \(P=1e10+19\),求\(\sum x[i] mod P = 0\)的方案数,其中\(0 \leq x[i] < ...
- 2019 牛客国庆集训派对day1-C Distinct Substrings(exkmp+概率)
链接:https://ac.nowcoder.com/acm/contest/1099/C来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言65536 ...
- 牛客国庆集训派对Day1 L-New Game!(最短路)
链接:https://www.nowcoder.com/acm/contest/201/L 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 1048576K,其他语言20 ...
随机推荐
- Best Free Hacking E-Books 2017 In PDF Format
1.Best Free Hacking E-Books 2017 In PDF Format: 电子书籍下载地址 后续我会更新在我的百度云资源 上,需要的留言Black Belt Hacking &a ...
- java----static关键字(包括final)
static修饰字段: 使用static关键字修饰一个字段:声明的static变量实际上就是一个全局变量 使用static关键字修饰一个方法:可以直接使用类调用方法,和对象没有关系了 使用static ...
- 远程执行shell脚本
ssh -p2016 apache@10.10.18.130 '/bin/sh /data/www/vhosts/WOStest3_ENV/update_env.sh' 需要设置shell远程免密码登 ...
- C++ Primer 笔记——OOP
1.基类通常都应该定义一个虚析构函数,即使该函数不执行任何实际操作也是如此. 2.任何构造函数之外的非静态函数都可以是虚函数,关键字virtual只能出现在类内部的声明语句之前而不能用于类外部的函数定 ...
- H: Dave的组合数组(二分法)
Dave的组合数组 Time Limit: C/C++ 1 s Java/Python 3 s Memory Limit: 128 MB Accepted: 3 ...
- phoenix表操作
phoenix表操作 进入命令行,这是sqlline.py 配置到path环境变量的情况下 sqlline.py localhost如果要退出命令行:!q 或者 !quit 3.4.1 创建表 ...
- What's news in Visual Studio 2017
文字总结: 1.高级智能提示 在属性列表中输入 M C即可查询属性中包含字母m\c的属性 2.更快的导航查询,在Go To All中输入任意查询的字符,可快速查到任何包含关键字的文件 3.代码智能分 ...
- 解决 Ionic 浏览器跨域问题
<system.webServer> <httpProtocol> <customHeaders> <add name="Access-Contro ...
- Codeforces 594D REQ 线段树
REQ 把询问离线, 我们从n 到 1遍历过去的时候, 把(1 - 1 / p)乘在最靠近当前位置的地方, 然后区间求乘积就好啦. #include<bits/stdc++.h> #def ...
- 在vue-cli项目中使用bootstrap的方法示例
在一个html页面中加入bootstrap是很方便,就是一般的将css和js文件通过Link和Script标签就行. 那么在一个用vue-cli生成的前端项目中如何加入?因为框架不一样了,略微要适应一 ...