组合数

求组合数$C_a^{b}modp$

小数据范围 $1\le a,b\le 2000$

$递推式:C_a^b=C_{a-1}^b+C_{a-1}^{b-1}$

初始化$C[maxn][maxn]$ 数组 $C[i][j]=C_i^j$

const int maxn = 2005; 
const int mod = 1e9+7;
int C[maxn][maxn];
void init()//init函数
{
	for (int i=0;i<maxn;++i)
		for (int j=0;j<=i;++j)
			if (!j) C[i][j] = 1;
			else C[i][j] = (C[i-1][j] + C[i-1][j-1]) % mod;
}

大数据 $1\le a,b\le 10^5$

$公式:C_a^b=\frac{a!}{b!(a-b)!}$

预处理$a!和b{!}^{-1}$

$fact[i]表示i!$

$infact[i]表示i{!}^{-1}(i的阶乘的逆元)$

由费马小定理 $a^{p-1}\equiv 1(modp)(p为质数且p\%a\ne 0)$

$两边同时除以a,a^{p-2}\equiv a^{-1}(modp)$

所以$a^{-1}(modp)是a^{p-2}(modp)$

typedef long long ll;
const int maxn = 1e5+5;
const int mod = 1e9+7;
int fact[maxn]; 
int infact[maxn];


//快速幂求逆元 inverse element
int ie(int a,int b,int p)
{
	int res = 1;
	a %= p;
	while(b)
	{
		if(b&1) res = (ll)res*a%p;
		a = (ll)a*a%mod;
		b>>=1; 
	}
	return res;
}
//初始化函数
void init()
{
    fact[0] = infact[0] = 1;
	for(int i=1;i<maxn;++i)
	{
		fact[i] = (ll)fact[i-1]*i%mod;//求阶乘
		infact[i] = (ll)infact[i-1]*ie(i,mod-2,mod) % mod; //求阶乘的逆元
	}
}
int main()
{
	int a,b;
    cin>>a>>b;
    int res = (ll)fact[a] * infact[b] % mod * infact[a-b] % mod;
    printf("%d\n",res);
	return 0;
}

超大数据 $1\le a,b\le 10^{18}$ ,$1\le p\le 10^5$

$卢卡斯定理:C_a^b=C_{a\%p}^{b\%p}\ast C_{\lfloor a/p\rfloor }^{\lfloor b/p\rfloor }$

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long ll;
int p; //模数 (全局变量)

//快速幂求逆元
int ie(int a,int b)
{
	int res = 1;
	a %= p;
	while(b)
	{
		if(b&1) res = (ll)res * a % p;
		a = (ll)a*a%p;
		b>>=1;
	}
	return res;
}
//求C[a,b]
int C(int a,int b)
{
	if(b>a) return 0;
	int res = 1;
	for(int i=1,j=a;i<=b;++i,--j)
	{
		res = (ll)res * j % p;			//分子 
		res = (ll)res * ie(i,p-2) % p; //分母 
	}
	return res;
}

int lucas(ll a,ll b)
{
	if(a<p && b<p) return C(a,b);
	else return (ll)C(a%p,b%p) * lucas(a/p,b/p) % p;
}

int main()
{
	int n;cin>>n;
	while(n--)
	{
		ll a,b;
		cin>>a>>b>>p;
		cout<<lucas(a,b)<<endl;
	}
	return 0;
}