欧几里得算法
一、 欧几里得算法
a).先假设x>y,x和y的最大公约数用f(x,y)表示
b).假设
x/y = a;
x%y = b;
=> a*y + b = x
将上面那个式子调换一下位置得到:b = x – a*y;
所以欧几里得的公式为
f(x,y) = f(y,x%y);
所以这个算法的实现也就是不停的迭代,直到找出了x%y等于0时,则停止迭代,那个时候最大公约数也就是y了(因为x%y都等于0了,所以x和y的最大公约数也就是y本身了)。
// An highlighted block
public static int gcd(int x, int y){
//不能输入为0
if(x == 0 || y == 0){return 0;}
//添加一个判断保证x > y
if(x < y){
int temp = x;
temp = y;
y = x;
}
//算法实现
if(x%y == 0){
return y;
}else
{
return gcd(y,x%y);
}
}
二、模运算基本性质
(1)若p|(a-b),则a≡b (% p)。例如 11 ≡ 4 (% 7), 18 ≡ 4(% 7)
(2)(a % p)=(b % p)意味a≡b (% p)
(3)对称性:a≡b (% p)等价于b≡a (% p)
(4)传递性:若a≡b (% p)且b≡c (% p) ,则a≡c (% p)
符合结合律,交换律,分配律
三、欧拉函数计算
欧拉定理:a^(φ(m)) ≡1(mod m) (a与m互质)
【欧拉函数的通式】
φ(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn)
p为素数
例如
6 = 2 * 3 ==> φ(6)=6 * ( 1 - 1/2 ) * ( 1 - 1/3) = 2
27=3 * 3 * 3 ==> φ(27)=27 * ( 1 - 1/3 ) = 18
技巧:
n为素数时φ(n) = n - 1
如果n为任意两个数a和b的积,那么φ(ab)=φ(a)φ(b)
φ(pi)=pφ(i) (当p%i==0时)
φ(p*i)=(p-1)*φ(i) (当p%i!=0时)