大整数求组合数取余(Lucas定理)
来源:互联网 发布:淘宝网开店需要什么 编辑:程序博客网 时间:2024/05/21 08:54
【卢卡斯(Lucas)定理】
Lucas定理用来求C(a,b)mod p的值,其中p为素数。
数学表达式为:
Lucas(a,b,q)=C(a%q,b%q)*Lucas(a/p,b/p,p);
Lucas(a,0,q)=1;
通过这个定理就可以很方便的把大数的组合转化成小数。但其中还是要求C(a%q,b%q)%p,所以这里引入逆元来求。
【定义】若整数a,b,p, 满足a·b≡1(mod p).则称a 为b 模p 的乘法逆元, 即a=b- 1mod p.其中, p 是模数。
应用到组合数中来就是:
a!/[b!*(a-b)!] % p == a! * [b!*(a-b)!]-1 %p
【逆元求法】:
应用费马小定理,ap-1=1 mod p ,即 a*ap-2=1 mod p
也就是说 ap-2就是a的逆元。
当然这里求出来的逆元是在取模p的逆元,对我们最终目标没有影响。这也是比较方便而且比较好的方法。
弄个模板来。
#include <cstdio> #include <iostream> #include <cmath> #include <cstring> #include <algorithm> using namespace std; #define maxn 100010 typedef long long LL; LL m,n,p; LL Pow(LL a,LL b,LL mod) { LL ans=1; while(b) { if(b&1) { b--; ans=(ans*a)%mod; } else { b/=2; a=(a*a)%mod; } } return ans; } LL C(LL n,LL m) { if(n<m) return 0; LL ans=1; for(int i=1;i<=m;i++) { ans=ans*(((n-m+i)%p)*Pow(i,p-2,p)%p)%p; } return ans; } LL Lucas(LL n,LL m) { if(m==0) return 1; return (Lucas(n/p,m/p)*C(n%p,m%p))%p; } int main() { int t; scanf("%d",&t); while(t--) { scanf("%lld%lld%lld",&n,&m,&p); printf("%lld\n",Lucas(n,m)); } return 0; }
- 大整数求组合数取余(Lucas定理)
- 大整数组合数求余(lucas定理)
- hdu-3037(lucas定理求大组合数取模)
- 求大组合数取模的模板(lucas定理)
- HDU3037(Lucas定理求大组合数取模)
- Lucas 定理 大组合数取余 【记录】
- Lucas 定理 大组合数取余 【学习】
- Lucas定理(大组合数取模)
- hdu3037 隔板法+Lucas定理求大组合取模
- Lucas定理求组合数
- Lucas定理 大组合数取模
- FZU 2020 组合 -- (大组合数取模 Lucas定理)
- 费马小定理+Lucas定理 大组合数取模
- Lucas定理and大组合数取余算法总结
- JOJ1595:Combinations(求解大组合数,Lucas定理)
- hdu 3037(Lucas定理,大组合数取模)
- Lucas定理(大组合数的取模)
- hdu 3037 Saving Beans (大组合数取模--Lucas定理)
- _bstr_t 类型,BSTR类型和CString类型到底有什么区别
- 转:java 基础知识
- ural 1029. Ministry
- 解决Thinkpad笔记本安装WIN8系统后Power Manager电源管理软件无法使用的问题
- 我的程序都是这样命名的:openeim001
- 大整数求组合数取余(Lucas定理)
- vs2012安全函数 error C4996: 'scanf':may be unsafe. use _CRT_SECURE_NO_WARNINGS
- VC++ MCI多媒体播放原理
- JDBC更新CLOB字段
- linux配置oracle客户端,sqlplus、sqlldr、exp、imp
- VS 2012 显示Link的参数
- java设计模式笔记之抽象工厂
- JNDI的两种书写方式的不同
- STL之list