ch3(数学)——约数

1.自然语言描述
试除法求约数：从1枚举到sqrt(n)（自然数n的算术平方根）若能够整除n则为n的约数。
约数个数：自然数n可以分解质因数为(p1^a1) * (p2^a2) * (p3^a3) * … * (pk^ak)（p为质数，a为指数部分）自然数n的约数个数为(a1+1) * (a2+1) * (a3+1) * … * (ak+1)
约数之和：自然数n的约数之和为 ( p1^0 + p1^1 + … + p1^a1 ) * ( p2^0 + p2^1 + … +p2^a2 ) * … * ( pk^0 + pk^1 + … + pk^ak )（展开后的每一项都是一个约数）
最大公约数：辗转相除法，除数为1时为递归边界。

2.代码描述
题目：Acwing.869 试除法求约数题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>

using namespace std;

vector<int> calc(int x)
{
    vector<int> res;
    for(int i=1;i<=x/i;i++){//从1开始枚举x的所有约数,最大约束为sqrt(x),像试除法判定质数枚举时一样,采用了等价写法
        if(x%i==0){
            res.push_back(i);
            if(i!=x/i)//sqrt(x)只需要一个,这里避免装入2个
                res.push_back(x/i);
        }
    }
    
    sort(res.begin(),res.end());
    
    return res;
}

int main(void)
{
    int n;
    cin>>n;
    
    while(n--){
        
        int x;
        cin>>x;
        
        vector<int> ans=calc(x);
        
        for(auto t:ans)
            cout<<t<<' ';
        cout<<endl;
    }
    
    return 0;
}

题目：Acwing.870 约数个数题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <unordered_map>

using namespace std;

const int MOD=1e9+7,MAXN=2000;
typedef long long LL;
int n;

int main(void)
{
    cin>>n;
    
    unordered_map<int,int> primes;//记录质数幂
    LL ans=1;
    while(n--){
        LL x;
        cin>>x;
        
        for(int i=2;i<=x/i;i++){//分解质因数,记录每个质数的幂
            if(x%i==0){
                while(x%i==0){
                    x/=i;
                    primes[i]++;
                }
            }
        }
        if(x>1)
            primes[x]++;
    }
    
    for(auto t:primes)//约数个数=(a1+1)*(a2+1)*(a3+1)*...*(ak+1),a为每个质数的次数
          ans=ans*(LL)(t.second+1)%MOD;
    
    cout<<ans<<endl;
    
    return 0;
}

题目：Acwing.871 约数之和题目链接

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <unordered_map>

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

int n;

int main(void)
{
    cin>>n;
    
    unordered_map<int,int> primes;
    while(n--){
        int x;
        cin>>x;
        
        for(int i=2;i<=x/i;i++)
            if(x%i==0){
                while(x%i==0){
                    x/=i;
                    primes[i]++;
                }
            }
        if(x>1)
            primes[x]++;
    }
    
    LL ans=1;
    for(auto it:primes){//约数之和=(p1^0+p1^1+...+p1^a1)*(p2^0+...+p2^a2)*...*(pk^0+..pk^ak)
        int p=it.first,a=it.second;
        LL t=1;
        while(a--)
            t=(t*p+1)%MOD;//(pi^0+pi^1+...+pi^ai)=pi*(pi*(...pi*(pi+1)))
        ans=ans*t%MOD;
    }
    
    cout<<ans<<endl;
    
    return 0;
}

题目：Acwing.872 最大公约数题目链接

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

using namespace std;

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

int main(void)
{
    int n;
    cin>>n;
    
    for(int i=0;i<n;i++){
        int x,y;
        cin>>x>>y;
        
        if(x<y)
            swap(x,y);
        cout<<gcd(x,y)<<endl;
    }
    
    return 0;
}