约数之和,约束个数

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

......

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值