洛谷 P3807 【模板】卢卡斯定理/Lucas 定理

最新推荐文章于 2024-09-26 22:04:22 发布

qq_38232157

最新推荐文章于 2024-09-26 22:04:22 发布

阅读量234

点赞数

分类专栏：洛谷数论文章标签：卢卡斯定理

本文链接：https://blog.csdn.net/qq_38232157/article/details/127433397

版权

本文介绍了卢卡斯定理在模运算中的应用，特别是在计算组合数c(m, n)模素数p时的角色。首先，回顾了费马小定理和逆元的概念，接着利用逆元计算组合数模p的值。接着，通过数组f[x]存储x!模p的值和g(x)表示x!的模p逆元，公式化c(m, n) mod p = f[n] * g[m] * g[n - m] (mod p)。最后，讨论了如何套用卢卡斯定理来解决更复杂的模运算问题。" 125581726,14574405,SpringBoot的@ConditionalOnBean注解详解,"['spring boot', 'java', '自动配置', 'bean']

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

1、如果 p 是素数， a, p 互素，(a, p) = 1, 则费马小定理 a^(p - 1) = 1 (mod p)
a^(p - 2) 是 a 模 p 的逆元。通过逆元求组合数 c(m, n) 模p 的值
组合公式：
c(m, n) = n! / (m! * (n - m)!), 1 <= n < p , 1 <= m < p
c(m, n) % p = (n! % p) * {(1 / (m! * (n - m)!)) % p}
2、数组 f[x] 表示 x! (mop p) 的值
g(x) 表示 (x!)^-1 (mop p) 的值
那么 c(m, n)(mod p) = f[n] * g[m] * g[n - m] (mod p)
3、套用卢卡斯定理

#include <cstdio>
#include <cstdlib>
#include <cstring>
typedef long long ll;
using namespace std;
const ll MaxN = 100010;
ll f[MaxN];
ll g[MaxN];	// 存逆元
ll t, n, m, p;

ll quick_pow(ll a, ll k)
{
   
	ll res = 1;