hdu--2669Romantic

对于欧几里得算法，仔细研究几天，反复看了几次百度百科和别人的博客详解，总算是明白了不少，始终相信，坚持就一定会有收获，有关于欧几里得算法的题目，推导出不定方程很重要，一旦要求的不定方程推导出来，对硬搬扩展欧几里得算法模板应该也能解决问题，关于深刻理解扩展gcd算法，由于之前写了青蛙的约会的经典扩展Gcd的题目，收藏下下面这个博客链接以后回头再理解理解：点击打开链接每一步都很详细

接下来就是关于这道牵手虐单身狗的2669，题目意思是：给a，b两个整数，求使得x*a+b*y=1；的x,和y，并要求x要是正数，自然一步了然是扩展gcd，所以死套模板，然而在杭电上交代码用long long int过不了，只好用--I64d,，这里不得不吐槽下杭电oj的一些bug，做KMP算法的时候，next数组用next不行，非要改成其他名称才能过，坑的一批

用完gcd算法，x和y值就会有了，然而这时候的x并不一定是正数，于是根据通解公式x=c*x+b*t;y=c*y-a*t;这里t是任意正数，题目中c又是1，所以只要x<0,就可以用循环x+b使其大于0，自然同时要y-a；这样求得的就是正确答案啦！

题目链接：点击打开链接

代码如下：

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
//typedef long long __int64;
__int64 kgcd(__int64 a,__int64 b,__int64 &x,__int64 &y){
	if(!b){
		x=1;
		y=0;
		return a;
	}
	__int64 ans=kgcd(b,a%b,y,x);
	y-=a/b*x;
	return ans;
}
int main(){
	__int64 a,b;
	while(cin>>a>>b){
		__int64 x,y;
		__int64 ans=kgcd(a,b,x,y);
		if(1%ans!=0){
			cout<<"sorry"<<endl;
			//printf("sorry\n");
			continue;
		}
		while(x<=0){
			x+=b;
			y-=a;
		}
		cout<<x<<" "<<y<<endl;
		//printf("%I64d %I64d\n",k1,k2);
	}
	return 0;
}