大意: 给定$n,a$, 求$n$个$3$的倍数, $or$和为$a$的方案数.

简单容斥题

可以求出$f_{x,y}$表示所有$3$的倍数中, 奇数位不超过$x$个$1$, 偶数位不超过$y$个$1$的个数.

假设$a$二进制奇数位$c_1$个$1$,偶数位$c_0$个$1$, 根据容斥就有

$ans=\sum\limits_{i=0}^{c_1}\sum\limits_{j=0}^{c_0}(-1)^{c_0+c_1-i-j}\binom{c_1}{i}\binom{c_0}{j}f_{i,j}^n$

#include <iostream>
#include <sstream>
#include <algorithm>
#include <cstdio>
#include <math.h>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <string.h>
#include <bitset>
#define REP(i,a,n) for(int i=a;i<=n;++i)
#define PER(i,a,n) for(int i=n;i>=a;--i)
#define hr putchar(10)
#define pb push_back
#define lc (o<<1)
#define rc (lc|1)
#define mid ((l+r)>>1)
#define ls lc,l,mid
#define rs rc,mid+1,r
#define x first
#define y second
#define io std::ios::sync_with_stdio(false)
#define endl '\n'
#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int P = 998244353, INF = 0x3f3f3f3f;
ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}
ll qpow(ll a,ll n) {ll r=1%P;for (a%=P;n;a=a*a%P,n>>=1)if(n&1)r=r*a%P;return r;}
ll inv(ll x){return x<=1?1:inv(P%x)*(P-P/x)%P;}
inline int rd() {int x=0;char p=getchar();while(p<'0'||p>'9')p=getchar();while(p>='0'&&p<='9')x=x*10+p-'0',p=getchar();return x;}
//headconst int N = 63;
int C[N][N],f[N][N];
int main() {REP(i,0,N-1) {C[i][0] = 1;REP(j,1,i) C[i][j]=(C[i-1][j]+C[i-1][j-1])%P;}REP(i,0,N-1) REP(j,0,N-1) REP(ii,0,i) REP(jj,0,j) if ((ii+2*jj)%3==0) {f[i][j] = (f[i][j]+(ll)C[i][ii]*C[j][jj])%P;}int t;scanf("%d", &t);while (t--) {ll n, a;scanf("%lld%lld", &n, &a);n %= P-1;int c[2]{};REP(i,0,N-1) if (a>>i&1) ++c[i&1];int ans = 0;REP(i,0,c[0]) REP(j,0,c[1]) {int ret = (ll)C[c[0]][i]*C[c[1]][j]%P*qpow(f[i][j],n)%P;if (c[0]+c[1]+i+j&1) ret = P-ret;ans = (ans+ret)%P;}printf("%d\n", ans);}
}

转载于:https://www.cnblogs.com/uid001/p/11258406.html

2019牛客多校四 E. triples II (容斥)相关推荐

  1. 2019牛客国庆集训派对day2 K 2018(容斥)

    链接:https://ac.nowcoder.com/acm/contest/1107/K 来源:2019牛客国庆集训派对day2 题目描述   Given a, b, c, d, find out ...

  2. 2019牛客多校第一场

    2019牛客多校第一场 题号 题目 知识点 A Monotonic Matrix B Symmetric Matrix C Fluorescent 2 D Two Graphs E Removal F ...

  3. 2019牛客多校训练第十场F Popping Balloons

    2019牛客多校训练第十场F Popping Balloons 题意:二维平面内给你若干个点,然后你可以在x轴和y轴分别射三枪(每一枪的间隔是R),问最多能射掉多少气球. 题解:贪心.这个应该只能算作 ...

  4. 2019牛客多校第九场AThe power of Fibonacci(广义BM)

    2019牛客多校第九场AThe power of Fibonacci(广义BM) 题目大意 求斐波那契数列m次方的前n项和 解题思路 显然,斐波那契的m次方前缀和依然是线性递推,因此考虑用exBM求解 ...

  5. 2019牛客多校第四场 I题 后缀自动机_后缀数组_求两个串de公共子串的种类数

    目录 求若干个串的公共子串个数相关变形题 对一个串建后缀自动机,另一个串在上面跑同时计数 广义后缀自动机 后缀数组 其他:POJ 3415 求两个串长度至少为k的公共子串数量 @(牛客多校第四场 I题 ...

  6. Knapsack Cryptosystem(2019牛客多校折半查询)

    链接:https://ac.nowcoder.com/acm/contest/889/D 来源:牛客网 Amy asks Mr. B problem D. Please help Mr. B to s ...

  7. Quadratic equation(二次剩余)2019牛客多校第九场

    链接:https://ac.nowcoder.com/acm/contest/889/B 来源:牛客网 题目描述 Amy asks Mr. B problem B. Please help Mr. B ...

  8. [2019牛客多校训练第3场]Median

    链接:https://ac.nowcoder.com/acm/contest/883/I 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言65536 ...

  9. 2019牛客多校训练营第一场 H题 HOR 题解

    题目描述: 输入描述: 输出描述: 示例1: 题解: 更多问题可关注牛客竞赛区,一个刷题.比赛.分享的社区. 传送门:https://ac.nowcoder.com/acm/contest/discu ...

最新文章

  1. Microbiome:山大杜宗军团队揭示捕食性细菌新类群-慢生单胞菌目细菌独特的生境适应性...
  2. mysql adodb_关于从MySQL转向ADODB的方法
  3. jzoj5123-diyiti【统计,容斥】
  4. “睡服”面试官系列第二十一篇之class基本语法(建议收藏学习)
  5. Eclipse快捷键 10个最有用的快捷键(转)
  6. 企业发展如何启动云的力量
  7. 7-3 奇偶分家 (10 分)
  8. 华为那道上机题 用?来区隔一句话的输入。
  9. 181007扇贝有道每日一句
  10. oracle可视化操作界面——plsql dev安装配置与使用
  11. 2019年软件评测师考试大纲
  12. 数据库基础笔试题(附带答案)
  13. 允许计算机usb调试,usb调试不弹出授权,电脑一直弹出无法识别USB
  14. matlab仿真没有synchr,Synchro交通仿真系统分析及应用
  15. 【内部流出,勿做商用】今日教大家如何抢聚划算
  16. Python 操作pdf文件-加密与解密 (二)
  17. [DAX] DIVIDE函数
  18. Android实现高德地图POI搜索
  19. 第二本书 第一章 Linux中的网络配置
  20. 计算 11+12+13+...+m python

热门文章

  1. 关于CleanMyMac常见问题与解答
  2. P2240 【深基12.例1】部分背包问题
  3. 项目管理中网络图的看法和相关参数阅读说明
  4. RTX5 | 内存池04 - 共享内存用于线程之间的通讯(阻塞方式)- 使用信号量
  5. STM32H743+CubeMX-ADC(16bit分辨率)+DMA采样三路模拟量,硬件过采样器实现1024倍过采样
  6. STM32H743+Cube-Keil上移植RTX5实时系统
  7. 操作系统课设之虚拟内存页面置换算法的模拟与实现
  8. input输入框匹配数字(带两位小数) js--金额匹配
  9. Linux笔记-iptables规则原理和组成
  10. 信息安全工程师笔记-10种端口扫描技术概念