质数+因数

昨天和今天的模拟赛总结就没有记录了。

数学（超级基础的）

一、质数：

1.试除法判断质数：

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

2.分解质因数：

void divide(int x) {
    for (int i = 2; i <= x / i; i++) {
        int p = 0;
        if (x % i == 0) {
            while (x % i == 0) {
                x /= i;
                p++;
            }
            printf("%d %d\n", i, p);
        }
    }
    if (x != 1) printf("%d %d\n", x, 1);
    printf("\n");
}

3.筛质数：

//埃氏筛
int get_primes() {
    int sum = 0;
    for (int i = 2; i <= n; i++) {
        if (!st[i]) {
            sum++;
            for (int j = i + i; j <= n / i; j += i)
                st[i * j] = 1;
        }
    }
    return sum;
}

//线性筛
int get_primes() {
    for (int i = 2; i <= n; i++) {
        if (!st[i]) primes[++cnt] = i;
        for (int j = 1; primes[j] <= n / i; i++) {
            st[primes[j] * i] = 1;
            if (i % primes[j] == 0) break;
        }
    }    
    return cnt;
}

二、约数

1.约数个数：

因为：
$N=p_1^{a_1}\ast p_2^{a_2}\ast p_3^{a_3}\ast p_4^{a_5}\ast ...\ast p_k^{a_k}$
所以：
$ans=(a_1+1)\ast (a_2+1)\ast (a_3+1)\ast ...\ast (a_k+1)$
证明：
因为根据算数基本定理：每个数的因式分解是唯一的，所以 $n$ 的约数就跟 $a_i$ 的选法相对应。由于 $a_i$ 不同， 那么约数 $d$ 就不同。
$a_1$ 一共有 $0-a_1$ 种取法
$a_2$ 一共有 $0-a_2$ 种取法
$a_3$ 一共有 $0-a_3$ 种取法
$a_k$ 一共有 $0-a_k$ 种取法
根据乘法原理，一共有 $ans$ 个约数

unordered_map<int, int> primes;
for (int i = 2; i <= x / i; i ++ )
    while (x % i == 0)
    {
        x /= i;
        primes[i] ++ ;
    }
if (x > 1) primes[x] ++ ;

这里使用unordered_map来存储是为了处理 $n$ 个数的乘积的约数个数。

2.约数之和：

$N=p_1^{a_1}\ast p_2^{a_2}\ast p_3^{a_3}\ast p_4^{a_5}\ast ...\ast p_k^{a_k}$
$ans=(p_1^0+p_1^1+...+p_1^{a_1})\ast (p_2^0+p_2^1+...+p_2^{a_2})\ast ...\ast (p_k^0+p_k^1+...+p_k^{a_k})$
证明：展开就好了。

3.最大公约数：

int gcd(int a, int b)
{
    return b ? gcd(b, a % b) : a;
}

三、小结

数学写题学思想，数学巩固推公式。