卢卡斯Lucas定理

最新推荐文章于 2021-08-17 00:05:16 发布

爱寂寞的时光

最新推荐文章于 2021-08-17 00:05:16 发布

阅读量271

点赞数

分类专栏： # 组合数学、概率论、博弈论与数理统计等基础数学

本文链接：https://blog.csdn.net/jiahonghao2002/article/details/117231285

版权

组合数学、概率论、博弈论与数理统计等基础数学专栏收录该内容

40 篇文章 2 订阅

订阅专栏

博客介绍了卢卡斯定理及其在计算组合数模p时的应用，通过生成函数和多项式展开证明了定理。同时，展示了如何利用费马小定理预处理阶乘和求逆元。在给定的C++代码示例中，演示了如何使用卢卡斯定理解决实际问题，例如在P3807题目中的计算。

摘要由CSDN通过智能技术生成

卢卡斯Lucas定理

$\mathscr{LUCAS}$

内容

$\binom{n}{m} \equiv \binom{\lfloor \frac{n}{p} \rfloor}{\lfloor \frac{m}{p} \rfloor} \cdot \binom{n \mod p}{m \mod p} \mod p$

其中 $p$ 为质数。

上式称为Lucas定理。

证明

生成函数 $1+x)^n$ 项 $x^{m}$ 前的系数为 $\binom{n}{m}$ ，因此

$(1+x)^n \equiv (1 + x)^{p \lfloor \frac{n}{p} \rfloor} (1 + x)^{n \mod p} \mod p$

多项式 $1 + x)^p \mod p$ ：

$x)^p \equiv \sum_{i = 0}^{p}\binom{p}{i}x^i = 1 + x^p \mod p$

此处考虑 $p$ 是质数，式子 $\binom{p}{i}$ ，当且仅当 $i = 0$ 或 $i = p$ 的时候模为 $1$ 。

故：

$x^p)^{\lfloor \frac{n}{p} \rfloor} (1 + x)^{n \mod p} \mod p$

那么卷积对 $x^m$ 项的贡献只有当前项取 $p$ 的倍数，后项取余数的时候才能取到，
因此：

$\binom{n}{m} \equiv \binom{\lfloor \frac{n}{p} \rfloor}{\lfloor \frac{m}{p} \rfloor} \cdot \binom{n \mod p}{m \mod p} \mod p$

证毕。

例题

P3807

预处理阶乘+费马小定理求逆元，Lucas定理的运用。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

#define FR freopen("in.txt", "r", stdin)
#define FW freopen("out1.txt", "w", stdout)

#define MAXT 100010
#define INV(x) fpow(x, p - 2, p)

ll fac[MAXT];

ll fpow(ll a, ll e, ll p)
{
    ll ans = 1;
    a %= p;
    for (; e; e >>= 1, a = (a * a) % p)
    {
        if (e & 1)
            ans = (ans * a) % p;
    }
    return ans;
}

void prework(ll p)
{
    fac[0] = fac[1] = 1;
    for (ll i = 2; i <= p; i++)
    {
        fac[i] = (fac[i - 1] * i) % p;
    }
}

ll C(ll n, ll m, ll p)
{
    if (m > n)
        return 0;
    if (m == 0)
        return 1;
    return (((fac[n] * INV(fac[m])) % p) * INV(fac[n - m])) % p;
}

ll lucas(ll n, ll m, ll p)
{
    if (m == 0)
        return 1;
    return (lucas(n / p, m / p, p) * C(n % p, m % p, p)) % p;
}

int main()
{
    int t;
    scanf("%d", &t);
    while (t--)
    {
        ll n, m, p;
        scanf("%lld %lld %lld", &n, &m, &p);
        prework(p);
        printf("%lld\n", lucas(n + m, n, p));
    }
    return 0;
}

爱寂寞的时光

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
卢卡斯Lucas定理

卢卡斯Lucas定理LUCAS\mathscr{LUCAS}LUCAS内容(nm)≡(⌊np⌋⌊mp⌋)⋅(nmod pmmod p)mod p\binom{n}{m} \equiv \binom{\lfloor \frac{n}{p} \rfloor}{\lfloor \frac{m}{p} \rfloor} \cdot \binom{n \mod p}{m \mod p} \mod p(mn)≡(⌊pm⌋⌊pn⌋)⋅(mmodpnmodp)modp其中ppp为质数。上式
复制链接

扫一扫

专栏目录