算术基本定理:任何一个大于 1 的整数都能唯一分解为有限个质数的乘积。
试除法:结合质数质数判定的“试除法”和质数筛选的“Eratosthenes筛法”,我们可以枚举 2~sqrt(n) 中的每个数 d ,若 d 能整除 N,则从 N 中除掉所有的因子 d ,同时累计除去因子 d 的个数。
因为一个合数的因子一定在扫描到这个合数之前就从 N 中被除掉了,所以在上述过程中能整除 N 的一定是质数。
特别地,若 N 没有被 2~sqrt(N) 的数整除,则 N 是质数,无需分解。
时间复杂度为O(sqrt(n))。Pollard's Rho算法的时间复杂度为 O(sqrt(sqrt(n))),但是不会0.0。
代码实现:
const int maxn=10005;
int c[maxn];
int p[maxn], m;
void divide(int n){
m = 0;
for(int i = 2; i <= sqrt(n); i++){
if(n % i == 0){// i 是质数
p[++m] = i, c[m] = 0;
while(n % i == 0){
n /= i, c[m]++;
}
}
}
if(n > 1){// n 是质数
p[++m]=n, c[m] = 1;
}
for(int i = 1; i <= m; i++)
printf("%d^%d\n", p[i], c[i]);
}