lucas定理（什么是Lucas定理）

本文目录

• 什么是Lucas定理
• 谁有Lucas定理+逆元解决大组合数求模的可运行matlab代码
• 概率c公式是什么 c表示什么
• 大数组合取模,Lucas定理,费马小定理的运用
• lucas定理中p为什么要为素数

什么是Lucas定理

• LUCAS就是光。光的定理。

• 陈冠希+张柏芝=Lucas

• 就是谢霆锋他儿子。

谁有Lucas定理+逆元解决大组合数求模的可运行matlab代码

#include 《stdio.h》#include 《math.h》#include 《string.h》#include 《stdlib.h》#include 《iostream》#include 《sstream》#include 《algorithm》#include 《set》#include 《queue》#include 《stack》#include 《map》#include 《bitset》#pragma comment(linker, "/STACK:102400000,102400000")using namespace std;typedef __int64 LL;c***t int inf=0x3f3f3f3f;c***t double pi= acos(-1.0);c***t double esp=1e-6;using namespace std;LL n,m,mod;LL modxp(LL a,LL b){ LL res=1; while(b》0) { if(b&1) res=res*a%mod; b=b》》1; a=a*a%mod; } return res;}LL C(LL n, LL m){ if(m》n) return 0; LL ans=1; for(int i=1; i《=m; i++) { LL a=(n+i-m)%mod; LL b=i%mod; ans=ans*(a*modxp(b, mod-2)%mod)%mod; } return ans;}LL Lucas(LL n,LL m){ if(m==0) return 1; return C(n%mod,m%mod)*Lucas(n/mod,m/mod)%mod;}int main(){ int T; scanf("%d",&T); while(T--) { scanf("%lld %lld %lld",&n,&m,&mod); printf("%lld\n",Lucas(n,m)); } return 0;}

概率c公式是什么 c表示什么

概率公式：C(n,k)=n(n-1)(n-2)...(n-k+1)/k!，其中k≤n。C表示组合数。

概率公式是什么 c表示什么

C表示组合数。

C(n,m) 表示n选m的组合数，其中n是下标 , m是上标 (C上面m,下面n)。

nCk是一个整体，是n个元素中，取k个元素的取法的个数，也叫n个元素中，取k

个k组合数，（C代表组合），算法是：

nCk＝n!/k!（n-k）!＝n（n-1）……（n-k+1）/k!

等于从n开始连续递减的m个自然数的积除以从1开始连续递增的m个自然数的积。

该概率公式的推导过程：

在这个证明中，表示n次实验中，成功的k次，取法的个数。

每次取定后，k次成功，n-k次失败，概率用乘法P＝p^k*(1-p)^(n-k)

总共有nCk个取法，即nCk个情况，概率用加法，每个情况的概率又相同，所以

成为nCk倍。

求组合数C的方法

1、当n,m都很小的时候可以利用杨辉三角直接求。

C(n,m)=C(n-1,m)+C(n-1,m-1)；

2、利用乘法逆元

乘法逆元：(a/b)%mod=a*(b^(mod-2)) mod为素数。

逆元可以利用扩展欧几里德或欧拉函数求得。

3、当n和m比较大，mod是素数且比较小的时候（10^5左右），通过Lucas定理计算

大数组合取模,Lucas定理,费马小定理的运用

从一个例题：【HDU 3037】 Saving Beans 来开始Lucas定理的应用。 题目大意为：松鼠要从n棵树上摘一共m个豆子，结果的方案数对素数p（不大于1e5）取模，求解。 思路： 可以理解为m个豆子分为n份，求分的方法个数。 由插板法来对m个数进行划分，由于可能某棵树没有摘豆子，可以理解为：x1+x2+x3+……+xn=m的解的个数，即为C(m+n-1,n-1)。（将m颗豆子加上n-1个板子的位置，得到的序列再从中取n-1个板子的位置）=C(m+n-1,m)。 由于m的值取0~m，那么就得sum=C(n-1,0)+C(n,1)+C(n+1,2)+C(n+2,3)+……+C(m+n-1,m)。 利用公式C（n，r）=C（n-1，r）+C（n-1，r-1）=C（n-1，r）+C（n-2，r-1）+C（n-3，r-2）…… sum=C（n+m，m）。 也就是说，接下来的算法变成了C（n+m，m）%p。 然后就是Lucas定理的运用： Lucas（m，n，p）=C（m%p，n%p，p）✲Lucas（m/p，n/p，p）。 Lucas（x，0，p）=1。 这里可以采用的方法是递归求解。 简单的理解就是： 以求解n! % p 为例，把n分段，每p个一段，每一段求得结果是一样的。但是需要单独处理每一段的末尾p,2p,...,把p提取出来，会发现剩下的数正好又是(n/p)! ，相当于划归了一个子问题，这样递归求解即可。 这个是单独处理n！的情况，当然C(n,m)就是n！/(m! *（n－m）！），每一个阶乘都用上面的方法处理的话，就是Lucas定理了. Lucas最大的数据处理能力是p在10^5左右。 而C(a,b) =a! / ( b! ✲ (a-b)! ) mod p 其实就是求 ( a! / (a-b)!) ✲ ( b! )^(p-2) mod p

上面这一步变换是根据费马小定理：假如p是质数，且a,p互质，那么a的（p-1）次方除以p的余数恒为1, 那么a和a^(p-2)互为乘法逆元，则（b / a) = (b * a^(p-2) ) mod p)

b!与b! (p-2)互为乘法逆元，即b!✲b! (p-2)=1，那么，

//快速幂a^b % k

//求C(n, m)%p p最大为10^5 n, m可以很大！

用下面的Lucas定理程序实现就能得出结果，实现过程中要注意乘法时的强制转换

lucas定理中p为什么要为素数

加快速幂解决。 #include #include typedef __int64 lld; c***t lld MAX=500005; bool ok,c=0; lld count(lld n,lld prime) { lld ret=0; while(n/prime) { ret+=n/prime; n/=prime; } return ret; ...