Lucas定理及组合数取模

首先给出这个Lucas定理：

A、B是非负整数，p是质数。AB写成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]) modp同余

即：Lucas(n,m,p)=c(n%p,m%p)*Lucas(n/p,m/p,p)

这个定理的证明不是很简单，我一直想找个很好的张明，但是，没找到，昨天看到了一个解题报告，基本上可以说明白这个Lucas定理是怎么回事了，具体的是说：

以求解n! % p为例，把n分段，每p个一段，每一段求的结果是一样的。但是需要单独处理每一段的末尾p, 2p, ...，把p提取出来，会发现剩下的数正好又是(n / p)!，相当于划归成了一个子问题，这样递归求解即可。

这个是单独处理n!的情况，当然C(n,m)就是n!/(m!*(n-m)!)，每一个阶乘都用上面的方法处理的话，就是Lucas定理了，注意这儿的p是素数是有必要的。

贴一个例题的代码，方便以后看。

hdu 3037

将不大于m颗种子存放在n颗树中，问有多少种存法。

首先是不大于m颗种子，我没可以认为少于m的那些种子存放在了第n+1颗树上，这样的话，问题就转化成了将m颗种子存放在n+1颗树上的方案数。ok这个是组合数学里面的公式，亦即插板法，也就是X1+X2+X3+……+Xn+1 = m;ok，答案是C(n+m,m);

然后就是上面说的Lucas定理解决大组合数问题了

#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;#define N 100010long long mod_pow(int a,int n,int p)
{long long ret=1;long long A=a;while(n){if (n & 1)ret=(ret*A)%p;A=(A*A)%p;n>>=1;}return ret;
}long long factorial[N];void init(long long p)
{factorial[0] = 1;for(int i = 1;i <= p;i++)factorial[i] = factorial[i-1]*i%p;//for(int i = 0;i < p;i++)//ni[i] = mod_pow(factorial[i],p-2,p);
}long long Lucas(long long a,long long k,long long p) //求C(n,m)%p p最大为10^5。a,b可以很大！
{long long re = 1;while(a && k){long long aa = a%p;long long bb = k%p;if(aa < bb) return 0; //这个是最后的改动！re = re*factorial[aa]*mod_pow(factorial[bb]*factorial[aa-bb]%p,p-2,p)%p;//这儿的求逆不可先处理a /= p;k /= p;}return re;
}int main()
{int t;cin >> t;while(t--){long long n,m,p;cin >> n >> m >> p;init(p);cout << Lucas(n+m,m,p) << "\n";}return 0;
}

//以上为前几届学长的lucas相关总结，在此基础上，做了少许的修改。

Lucas定理及组合数取模相关推荐

2015 ICL, Finals, Div. 1 Ceizenpok’s formula（组合数取模，扩展lucas定理）
J. Ceizenpok's formula time limit per test 2 seconds memory limit per test 256 megabytes input stand ...
组合数学 —— 组合数取模
[概述] 组合数取模,即计算组合数 ,由于 ,同余定理对除法不适用,因此需要使用别的方法来解决这个问题常见的方法有:使用逆元对组合数取模.递推打表取模.卢卡斯定理.扩展卢卡斯定理等,这些方法应用的场 ...
Codeforces 869C The Intriguing Obsession 组合数取模
Codeforces 869C The Intriguing Obsession 思考一下人生. 这是一场物语场,而且A题直接puts("Karen")能过,我对此印象非常深.我不 ...
组合数取模 Lucas定理
对于C(n, m) mod p.这里的n,m,p(p为素数)都很大的情况.就不能再用C(n, m) = C(n - 1,m) + C(n - 1, m - 1)的公式递推了. 这里用到Lusac定理 ...
组合数学 —— 组合数取模 —— 卢卡斯定理与扩展卢卡斯定理
[卢卡斯定理] 1.要求:p 是质数,m.n 很大但 p 很小或者 n.m 不大但大于 p 2.定理内容其中, 3.推论当将 n 写成 p 进制:,将 m 写成 p 进制: 时,有: 4.实现 ...
bzoj1951 组合数取模中国剩余定理
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int a[4]={2,3,4679,35 ...
组合数取模 - Lucas/exLucas - LibreOJ #181. 二项式系数
文章目录 Lucas定理证明拓展问题一些优化(LibreOJ #181. 二项式系数) 预处理部分复杂度分析大佬的玄学优化(取模) AC代码证明和代码分开,可以根据自己的需要跳转. 去我的 ...
大数组合数取模(逆元+打表)
将阶乘O(n)打表之后C(n,m)便可O(1)求出,除法取模用逆元解决 hdu5698瞬间移动 #include<bits/stdc++.h> using namespace std; c ...
BZOJ-1951 古代猪文（组合数取模Lucas+中国剩余定理+拓展欧几里得+快速幂）...
数论神题了吧算是 1951: [Sdoi2010]古代猪文 Time Limit: 1 Sec Memory Limit: 64 MB Submit: 1573 Solved: 650 [Submit ...

Lucas定理及组合数取模

Lucas定理及组合数取模相关推荐

最新文章

热门文章