java1e6,大组合数取模之lucas定理模板,1<=n<=m<=1e9,1<p<=1e6，p必须为素数

typedef long long ll;

/**********************************

大组合数取模之lucas定理模板,1<=n<=m<=1e9,1

输入：C(n,m)%p 调用lucas(n,m,p)

复杂度：min(m,p)*log(m)

***********************************/

//ax + by = gcd(a,b)

//传入固定值a,b.放回 d=gcd(a,b), x , y

void extendgcd(ll a,ll b,ll &d,ll &x,ll &y)

{

if(b==){d=a;x=;y=;return;}

extendgcd(b,a%b,d,y,x);

y-=x*(a/b);

}

//Ax=1(mod M)，gcd(A,M)==1

//输入：10^18>=A,M>=1

//输出：返回x的范围是[1,M-1]

ll GetNi(ll A,ll M)

{

ll rex=,rey=;

ll td=;

extendgcd(A,M,td,rex,rey);

return (rex%M+M)%M;

}

ll C(ll n,ll m,ll p)

{

if(m>n) return ;

ll up=,dn=;

for(int i=;i

{

up = up*(n-i)%p;

dn = dn*(i+)%p;

}

return up*GetNi(dn, p)%p;

}

ll lucas(ll n,ll m,ll p)

{

if(m==) return ;

return C(n%p,m%p,p)*lucas(n/p,m/p,p) % p;

}

组合数取模(lucas定理&plus;CRT合并)(AC)

#include #define re register #define int long long using namespace std; ; inlin ...

组合数取模介绍----Lucas定理介绍

转载https://www.cnblogs.com/fzl194/p/9095177.html 组合数取模方法总结(Lucas定理介绍) 1.当n,m都很小的时候可以利用杨辉三角直接求. C(n,m) ...

lucas定理解决大组合数取模

LL MyPow(LL a, LL b) { LL ret = ; while (b) { ) ret = ret * a % MOD; a = a * a % MOD; b >>= ; ...

HDU 5698 大组合数取模&lpar;逆元&rpar;

瞬间移动 Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submis ...

bzoj1951 组合数取模 中国剩余定理

#include using namespace std; typedef long long ll; const int a[4]={2,3,4679,35 ...

Codeforces 57C &lpar;1-n递增方案数，组合数取模，lucas&rpar;

这个题相当于求从1-n的递增方案数,为C(2*n-1,n); 取模要用lucas定理,附上代码: #include using namespace std; ...

组合数取模&amp&semi;&amp&semi;Lucas定理题集

题集链接: https://cn.vjudge.net/contest/231988 解题之前请先了解组合数取模和Lucas定理 A : FZU-2020  输出组合数C(n, m) mod p (1 ...

Uva12034 (组合数取模)

题意:两匹马比赛有三种比赛结果,n匹马比赛的所有可能结果总数 解法: 设答案是f[n],则假设第一名有i个人,有C(n,i)种可能,接下来还有f(n-i)种可能性,因此答案为 ΣC(n,i)f(n-i ...

组合数取模Lucas定理及快速幂取模

组合数取模就是求的值,根据,和的取值范围不同,采取的方法也不一样. 下面,我们来看常见的两种取值情况(m.n在64位整数型范围内) (1)  , 此时较简单,在O(n2)可承受的情况下组合数的计算可以 ...

随机推荐

开源服务专题之------sshd服务安装管理及配置文件理解和安全调优

本专题我将讨论一下开源服务,随着开源社区的日趋丰富,开源软件.开源服务,已经成为人类的一种公共资源,发展势头可谓一日千里,所以不可不知.SSHD服务,在我们的linux服务器上经常用到,很重要,涉及到 ...

springmvc restful配置有一个小小的坑坑

首先web.xml配置 springServlet

拖放 js

之前被小伙伴问自己能不能写一个简单的原生的 我稍微犹豫了下  这次重新学习下拖拽的过程  分享下  参考 JavaScript高级程序设计 必要的准备 自定义事件(实现事件模型)  简单来说事件模型就 ...

elixir 高可用系列(三) GenEvent

概述 GenEvent 是事件处理的通用部分的抽象. 通过 GenEvent ,我们给已有的服务 动态 的添加 事件处理. GenEevent 和 GenServer 的区别 之前已经介绍了 GenS ...

Android Sutido 编译速度优化

虽然Android Studio 此时已经更新到了Android Studio 2.1版本,build 版本android-studio-bundle-143.2739321.但是在安装该版本都是根据 ...

Linux之shell编程基础

一.变量 变量在shell中分为:本地变量.环境变量.位置参数: 本地变量:仅可在用户当前shell生命期的脚本中使用的变量,本地变量随着shell进程的消亡而无效,本地变量在新启动的shell中依旧 ...

JVM学习笔记二：垃圾收集算法

垃圾回收要解决的问题: 哪些内存需要回收? 线程私有区域不需要回收,如PC.Stack.Native Stack:Java 堆和方法区需要 什么时候回收? 以后的文章解答 如何回收? 首先进行对象存活 ...

C&plus;&plus;笔记004：C&plus;&plus;类通俗点说

核心: C++的类就是对C语言的结构体进行了扩展,C++的结构体可以包含函数! ------------------------------------------------------ 我们学习C ...

Linux查看cpu个数

[root@lidongbo~]# cat /proc/cpuinfo processor       : 0 vendor_id       : GenuineIntel cpu family    ...

css3新增的属性

由于CSS5标准还未完全订下来,所以各种内核的浏览器都有自己的标准,为了不使属性混淆,所以各家在各自标准前加了一个前缀, 如:-moz-       firefox火狐 -ms-         IE ...