E-triples II_2019牛客暑期多校训练营(第四场)
求用n个3的倍数的数按位或出数字a的方案数有多少种(0也算3的倍数)
题解
- 若数b的每个二进制位上的1,在a中也为1,则称b为a的子集
- 容易知道任意个a的子集按位或出来的结果还是a的子集
- 若问题改为按位或出来的结果是a的子集的方案数,那么答案就是a的子集中是3的倍数的子集个数的n次方
接着我们对子集按二进制上的1 mod 3的个数划分,例如1101有两个1mod3=1, 一个1mod3 = 2,设\(S[i][j]\)表示a的子集中有i个mod3=1,j个mod3=2的子集的子集 中是3的倍数的个数,例如a = 1101的一个子集1001表示的状态为\(S[1][1]\), 1001的子集中是3的倍数的有1001和0000所以\(S[1][1] = 2\),那么\(S[i][j]\)的n次方就可以表示为用n个3的倍数的数按位或出来的结果的状态是S[i][j]的子集方案数
那么\(\sum_{i=1}^kS[i][k-i]\)就表示或出来的结果最多匹配上a中K个1的方案数,那么我们就可以用最多匹配上a中K个1的方案数,减去匹配上a中K-1个1的方案数得出答案,但是这样简单的相减是不行的因为\(S[i][k-i]\)的子集是会有重叠的,会多扣掉最多匹配k-2个1的方案数,根据容斥原理应当减去最多匹配K-1的方案数,加上最多匹配K-2的方案数,扣掉K-3加上K-4...
代码
#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int mx = 65;
const ll mod = 998244353;
int C[mx][mx], S[mx][mx];ll pow_mod(ll a, ll b) {ll ans = 1;while (b > 0) {if (b & 1) ans = ans * a % mod;a = a * a % mod;b /= 2;}return ans;
}int main() {C[0][0] = 1;for (int i = 1; i < mx; i++) {C[i][0] = 1;for (int j = 1; j <= i; j++) {C[i][j] = (C[i-1][j-1] + C[i-1][j]) % mod;}}for (int i = 0; i < mx; i++) {for (int j = 0; j < mx; j++) {for (int p = 0; p <= i; p++) {for (int q = 0; q <= j; q++) {if ((p + 2*q) % 3 != 0) continue;S[i][j] += C[i][p] * C[j][q] % mod;S[i][j] %= mod;}}}}S[0][0] = 1;int T;scanf("%d", &T);while (T--) {ll n, a, x = 0, y = 0;scanf("%lld%lld", &n, &a);for (int i = 0; i < 64; i++) {if (a & (1LL<<i)) {if (i % 2 == 0) x++;else y++;}}ll ans = 0;for (int i = 0; i <= x; i++) {for (int j = 0; j <= y; j++) {ll tmp = C[x][i] * C[y][j] % mod * pow_mod(S[i][j], n) % mod;if ((x+y-i-j) % 2) tmp *= -1;ans = (ans + tmp) % mod;}}ans = (ans + mod) % mod;printf("%lld\n", ans);}return 0;
}
转载于:https://www.cnblogs.com/bpdwn-cnblogs/p/11289740.html
E-triples II_2019牛客暑期多校训练营(第四场)相关推荐
- 2020牛客暑期多校训练营(第一场)
文章目录 A B-Suffix Array B Infinite Tree C Domino D Quadratic Form E Counting Spanning Trees F Infinite ...
- 2020牛客暑期多校训练营(第二场)
2020牛客暑期多校训练营(第二场) 最烦英语题 文章目录 A All with Pairs B Boundary C Cover the Tree D Duration E Exclusive OR ...
- E Groundhog Chasing Death(2020牛客暑期多校训练营(第九场))(思维+费马小定理+质因子分解)
E Groundhog Chasing Death(2020牛客暑期多校训练营(第九场))(思维+费马小定理+质因子分解) 链接:https://ac.nowcoder.com/acm/contest ...
- 【2019牛客暑期多校训练营(第二场) - H】Second Large Rectangle(单调栈,全1子矩阵变形)
题干: 链接:https://ac.nowcoder.com/acm/contest/882/H 来源:牛客网 题目描述 Given a N×MN \times MN×M binary matrix. ...
- 2019牛客暑期多校训练营(第一场) A Equivalent Prefixes ( st 表 + 二分+分治)
链接:https://ac.nowcoder.com/acm/contest/881/A 来源:牛客网 Equivalent Prefixes 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/ ...
- 【2019牛客暑期多校训练营(第二场)- E】MAZE(线段树优化dp,dp转矩阵乘法,线段树维护矩阵乘法)
题干: 链接:https://ac.nowcoder.com/acm/contest/882/E?&headNav=acm 来源:牛客网 Given a maze with N rows an ...
- 【2019牛客暑期多校训练营(第二场)- F】Partition problem(dfs,均摊时间优化)
题干: 链接:https://ac.nowcoder.com/acm/contest/882/F 来源:牛客网 Given 2N people, you need to assign each of ...
- 【2019牛客暑期多校训练营(第二场) - D】Kth Minimum Clique(bfs,tricks)
题干: 链接:https://ac.nowcoder.com/acm/contest/882/D 来源:牛客网 Given a vertex-weighted graph with N vertice ...
- 【2019牛客暑期多校训练营(第一场) - A】Equivalent Prefixes(单调栈,tricks)
题干: 链接:https://ac.nowcoder.com/acm/contest/881/A 来源:牛客网 Two arrays u and v each with m distinct elem ...
- 【2019牛客暑期多校训练营(第一场) - H】XOR(线性基,期望的线性性)
题干: 链接:https://ac.nowcoder.com/acm/contest/881/H 来源:牛客网 Bobo has a set A of n integers a1,a2,-,ana1, ...
最新文章
- Go语言调度器之调度main goroutine(14)
- linux redis ruby,redis requires ruby version 2.2.2的解决方案
- POJ 3580 SuperMemo
- Codeigniter 用户登录注册模块
- STC12C5A60S2笔记8(串口)
- ajax 赋值return
- 信息学奥赛一本通 2050:【例5.20】字串包含 | OpenJudge NOI 1.17 19:字符串移位包含问题
- 柴犬为什么总是被卡住狗头?
- 【Elasticsearch】分片未分配 (UNASSIGNED) failed to obtain in-memory shard lock
- python spark2.0_Python+Spark2.0+hadoop学习笔记——Python Spark MLlib决策树二分类
- MySQL数据库 -- Navicat、pycharm连接数据库
- CoffeeScript、TypeScript 、 JavaScript 之间
- 阻塞非阻塞、同步异步
- 程序员简洁简历模板分享
- 艾迈斯半导体推出新款高性能读取IC,推动医疗和工业数字化X射线设备制造商降本增效
- c++::举例_举例说明:网络
- 前端工程师的摸鱼日常(2)
- 2017计科01-08编译原理模拟测试2--chap03
- 拼多多第二季度表现亮眼,“贪便宜”的你贡献了多少?
- Java工具使用(导读)
热门文章
- 纯虚函数--抽象类中的this到底指的是哪个子类实例呢
- [react] 经常用React,你知道React的核心思想是什么吗?
- React开发(229):react删除的实现
- React开发(189):ant design textarea允许清除allowClear
- React开发(155):请求方式得问题
- 前端学习(3071):vue+element今日头条管理-删除文章失败
- [html] 举例说明当我们在写布局时,都有哪些边界的情况需要关注的?
- [html] 你知道什么是粘性布局吗?
- 工作377-处理url拼接里面的参数方法
- [vue] 在vue事件中传入$event,使用e.target和e.currentTarget有什么区别?