浅谈乘法逆元

浅谈乘法逆元

前言

乘法逆元……难以言表，一直觉得没有什么用，但是面对现实又不得……

正题

定义

若 $a*x\equiv1 (\bmod {b})$，且a与b互质，那么我们就能定义: x为a的逆元，记为 $a^{-1}$
，所以我们也可以称x为a的倒数(我的理解是不在模P意义下)
在模P意义下
所以对于 $\frac{a}{b} (\bmod {p})$
我们就可以求出b在 $\bmod {p}$下的逆元，然后乘上 $a$，再 $\bmod {p}$，就是这个乘法逆元的值了。

求法

扩欧的乘法逆元求法我已经在我的上篇博客介绍过了，
乘法逆元的扩展欧几里得求法

我现在来介绍一点玄妙的东西

快速幂

费马小定理大家应该会的
$a^{p−1}≡1(modp)$
不会自行百度

那么
$a*x\equiv1 (\bmod {b})$
可得
$a*x\equiv a^{p−1} (\bmod {b})$
可得
$x\equiv a^{p-2} (\bmod {b})$

ll fpm(ll x, ll power, ll mod) 
{
    x%=mod;
    ll ans=1;
    while (power) {
        if (power & 1) ans = (ans * x) % mod;
        x = (x * x) % mod;
        power >>= 1;
    }
    return ans;
}
printf ("%lld\n",fpm(i,p-2,p));

线性求一连串数字模P的乘法逆元

这个真的好玄啊
首先
很容易知道
$1^{-1}\equiv 1(mod p)$
设
$p=k*i+r$
$k*i+r\equiv 0(mod p)$
乘上 $i^{-1}r^{-1}$
得
$k*r^{-1}+i^{-1}\equiv 0(mod p)$
得
$i^{-1}\equiv -k*r^{-1}(mod p)$
得
$i^{−1}≡−⌊\frac{p}{i}⌋∗(pmodi)^{−1}(modp)$

好神奇啊（雾）

#include<bits/stdc++.h>
using namespace std;

long long a[3000010],n,m,p;

int main()
{
	scanf("%lld%lld",&n,&p);
	a[1]=1;
	printf("%lld\n",1);
	for (int i=2; i<=n; i++)
	  {
	  	a[i]=((-(p/i)*a[p%i])%p+p)%p;
	  	printf("%lld\n",a[i]);
	  }
}