数论初步———扩展欧几里德算法

问题引入：直线上的点
求直线$ax$ + $by$ + c = 0上有多少个整数点(x，y)满足$x$ $\epsilon$ [$x_{1}$，$x_{2}$]，$y$ $\epsilon$ [$y_{1}$，$y_{2}$]。
分析：
在解决这个问题之前，首先学习扩展欧几里德算法——找一对整数$（x，y）$，使得$ax + by$ = $gcd（a，b）$。
注意这里的 $x 和 y$ 不一定是正数，也可能是负数或者是0,。例如，$gcd（6，15）$ = 3，6*3 - 15*1 = 3，其中$x$ = 3，$y$ = 1。这个方程还有其他解，如$x$ = - 2，$y$ = 1。
下面是扩展欧几里德算法的程序：

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

证明

(证明过程参考自百度百科)

原式: $ax+by=gcd(a,b)$(假设a≥b)

当$b=0$时有$gcd(a,b)=a$，此时$x=1，y=0$；
当$b\neq0$时,根据欧几里得定理$gcd(a,b)=gcd(b,a$%$b)$可得$ax+by=gcd(a,b)=gcd(b,a$%$b)=bx′+(a$%$b)y′$,即$ax+by=bx′+(a$%$b)y′=bx′+(a−b∗⌊a/b⌋)y′$

移项得
$ax+by=bx′+(a$%$b)y′=ay′+b(x′−⌊a/b⌋y′)$

根据恒等定理,有

$$\begin{cases} x=y′ \\ y=x′−⌊a/b⌋y′ \end{cases}$$
证毕
上面求出 $ax+by=gcd(a,b)$的一组解 $（x_{0}，y_{0}）$，那么其他解呢？
任取另外一组解 $（x，y）$，则 $ax_{0}+by_{0}=ax+by$(他们都等于 $gcd（a，b）$)；
变形得： $a（x_{0}-x）= b（y_{0}-y）$
假设 $gcd（a，b）= g$，
方程左右同时除以g，得： $a^{'}（x_{0}-x）= b^{'}（y_{0}-y）$
其中 $a^{'} = \frac{a}{g}$， $b^{'} = \frac{b}{g}$
这里 $a^{'}$， $b^{'}$互素，因此 $x_{0}-x$一定是 $b^{'}$的倍数， $y_{0}-y$一定是 $a^{'}$的倍数；
$x_{0}-x=kb^{'}$；
$y - y_{0}=ka^{'}$；
所以得出结论：

假设a，b，c为任意整数。若方程$ax+by=c$的一组解为$（x_{0}，y_{0}）$，则它的任意常数解都可以写成$（x_{0} + kb^{'}，y_{0} - ka^{'}）$，其中$a^{'} = \frac{a}{gcd（a，b）}$，$b^{'} = \frac{b}{gcd（a，b）}$，$k$取任意整数。
当$c$是$g$的倍数时$ax+by=c$的一组解是$（\frac{x_{0}c}{g}，\frac{y_{0}c}{g}）$;
当$c$不是是$g$的倍数时无整数解。