卢卡斯定理

#include<bits/stdc++.h>
using namespace std;
#define int long long 
const int N=1e5+10;
int f[N],inv[N];

int C(int n,int m,int p)
{
    if(n<m) return 0;
    return f[n]*inv[m]%p*inv[n-m]%p;
}

int Lucas(int n,int m,int p)
{
    if(!m) return 1;
    return Lucas(n/p,m/p,p)*C(n%p,m%p,p)%p;
}

signed main()
{
	ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    int t; cin>>t;
    while(t--)
    {
        int n,m,p; cin>>n>>m>>p; n+=m;
        f[1]=f[0]=inv[1]=inv[0]=1;
        for(int i=2;i<=p;i++)
        f[i]=f[i-1]*i%p,inv[i]=(p-p/i)*inv[p%i]%p;
        for(int i=2;i<=p;i++)
        inv[i]=inv[i]*inv[i-1]%p;
        cout<<(Lucas(n,m,p)+p)%p<<'\n';
    }
}

24/8/27