组合数的三种模板

常规组合数

#include<bits/stdc++.h>
using namespace std; 
typedef long long ll;
ll C[2333][2333];
const int MOD=1e9+7;
int main() {
    for(int i=0; i<=2000; ++i) {
        C[i][0]=1;
        for(int j=1; j<=i; ++j)
            C[i][j]=(C[i-1][j-1]+C[i-1][j])%MOD;
    }
    int a,b,c,d;
    cin>>a>>b>>c>>d;
    ll t=C[a+b+c+d][a+b+c]*C[a+b+c][a+b]%MOD*C[a+b][a]%MOD;
    cout<<t<<"\n";
}

卢卡斯组合数

对于模数不是很大情况

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+5;
ll n,m,p;
ll fac[maxn];
void getfac(ll p) { //预处理阶层
	fac[0]=1;
	for(int i=1; i<=p; i++)	fac[i]=fac[i-1]*i%p;
}
ll power(ll a,ll n,ll p) { //快速幂运算
	ll ans=1;
	while(n) {
		if(n&1)	ans=ans*a%p;
		a=a*a%p;
		n/=2;
	}
	return ans;
}
ll lucas(ll n,ll m,ll p) {
	ll ans=1;
	while(n&&m) {
		ll a=n%p,b=m%p;
		if(a<b) return 0;
		ans=(ans*fac[a]*power(fac[b]*fac[a-b]%p,p-2,p))%p;
		n/=p,m/=p;
	}
	return ans;
}
int main() {
	int t;
	cin>>t;
	while(t--) {
		cin>>n>>m>>p;
		getfac(p);
		cout<<lucas(n,m,p)<<endl;
	}
	return 0;
}

组合数（忘记从哪里瞟来的了）

对于n，m不是很大但是模数很大情况

#include<bits/stdc++.h> 
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int maxn=2555;
ll fac[maxn],inv[maxn];
ll fpow(ll x,ll y) {
    ll ans=1;
    while(y) {
        if(y&1) ans=(ans*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return ans;
}
ll C(int n,int m) {
    if(m>n || n<0 || m<0) return 0;
    ll ans=1;
    ans=(ans*fac[n])%mod;
    ans=(ans*inv[m])%mod;
    ans=(ans*inv[n-m])%mod;
    return ans;
}
void init(int N) {
	fac[0]=inv[0]=1;
    for(int i=1;i<=N;i++) fac[i]=(1LL*fac[i-1]*i)%mod;
    inv[N]=fpow(fac[N],mod-2);
    for(int i=N-1;i>=1;i--) inv[i]=(1LL*inv[i+1]*(i+1))%mod;
}
int main() {
	int a,b,c,d;
	init(maxn-1);
	scanf("%d%d%d%d",&a,&b,&c,&d);
	printf("%lld\n",1LL*C(a+b+c+d,a)%mod*C(b+c+d,b)%mod*C(c+d,c)%mod*C(d,d)%mod);
	return 0;
}