求直线上的整数点算法

问题描述：

求直线上有多少个整数点(x,y)满足。

问题分析：

把直线方程进行移项，得到。

假设已经得到直线上一个点的坐标，假设直线存在另外一个点坐标为，那么把两点代入直线得到方程，整理得到 。

假设,上上式两边同除g，可以得到，其中。把等式再变形得到，此时的和互素，必定是的整数倍，即k可以为任意整数。整理得到，。

由上面的分析知道，只需要就出方程的任意一组解，就可以求出其他解。

扩展欧几里德算法：

找出一对整数,使得。

对应的算法程序如下：

void gcd(int a,int b, int& d, int& x, int& y){

      if (!b){ d = a; x = 1; y = 0; }

      else{ gcd(b, a%b, d, y, x); y -= x*(a /b); }

}//d为gcd(a,b),x和y为解

对于a和b最大公约数的求解在欧几里德算法中就介绍过，唯一需要记忆的就是x和y的初始条件x=1，y=0，x和y的递归关系y -= x*(a /b);。

假设,的一组解是。如果-c是g的整数倍，则的一组解是；当c不是g的整数倍是无整数解。

综上，先求的解，再求的解，最后求的其他解，分三步走。程序实现如下：

#include<iostream>
#include<cstdio>
using namespace std;

void gcd(int a, int b, int& g, int& x, int& y){
	if (!b){ g = a; x = 1; y = 0; }
	else{ gcd(b, a%b, g, y, x); y -= x*(a / b); }
}

int main()
{
	int a, b, c, x1, x2, y1, y2, g, x, y;
	cin >> a >> b >> c >> x1 >> x2 >> y1 >> y2;
	gcd(a, b, g, x, y);//ax+by=gcd(x,y)
	if (c%g == 0){
		int x0 = -x*(c / g), y0 = -y*(c / g);//ax+by=-c
		printf("(%d,%d)\n",x0,y0);
		bool kup = true, kdown = true;
		for (int k = 1;; k++){
			int xn = x0 - k*(b / g), yn = y0 + k*(a / g);
			if (kup&&xn >= x1&&xn <= x2&&yn >= y1&&yn <= y2){
				printf("(%d,%d)\n",xn,yn);
			}
			else kup = false;
			xn = x0 + k*(b / g);
			yn = y0 - k*(a / g);
			if (kdown&&xn >= x1&&xn <= x2&&yn >= y1&&yn <= y2){
				printf("(%d,%d)\n", xn, yn);
			}
			else kdown = false;
			if (!(kup || kdown))break;
		}
	}
	else printf("No point.\n"); 
	system("pause");
	return 0;
}