关于数论的一些知识

1.组合数取膜（Lucas定理）

求 C n m % p ， 其 中 1 ≤ p ≤ 1 0 9 ， 并 且 p 是 素 数 求 C_{n}^{m} \% p，其中1\leq p \leq 10^{9}， 并且p是素数 求Cnm​%p，其中1≤p≤109，并且p是素数

#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long
ll n,m,p;
ll pow_m(ll a,ll k,ll p)
{
    ll ans=1;
    ll tmp=a%p;
    while(k)
    {
        if(k&1)ans=ans*tmp%p;
        tmp=tmp*tmp%p;
        k>>=1;
    }
    return ans;
}
ll C(ll n,ll m,ll p)
{
    if(m>n)return 0;
    ll a=1,b=1;
    for(int i=1;i<=m;i++)
    {
        a=a*(n+i-m)%p;
        b=b*i%p;
    }
    return a*pow_m(b,p-2,p)%p;
}
ll Lucas(ll n,ll m,ll p)
{
    ll ans=1;
    while(n&&m)
    {
        ans=ans*C(n%p,m%p,p)%p;
        n/=p;
        m/=p;
    }
    return ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&m,&p);
        printf("%lld\n",Lucas(n,m,p));
    }
    return 0;
}

标签：tmp,return,Lucas,数论,ll,知识,int,关于,ans
来源： https://blog.csdn.net/qq_42752207/article/details/120614891