数论 - 约数 - 试除法求约数 + 约数个数 + 约数之和

数论 - 约数 - 试除法求约数 + 约数个数 + 约数之和 + 最大公约数

1、 试出法求约数

给定n个正整数ai，对于每个整数ai,请你按照从小到大的顺序输出它的所有约数。

输入格式
第一行包含整数n。

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

输出格式
输出共n行，其中第 i 行输出第 i 个整数ai的所有约数。

数据范围
1≤n≤100,
2≤ai≤2∗10⁹

输入样例：
2
6
8
输出样例：
1 2 3 6 
1 2 4 8 

---

分析：

$O(\sqrt{n})的做法，枚举一个较小的因子，另一个直接用原数去除。$

$注意完全平方数只要加入一个因子。$

代码：

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

using namespace std;

const int N=110;

int n;

void solve(int x)
{
    vector<int> p;
    for(int i=1;i<=x/i;i++)
        if(x%i==0) 
        {
            p.push_back(i);
            if(i!=x/i) p.push_back(x/i);
        }
    
    
    sort(p.begin(),p.end());
    
    for(int i=0;i<p.size();i++) cout<<p[i]<<' ';
    cout<<endl;
}

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

---

2、约数个数

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

输入格式
第一行包含整数n。

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

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

数据范围
1≤n≤100,
1≤ai≤2∗10⁹

输入样例：
3
2
6
8
输出样例：
12

---

$$\begin{gathered} 正整数a的约数个数： \\ 由算术基本定理，设A=P_1^{a_1}{P}_2^{a_2}...P_k^{a_k}，则A的约数个数为(a_1+1)(a_2+1)...(a_k+1)。 \end{gathered}$$

$$\begin{gathered} 说明： \\ 设A的一个约数A^{\prime }=P_1^{b_1}{P}_2^{b_2}...P_k^{b_k}，其中b_i\in [0,a_i]，i\in [1,k]。 \end{gathered}$$

$每一个k元组(b_1,b_2,...,b_k)均能对应A的一个约数，对于每个b_i都有a_i+1种取值，$

$根据乘法原理，不同的k元组的个数即A的所有约数的个数为C_{a_1+1}^1{C}_{a_2+1}^1...C_{a_k+1}^1=(a_1+1)(a_2+1)...(a_k+1)。$

具体落实:

$问题转化分解A的所有质因数并保存每个质因数的指数，再用公式。$

代码：

#include<iostream>
#include<unordered_map>

#define ll long long
#define x first
#define y second

using namespace std;

const int mod=1e9+7;

int n;
unordered_map<int,int> P;

int main()
{
    cin>>n;
    while(n--)
    {
        int x;
        cin>>x;
        for(int i=2;i<=x/i;i++)
            while(x%i==0)
            {
                x/=i;
                P[i]++;
            }
        if(x>1) P[x]++;
    }
    
    ll res=1;
        for(auto p: P)
            res=res*(p.y+1)%mod;
        
    cout<<res<<endl;
    
    return 0;
}

---

3、约数之和

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

输入格式
第一行包含整数n。

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

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

数据范围
1≤n≤100,
1≤ai≤2∗10⁹

输入样例：
3
2
6
8
输出样例：
252

---

$$\begin{gathered} 正整数a的约数个数： \\ 由算术基本定理，设A=P_1^{a_1}{P}_2^{a_2}...P_k^{a_k}，则A的约数之和为{\prod }_{i=1}^k{\sum }_{j=0}^{a_i}{P}_i^j。 \end{gathered}$$

$即(P_1^0+P_1^1+...+P_1^{a_1})(P_2^0+P_2^1+...+P_2^{a_2})...(P_k^0+P_k^1+...+P_k^{a_k})$

$说明：$
$设A的一个约数A^{\prime }=P_1^{b_1}{P}_2^{b_2}...P_k^{b_k}，其中b_i\in [0,a_i]，i\in [1,k]。$

$每一个k元组(b_1,b_2,...,b_k)均能对应A的一个约数，对于每个质因子P_i^x，x均有a_i+1种选法，$

现 需 在 k 个 形 如 { P i 0 ， P i 1 ， . . . ， P i a i } 的 集 合 中 选 一 个 p i x ， i ∈ [ 1 , k ] , x ∈ [ 0 , a i ] ， 现需在k 个形如\{P_i^0，P_i^1，...，P_i^{a_i}\}的集合中选一个p_i^{x}，i∈[1,k],x∈[0,a_i]， 现需在k个形如{Pi0​，Pi1​，...，Piai​​}的集合中选一个pix​，i∈[1,k],x∈[0,ai​]，

$选k次，每次对应一个约数P_1^{x_1}{P}_2^{x_2}...P_k^{x_k}，最终求和结果为:$ $$(P_1^0+P_1^1+...+P_1^{a_1})(P_2^0+P_2^1+...+P_2^{a_2})...(P_k^0+P_k^1+...+P_k^{a_k})$$

具体落实：

$①、首先对A分解质因数，并记录每个质因数的指数，用一个哈希表存储。$

$$\begin{gathered} ②、求P_i^0+P_i^1+P_i^2+...+P_i^{a_i}可以用迭代来做，t=t\ast P_i+1，迭代a_i+1次。 \\ 当然也可以用等比数列求和公式、分治等方法。 \end{gathered}$$

代码：

#include<iostream>
#include<unordered_map>

#define ll long long
#define x first
#define y second

using namespace std;

const int mod=1e9+7;

int n;
unordered_map<int,int> prime;

int main()
{
    cin>>n;
    while(n--)
    {
        int x;
        cin>>x;
        for(int i=2;i<=x/i;i++)
            while(x%i==0)
            {
                x/=i;
                prime[i]++;
            }
        if(x>1) prime[x]++;
    }
    
    ll res=1;
    for(auto p : prime)
    {
        ll t=0;
        for(int i=0;i<=p.y;i++)
            t=(t*p.x+1)%mod;
        res=res*t%mod;
    }
    
    cout<<res<<endl;
    
    return 0;
    
}

---