数学问题-欧几里得算法扩展

欧几里得算法-计算a,b最大公约数

计算方法

1.当b≠0时：a=b,b=a%b;
2.当b=0时，最大公约数为a;

算法

int gcd(int a,int b){
	return !b?a:gcd(b,a%b);
}

扩展-计算ax+by=gcd(a,b)

方法(其中a,b是在不断变化的）

1. 首先计算出gcd=gcd(a,b);
2. $a{x}_1+b{y}_1=gcd$且$b{x}_2+(a\%b)y_2=gcd$,所以$a{x}_1+b{y}_1=b{x}_2+(a\%b)y_2$，又$a\%b=a-(a/b)\ast b$,所以整理得到$a{x}_1+b{y}_1=a{y}_2+b(x_2-(a/b)y_2)$,所以$x_1=y_2,y_1=x_2-(a/b)\ast y_2$;
3. 这样通过递归计算就可以得到一组解$x,y$,此时显然$a\ast (x+s_1)+b\ast (y-s_2)=gcd$,化简得到$s_1/s_2=b/a=(b/gcd)/(a/gcd)$,所以$x_{all}=x+b/gcd\ast k_i,y_{all}=y-a/gcd\ast k_i,i=0,1,2,3,...$，可以看出$x_{all}，y_{all}$分别以$b/gcd,a/gcd$为周期，所以最小的正整数$x=(x\%(b/gcd)+(b/gcd))\%(b/gcd)$,(这里是为了保证x为正整数,因为等式右侧的$x$可能为负数)；

代码

int getgcd(int a,int b,int &x,int &y){//这里找出x,y与a和b的最大公约数
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	int g = getgcd(b,a%b,x,y);//获得最大公约数
	//递归找到最初的x与y
	int temp=x;
	x=y;
	y=temp-(a/b)*y;
	return g;
}

进一步扩展-计算ax+by=c

方法

1. 首先我们计算$ax+by=gcd(a,b)=gcd$得$x$;
2. 然后两边同时乘以$c/gcd$得$ax\ast (b/gcd)+by\ast (b/gcd)=gcd\ast (c/gcd)$,（注意到$c\%gcd==0$才成立）；
3. 所以我们由上一个问题可以得到该问题的解$x_{all}=x\ast (b/gcd)+b/gcd\ast k_i,y_{all}=y\ast (b/gcd)-a/gcd\ast k_i,i=0,1,2,3,...$;

再次扩展-同余式计算

概念：

a与c模m同余：即（a-b）%m=0;可以写成：ax=c(mod m);

方法

1. $(ax-m)\%c=0$,可以转化为$ax+ck=m$,这样就转化成为了上一个问题同样有$m/gcd(a,m)=0$。此时$x_{all}=x+(m/gcd)\ast k$其中，k的取值只能在0,1，…，m-1中取，否则会被在这个范围内的数取代.

再次扩展-逆元

概念

若a*b=1(mod m),即$(a\ast b)\%m=1$,这样即是同余式中的余数1；这里注意要求gcd(a,m)=1,否则无解，因为此时c/gcd(a,m)=1/gcd(a,m),由第二个问题可得，必然要求gcd(a,m)=1.

扩展

当模m为素数时，所要求的a*b=1(mod m),已知a,且a≠0,则$b=a^{m-1}$。