扩展欧几里得算法求解不定方程||中国剩余定理

对于
a x + b y = g c d ( a , b ) ax+by=gcd(a,b) ax+by=gcd(a,b)
这样的方程,可以用扩展欧几里得算法exgcd求出一组通解。

根据欧几里得求gcd:
g c d ( a , b ) = g c d ( b , a % b ) gcd(a,b)=gcd(b,a\%b) gcd(a,b)=gcd(b,a%b)
可得
b x + ( a % b ) y = g c d ( b , a % b ) bx+(a\%b)y=gcd(b,a\%b) bx+(a%b)y=gcd(b,a%b)
根据
a % b = a − ( a / b ) ∗ b a\%b=a-(a/b)*b a%b=a(a/b)b
可得
b x + a y − ( a / b ) b ∗ y = g c d ( b , a % b ) bx+ay-(a/b)b*y=gcd(b,a\%b) bx+ay(a/b)by=gcd(b,a%b)
化简得
a y + b ( x − ( a / b ) y ) = g c d ( b , a % b ) ay+b(x-(a/b)y)=gcd(b,a\%b) ay+b(x(a/b)y)=gcd(b,a%b)

x ′ = y , y ′ = ( x − ( a / b ) y ) x' = y , y' = (x-(a/b)y) x=y,y=(x(a/b)y)
可得
a x ′ + b y ′ = g c d ( b , a % b ) &lt; = &gt; a x + b y = g c d ( a , b ) ax&#x27;+by&#x27;=gcd(b,a\%b) &lt;=&gt; ax+by=gcd(a,b) ax+by=gcd(b,a%b)<=>ax+by=gcd(a,b)
根据
g c d ( a , 0 ) = a gcd(a,0)=a gcd(a,0)=a a x + b y = a ax+by=a ax+by=a
可以得出一组平凡解
x = 1 , y = 0 x=1,y = 0 x=1,y=0

所以一直递归下去可以得出一组平凡解,然后再往回带得出 a x + b y = g c d ( a , b ) ax+by=gcd(a,b) ax+by=gcd(a,b)的一组解。

泛化来看不定方程
a x + b y = c ax+by=c ax+by=c
只有满足 c % g c d ( a , b ) = = 0 c\%gcd(a,b) == 0 c%gcd(a,b)==0才有解。

求解同余方程
a x ≡ b   m o d   n ax \equiv b\ mod\ n axb mod n
a x = b + n y ax = b+ny ax=b+ny
令 y = -y
a x + n y = b ax+ny = b ax+ny=b
就转变成了上述形式

中国剩余定理
https://blog.csdn.net/niiick/article/details/80229217
https://www.cnblogs.com/freinds/p/6388992.html

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值