关于卢卡斯定理的粗略理解

今天突然想起了卢卡斯定理还不会就简单学习了一下，在这里进行分享。

在学习卢卡斯定理之前必须要掌握逆元的概念，逆元详解
在掌握了逆元之后，开始介绍卢卡斯定理。
卢卡斯定理应用于大组合数取模，通常情况下，组合数可以用下面方法进行求取：
（1）.C(n, m) = C(n - 1，m) + C(n - 1, m - 1)
（2）.C(n, m) mod p = n!/(m!(n - m)!)；
（3）.C(n, m) =(n-i+1) / i; (i从1开始取值)

但是当n，m很大时运算会超时，所以卢卡斯定理，应运而生，但是为什么和逆元有关系呢，因为这时候组合数的数值是非常大的，所以需要对这个数进行取模操作，所以就出现了除法取余，这就需要用到费马小定理求取逆元啦。
————_
注意：p最大为1e5且为素数，再多的话可能卢卡斯定理也无法解决。
————-_

Lucas 定理：A、B是非负整数，p是质数。A、B写成p进制：A=a[n]a[n-1]…a[0]，B=b[n]b[n-1]…b[0] 这里的每一个数组元素表示其p进制的每一位。

表达式C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p ；

卢卡斯定理的证明一定要好好看这篇博客中的最后那张图，大神写的真心流批。

卢卡斯定理模板：

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

const int maxn=1e5+10;
const int p=10007;
ll f[maxn];

void  init()
{
    f[0]=f[1]=1;
    for(int i=2;i<=p;i++) f[i]=f[i-1]*i%p;
}

ll  quick_mod(ll x,ll n)
{
    ll res=1;
    while(n)
    {
        if(n&1) res=res*x&p;
        x=x*x%p;
        n>>=1;
    }
    return res;
}

ll lucas(ll n,ll m)
{
    ll ret=1;
    while(n&&m)//求取p进制
    {
        ll nn=n%p,mm=m%p;
        if(nn<mm) return 0;
        ret=ret*f[nn]*quick_mod(f[nn-mm]*f[mm]%p,p-2)%p;//求取出除法逆元
        n/=p,m/=p;//求取p进制的各个位数
    }

    return ret;
}

int main()
{
    init();
    printf("%I64d\n",lucas(n,m));
    return 0;
}