算法：约数，约数个数，约数和，最大公约数(b ? gcd(b, a%b) : a；)

约数

计算约数时，只需计算到$\sqrt{n}$就可以，大于$\sqrt{n}$的约数可以通过用n除以小于$\sqrt{n}$的约数得到。并且需要特判一下边界n / i != i。

void get_divosors(int n){
    vector<int> ans;
    for(int i = 1; i <= n / i; ++i){
        if (n % i == 0) {
            ans.push_back(i);
            if(n / i != i){
                ans.push_back(n / i);
            }
        }
    }
    sort(ans.begin(), ans.end());
    for(int i = 0; i < ans.size(); ++i){
        cout << ans[i] << " ";
    }
    cout << endl;
    return;
}

约数个数

约数个数定理

则n的正约数的个数就是

其中a1、a2、a3…ak是p1、p2、p3，…pk的指数。

首先同上，n可以分解质因数：$n=p_1^{a_1}\times p_2^{a_2}\times p_3^{a_3}\ast \dots \ast p_k^{a_k}$，
由约数定义可知$p_1^{a_1}$的约数有:$p_1^0,p_1^1,p_1^2......p_1^{a1}$ ，共$(a_1+1)$个;同理$p_2^{a_2}$的约数有$(a_2+1)$个…$p_k^{a_k}$的约数有$(a_k+1)$个。
故根据乘法原理：n的约数的个数就是$(a_1+1)(a_2+1)(a_3+1)\dots (a_k+1)$。
用unordered_map存储每一个质因数的指数。
问题描述

给定n个正整数ai，请你输出这些数的乘积的约数之和，答案对10^9+7取模。

#include <iostream>
#include <unordered_map>
using namespace std;
typedef long long LL;

const int N = 110, mod = 1e9 + 7;

int main(){
    int n;
    cin >> n;
    unordered_map<int, int> primes;
    
    while(n--){
        int x;
        cin >> x;
        
        for(int i = 2; i <= x / i;  ++ i)
            while(x % i == 0){
                x /= i;
                primes[i] ++;
            }
            
        if (x > 1) primes[x] ++;
    }
    LL res = 1;
    for (auto p : primes) res = res * (p.second + 1) % mod;
    cout << res << endl;
    return 0;
}

约数之和

约数之和
其中p为质因子，a为质因子的指数
$$f(n)=(p_1^0+p_1^1+p_1^2+\dots p_1^{a1})(p_2^0+p_2^1+p_2^2+\dots p_2^{a2})\dots (p_k^0+p_k^1+p_k^2+\dots p_k^{ak})$$

#include<iostream>
#include<unordered_map>
using namespace std;
const int N = 100;
typedef long long LL;
LL mod = 1e9 + 7;

int main(){
    int n;
    cin >> n;
    unordered_map<int, LL> primes;
    while(n--){
        LL a;
        cin >> a;
        for(int i = 2; i <= a / i; ++i)
            while(a % i == 0){
                a /= i;
                primes[i] ++;
            }
            if(a > 1) primes[a] ++;
    }
    LL ans = 1;
    for (auto p : primes){
        LL cnt = 0;
        LL c = 1;
            LL a = p.first, b = p.second;
            LL t= 1;
            while(b --) t = (t * a + 1) % mod;
            ans = ans * t % mod;
    }
    cout << ans << endl;
    return 0;
}

最大公约数

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