组合(Lucas)
2014-07-14 10:04
357 查看
Description
给出组合数C(n,m), 表示从n个元素中选出m个元素的方案数。例如C(5,2) = 10, C(4,2) = 6.可是当n,m比较大的时候,C(n,m)很大!于是xiaobo希望你输出 C(n,m) mod p的值!
Input
输入数据第一行是一个正整数T,表示数据组数 (T <= 100) 接下来是T组数据,每组数据有3个正整数 n, m, p (1 <= m <= n <= 10^9, m <= 10^4, m < p < 10^9, p是素数)
Output
对于每组数据,输出一个正整数,表示C(n,m) mod p的结果。
Sample Input
Sample Output
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
using namespace std;
__int64 power(__int64 a,__int64 b,__int64 n)
{
__int64 res=1;
while(b)
{
if(b&1)
res=res*a%n;
b>>=1;
a=a*a%n;
}
return res;
}
__int64 cal(__int64 n,__int64 r,__int64 p)
{
__int64 i,res=1,re;
for(i=1;i<=r;i++)
{
res=res*(n-i+1)%p;
re=power(i,p-2,p);
res=res*re%p;
}
return res;
}
__int64 lucas(__int64 n,__int64 m,__int64 p)
{
if(n<m)
return 0;
else
return cal(n,m,p);
}
int main()
{
__int64 n,m,p,res;
int t;
scanf("%d",&t);
while(t--)
{
scanf("%I64d%I64d%I64d",&n,&m,&p);
res=1;
while(n&&m)
{
res=res*lucas(n%p,m%p,p);
if(res==0)
break;
n=n/p;
m=m/p;
}
printf("%I64d\n",res);
}
return 0;
}
给出组合数C(n,m), 表示从n个元素中选出m个元素的方案数。例如C(5,2) = 10, C(4,2) = 6.可是当n,m比较大的时候,C(n,m)很大!于是xiaobo希望你输出 C(n,m) mod p的值!
Input
输入数据第一行是一个正整数T,表示数据组数 (T <= 100) 接下来是T组数据,每组数据有3个正整数 n, m, p (1 <= m <= n <= 10^9, m <= 10^4, m < p < 10^9, p是素数)
Output
对于每组数据,输出一个正整数,表示C(n,m) mod p的结果。
Sample Input
2 5 2 3 5 2 61
Sample Output
110
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
using namespace std;
__int64 power(__int64 a,__int64 b,__int64 n)
{
__int64 res=1;
while(b)
{
if(b&1)
res=res*a%n;
b>>=1;
a=a*a%n;
}
return res;
}
__int64 cal(__int64 n,__int64 r,__int64 p)
{
__int64 i,res=1,re;
for(i=1;i<=r;i++)
{
res=res*(n-i+1)%p;
re=power(i,p-2,p);
res=res*re%p;
}
return res;
}
__int64 lucas(__int64 n,__int64 m,__int64 p)
{
if(n<m)
return 0;
else
return cal(n,m,p);
}
int main()
{
__int64 n,m,p,res;
int t;
scanf("%d",&t);
while(t--)
{
scanf("%I64d%I64d%I64d",&n,&m,&p);
res=1;
while(n&&m)
{
res=res*lucas(n%p,m%p,p);
if(res==0)
break;
n=n/p;
m=m/p;
}
printf("%I64d\n",res);
}
return 0;
}
相关文章推荐
- HDU5894分位置(组合数,lucas,乘法逆元)
- Codeforces Round #313 (Div. 2) E. Gerald and Giant Chess 定理Lucas求大组合数
- fzu 2020 Lucas 定理,组合数求模模板
- 【组合数】【Lucas】Number Theory Special Training T3 comb 题解
- 4403: 序列统计 组合数学+Lucas定理
- hdu3037 隔板法+Lucas定理求大组合取模
- BZOJ 4403: 序列统计 (组合数 Lucas 数论推导)
- lucas定理,组合数学问题
- 洛谷 1869 愚蠢的组合数 Lucas定理 解题报告
- HDU 3037 Saving Beans(组合数学+Lucas定理)
- FZU 2020 组合(Lucas定理)
- FZU 2020-组合(Lucas定理+逆元解决大组合数求模)
- FOJ 2020 组合(组合数取素数摸模板:Lucas)
- hdu5698 Lucas 大组合数
- 【hdu 5894】【组合数学 lucas 费马小求逆元】 【m个人坐n个围成一周的位置每个人距离至少为k,求方法数】
- 51nod 1120 机器人走方格V3(组合数学+卡特兰数+Lucas)
- Lucas定理 大组合数取模
- poj 1430 Binary Stirling Number 求斯特林数奇偶性 数形结合| 斯特林数奇偶性与组合数的关系+lucas定理 好题
- UOJ#86:mx的组合数 (Lucas定理+原根+NTT+高精度)
- 51Nod1778 小Q的集合 【组合数】【Lucas定理】