卢卡斯定理的详解和线性求解逆元（含题目）

【卢卡斯定理】

卢卡斯定理是用来求解：$C_{n+m}^m$%mod的值，而mod是一个质数。

我也是看了很长时间的讲解才看懂的（感觉好弱啊~~~）

这张图是我认为讲解最全的，（可能一时看不懂，我再补充一下）。

首先，你要看懂

$$(1+x)^p≡1+x^p\%p$$
并推出 $$(1+x)^a≡(1+x)^{a_0}*(1+x^p)^{a_1}*•••*((1+x)^{p^k})^{a_k}\%p$$
以上都不难，最难的是最后一步，我一直没有想明白，想了很久才明白的（太菜了）。
就是左右两边的 $x^b$的项系数是相同的，
左边是 $$C_a^b*x^b$$
右边是,每一项提取一个式子 $$C_{a_0}^{b_0}x^{b_0}*C_{a_1}^{b_1}x^{pb_1}*C_{a_2}^{b_2}x^{p^2b_2}*•••*C_{a_k}^{b_k}x^{p^kb_k}$$ $$=C_{a_0}^{b_0}*C_{a_1}^{b_1}*C_{a_2}^{b_2}*•••*C_{a_k}^{b_k}*x^{b_0+pb_1+p^2b_2+•••+p^kb_k}$$
而 $b_0+pb_1+p^2b_2+•••+p^kb_k$是p进制的b，所以是 $x^b$所以
$$C_a^b=C_{a_0}^{b_0}*C_{a_1}^{b_1}*C_{a_2}^{b_2}*•••*C_{a_k}^{b_k}$$
证明完毕了！

也可以用递归的方法来写：Lucas(a,b)=Lucas(n%p,m%p)*Lucas(n/p,m/p)%p;
继续对Lucas(n/p,m/p)递归即可。

【线性求解逆元】

这个博客讲解的很清楚了：线性求解逆元
求解n以内的逆元，时间复杂度O(n)
设 t=M/i , k=M%i;

inv数组记录的就是i的逆元。

下面是一道例题：
卢卡斯定理

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;

typedef long long ll;
const int maxn=1e5+9;

ll factorial[maxn];//记录阶乘
ll inelement[maxn];//记录，先记录每个数的逆元，再记录阶乘的逆元

ll lucas(int n,int m,int mod)
{
    if(n<m)
        return 0;
    else if(n<mod)
        return factorial[n]*inelement[m]*inelement[n-m]%mod;
    else return lucas(n%mod,m%mod,mod)*lucas(n/mod,m/mod,mod)%mod;
}
int main()
{
    int t,n,m,mod;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d",&n,&m,&mod);
        factorial[0]=factorial[1]=inelement[0]=inelement[1]=1;
        //阶乘
        for(int i=2;i<=m+n;i++)
        {
            factorial[i]=factorial[i-1]*i%mod;
        }

        //逆元，先先求每一个数的逆元，线性求解

        for(int i=2;i<=m+n;i++)
        {
            inelement[i]=(mod-mod/i)*inelement[mod%i]%mod;
        }

        //求阶乘的逆元
        for(int i=2;i<=m+n;i++)
        {
            inelement[i]=inelement[i-1]*inelement[i]%mod;
        }
        printf("%lld\n",lucas(n+m,m,mod));

    }
    return 0;
}