小菜想破茧成蝶为中佬——算法基础笔记（AcWing算法基础课第四章数学知识）

前言：由于方便记下我学习算法的一路历程，特编辑在CSDN平台，如有问题，欢迎大家指正！

今日2022.10.18，日后我会多在此篇文章进行更新，加油，希望来年我能顺利到达我的彼岸！

终极目标：2023.12.23考研——加油，星光不问赶路人,时光不负有心人！

1.试除法判定质数

bool check(int a)
{
    if(a < 2) return false;
    for(int i = 2; i <= a / i; i ++)
    {
        if(a % i == 0)
            return false;
    }
    return true;
}

补充

• 素数即质数，质数是指在大于1的自然数中，除了1和它本身以外不再有其他因数的自然数。

• 约数即因数，是指整数a除以整数b(b≠0) 的商正好是整数而没有余数，我们就说b是a的因数。a称为b的倍数，b称为a的约数

• “1”不是质数。

2.分解质因数

void check(int x)
{
    for(int i = 2; i <= x / i; i ++)
    {
        if(x % i == 0)
        {
            int cnt = 0;
            while(x % i == 0)
            {
                x /= i;
                cnt ++;
            }
            printf("%d %d\n", i, cnt);
        }
    }
    if(x > 1) printf("%d %d\n", x, 1);
    puts("");
}

3.1埃氏筛法

void prime(int n)
{
    for(int i = 2; i <= n; i ++)
    {
        if(!st[i])
        {
            a[cnt ++] = i;
            for(int j = i + i; j <= n; j += i)
                st[j] = true;
        }
    }
}

补充

• 比较直接，2是最小的质数，并在自然数轴上将 i 的 倍数全部进行标记，因此if判断则为质数

3.2线性筛法

void prime(int n)
{
    for(int i = 2; i <= n; i ++)
    {
        if(!st[i]) primes[cnt ++] = i; // 如果st[]数组没有标记，那说明此刻 i 为质数，则加入primes数组，并cnt计数
       
        for(int j = 0; primes[j] <= n / i; j ++) 
        {
            st[primes[j] * i] = true; //如果满足for循环，则首先质数primes[0] = 2 * i。标记为因数，也称约数。注意，这里就是线性筛的核心，即用最小质数不断枚举 i 去筛
            if(i % primes[j] == 0) break; //如果该 i 可以整除primes数组的质数则结束内循环
        }
    }
}

3.3试除法求约数

vector<int> get_divisors(int n)
{
    vector<int> res;
    
    for(int i = 1; i <= n / i; i ++)
    {
        if(n % i == 0)
        {
            res.push_back(i);
            if(i != n / i) res.push_back(n / i);
        }
    }
    sort(res.begin(), res.end());
    return res;
}

补充

• 利用vector变长数组方便排序以及存储。