【数论】Lucas定理推导过程+模板

Lucas定理主要用于组合数取模。

$$C_n^m \equiv C_{n / p} ^ {m / p} \cdot C_{n \% p}^{m \% p} \pmod p$$

$$p为素数，m和n为非负整数$$

$$\binom{a}{b}\quad 即 \quad C^b_a$$

---

证明过程:

已知p为素数，将非负整数a转化为p进制表示

$$a = a_kp^k + a_{k - 1} p^{k - 1} + \cdots + a_1k + a_0$$

因为p是素数，所以对于

$$1 \le j \le p - 1$$

都有

$$C^j_p = {p\over j} \cdot C_{p-1}^{j-1} \equiv 0 \pmod p$$

于是有

$$(1 + x)^p = 1 + C_p^1 + C_p^2 + \cdots + C_p^{p - 1}x^{p-1} + C_p^px^p = 1 + \sum_{i=1}^{p-1}C_p^j + x^p \equiv 1 + x^p \pmod p\qquad①$$

设

$$n = sp + q$$ $$m = tp + r$$

即

$$s = n \div p$$ $$q = n \% p$$

$$t = m \div p$$ $$r = m \% p$$

则

$$(1 + x)^n = (1 + x)^{sp +q} = (1 + x)^{sp} \cdot (1 + x)^q = \big [(1 + x)^s\big]^p \cdot (1 + x) ^q$$

代入上式得

$$(1 + x) ^ n \equiv (1 + x^p)^s \cdot (1 + x) ^q \pmod p \qquad②$$

对$(1 + x^p) ^ s$ 和 $(1 + x^p) \cdot (1 + x) ^q$ 分别进行二项式展开得

$$(1 + x ^ p) ^ s \cdot (1 + x) ^q \equiv \sum^s_{i = 0} \binom{s}{i} x^{ip} \cdot \sum^q_{j=0}\binom{q}{j}x^j \pmod p$$

即

$$(1 + x) ^ n \equiv \sum_{i = 0}^s \binom {s}{i}x^ip \cdot \sum_{j = 0}^q \binom{q}{j}x^j mod \pmod p\qquad ③$$

对$(1 + x) ^ n$进行普通二项式展开得：

$$(1 + x)^n = \sum _{i = 0}^{sp + q} \binom {sp + q}{k}x^k$$

③式中$x^{tp + r}$的系数为$\binom {sp + q}{tp+ r}$

寻找②式中$x ^{tp + r}$这一项，当且仅当$i = t$，$j = r$时，能得到$x^{tp + r}$的系数，即$\binom {s}{t}$ $\binom {q}{r}$

即

$$\binom {sp + q}{tp + r} \equiv \binom {s}{t}\binom{q}{r} \pmod p$$

即

$$\binom{n}{m} \equiv \binom{n / p}{m / p} \quad \binom{m \% p}{ n \% p} \pmod p$$

证明完毕。

---

模板：

/***********************
author : Meng Xiangcong
***********************/
//#include <bits/stdc++.h>
#include <map>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <string>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <utility>
#include <climits>
#include <ctype.h>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;

#define PI         acos(-1)
#define INF        (1<<30)
#define EPS        1e-6
#define SETDATA    memset(data, 0, sizeof(data))
#define SET(a, b)  memset(a, b, sizeof(a))
typedef long long LL;
//template:
#define MOD 1000000007
LL quickPower(LL a, LL b)
{
    LL ans = 1;
    a %= MOD;
    while (b)
    {
        if (b & 1)
        {
            ans = ans * a %MOD;
        }
        b >>= 1;
        a = a * a % MOD;
    }
    return ans;
}
LL c(LL n, LL m)
{
    if (m > n)
    {
        return 0;
    }
    LL ans = 1;
    for (int i = 1; i <= m; i++)
    {
        LL a = (n + i - m) % MOD;
        LL b = i % MOD;
        ans = ans * (a * quickPower(b, MOD - 2) % MOD) % MOD;
    }
    return ans;
}
LL lucas(LL n, LL m)
{
    if (m == 0)
    {
        return 1;
    }
    return c(n % MOD, m % MOD) * lucas(n / MOD, m / MOD) % MOD;
}
int main()
{
    LL a, b;
    while (~scanf("%lld %lld", &a, &b))
    {
        printf("%lld\n", lucas(a, b)); //输出Cab
    }
    return 0;
}

相关题目

HDU3037解法

数学真是太有趣辣~