扩展欧几里得算法(求解二元一次方程)

假设方程为:aX+bY=c;
首先判断是否有整数解:c%gcd(a,b)==0则有解
代码:

#include "iostream"
#include
using namespace std;
void exgcd(int a,int b,int &d,int &x,int &y)
{
  if(b==0)
  {
    x=1;
    y=0;
    d=a;//d==gcd(a,b);
    return ;
  }
  exgcd(b,a%b,d,x,y);
  int tmp=x;
  x=y;
  y=tmp-(a/b)*y;
}
void fun(int a,int b,int c,int &x,int &y)
{
  int d;
  exgcd(a,b,d,x,y);
  if(c%d) return ;
  while(x<0)//求最小的非负数
  {
    x+=b;
    y-=a;
  }
  cout<<x%b<<endl;
  // return true;
}
int gcd(int a,int b)
{
  return b==0?a:gcd(b,a%b);
}
int main(int argc, char const *argv[]) {
  int a,b,x,y,t;
  scanf("%d",&t);
  while(t--)
  {


  scanf("%d %d",&a, &b);

    // cout<<gcd(a,b)<<endl;
    if(gcd(a,b)==1)
  {
    fun(a,b,1,x,y);
  }
  else
  {
    printf("-1\n" );
  }
  // exgcd(a,b,x,y);
}
return 0;
}
//x,y则为解

你可能感兴趣的:(知识整理)