LightOj-1067-Lucas定理

现在才看Lucas，真是菜QAQ。

Lucas首先是解决大组合数取模问题的，C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p，循环求解这个公式就是Lucas。

#include <iostream>
#include <cmath>
#include <vector>
#include <cstring>
#include <cstdio>
#include <algorithm>

using namespace std;

typedef long long int ll;
const int maxn=1e6+5;
const int mod=1000003;
int t;
ll n,m;
ll f[maxn];

void fac(ll x)//处理出0-mod内的阶乘
{
    f[0]=1;//一定要处理0，0的阶乘是1，容易忘
    for(ll i=1;i<=x;i++)
        f[i]=(f[i-1]*i)%x;
}

ll quick_pow(ll a,ll b)//快速幂
{
    ll res=1;
    while(b>0)
    {
        if(b%2) res=(res*a)%mod;
        a=(a*a)%mod;
        b/=2;
    }
    return res;
}

ll Lucas(ll n,ll k)//Lucas
{
    ll res=1;
    while(n && k)
    {
        ll nn=n%mod,kk=k%mod;
        if(nn<kk) return 0;
        res=(res*f[nn]*quick_pow(f[kk]*f[nn-kk]%mod,mod-2))%mod;//注意最后取模
        n/=mod;
        k/=mod;
    }
    return res;
}

int main()
{
    fac(mod);
    scanf("%d",&t);
    for(int i=1;i<=t;i++)
    {
        scanf("%lld%lld",&n,&m);
        ll ans=Lucas(n,m,mod);
        printf("Case %d: %lld\n",i,ans);
    }
    return 0;
}

 

C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p