素数分解
每一个数都可以分解成素数的乘积,例如 84 = 22 * 31 * 50 * 71 * 110 * 130 * 170 * …
整除
令 x = 2m0 * 3m1 * 5m2 * 7m3 * 11m4 * …
令 y = 2n0 * 3n1 * 5n2 * 7n3 * 11n4 * …
如果 x 整除 y(y mod x == 0),则对于所有 i,mi <= ni。
最大公约数最小公倍数
x 和 y 的最大公约数为:gcd(x,y) = 2min(m0,n0) * 3min(m1,n1) * 5min(m2,n2) * …
x 和 y 的最小公倍数为:lcm(x,y) = 2max(m0,n0) * 3max(m1,n1) * 5max(m2,n2) * …
生成素数序列
埃拉托斯特尼筛法在每次找到一个素数时,将能被素数整除的数排除掉。
function countPrimes(n) {
var notPrimes = [];
for (var i = 0; i < n+1; i++) {
notPrimes.push(false);
}
var count = 0;
for(var j = 2; j < n; j++) {
if(notPrimes[j]) {
continue;
}
count++;
for(var k = j*j; k < n; k += j) {
notPrimes[k] = true;
}
}
return count;
}
def notPrimes(n):
notPrimes = [False for _ in range(n+1)]
count = 0
for i in range(2, n+1):
if notPrimes[i]:
continue
count += 1
for j in range(i*i, n+1, i):
notPrimes[j] = True
return count
最大公约数
function gcd(a, b) {
return b == 0? a :gcd(b, a%b);
}
def gcd(a, b):
return a if b == 0 else gcd(b, a%b)
最小公倍数为两数的乘积除以最大公约数。
function lcm(a, b){
return a*b/gcd(a, b)
}
def lcm(a, b):
return a*b//gcd(a, b)
使用位操作和减法求解最大公约数
对于 a 和 b 的最大公约数 f(a, b),有:
- 如果 a 和 b 均为偶数,f(a, b) = 2*f(a/2, b/2);
- 如果 a 是偶数 b 是奇数,f(a, b) = f(a/2, b);
- 如果 b 是偶数 a 是奇数,f(a, b) = f(a, b/2);
- 如果 a 和 b 均为奇数,f(a, b) = f(b, a-b);
乘 2 和除 2 都可以转换为移位操作。
function gcd(a,b) {
//使用位操作和减法求最大公约数
if(a<b) {
return gcd(b,a);
}
if(b == 0) {
return a;
}
function isEven(n){
if(typeof n !== 'number'