卢卡斯定理/Lucas 定理【模板】

最新推荐文章于 2024-08-01 16:31:59 发布

ssl_fuyang

最新推荐文章于 2024-08-01 16:31:59 发布

阅读量175

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/qq_43010386/article/details/119676900

版权

数论专栏收录该内容

62 篇文章 0 订阅

订阅专栏

>Link

luogu P3807

>解题思路

（早上打比赛写不出来的时候去洛谷看了下lucas

求 $C^n_mmodp$ （ $1≤n,m,p≤10^5$ ）

因为 $p$ 是质数（太爱质数啦），所以 $\binom{p}{j}\equiv 0(modp)$
通过化简可以得到 $(1+j)^p=\sum_{i=0}^{p}\binom{p}{i}*j^i$
$(1+j)^p=1+\sum_{i=1}^{p-1}\binom{p}{i}*j^i+j^p\equiv 1+j^p(modp)$

设 $n = s p + q, m = l p + r$ $(q < r)$
所以 $(1+x)^{sp+q}=(1+x)^{sp}*(1+x)^q\equiv (1+x^p)^s*(1+x)^q$
$(1+x)^{sp+q}\equiv (\sum _{i=0}^{s}\binom{s}{i}x^{pi})*(\sum _{i=0}^{q}\binom{q}{i}x^{i})$

因为 $p$ 是质数（太爱质数啦），所以两边可以同时约去 $x^m$
$\binom{n}{m}*x^m\equiv \binom{s}{l}*x^{lp}*\binom{q}{r}*x^r\equiv \binom{s}{l}*\binom{q}{r}$

所以
${\color{Red} C\binom{n}{m}\equiv C\binom{\left \lfloor n/p \right \rfloor}{\left \lfloor m/p \right \rfloor}*C\binom{nmodp}{m modp}(modp)}$

PS：组合数里面有除法，所以要用逆元

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;

int T;
LL p;

LL power (LL a, LL b)
{
	a %= p;
	LL ret = 1;
	for (; b; b >>= 1, a = a * a % p)
	  if (b & 1) ret = ret * a % p;
	return ret;
}
LL C (LL a, LL b)
{
	if (a < b) return 0; //一定要加!!!
	LL ret = 1;
	for (LL i = b + 1; i <= a; i++) ret = ret * (i % p) % p;
	for (LL i = 1; i <= a - b; i++) ret = ret * power (i, p - 2) % p;
	return ret;
}
LL lucas (LL n, LL m)
{
	return C (n % p, m % p) * C (n / p, m / p) % p;
}

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