【算法基础 & 数学】快速幂求逆元(逆元、扩展欧几里得定理、小费马定理)


原文链接

为什么需要逆元

首先,在算法竞赛中,很多情况下会遇到数值很大的数据,这个时候,题目往往会让我们对某个数去摸,来控制数据范围。

在±*运算中,我们可以对每个数单独取模,然后再对运算之后的数取模。

但是除法比较特殊,例如: ( 40 ÷ 5 ) m o d 10 ≠ ( ( 40 m o d 10 ) ÷ ( 5 m o d 10 ) ) ) m o d 10 (40\div5)mod 10\neq((40mod 10)\div(5mod 10)))mod10 (40÷5)mod10=((40mod10)÷(5mod10)))mod10

那我们可以将其转化为乘法,就相当于:除以一个数,就是乘以一个数的倒数

这个时候就要用到逆元

逆元的概念

1.单位元

在一个集合中,对于某种运算,如果对于任何的集合元素a和元素e运算,得到还是集合元素a本身,就称e为这个运算下的单位元

  • 例如

在加法中:对任意实数a: a + e = e + a = a a + e = e + a = a a+e=e+a=a, 则e为单位元, e = 0 e = 0 e=0
在乘法中:对任意实数a: a × e = e × a = a a \times e = e \times a = a a×e=e×a=a, 则e为单位元, e = 1 e = 1 e=1
(每种运算都有自己的单位元)

2.逆元

在一个集合中,对于某种运算,如果任意两个元素的运算结果等于单位元,则称这两个元素互为逆元

  • 例如

在加法中:对任意实数a: a + − a = 0 a + -a = 0 a+a=0, 而0是加法的单位元,则-a为加法的逆元
在乘法中:对任意实数a: a × a − 1 = 1 a \times a^{-1} = 1 a×a1=1,则 a − 1 a^{-1} a1为乘法的逆元

3.模乘的单位元

对于模 n n n乘法,所有模 n n n a a a同余的数都可以表示成:
a ( m o d   n ) = k n + a    ( k ∈ Z ) a(mod~n)=kn+a~~(k\in Z) a(mod n)=kn+a  (kZ)
令单位元为 e ( m o d   n ) e(mod~n) e(mod n),将 a ( m o d   n ) a(mod~n) a(mod n) e ( m o d   n ) e(mod~n) e(mod n)进行模乘运算,得到:
a ( m o d   n ) × e ( m o d   n ) = ( k 1 n + a ) ( k 2 n + e ) = ( k 1 k 2 n 2 + k 1 e n + k 2 a n + a e ) = ( k 1 k 2 n + k 1 e + k 2 a ) n + a e \begin{equation*} %加*表示不对公式编号 \begin{split} & a(mod~n)\times e(mod~n)\\ & = (k_1n+a)(k_2n+e)\\ & =(k_1k_2n^2+k_1en+k_2an+ae)\\ & =(k_1k_2n+k_1e+k_2a)n+ae\\ \end{split} \end{equation*} a(mod n)×e(mod n)=(k1n+a)(k2n+e)=(k1k2n2+k1en+k2an+ae)=(k1k2n+k1e+k2a)n+ae

根据单位元定义,
a ( m o d   n ) × e ( m o d   n ) = a ( m o d   n ) a(mod~n)\times e(mod~n) = a(mod~n) a(mod n)×e(mod n)=a(mod n)

带入定义式:
( k 1 k 2 n + k 1 e + k 2 a ) n + a e = k n + a (k_1k_2n+k_1e+k_2a)n+ae = kn+a (k1k2n+k1e+k2a)n+ae=kn+a


{ k = k 1 k 2 n + k 1 e + k 2 a e = 1 \begin{cases} k = k_1k_2n+k_1e+k_2a\\ e = 1\\ \end{cases} {k=k1k2n+k1e+k2ae=1

也就是说,模乘的单位元是1

4.模乘的逆元

模乘运算中,任意整数 a ( m o d   n ) a(mod~n) a(mod n)的逆元表示为:
a − 1 ( m o d   n ) a^{-1}(mod~n) a1(mod n)
根据定义,满足:
a a − 1 ≡ 1 ( m o d   n ) aa^{-1}\equiv 1(mod~n) aa11(mod n)

可以理解为 a a a a − 1 a^{-1} a1 n n n的作用下发生了反应,变成了1

但是,不像加减法和乘除法每个数(乘除法除0以外)都有逆元

对于每个数 a a a均有唯一的与之对应的乘法逆元 x x x,使得 a x ≡ 1 ( m o d   n ) ax\equiv 1(mod~n) ax1(mod n)
逆元存在的充要条件: g c d ( a , n ) = 1 gcd(a,n)=1 gcd(a,n)=1(这个数和模数互质)

此时,在模 n n n的意义下,一个数 a a a如果有逆元 x x x,那么除以 a a a就相当于乘以 x x x

开始求逆元

如何在给定 a a a n n n的前提下给出逆元?

先在这里附上欧几里得定理

1.扩展欧几里得定理

已知整数a、b,扩展欧几里得算法可以在求得a、b的最大公约数的同时,能找到整数x、y(其中一个很可能是负数),使它们满足贝祖等式 a x + b y = g c d ( a , b ) ax + by = gcd(a,b) ax+by=gcd(a,b)

扩展欧几里得常用来求形如方程:ax + by = c 的整数通解或者特解。

等式方程 ax+by = c 是不一定有整数解x、y的
但是已知:若 c c%gcd(a,b)=0 c,则方程 a x + b y = c ax + by = c ax+by=c必定存在整数解,否则必定无解(其推导如下)。
对于最简单的情况: 对于不完全为0的非负整数 a , b , g c d ( a , b ) a,b,gcd(a, b) a,b,gcd(a,b) 表示 a , b a,b a,b的最大公约数,必定存在整数对 x , y x, y x,y满足 a × x + b × y = = g c d ( a , b ) a\times x+b\times y==gcd(a, b) a×x+b×y==gcd(a,b)。 我们一般根据欧几里算法与最大公约数的关系由最简单的情况来拓展推导方程的通解。

已知:a % gcd(a,b) == 0,b % gcd(a,b) == 0;若整数x、y为方程 ax + by = c 的一组解
则 ax % gcd(a,b) == 0,by % gcd(a,b) == 0
则 ax + by % gcd(a,b) == 0
即 c % gcd(a,b) == 0
因此 方程 ax + by = c 有解的充要条件是 c % gcd(a,b) == 0
————————————————
版权声明:本文为CSDN博主「阿阿阿安」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/qq_40772692/article/details/81183174

  • 下面举一个例子

【例题】给定正整数 a , b a,b a,b,求满足等式 a x + b y = 1 ax+by=1 ax+by=1 x x x的最小正整数解。如果不存在,返回-1。

首先找出 a , b a,b a,b的最大公约数,令 g = g c d ( a , b ) g = gcd(a,b) g=gcd(a,b)
则,原式可以转化为:
g ( a g x + b g y ) = 1 g(\frac{a}{g}x+\frac{b}{g}y) = 1 g(gax+gby)=1
如果 a b = 1    ( a , b ∈ Z + ) ab = 1~~(a,b\in Z^+) ab=1  (a,bZ+),则 a = 1 , b = 1 a=1, b=1 a=1,b=1

所以 g = 1 g = 1 g=1

所以,只需要考虑 a , b a,b a,b互素的情况

下面给出扩展欧几里得定理的代码:

int extend_gcd(int a,int b,int &x,int & y)
{
	if(b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	// x,y调换传给下一次递归等价于x1 = y2
	int t = extend_gcd(b,a%b,y,x);
	//等价y1 = x2 -(a/b) * y2
	y -= a / b * x;
	return t;
}

那上面这个跟逆元有什么联系呢?
来看这道题

【例题】给定正整数 p p p a a a,求满足 a x ≡ 1 ( m o d   p ) ax\equiv 1(mod~p) ax1(mod p)的最小正整数 x x x,如果不存在,返回-1

乍看上去跟上面一点关系都没有
将原式换成:
a x = k n + 1    ( k ∈ Z ) ax = kn+1~~(k\in Z) ax=kn+1  (kZ)
移项后得到
a x − k n = 1 ax-kn = 1 axkn=1
由于k是个整数,不分正负,所以可以表示为
a x + k n = 1 ax+kn = 1 ax+kn=1
于是,可表示为
a x + b y = 1 ax+by = 1 ax+by=1

ll inv(ll a, ll n){
    ll x, y;
    ExpGcd(a, n, x, y);
    x = (x % n + n) % n;
    return x;
}

2.费马小定理

b b b存在乘法逆元的充要条件是 b b b与模数 m m m互质。当模数 m m m为质数时, b m − 1 ≡ 1 ( m o d   m ) b^{m-1}\equiv 1(mod~m) bm11(mod m), b m − 2 b^{m-2} bm2为b的逆元

【例题】给定素数 p p p和正整数 a a a,求满足 a x ≡ 1 ( m o d   p ) ax\equiv 1(mod~p) ax1(mod p)的最小正整数 x x x,如果不存在,返回-1

a a a p p p的倍数时, a x ≡ 0 ( m o d   p ) ax\equiv 0(mod~p) ax0(mod p),所以一定不存在
当不是倍数(a和p互质), a p − 1 ≡ 1 ( m o d   p ) ⇒ a × a p − 2 ≡ 1 ( m o d   p ) a^{p-1}\equiv 1(mod~p) \Rightarrow a\times a^{p-2}\equiv 1(mod~p) ap11(mod p)a×ap21(mod p)

此时 a p − 2 a^{p-2} ap2就是a的逆元

  • 这道题就是求的 a p − 2 a^{p-2} ap2,考的费马定理和快速幂

题目描述:给定 n n n a i , p i a_i,p_i ai,pi,其中 p i p_i pi 是质数,求 a i a_i ai p i p_i pi的乘法逆元,若逆元不存在则输出 impossible

注意:请返回在 0 ∼ p − 1 0∼p−1 0p1 之间的逆元。

#include<iostream>
#include<algorithm>
using namespace std;

typedef long long LL;

LL qmi(int a, int b, int p){
    LL res = 1 % p;
    while(b){
        if(b & 1) res = res * a % p;
        a = a * (LL)a % p;
        b >>= 1;
    }
    return res;
}

int main(){
    int n, a, p;
    cin >> n;
    
    while(n--){
        cin >> a >> p;
        int res = qmi(a, p - 2, p);
        if(a % p) cout << res << endl;      //保证a不是p的倍数
        else cout << "impossible" << endl;
    }
    
    return 0;
}
  • 26
    点赞
  • 29
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
扩展欧几里得算法解一元一次不定方程 ax + by = gcd(a,b) 的一种方法,其中 a 和 b 是整数,gcd(a,b) 是它们的最大公约数,x 和 y 是整数解。逆元是指在模运算下,一个数的乘法逆元是指与它相乘后模运算得到 1 的数。在数论中,常常需要一个数在模意义下的逆元,即一个数 k 满足 (k * x) % m = 1,其中 m 是模数。 下面是扩展欧几里得算法逆元的 C 语言实现: ```c #include <stdio.h> // 扩展欧几里得算法 int exgcd(int a, int b, int *x, int *y) { if (b == 0) { *x = 1; *y = 0; return a; } int gcd = exgcd(b, a % b, y, x); *y -= a / b * (*x); return gcd; } // 逆元 int modinv(int a, int m) { int x, y; int gcd = exgcd(a, m, &x, &y); if (gcd != 1) { return -1; // a 和 m 不互质,不存在逆元 } else { return (x % m + m) % m; // 转化为正整数 } } int main() { int a = 3, m = 11; int inv = modinv(a, m); if (inv == -1) { printf("%d 在模 %d 意义下不存在逆元\n", a, m); } else { printf("%d 在模 %d 意义下的逆元是 %d\n", a, m, inv); } return 0; } ``` 这个程序中,exgcd 函数通过递归实现扩展欧几里得算法,返回 a 和 b 的最大公约数,并且出 x 和 y 的值。在 modinv 函数中,我们调用 exgcd 函数出 a 和 m 的最大公约数,并且判断 a 和 m 是否互质,如果不互质则不存在逆元。否则,根据扩展欧几里得算法的结果,出 x 的值作为 a 在模 m 意义下的逆元。注意,由于 x 可能是负数,所以要将其转化为正整数。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

为梦而生~

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值