Number Theory 4.1 Divisibility

Number Theory 4.1 Divisibility

哈哈，第一篇，试写....

最熟悉的一个概念，最大公约数gcd，k/m 表示k能除尽m
注意“k能除尽m”和“m是k的倍数”并不完全一样，如k=0
gcd(m, n) = max{k| k/m and k/n}; ①
欧几里德算法的递归形式：
gcd(0, n) = n;
gcd(m, n) = gcd(n%m, m); m > 0 ②

还有一种扩展形式：
xm + yn = gcd(m, n);

这个x，y得出有个很美妙的递归求法。假定 n>m≥0，x₁m + y₁n = gcd(m, n);
1）当m=0时，x₁=0，y₁=1是它其中的一组可能解。
2）假设 r=n%m，则有 gcd(r, m) = x₂r + y₂m;
根据②可知 gcd(r, m) = gcd(m, n);
则有 x₁m + y₁n = x₂r + y₂m; ③
因为 r = n%m = n - [n/m]*m;
③可得 x₁m + y₁n = x₂(n - [n/m]*m) + y₂m = x₂n + (y₂ - [n/m]*x₂)m;
则有解 y₁ = x₂;  x₁ = y₂ - [n/m]*x₂;
可以看出x₂, y₂交叉赋值给y₁, x₁。

这样就得到了扩展欧几里德算法。

//  extended_euclid(m, n) = mx + ny
int  extended_euclid( int  m,  int  n,  int   & x,  int   & y)
... {
    if(m == 0) ...{
        x = 0;
        y = 1;
        return n;
    }
    int gcd = extended_euclid(n%m ,m ,y ,x);
    x -= n/m * y;
    return gcd;
}

 

还有个公式就是

∑a_m = ∑a_n/m (n > 0, m/n)

还不大清楚它的用途，下次碰到了再看下吧。