题意: 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;
}