【4.数学知识-模板】2.约数（模板）

一. 试除法求约数

题目来源：试除法求约数

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

using namespace std;

int n;

void 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)
            {  
                //避免i==n/i重复放入 
                res.push_back(n / i);
            }
        }
    }

    sort(res.begin(), res.end());
    
    
    for (auto item : res)
    {
        cout << item << " ";
    }
    puts("");
}



int main()
{
    cin >> n;
    while (n--)
    {
        int x;
        cin >> x;
        get_divisors(x);
    }
    return 0;
}

二.约数个数

题目来源：约数个数

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

using namespace std;

const int mod=1e9+7;

unordered_map<int,int> primes;//int 质数;int 指数

//分解质因数
void divide(int n)
{   
    for (int i = 2; i <= n / i; i ++ )
    {
        if(n%i==0)
        {
            while (n % i == 0)
            {
                n /= i;
                primes[i] ++ ;
            }
        }
    }
    if (n > 1) primes[n] ++ ;
}



int main()
{
    int n;
    cin>>n;

    while(n--)
    {
        int x;
        scanf("%d",&x);
        divide(x);
    }


    long long res=1;
    
    for(auto p:primes) res=res*(p.second+1)%mod;//将统计出来的数按照由图中公式所得出来的结论得出答案

    printf("%lld\n",res);

    return 0;
}

三.约数之和

题目来源：约数之和

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

using namespace std;

typedef long long LL;

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

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)
            {
                n /= i;
                primes[i] ++ ;
            }
        }
        
        if (n > 1) primes[n] ++ ;
}



int main()
{
    int n;
    cin >> n;

    

    while (n -- )
    {
        int x;
        cin >> x;
        divide(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<bits/stdc++.h>
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;
}