质数与约数

质数与约数

1. 算法分析

1.1 基本概念

质数: 作为最小单位，无法再进行拆分
约数: 8的约数有1，2，4，8

1.2 常用公式和性质

1.2.1 质数

重要结论

1. int范围内最多只有10个不重复质因子，且所有质数的次数总和不能超过30
2. 素数的大小:x/lnx,素数的个数:xlnx
3. 任何一个合数n一定存在一个<=sqrt(n)的质因子

1.2.2 约数

常用公式

1. n的约数个数: 对于一个大于1正整数可以分解质因数：n=p1^a1 * p2^a2 * … *pk^ak,则n的正约数个数为:f(n) = (a1+1)(a2+1)…(ak+1)
2. n的约数和：对于一个大于1正整数n可以分解质因数：n=p1^a1 * p2^a2 * … *pk^ak,那么n的(a₁+1)(a₂+1)(a₃+1)…(ak+1)个正约数的和为：f(n)=(p1⁰ + p1¹ + p1² + … + p1^a1)(p2⁰ + p2¹ + p2²+ … + p2^a2)…(pk⁰ + pk¹ + pk² +…+ pk^ak），其中：计算每个(pi⁰ + pi¹ + pi² + … + pi^ai)时可以借助秦九韶公式优化运算
3. 1~N中所有数的约数个数和
   N/1+N/2+…+N/N = N(lnN + C) ,C= 0.57721566490153286060651209
4. 1~N中所有数的约数和:
   1 * N / 1 + 2 * N / 2 + … + N * N / N

重要结论

1. int范围内拥有最多约数的那个数字拥有1600个约数
2. 1~N中任何数的不同质因子不会超过10个，且所有质因子的指数总和不超过30

2. 板子

2.1 质数

2.1.1 素数判定

1. 试除法判定素数(小素数)

// 判断是否为素数
bool is_prime(int x) {
   
    if (x < 2) return false;  // 1不是素数
    for (int i = 2; i <= x / i; ++i)  // 从2循环到sqrt(x),但是循环进行条件不写i <= sqet(x),因为太慢；不写i * i <= x，因为可能有溢出的风险
        if (x % i == 0) return false;
    return true;
}

2. miller-rabin判断素数(大素数)

#include <bits/stdc++.h>
using namespace std ;
typedef unsigned long long ll;

ll multi(ll a, ll b, ll mod) {
   
    ll ret = 0;
    while(b) {
   
        if(b & 1) ret = ret + a;
        if(ret >= mod) ret -= mod;

        a = a + a;
        if(a >= mod) a -= mod;
        b >>= 1;
    }
    return ret;
}

ll quick_pow(ll a, ll b, ll mod) {
   
    ll ret = 1;
    while(b) {
   
        if(b & 1) ret = multi(ret, a, mod);
        a = multi(a, a, mod);
        b >>= 1;
    }
    return ret;
}

bool Miller_Rabin(ll n) {
   
    ll u = n - 1, pre, x;
    int i, j, k = 0;
    if(n == 2 || n == 3 || n == 5 || n == 7 || n  == 11) return true;
    if(n == 1 || (!(n % 2)) || (!(n % 3)) || (!(n % 5)) || (!(n % 7)) || (!(n % 11)))
        return false;
    for(; !(u & 1); k++, u >>= 1);
    srand(time(NULL));
    for(i = 0; i < 5; i++) {
   
        x = rand() % (n - 2) + 2;
        x = quick_pow(x, u, n);
        pre = x;
        for(j = 0; j < k; j++) {
   
            x = multi(x, x, n);
            if(x == 1 && pre != 1 && pre != (n - 1))
                return false;
            pre = x;
        }
        if(x != 1) return false;
    }
    return true;
}

2.1.2 质因数分解

1. 试除法分解质因数(小整数) : O(sqrt(x))

int x;
cin >> x;
for (int i = 2; i <= x / i; ++i) {
     // 因为一个数x的所有因子内只会有一个大于sqrt(x), 所以只要枚举到sqrt(x),而后单独判断大于sqrt(x)的那个即可
    if (x % i == 0) {
     // 如果能够整除
        int s = 0;  // 记录能够整除几次
        while (x % i == 0) {
   
            s++;
            x /= i;
        }
        printf("%d %d\n", i, s);  // 输出因子和次数
    }
}
if (x > 1) printf("%d %d\n", x, 1); // 如果有剩下的数字必然比sqrt(x)大，因为如果比sqrt(x)小，那么在做上一步的枚举时就已经被除掉了，一旦有，就必然次数为1 

2. 更快的质因数分解:先筛素数，然后拿素数判断是不是质因数 O(sqrt(x/lnx))

void get_primes(int n) {
   
    for (int i = 2; i <= n; i ++ ) {
   
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= n / i; j ++ ) {
   
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

int main() {
   
    ...
    for (int i = 0; prime[</