lucas定理求组合数

#include<bits/stdc++.h> 
using namespace std;
const int N=100002;
long long a,b;
int p;
long long qmi(long long a,long long b,int p)
{
    long long ans=1;
    while(b)
    {
        if(b&1)
            ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}
long long C(long long a,long long b,int p)
{
    if(b>a) return 0;
    long long ans=1;
    for(int i=1,j=a;i<=b;i++,j--)
    {
        ans=ans*j%p;
        ans=ans*qmi(i,p-2,p)%p;
    }
    return ans;
}
long long lucas(long long a,long long b,int p)
{
    if(a<p&&b<p) return C(a,b,p);
    return  C(a%p,b%p,p)*lucas(a/p,b/p,p)%p;
}
int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        cin>>a>>b>>p;
        printf("%lld\n",lucas(a,b,p));
    }
    return 0;
}