约数之和，约束个数

1.求所有约数

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

 

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
void get_prime(int x)
{
    vector<int> res;
    for(int i=1;i <= x / i;i++)
    {
        if(x % i == 0)
        {
            res.push_back(i);
            if(i != x/i) res.push_back(x/i);
        }
    }
    sort(res.begin(),res.end());
    for(auto re :res) cout<<re<<" ";
    cout<<endl;
}

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

2. 约数个数

公式：所有质因数个数的乘积。[(a1+1)*(a2+1)*(a3+1)*......*(ak+1)]

一个正整数N = p1^a1 * p2^a2 * ...... pk^ak 

p1,p2,p3,......,pk 是正整数的所有质因子

a1,a2,a3,......,ak 是对应的质因数存在的个数 （6 的质因子是 2 和3, 所以它有1个2，1个3，也就是2^1 * 3^1)

公式证明： 一个约数 X = p1^b1*p2^b2*......pk^bk (0<= bi <= ai)

(b1,b2,....,bk)  bi 只要在每一个【0，a】区间内任选一个值，都可以得到一个值，这个值就是N的约数。

6的约数 1， 2 ，3 ，6；

1=2^0 * 3^0

2= 2^1 * 3^0

3=2^0 * 3^1

6=2^1 * 3^1

每一个质因子有a+1（[0,a] 闭区间）种选法，所以k个质因子有(a1+1)*(a2+1)*(a3+1)*......*(ak+1)种选法，得到的选法种数就是N的约数个数

#include<iostream>
#include<unordered_map>
#include<algorithm>
using namespace std;

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

int main()
{
    int n;
    cin>>n;
    
    unordered_map<int,int> prime;
    
    while(n--)
    {
        int a;
        cin>>a;
        
        for(int i=2; i <= a / i; i++)
        {
            while(a % i == 0)
            {
                a /= i;
                prime[i] ++;
            }
        }
        if(a > 1) prime[a] ++;
    }
    
    LL res = 1;
    for(auto p : prime) 
    {
        res = res * (p.second + 1) % mod;
        cout<<p.first<<" * "<<p.second<<endl;
    }
    cout<<res<<endl;
    return 0;
}

 

3.约数之和

公式：（p1^0 + p1^1 + ..... +p1^k)(p2^0 + p2^1 + ..... +p2^k)......(pk^0 + pk^1 + ..... +pk^k)

证明：把公式按乘法分配律展开 （从上述每一个括号中选一个数相乘）+（）+（） + .....+ ()

可以看出展开之后每个括号都是一个约数，共有约数个（上面求过的约数个数）括号。

#include<iostream>
#include<unordered_map>
#include<algorithm>
using namespace std;
const int mod = 1e9 + 7;
typedef long long LL;

int main()
{
    int n;
    cin>>n;
    
    unordered_map<int,int> prime;
    
    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)
    {
        int a = p.first  ,b = p.second;
        LL t = 1;
        while(b--) t = (t * a + 1) % mod;//求每一项p的乘积
        res = res * t % mod;
    }
    cout<<res<<endl;
    return 0;
}

p^0+p^1+p^2+.....+p^k =

while循环

第一次循环 t = p+1

第二次 t = (p+1) *p + 1

......