大组合数取模 hdu 3037 Saving Beans lucas定理
来源:互联网 发布:晋业进销存软件8.25 编辑:程序博客网 时间:2024/06/12 01:40
from yangboy
Lucas定理
A、B是非负整数,p是质数。A B写成p进制:A=a[n]a[n-1]...a[0],B=b[n]b[n-1]...b[0]。
则组合数C(A,B)与C(a[n],b[n])*C(a[n-1],b[n-1])*...*C(a[0],b[0]) mod p同余
即:Lucas(n,m,p)=C(n%p,m%p)*Lucas(n/p,m/p,p)
//快速幂a^b % k
- ll PowerMod(ll a, ll b, ll k) {
- ll tmp = a, ret = 1;
- while (b) {
- if (b & 1) ret = ret * tmp % k;
- tmp = tmp * tmp % k;
- b >>= 1;
- }
- return ret;
- }
//求C(n, m)%p p最大为10^5 n, m可以很大!
- ll Lucas(ll n, ll m, ll p) {
- ll ret = 1;
- while (n && m) {
- ll nn = n%p, mm = m%p;
- if (nn < mm) return 0;
- //fac[nn]为预处理的 fac[n] = n!%p
- ret = ret*fac[nn]*PowerMod(fac[mm]*fac[nn-mm]%p, p-2, p)%p;
- n /= p;
- m /= p;
- }
- return ret;
- }
- //C(n, m) % p
- Lucas(n, m, p);
AC代码:
- #include <cstdio>
- #include <algorithm>
- #include <cmath>
- #include <iostream>
- using namespace std;
- typedef long long ll;
- ll fac[100003];
- void init(ll p) {
- fac[0] = 1;
- for (int i=1; i<=p; i++) fac[i] = fac[i-1]*i%p;
- }
- ll PowerMod(ll a, ll b, ll k) {
- ll tmp = a, ret = 1;
- while (b) {
- if (b & 1) ret = ret * tmp % k;
- tmp = tmp * tmp % k;
- b >>= 1;
- }
- return ret;
- }
- ll Lucas(ll n, ll m, ll p) {
- ll ret = 1;
- while (n && m) {
- ll nn = n%p, mm = m%p;
- if (nn < mm) return 0;
- ret = ret*fac[nn]*PowerMod(fac[mm]*fac[nn-mm]%p, p-2, p)%p;
- n /= p;
- m /= p;
- }
- return ret;
- }
- int main() {
- int T;
- ll n, m, p;
- cin >> T;
- while (T--) {
- cin >> n >> m >> p;
- init(p);
- cout << Lucas(n+m, m, p) << endl;
- }
- return 0;
- }
- hdu 3037 Saving Beans 【大组合数取模-Lucas定理+逆元取模】
- 大组合数取模 hdu 3037 Saving Beans lucas定理
- hdu 3037 Saving Beans (大组合数取模--Lucas定理)
- hdu-3037-Saving Beans(Lucas定理+大组合数取模)
- HDU 3037 Saving Beans 大组合数 lucas定理
- [ACM] hdu 3037 Saving Beans (Lucas定理,组合数取模)
- HDU 3037 Saving Beans 组合数取模 Lucas定理
- HDU - 3037 Saving Beans (数论,组合数取模,lucas定理)
- hdu 3037 Saving Beans(组合数学+lucas定理)
- HDU 3037 Saving Beans(组合数学+Lucas定理)
- hdu 3037 Saving Beans(Lucas定理)
- HDU 3037 Saving Beans (Lucas定理)
- HDU 3037 Saving Beans(lucas定理)
- HDU 3037 Saving Beans (Lucas定理)
- [Lucas定理] hdu 3037 Saving Beans
- hdu 3037 saving beans (lucas定理)
- HDU 3037 Saving Beans (Lucas定理)
- hdu 3037 Saving Beans lucas定理
- .在二元树中找出和为某一值的所有路径
- SQL SERVER 创建GHUID命令:select newid();
- 浅析php中常量,变量的作用域和生存周期
- php动态生成JavaScript代码
- /etc/sysctl.conf 参数说明
- 大组合数取模 hdu 3037 Saving Beans lucas定理
- OpenCV2计算机视觉编程手册(1)
- Eclipse/MyEclipse快捷键失效
- Ubuntu 12.04 成功编译安装 Erlang R16B01
- 数据库建表的各徙关系-- 一对多/多对一/一对一/多对多 关系分析
- H3C交换机802.1x用户的RADIUS AAA方案配置示例
- java基础12:线程间通信----等待唤醒机制
- MyEclipse 启动报错:'Building workspace' has encountered a problem解决方法
- WPF Application启动界面设置