【Gym 100633J】Ceizenpok’s formula 扩展Lucas详解

最新推荐文章于 2022-05-05 21:48:00 发布

我是一只计算鸡

最新推荐文章于 2022-05-05 21:48:00 发布

阅读量292

点赞数

分类专栏： Gym 数论与数论函数文章标签：扩展Lucas Gym 100633J

本文链接：https://blog.csdn.net/giftedpanda/article/details/100547866

版权

数论与数论函数同时被 2 个专栏收录

16 篇文章 0 订阅

订阅专栏

Gym

1 篇文章 0 订阅

订阅专栏

$C_n^m \ mod \ p$ , 不为素数

首先对p进行因式分解， $p = p_1^{a_1}p_2^{a_2}p_3^{a_3}...p_n^{a_n}$

然后用中国剩余定理合并

$\left \{ \begin{matrix} C_n^m\ mod \ p_1^{a_1}\\ C_n^m\ mod \ p_2^{a_2}\\ ...\\ C_n^m\ mod \ p_n^{a_n}\\ \end{matrix}$

现在的问题是怎么求出 $C_n^m\ mod\ p^t$

$C_n^m\ mod\ p^t = \frac{n!}{m!(n-m)!}\ mod \ p^t$

因为、 (n - m)! 不一定和 p^t 互质, 所以、的逆元不一定存在

因此我们可以化简一下

$\frac{\frac{n!}{p^x}} {\frac{m!}{p^y} \frac{(n-m)!}{p^z}}\ p^{x-y-z}\ mod \ p^t$ ,这样形如 $\frac{n!}{p^x}$ 与 p^t 一定互质，逆元一定存在

那我们怎么计算 $\frac{n!}{p^x}\ mod \ p^t$ , 我们可以先计算 $n!\ mod\ p^t$

为了方便理解，我们先假设 n = 22, p = 3, t = 2

22!

$=(1\cdot2\cdot3\cdot4\cdot5\cdot6\cdot7\cdot8\cdot9\cdot10\cdot11\cdot12\cdot13\cdot14\cdot15\cdot16\cdot17\cdot18\cdot19\cdot20\cdot21) \ mod \ 3^2$

$=(1\cdot2\cdot4\cdot5\cdot7\cdot8)(10\cdot11\cdot13\cdot14\cdot16\cdot17)(19\cdot20\cdot22)(3\cdot6\cdot9\cdot12\cdot15\cdot18\cdot21)$

$=3^7\cdot7!(1\cdot2\cdot4\cdot5\cdot7\cdot8)^2(19\cdot20\cdot22)\ mod \ 3^2$

$n! \ mod \ p^t = p^{\left \lfloor \frac{n}{p} \right \rfloor}\left \lfloor \frac{n}{p} \right \rfloor!{(\sum_{i=1,i \not\equiv 0\ mod\ p}^{p^t} i)}^{\left \lfloor \frac{n}{p^t} \right \rfloor}({\sum_{i = \left \lfloor \frac{n}{p^t} \right \rfloor p^t , i \not\equiv 0\ mod\ p^t }^n i})$

因为我们要保证互质，逆元才存在，所以 $p^{\left \lfloor \frac{n}{p} \right \rfloor}$ 要被除掉

我们定义 $f(n) = \frac{n!}{p^x}$

$f(n) = f(\left \lfloor \frac{n}{p} \right \rfloor){(\sum_{i=1,i \not\equiv 0\ mod\ p}^{p^t} i)}^{\left \lfloor \frac{n}{p^t} \right \rfloor}({\sum_{i = \left \lfloor \frac{n}{p^t} \right \rfloor p^t , i \not\equiv 0\ mod\ p^t }^n i})$

$\frac{\frac{n!}{p^x}} {\frac{m!}{p^y} \frac{(n-m)!}{p^z}}\ p^{x-y-z}\ mod \ p^t$

$= \frac{f(n)}{f(m)f(n-m)}p^{x-y-z} \ mod\ p^t$ ，那么现在就可以求逆元了

接下来讲解如何求解指数

我们令 g(n) = x ，显然 $p^{\left \lfloor \frac{n}{p} \right \rfloor}$ 是我们需要的指数，但是 $\left \lfloor \frac{n}{p} \right \rfloor!$ 可能还有的倍数

所以, $g(n) = \left \lfloor \frac{n}{p} \right \rfloor + g(\left \lfloor \frac{n}{p} \right \rfloor )$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll power(ll a, ll b, ll n)  // 快速幂
{
	a %= n;
	ll ans = 1;
	while(b) {
		if(b & 1) ans = (ans * a) % n;
		a = (a * a) % n;
		b >>= 1;
	}
	return ans;
}
ll exgcd(ll a, ll b, ll &x, ll &y) // 扩展欧几里得
{
	if(b == 0) {
		x = 1;
		y = 0;
		return a;
	}
	ll d = exgcd(b, a % b, x, y);
	ll t = x;
	x = y;
	y = t - a / b * y;
	return d;
}
ll inverse(ll a, ll p)  // 逆元
{
	ll x, y;
	exgcd(a, p, x, y);
	x = (x % p + p) % p;
	return x;
}
ll CRT(ll *a, ll *m, ll n) // 中国剩余定理
{
	ll M = 1, ans = 0, x, y;
	for(ll i = 1; i <= n; i++) M *= m[i];
	for(ll i = 1; i <= n; i++) {
		ll w = M / m[i];
		exgcd(w, m[i], x, y);
		x = (x % m[i] + m[i]) % m[i];
		ans = (ans + (a[i] % M + M) % M * w % M * x % M) % M;
	}
	return (ans % M + M) % M;
}
ll calc(ll n, ll x, ll p)  // n! mod p
{
	if(n == 0) return 1;
	ll s = 1;
	for(ll i = 2; i <= p; i++) if(i % x) s = s * i % p;
	s = power(s, n / p, p);
	for(ll i = 2; i <= n % p; i++) if(i % x) s = s * i % p;
	return s * calc(n / x, x , p) % p;
}
ll multilucas(ll n, ll m, ll x, ll p)  //
{
	ll cnt = 0;
	for(ll i = n; i; i /= x) cnt += i / x;
	for(ll i = m; i; i /= x) cnt -= i / x;
	for(ll i = n-m; i; i /= x) cnt -= i / x;
	return power(x, cnt, p) % p * calc(n, x, p) % p
		   * inverse(calc(m, x, p), p) % p * inverse(calc(n-m, x, p), p) % p;
}
ll exlucas(ll n, ll m, ll p)  //扩展Lucas
{
	ll cnt = 0;
	ll prime[100], a[100];
	for(ll i = 2; i * i <= p; i++) {
		if(p % i == 0) {
			prime[++cnt] = 1;
			while(p % i == 0) prime[cnt] *= i, p /= i;
			a[cnt] = multilucas(n, m, i, prime[cnt]);
		}
	}
	if(p > 1) prime[++cnt] = p, a[cnt] = multilucas(n, m, p, p);
	return CRT(a, prime, cnt);
}
int main()
{
	ll n, m, p;
	while(scanf("%lld %lld %lld", &n, &m, &p) == 3) {
		printf("%lld\n", exlucas(n, m, p));
	}
	return 0;
}

我是一只计算鸡

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
【Gym 100633J】Ceizenpok’s formula 扩展Lucas详解

,不为素数首先对p进行因式分解，然后用中国剩余定理合并现在的问题是怎么求出因为、不一定和互质, 所以、的逆元不一定存在因此我们可以化简一下,这样形如与一定互质，逆元一定存在那我们怎么计算, 我们可以先计算为了方便理解，我们先假设因为我们要保证互质，逆元才存在，所以要被除掉我们定义，那么现在就可以求逆...
复制链接

扫一扫

专栏目录