扩展欧几里得算法

1、

找出一对整数（x，y），使得ax+by=gcd(a,b).

若方程ax+by=gcd(a,b)的一组解为(x0,y0)，则其他解可写为（x0+k*b/gcd(a,b）,y0-k*a/gcd(a,b))

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int exgcd(int a, int b, int &x, int &y)
{
    if (b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int r = exgcd(b, a%b, x, y);
    int t = x;
    x = y;
    y = t - (a/b) * y;
    return r;
}

int main()
{
    int a, b, x, y;
    while (cin >> a>> b)
    {
        int r = exgcd(a, b, x, y);
        cout << "最大公约数为"<< r<< "   "<< "x、y的值分别为" << x << " "<< y<< endl;
        cout << "方程的每一个解都可以由 "<< x <<"+ k*"<< b/r<< " "<< y << "- k*"<< a/r<< " 得到！"<< endl<< endl;
    }
    return 0;
}

注：设a,b,c为任意整数，g=gcd(a,b),方程ax+by=g的一组解是（x0,y0）,则当c是g的倍数时，ax+by=c的一组解是（x0*x/g,y0*c/g）;当c不是g的倍数时无整数解。

2、

求直线ax+by+c=0上有多少个整点（x,y）满足x属于【x1,x2】，y属于【y1,y2】。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int gcd(int x,int y)
{
    return!y?x:gcd(y,x%y);   //欧几里得算法
}
void gcd(int a,int b,int &x,int &y)      //扩展欧几里得算法
{
    if(!b)
    {
        x=1;
        y=0;
    }
    else
    {
        gcd(b,a%b,y,x);
        y-=x*(a/b);
    }
}
int max(int x,int y)
{
    return x>y?x:y;   //最大值算法
}
int min(int x,int y)
{
    return x<y?x:y;   //最小值算法
}
int main()
{
    int T,a,b,c,x1,x2,y1,y2;
    int aa,bb,g,x0,y0;
    int lx,rx,ly,ry;
    while(~scanf("%d",&T))
    {
        while(T--)
        {
            scanf("%d%d%d",&a,&b,&c);
            scanf("%d%d%d%d",&x1,&x2,&y1,&y2);
            g=gcd(a,b);
            aa=a/g;
            bb=b/g;
            gcd(a,b,x0,y0);
            if(c%g==0)
            {
                x0*=c/g;
                y0*=c/g;
                lx=(x1-x0)/bb;
                rx=(x2-x0)*1.0/bb-0.99999999;   //对于下界，向上取整
                ly=(y0-y2)/aa;
                ry=(y0-y1)*1.0/aa-0.99999999;   //对于下界，向上取整
                if(rx<ly||ry<lx) printf("0\n");
                else
                {
                    lx=max(lx,ly);
                    rx=min(rx,ry);
                    printf("%d\n",rx-lx+1);
                }
            }
            else printf("0\n");
        }
   }
return 0;
}