扩展欧几里得算法与唯一分解定理

欧几里得算法

  欧几里得算法又称辗转相除法。

功能

  首先说一下这个算法是用来干啥的。
  欧几里得算法的作用是求两个数$a$和$b$的最大公约数，即$gcd(a,b)$.

算法过程

  这里再说一下怎么用这个算法求两个数$a$和$b$的最大公约数。
  我们不妨设$a>b$。
  我们发现如果$b$是$a$的约数，那么$b$就是二者的最大公约数。我们讨论不能整除的情况，即$a=b\ast q_1+r_1$,其中$a,b,q_1,r_1$都是整数、$r_1=a\%b$，我们通过证明可以得到$gcd(a,b)=gcd(b,a\%b)$.

证明

  我们只需要证明 $a$和$b$ 与 $b$和$a\%b$ 有相同的公因子。
  我们设$d$是 $a$和$b$ 的最大公因子，即$d|a$(表示d是a的因子)且$d|b$.注意到$r=a-b\ast q$，所以$d|r$。又因为$d|b$，所以$d$是$b$和$r$的公因子。
  反之，如果设$d|b$、$d|(a\%b)$，我们知道$a=b\ast q+r$，所以$d|a$，又因为$d|b$，所以$d$是$a$和$b$的公因子。
  既然两式的公因子相同，那么两式的最大公因子也相同。即$gcd(a,b)=gcd(b,a\%b)$.
  由此，我们可以持续递归操作，等到$b==0$，即上一步$a\%b==0$的情况下，求得$a$和$b$的最大公约数。
代码：

int gcd(int a, int b)
{
	if(b==0)
		return a;
	else
		return gcd(b,a%b);
}

扩展欧几里得算法

  扩展欧几里得算法又称裴蜀定理或贝祖定理，可以用来求一元线性同余方程的解，也即是求一元二次方程的解。

内容

  定理：设$a$,$b$是不全为零的整数，则存在整数$x$,$y$，使得$a\ast x+b\ast y=gcd(a,b)$.
  这个算法是用来判断一元二次方程$ax+by=c$是否存在整数解，如果存在求出其的通解。

求解过程

  这里给出结论，不做证明，当$gcd(a,b)|c$时，一元二次方程有整数解，否则，无整数解。
  求解过程：我们假设$a\ge b$可知$a=b\ast q+r$，$r=a-b\ast q$,由$gcd(a,b)=gcd(b,a\%b)$知，$a{x}_1+b{y}_1=c$与$b{x}_2+r{y}_2=c$有相同的解。即$a{x}_1+b{y}_1=b{x}_2+r{y}_2$，如果我们知道$x_2$与$y_2$，我们就可以求$x_1$与$y_1$。
$$a{x}_1+b{y}_1=b{x}_2+r{y}_2$$$$\because r=a-b\ast q$$$$\because q=\lfloor \frac{a}{b}\rfloor$$$$\therefore a{x}_1+b{y}_1=b{x}_2+(a-\lfloor \frac{a}{b}\rfloor \ast b)y_2$$$$\therefore a{x}_1+b{y}_1=b{x}_2+a{y}_2-\lfloor \frac{a}{b}\rfloor \ast b{y}_2$$$$\therefore a{x}_1+b{y}_1=a{y}_2+b(x_2-\lfloor \frac{a}{b}\rfloor y_2)$$$$\therefore x_1=y_2$$$$\therefore y_1=x_2-\lfloor \frac{a}{b}\rfloor y_2$$
  我们知道一直递归当$b==0$时开始返回，那么当$b==0$时，$x_n$和$y_n$应该是多少呢？首先我们来看一下当$b==0$时的一元二次方程是什么？答案是：$a{x}_n=c$，对于这个结果我们分开讨论。
  当$gcd(a,b)==c$,所以$x_n=1$，为了方便计算，我们令$y_n=0$，到此我们求出$x_n$和$y_n$，当我们求出递归本层的$x_n$和$y_n$，我们就可以计算出上一层的$x_{n-1}$和$y_{n-1}$，依次最终可以求出$x_1$和$y_1$，这既是一元二次方程的整数解。
  当$gcd(a,b)|c$,所以$x_n=\frac{c}{gcd(a,b)}$,$x_n$是一个整数，相当于对最终的解扩大了$\frac{c}{gcd(a,b)}$，所以我们可以先求出当$gcd(a,b)==c$的解，在此基础上扩大$\frac{c}{gcd(a,b)}$倍。
  当$gcd(a,b)$不是$c$的一个因子，那么$x_n=\frac{c}{gcd(a,b)}$,$x_n$不是一个整数，那么最终解非整数，即一元二次方程无整数解。
  到此我们给出了一元二次方程解的所有情况，我们知道当存在整数解时，求得的一元二次方程的解不止一个。我们可以求出其通解为$x=x+b/d\ast k$,$y=y-a/d\ast k$.
下面给出代码：

int exgcd(int a,int b, int &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	int d = exgcd(b,a%b,x,y);
	int x1 = x;
	x = y;
	y = x1-a/b*y;
	return d;
}
int main()
{
	int a,b,c,x,y;
	scanf("%d%d%d",&a,&b,&c);
	int d = exgcd(a,b,x,y);
	if((a==0&&b==0)||c%d!=0)
	{
		printf("无整数解");
	}
	else
	{	
		int n = c/d;
		x *= n;
		y *= n;
		printf("一个整数解为x=%d y=%d\n",x,y);
		printf("通解为x=x+b/d*k,y=y-a/d*k\n");
	}
	return 0;
}

唯一分解定理

定理

  每个大于1的正整数n都可以被唯一地写成素数的乘积，在乘积中的素因子按照非降序排列。正整数n的分解式$n=p_1^{e_1}\ast p_2^{e_2}\ast \cdot \cdot \cdot \cdot \cdot \cdot \ast p_k^{e_k}$，称为n的标准分解式，其中$p_1$,$p_2$,$\cdot \cdot \cdot$,$p_k$是素数，$p_1\le p_2\le \cdot \cdot \cdot \le p_k$，且$e_1$,$e_2$,$\cdot \cdot \cdot$,$e_k$是正整数。

性质

1. 若n的标准素因子分解表达式为$n=p_1^{e_1}\ast p_2^{e_2}\ast \cdot \cdot \cdot \cdot \cdot \cdot \ast p_k^{e_k}$，设d(n)为n的正因子的个数。s(n)为n的所有因子之和，则有$$d(n)=(e_1+1)\ast (e_2+1)\ast \cdot \cdot \cdot \ast (e_k+1)$$$$s(n)=\frac{p_1^{e_1+1}-1}{p_1-1}\ast \frac{p_2^{e_2+1}-1}{p_2-1}\ast \cdot \cdot \cdot \ast \frac{p_k^{e_k+1}-1}{p_k-1}$$
   上边公式不做证明，给出解释。例如$n=12$，$n$的因子有$1,2,3,4,6,12$，通过暴力计算，因子个数$d_1(12)=6$,因子和$s_1(12)=28$，唯一分解式$12=2^2\ast 3^1$，通过公式计算，因子个数$d(12)=(2+1)\ast (1+1)=6$,因子和$s(12)=\frac{2^{2+1}-1}{2-1}\ast \frac{3^{1+1}-1}{3-1}=7\ast 4=28$
2. 设$a=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdot \cdot \cdot \cdot \cdot \cdot \ast p_k^{{\alpha }_k}$,$b=p_1^{{\beta }_1}\ast p_2^{{\beta }_2}\ast \cdot \cdot \cdot \cdot \cdot \cdot \ast p_k^{{\beta }_k}$，则有：$$gcd(a,b)=p_1^{min({\alpha }_1,{\beta }_1)}\ast p_2^{min({\alpha }_2,{\beta }_2)}\ast \cdot \cdot \cdot \ast p_k^{min({\alpha }_k,{\beta }_k)}$$$$lcm(a,b)=p_1^{max({\alpha }_1,{\beta }_1)}\ast p_2^{max({\alpha }_2,{\beta }_2)}\ast \cdot \cdot \cdot \ast p_k^{max({\alpha }_k,{\beta }_k)}$$
3. 如果$a$和$b$是正整数，则:$$lcm(a,b)=a\ast b/gcd(a,b)$$
4. $n!$的素因子分解中的素数p的幂为：$$\lfloor \frac{n}{p}\rfloor +\lfloor \frac{n}{p^2}\rfloor +\lfloor \frac{n}{p^3}\rfloor +\cdot \cdot \cdot$$