[数学] Lucas定理模板
2015-08-11 21:43
274 查看
今天做多校联合的时候有一个问题需要用到大组合数计算,看起来以后可能会用到卢卡斯定理,所以把这部分的内容记录一下。
long long F[100010];
void init(long long p)
{
F[0] = 1;
for(int i = 1;i <= p;i++)
F[i] = F[i-1]*i % (1000000007);
}
long long inv(long long a,long long m)
{
if(a == 1)return 1;
return inv(m%a,m)*(m-m/a)%m;
}
long long Lucas(long long n,long long m,long long p)
{
long long ans = 1;
while(n&&m)
{
long long a = n%p;
long long b = m%p;
if(a < b)return 0;
ans = ans*F[a]%p*inv(F[b]*F[a-b]%p,p)%p;
n /= p;
m /= p;
}
return ans;
}
Cmn%p
C(n,m)=C([n/p],[m/p])∗C(n%p,m%p)%(modp)
long long F[100010];
void init(long long p)
{
F[0] = 1;
for(int i = 1;i <= p;i++)
F[i] = F[i-1]*i % (1000000007);
}
long long inv(long long a,long long m)
{
if(a == 1)return 1;
return inv(m%a,m)*(m-m/a)%m;
}
long long Lucas(long long n,long long m,long long p)
{
long long ans = 1;
while(n&&m)
{
long long a = n%p;
long long b = m%p;
if(a < b)return 0;
ans = ans*F[a]%p*inv(F[b]*F[a-b]%p,p)%p;
n /= p;
m /= p;
}
return ans;
}
C(n,m)=C([n/p],[m/p])∗C(n%p,m%p)%(modp)
相关文章推荐
- iOS自动打包流程
- Android 编程 权威 指南 阅读 笔记
- Hibernate之Session对象的相关方法以及持久化对象的状态
- Hibernate之Session对象的相关方法以及持久化对象的状态
- NBU基础知识
- atitit。流程图的设计与制作 attilax 总结
- 克鲁斯卡尔算法 畅通工程再续 Truck History 畅通工程
- rsync 安装使用详解
- assert的用法
- atitit。流程图的设计与制作 attilax 总结
- atitit。流程图的设计与制作 attilax 总结
- Windows10
- hdu 4061 福州赛区网络赛A 数学 ***
- go-内存分配
- java中String对象的常用方法
- C++ 文件结束符
- iOS常用第三方类库
- pwn中的一个小问题
- POJ 1789 Truck History
- Java学习整理系列之ThreadLocal的理解