博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷 P4317 花神的数论题(组合数)
阅读量:5224 次
发布时间:2019-06-14

本文共 1512 字,大约阅读时间需要 5 分钟。

题面

题解

组合数

枚举有多少个\(1\),求出有多少种数

扫描\(n\)的每一位\(1\), 强制选\(0\)然后组合数算一下有多少种方案

Code

#include
#define LL long long#define RG registerusing namespace std;template
inline void read(T &x) { x = 0; RG char c = getchar(); bool f = 0; while (c != '-' && (c < '0' || c > '9')) c = getchar(); if (c == '-') c = getchar(), f = 1; while (c >= '0' && c <= '9') x = x*10+c-48, c = getchar(); x = f ? -x : x; return ;}template
inline void write(T x) { if (!x) {putchar(48);return ;} if (x < 0) x = -x, putchar('-'); int len = -1, z[20]; while (x > 0) z[++len] = x%10, x /= 10; for (RG int i = len; i >= 0; i--) putchar(z[i]+48);return ;}const int N = 55, Mod = 10000007;LL C[N][N], cnt, a[N];inline LL Pow(int a, LL b) { LL s = 1; for (LL x = a; b; b >>= 1, x = x*x%Mod) if (b&1) s = s*x%Mod; return s;}int main() { //freopen(".in", "r", stdin); //freopen(".out", "w", stdout); for (int i = 0; i <= 50; i++) C[i][i] = C[i][0] = 1; for (int i = 2; i <= 50; i++) for (int j = 1; j < i; j++) C[i][j] = C[i-1][j-1]+C[i-1][j]; LL n; read(n); for (int i = 50; i >= 0; i--) { if ((n>>i)&1) { for (int j = 1; j <= i; j++) a[j+cnt] += C[i][j];//至少选1个的方案 a[++cnt]++;//不选的方案(必须分开算,不然会有重复计算) } } LL ans = 1; for (int i = 1; i <= 50; i++) ans = ans*Pow(i, a[i])%Mod; printf("%lld\n", ans); return 0;}

转载于:https://www.cnblogs.com/zzy2005/p/10223743.html

你可能感兴趣的文章
手机网页支付
查看>>
hdu_5862_Counting Intersections(扫描线)
查看>>
1442: Neo 的简单字符串(字符串)
查看>>
241. Different Ways to Add Parentheses
查看>>
android 线程通信Handler Bundle
查看>>
JVM内存管理 + GC垃圾回收机制
查看>>
HTML5快速写页面的方法
查看>>
js短路表达式
查看>>
storm1.1运行时问题
查看>>
storm排错
查看>>
如何保持工作热情
查看>>
css代码实现
查看>>
[RPi浙大yq向]设置RPi的mac地址,静态ip地址,dns服务器地址等。比较普通的上内网的方法。...
查看>>
struts2(七)之s标签和#、$、%d的使用
查看>>
Swift中类的两段式构造(类的构造过程)
查看>>
bzoj1013 [JSOI2008]球形空间产生器sphere
查看>>
linux 得到内网外网ip
查看>>
简单的移动端页面
查看>>
DTO的一些理解(转载)
查看>>
mysql left join,right join,inner join的区别
查看>>