数学专题（一）

数学问题总结

1. 数论
2. 组合计数
3. 高斯消元
4. 简单博弈论
   一、判素数：
5. 试除法
6. 埃氏筛法：O(nlognlogn)
7. 线性筛法：O(n)
   二、分解质因数
8. 试除法求约数
9. 约数个数：INT_MAX范围内约数个数1600
10. 约数之和
    三、最大公约数（辗转相除法）

870. 约数个数

#include <iostream>
#include <algorithm>
#include <unordered_map>
#include <vector>

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;
}

AcWing 871. 约数之和

#include <iostream>
#include <algorithm>
#include <unordered_map>
#include <vector>

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 prime : primes) 
    {
        LL a = prime.first, b = prime.second;
        LL t = 1;
        while (b -- ) t = (t * a + 1) % mod;
        res = (res * t) % mod;
    }

    cout << res << endl;

    return 0;
}

872. 最大公约数

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

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

int main()
{
    int n;
    cin >> n;
    
    while (n -- )
    {
        int a, b;
        scanf("%d %d", &a, &b);
        
        cout << gcd(a, b) << endl;
    }
    return 0;
}