Lucas定理

#include<bits/stdc++.h>
#define N 100005
#define ll long long
using namespace std;
ll a[N],b[N];int p;
ll lucas(int x,int y){//求CXY
	if(x<y)return 0;//X比Y小返0，即N是负值不用计算（觉得没什么意义）
	else
		if(x<p)return b[x]*a[y]*a[x-y]%p;
		else return lucas(x/p,y/p)*lucas(x%p,y%p)%p;
//b[x]*a[y]*a[x-y]就是X的阶乘除以Y的阶乘再除以(x-y)的阶乘，X是下标
}
int main(){
	int T;scanf("%d",&T);
	while(T--){//T组数据
		int m,n;//m是下标，N是上标
		scanf("%d%d%d",&m,&n,&p);//求CMN模P
		a[0]=b[0]=a[1]=b[1]=1;
		for(int i=2;i<p;i++)b[i]=b[i-1]*i%p;//这是直接阶乘模P
		for(int i=2;i<p;i++)a[i]=(p-p/i)*a[p%i]%p;//此时A[I]指任意数X求(X/i)mod P就等于乘上A[I]mod p
		for(int i=2;i<p;i++)a[i]=a[i-1]*a[i]%p;//此时A[i]指任意数X求(X/i!)mod P就等于乘上A[I]mod p
		printf("%lld\n",lucas(m,n));//求CMN
	}
}

定义：