博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
FZU 2020-组合(Lucas定理+逆元解决大组合数求模)
阅读量:2029 次
发布时间:2019-04-28

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

题目地址:

题意:求C(n,m)%p的值(1 <= m <= n <= 10^9, m <= 10^4, m < p < 10^9, p是素数)。
思路:
对于这里写图片描述这里写图片描述并且p是素数,我们一般采用Lucas定理来解。
1).Lucas定理是用来求 C(n,m) mod p的值,p是素数。其描述为:
如果
这里写图片描述
那么得到
这里写图片描述
Lucas(n,m,p)=C(n%p,m%p)* Lucas(n/p,m/p,p)
Lucas(n,0,p)=1;
2).对于大组合数求模C(N,M)%P=N! / (M! * (N-M)! ) % mod
=( N-M+i )! / M!*(i>=1&&i<=M)。然后根据乘法逆元将除法变成乘法即可。

#include 
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#pragma comment(linker, "/STACK:102400000,102400000")using namespace std;typedef __int64 LL;const int inf=0x3f3f3f3f;const double pi= acos(-1.0);const 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;}
你可能感兴趣的文章
valgrind 的使用简介
查看>>
C++中虚继承的作用
查看>>
构建高并发高可用的电商平台架构实践
查看>>
服务器架构设计,常见问题分析
查看>>
QQ游戏百万人同时在线服务器架构实现
查看>>
几种经典的网络服务器架构模型的分析与比较
查看>>
Hadoop学习路线图
查看>>
谷歌技术"三宝"之MapReduce
查看>>
HADOOP中HDFS工作原理
查看>>
Design TinyURL 设计短网址系统
查看>>
Boost智能指针——shared_ptr
查看>>
Java HashMap工作原理及实现
查看>>
干货|Tomcat 连接数与线程池详解
查看>>
Spring @Configuration 和 @Component 区别
查看>>
MySQL常见的三种存储引擎
查看>>
Executors与ThreadPoolExecutor(阿里发布的 Java开发手册中强制线程池不允许使用 Executors 去创建)
查看>>
Java 并发编程:线程池的使用
查看>>
Java并发编程:CountDownLatch、CyclicBarrier和 Semaphore
查看>>
Spring线程池ThreadPoolTaskExecutor配置及详情
查看>>
【原创】分布式之redis复习精讲
查看>>