约数个数 + 约数之和（ 数学知识 + 约数 + 试除法 + 因式分解 ）

约数个数 

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

输入格式

第一行包含整数 n。

接下来 n 行，每行包含一个整数 ai。

输出格式

输出一个整数，表示所给正整数的乘积的约数个数，答案需对 109+7 取模。

数据范围

1≤n≤100,
1≤ai≤2×109

输入样例：

3
2
6
8

输出样例：

12

 

 

 

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 110, mod = 1e9 + 7;
int n;
int main(){
    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] ++ ; // 最后 x > 1 是是特判大于 sqrt(x) 的那个唯一的质因子
    }
    LL res = 1;
    for(auto p : primes)
        res = res * (p.second + 1) % mod;
    cout << res << endl;
    return 0;
}

---

约数之和(基础版）

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

输入格式

第一行包含整数 n。

接下来 n 行，每行包含一个整数 ai。

输出格式

输出一个整数，表示所给正整数的乘积的约数之和，答案需对 109+7 取模。

数据范围

1≤n≤100,
1≤ai≤2×109

输入样例：

3
2
6
8

输出样例：

252

#include<bits/stdc++.h>
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){
        LL a = p.first, b = p.second;
        LL t = 1;
        while(b -- ) t = (t * a + 1) % mod;
        res = res * t % mod;
    }
    cout << res << endl;
    return 0;
}

---

约数之和（提高版）

假设现在有两个自然数 A 和 B，S 是 A^B 的所有约数之和。

请你求出 Smod9901 的值是多少。

输入格式

在一行中输入用空格隔开的两个整数 A 和 B。

输出格式

输出一个整数，代表 Smod9901 的值。

数据范围

0≤A,B≤5×107

输入样例：

2 3

输出样例：

15

注意: A 和 B 不会同时为 0。

 再次预习一下哈哈哈哈

 

#include<bits/stdc++.h>
#define x first
#define y second
using namespace std;
typedef long long LL;
const int mod = 9901;
int a, b;
unordered_map<int, int> primes;
void divide(int n){
    for(int i = 2; i <= n / i; i ++ ){
        if(n % i == 0){
            while(n % i == 0){
                primes[i] ++ ;
                n /= i;
            }
        }
    }
    if(n > 1) primes[n] ++ ;
}
int qmi(int a, int b){
    int res = 1;
    while(b){
        if(b & 1) res = (LL) res * a % mod;
        b >>= 1;
        a = (LL) a * a % mod;
    }
    return res;
}
int sum(int p, int k){
    if(k == 1) return 1;
    if(k % 2 == 0){
        return (LL) (qmi(p, k / 2) + 1) * sum(p, k / 2) % mod;
    }
    return (LL) (qmi(p, k - 1) + sum(p, k - 1)) % mod;
}
int main(){
    cin >> a >>b;
    divide(a);
    int res = 1;
    for(auto it : primes){
        int p = it.x, k = it.y * b;
        res = (LL) res * sum(p, k + 1) % mod;
    }
    if(!a) res = 0;
    cout << res << endl;
    return 0;
}