扩展欧几里得模板

扩展欧几里得：对于不完全为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;
}