【数论】扩展欧几里得

大家一定都知道欧几里得算法gcd(a,b),它能够在O(lg(ab))的时间内求出整数a和b的最大公约数,是一个非常优秀的算法。这么一个优秀的算法,能不能在不增加时间复杂度的情况下,让它能处理更多的信息,以便解决更多的问题呢?今天,我们来介绍一下 扩展欧几里得算法。

算法解决对象

扩展欧几里得算法可以用来求形如ax+by=GCD(a,b)的方程的 整数解。(当然,a,b是整数)

算法证明

首先,一定存在解符合题意。(我也不知道为什么?)
其次,因为GCD(a,b)=GCD(b,a%b),所以x*a+y*b=GCD(a,b)=GCD(b,a%b)=p*b+q*a%b=p*b+q*(a-a/b*b)=q*a+(p-a/b*q)*b. 说成人话就是: 假设方程 b*x+a%b*y=GCD(b,a%b) 的解为x0,y0;那么原方程的解就是 x=y0,y=x0-(a/b)*y0.
基本情况(Base Case):当b=0时,x=1,y=0.

算法实现

#include<bits/stdc++.h>
using namespace std;
int extended_gcd(int a, int b, int &x, int &y){
 int res,x0,y0;
 if (!b){
  x=1;
  y=0;
  return a;
 }
 res=extended_gcd(b,a%b,x0,y0);
 x=y0; y=x0-a/b*y0;
 return res;
} 
int main(){
 int a,b,x,y,z;
 scanf("%d%d", &a, &b);
 z=extended_gcd(a,b,x,y);
 printf("%d %d %d\n",z,x,y);
 return 0;
} 
/*
模板名称:扩展欧几里得
模板作用:求方程ax+by=GCD(a,b)的一组整数解
*/

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值