扩展欧几里得:对于不完全为0的非负整数a和b,gcd(a,b)表示a,b的最大公约数,必然存在整数x,y,使得ax+by=gcd(a,b)。
因为:gcd(a,b) = gcd(b,a%b)
所以:bx1+(a%b)y1 = ax+by
bx1+(a-a/b*b)y1 = ax+by
ay1+b(x1-a/b*y1) = ax+by
令:x=y1,则:y=x1-a/b*y1
那么:ax+by=c又怎么求呢?
分两种情况:
假设d=gcd(a,b)
1、如果c%d != 0,即c不是gcd(a,b)的倍数,那么无解。
2、如果c%d == 0,且t = c/d。
ax0+by0 = d
ax0+ab/d + by0-ab/d = d
a(x0+b/d) + b(y0-a/d) = d
所以:x=x0+b/d,y=y0-a/d
ax'+by' = c = t*d
所以求出x,y以后,分别乘 t 就求出了x',y'
//扩展欧几里得模板
#include<iostream>
#include<cstdio>
#include<string.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y)
{
if(a==0&&b==0) return -1;
if(b==0)
{
x=1;
y=0;
return a;
}
int d=exgcd(b,a%b,x,y);
int x1=x,y1=y;
x=y1;
y=x1-(a/b)*y1;
return d;
}
int main()
{
int a,b,c,x0,y0,x,y;
scanf("%d%d%d",&a,&b,&c);
int d=exgcd(a,b,x0,y0);
if(c%d!=0)
printf("no solution\n");
else
{
int t=c/d;
x=(x0+b/d)*t;
y=(y0-a/d)*t;
printf("x0=%d y0=%d\n",x0,y0);//x0,y0是ax+by=gcd(a,b)的解
printf("x=%d y=%d\n",x,y);//x,y是ax+by=c的解
}
return 0;
}