扩展欧几里得

View Code

#include<stdio.h>
#include<stdlib.h>
__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    if(b==0) {x=1;y=0;return a;}
    else 
    {
        __int64 d=exgcd(b,a%b,x,y);
        __int64 t=x;
        x=y;
        y=t-a/b*y;
        return d;
    }
}
int main()
{
    __int64 a,b,x,y;
    while(scanf("%I64d%I64d",&a,&b)!=EOF)
    {
        __int64 d=exgcd(a,b,x,y);
        if(d!=1) {printf("sorry\n");continue;}
        while(x<0) x+=b,y-=a;
        printf("%I64d %I64d\n",x,y);
    }
    return 0;
}