hdu 2669 C - Romantic

题意： a*x+ by= 1; 给出a,b;求x,y;  （0<a,b<2^31）

  

这道题摆明的扩展GCD，但是我写了很多很多次都没写对，最后发现是求出X',Y'（ax'+by'= gcd(a,b)）;  然后求最小非负整数解的时候出现了问题

设 at+bp= gcd(a,b);

通过这题我证明了

ans1= t;
  ans1= (ans1%b + b)%b;// b为方程的循环节

这种方法求最小非负整数是正确的（不用枚举，尽管这题枚举也能过）

注意定义的时候用__int64， 不然会wa的

代码：

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
#include<limits.h>
using namespace std;
__int64 t,p;
__int64 gcd(__int64 a,__int64 b)
{
	return b==0? a: gcd(b,a%b);
}
void exgcd(__int64 a,__int64 b)
{
	if(b==0)
	{
		t= 1;
		p= 0;
	}
	else
	{
		exgcd(b,a%b);
		int temp= t;
		t= p;
		p= temp- a/b*p;
	}
}
int main()
{
	__int64 a,b,c;
	while(scanf("%I64d%I64d",&a,&b)!=EOF)
	{
		__int64 d= gcd(a,b);
		c= 1;
		if(d!=1)
		{
			printf("sorry\n");
		}
		else
		{
			exgcd(a,b);
			__int64 ans1,ans2;
			ans1= t;
			ans1= (ans1%b + b)%b;
			ans2= (1- ans1*a)/ b;
			printf("%I64d %I64d\n",ans1,ans2);
		}
	}
return 0;
}