数学知识：约数

试除法求约数：

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

using namespace std;
//时间复杂度:O(n * sqrt(a))

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

int main()
{
    int n;
    scanf("%d", &n);
    while(n--)
    {
        int x;
        cin >> x;
        auto res = get_divisors(x);
        for(auto t : res)
            cout << t << ' ';
        cout << endl;
    }
    return 0;
}

求约数个数：

一个数

N = 

则约数个数ans =

推导：

N的任何一个约数d = 

每一项的不同，那么每一项的d就不同，每一个都有0 ~ 种方式

由乘法原理得到由ans个约数

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

using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;

//int范围内约数最多是1500多个
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 ans : primes)
        res = res * (ans.second + 1) % mod;
    cout << res << endl;
    return 0;
}

约数之和：

若： N =  

则约数之和

公式：

while (b -- ) t = (t * a + 1) % mod;

t = t∗p+1
t = 
t = 
t = 
 ……
t = 

#include <iostream>
#include <unordered_map>

using namespace std;
typedef long long LL;
const int 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 的 最大公约数 d
则有 gcd(a,b) = gcd(b,a%b)
证明：
    设a%b = a - k*b 其中k = a/b(向下取整)
    若d是(a,b)的公约数 则知 d|a 且 d|b 则易知 d|a-k*b 故d也是(b,a%b) 的公约数
    若d是(b,a%b)的公约数 则知 d|b 且 d|a-k*b 则 d|a-k*b+k*b = d|a 故而d|b 故而 d也是(a,b)的公约数
    因此(a,b)的公约数集合和(b,a%b)的公约数集合相同 所以他们的最大公约数也相同

#include <iostream>

using namespace std;

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

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