Acwing-基础算法课笔记之数学知识(扩展欧几里得算法)

一、扩展欧几里得算法

1、裴蜀定理

对于任意正整数 a a a b b b,那么一定存在非零整数 x x x y y y,使得 a x + b y = e x g c d ( a , b ) ax+by=exgcd(a,b) ax+by=exgcd(a,b)

2、过程模拟

例如求 g c d ( a , b ) gcd(a,b) gcd(a,b)

∙ \bullet b = 0 b=0 b=0 时,则可以直接返回 a a a 的值,即最大公约数,推理如下:

根据公式 a x + b y = e x g c d ( a , b ) ax+by=exgcd(a,b) ax+by=exgcd(a,b),得:

a × x + 0 × y = a a\times x+0\times y=a a×x+0×y=a

⇔ a × x = a \Lrarr a\times x=a a×x=a

⇔ x = 1 \Lrarr x=1 x=1

⇒ y = 0 \Rarr y=0 y=0

∙ \bullet b ≠ 0 b\not =0 b=0 时,求得扩展欧几里得算法 e x g c d ( b , a % b , y , x ) exgcd(b,a\%b,y,x) exgcd(b,a%b,y,x),推理如下:

b y + ( a by+(a by+(a m o d mod mod b ) x = e x g c d ( a , b ) b)x=exgcd(a,b) b)x=exgcd(a,b)

⇒ a \rArr a a m o d mod mod b = a − ⌊ a b ⌋ b b=a-⌊\frac{a}{b}⌋b b=abab

⇒ b y + ( a − ⌊ a b ⌋ b ) x = e x g c d ( a , b ) \rArr by+(a-⌊\frac{a}{b}⌋b)x=exgcd(a,b) by+(abab)x=exgcd(a,b)

⇒ a x + b ( y − ⌊ a b ⌋ x ) = e x g c d ( a , b ) \rArr ax+b(y-⌊\frac{a}{b}⌋x)=exgcd(a,b) ax+b(ybax)=exgcd(a,b)

3、代码模板

// 求x, y,使得ax + by = gcd(a, b)
int exgcd(int a, int b, int &x, int &y)
{
    if (!b)
    {
        x = 1; y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= (a/b) * x;
    return d;
}

Acwing-扩展欧几里得算法

二、线性同余方程

1、定义

给定 n n n 组数据 a i a_i ai b i b_i bi m i m_i mi,对于每组数求出一个 x i x_i xi,使其满足 a i × x i ≡ b i ( m o d a_i\times x_i\equiv b_i(mod ai×xibi(mod m i ) m_i) mi)

2、模拟过程

a = 2 a=2 a=2 b = 3 b=3 b=3 m = 6 m=6 m=6,此时以上并不满足条件。

a = 4 a=4 a=4 b = 3 b=3 b=3 m = 5 m=5 m=5,要使 4 x % 5 = 3 4x\%5=3 4x%5=3,所以 x = − 3 x=-3 x=3 x = 2 x=2 x=2

3、结论证明

a × x ≡ b ( m o d a\times x\equiv b(mod a×xb(mod m ) m) m)

⇔ \Lrarr ∃ y ∈ Z \exist y\in Z yZ,使得 a x = m y + b ax=my+b ax=my+b

⇔ a x − m y = b \Lrarr ax-my=b axmy=b

⇔ y ′ = − y \Lrarr y'=-y y=y

⇔ a x + m y ′ = g c d ( a , m ) \Lrarr ax+my'=gcd(a,m) ax+my=gcd(a,m)(条件: g c d ( a , m ) ∣ b gcd(a,m)|b gcd(a,m)b

  • 49
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

不会敲代码的狗

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

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

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

打赏作者

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

抵扣说明:

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

余额充值