1. 程式人生 > >hdu2669 擴充套件歐幾里德 二元一次解不定方程

hdu2669 擴充套件歐幾里德 二元一次解不定方程


題意:
給出兩個非負整數a,b 求x,y 使ax+by=1,而且x非負並最小的答案


題解:

樸素的歐幾里德原理: gcd(a,b)=gcd(b,a mod b)
擴充套件歐幾里德定理:
 對於不完全為 0 的非負整數 a,b,gcd(a,b)表示 a,b的最大公約數,
必然存在整數對 x,y ,使得 gcd(a,b)=ax+by。

令a>b
當b==0 ,x=1,y=0;
否則,根據樸素的歐幾里德原理,有
a*x1+b*y1=b*x2+(a%b)*y2
=b*x2+(a-a/b*b)*y2=b*x2+a*y2- a/b*b*y2
所以,
x1=y2,y1=x2-(a/b)*y2;
x1,y1 就是所求值之一。

若x不是非負數,可讓x+=b,y-=a( 左邊加上a*b再減a*b,等式依然成立 ),直到符合;

對題目中的gcd(a,b)==1才有解

程式碼如下:

#include<iostream>
#include<cstdio>
#include<cmath>

using namespace std;
__int64 a,b,c,d,x,y;

__int64 exgcd(__int64 a,__int64 b)
{
	if(b==0)
	{
		x=1;y=0;
		return a;
	}
	d=exgcd(b,a%b);
	__int64 xx=y,yy=x-(a/b)*y;
	x=xx;y=yy;
	return d;
}

int main()
{
	c=1;
	while(scanf("%I64d%I64d",&a,&b)!=EOF)
	{
		if(a%2==0&&b%2==0){printf("sorry\n");continue;}
		d=exgcd(a,b);
		if(d !=1){printf("sorry\n");continue;}
		while(x<0)
		{
			x+=b;
			y-=a;
		}
		printf("%I64d %I64d\n",x,y);
	}
	return 0;
}