分数取模

最新推荐文章于 2024-10-14 16:36:25 发布

Y390d

最新推荐文章于 2024-10-14 16:36:25 发布

阅读量2w

点赞数 12

分类专栏：数学扩展欧几里得

本文链接：https://blog.csdn.net/godleaf/article/details/79844074

版权

数学同时被 2 个专栏收录

26 篇文章 0 订阅

订阅专栏

扩展欧几里得

2 篇文章 0 订阅

订阅专栏

这里介绍一种分数取模的代码：

假设要求 $\frac{1}{a}\quad mod \quad p$

这里要引用小费马定理 $a^{p-1}\quad mod \quad p = 1 \quad mod \quad p$ (有兴趣的可以百度查下证明过程) ，这里对这个公式做点改动

$\because \quad a^{p-1} \quad = a^{p-2}\times a$

$\therefore \quad (a^{p-2}\times a)\quad mod \quad p = 1\quad mod\quad p$

$\therefore \quad a^{p-2}\quad mod \quad p = a^{-1}\quad mod\quad p$

在代码里面，我们可以通过扩展欧几里得定理来求一个数的模，但是扩展欧几里得暂时还不能求分数的模，那么这里可以将这个分数的取模换成对整数求模来得出相同的结果，也就是说把 $\quad \frac{1}{a}\quad mod \quad p \quad$ 这个问题转化求 $\quad a^{p-2}\quad mod \quad p$

现在把代码放上来：

long long fast_mod(long long a,long long b) {
	long long r = 1;
	a %= mod;
	while (b) {
		if (b & 1) r = (r*a) % mod;
		a = (a*a) % mod;
		b >>= 1;
	}
	//cout << r << endl;
	return r;
}

同样的，如果分子不为1，那么可以用 $(m\times a^{-1}) \quad mod \quad p = ((m\quad mod\quad p)\times (a^{-1}\quad mod\quad p))\quad mod \quad p$ (这是取模的运算规则，可以在百度查找一下，对取模十分有用)