UVA 11768 - Lattice Point or Not

首先本题需要用到扩展欧几里得算法……

关于exgcd算法的一点简略证明：

那么，对于函数exgcd(a,b)=(d,x,y)，其中d满足d=gcd(a,b); (x,y)满足ax+by=d;

则exgcd(b,a mod b)=(d,x',y')，而此时，使用 x = y' ;  y = x' - floor(a/b) * y' = x' - floor(a/b) * x 就能得到exgcd(a,b)的值。

故我们可以有扩展欧几里得算法如下：

void exgcd(int a,int b,int &d,int &x,int &y){
    if(!b){d=a;x=;y=;}
    else {exgcd(b,a%b,d,y,x);y-=(a/b)*x;}
}

此处的边界条件是exgcd(a,0)=(a,1,0)

那么此时，若c mod gcd(a,b) = 0，就容易求得一个整点(x,y)，那么如何得到所有的点？

然后是题目的代码：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 typedef long long ll;
 using namespace std;
 void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
     if(!b){d=a;x=;y=;}
     else {exgcd(b,a%b,d,y,x);y-=(a/b)*x;}
 }
 int main() {
     double X1,Y1,X2,Y2;
     int t;
     scanf("%d",&t);
     while(t--){
         scanf("%lf%lf%lf%lf",&X1,&Y1,&X2,&Y2);//获得初始坐标(x1,y1),(x2,y2)
         ll x1=X1*,x2=X2*,y1=Y1*,y2=Y2*;//将坐标都换成整型
         if(x1==x2)//如果这条线是一条竖直的线
         {
             if(x1%!=){//如果此时的直线l:x=c，c不是整数，那么这条线上显然没有整点
                 printf("0\n");
                 continue;
             }
             y2=floor(max(Y1,Y2));
             y1=ceil(min(Y1,Y2));
             printf("%lld\n",y2-y1+);
             continue;
         }
         if(y1==y2)//如果这条线是一条水平的线
         {
             if(y1%!=){//如果此时的直线l:y=c，c不是整数，那么这条线上显然也没有整点
                 printf("0\n");
                 continue;
             }
             x2=floor(max(X1,X2));
             x1=ceil(min(X1,X2));
             printf("%lld\n",x2-x1+);
             continue;
         }
         ll a=(y2-y1)*,b=(x1-x2)*,c=x1*y2-x2*y1;//计算出l:ax+by=c
         ll x,y,d;
         if(X2<X1) swap(X2,X1);//确保x2>x1
         x1=ceil(X1),x2=floor(X2);
         if(x1>x2){//如果一条线的斜率过大，使得x1,x2有m<x1<x2<m+1 (m为整数)，那么显然这个范围内没有整点
             printf("0\n");
             continue;
         }
         ll k=;
         exgcd(a,b,d,x,y);//得到满足ax+by=gcd(a,b)的(x0,y0)
         if (c%d==)//如果c是gcd(a,b)的整数倍，即c = k * gcd(a,b)，那么显然点( k * x0 , k * y0 )是ax+by=c上的一个整点
         {
             a/=d,b/=d;
             x*=c/d,y*=c/d;
             if(b<) b=-b;
             x=x-(x-x1)/b*b;
             while(x<x1) x+=b;
             while(x+k*b<=x2) k++;
             printf("%lld\n",k);
         }
         else printf("0\n");
     }
     return ;
 }