拓展欧几里得解方程ax+by=gcd(a,b)

好吧，过年在家还是看了点点书，还自己给自己出了两道坑爹题。

这道题目一开始想了半天敲出来了，但是只对一组数据得出的答案是对的。。。。后来又写了一遍就对了。。。。

描述：

  求出方程：ax+by=gcd(a,b)的初始解,输入a,b（int），输出x,y（空格隔开）;

输入：

60 22

22 60

12453 2467

输出：

-4  11

11  -4

304  -1613

代码：

#include<stdio.h>
void swap(__int64 &a,__int64 &b)
{
     a^=b;
     b^=a;
     a^=b;
}
__int64 gcd(__int64 a,__int64 b)
{
     return b?gcd(b,a%b):a;
}
void func(__int64 a,__int64 b)
{
     int f;
     f=0;
     if(a<b)
    {
        swap(a,b);
        f=1;
     }
     __int64 x1=1,x2,x,y1=1,y2,y,p,q,r1,r2,r;
     r1=a%b;r2=b%r1;
    q=a/b;p=b/r1;r=r1%r2;
     x1=1;y1=-q;
     x2=-p;y2=1+p*q;
 //    printf("1\n x1:%d y1:%d\n x2:%d y2:%d\nq:%d\n",x1,