约数之和 (简单+普通）

胖川怎么可能不胖

已于 2023-11-03 17:27:49 修改

阅读量419

点赞数

分类专栏：算法整理文章标签：算法数据结构 1024程序员节 c++

于 2023-11-03 17:27:31 首次发布

本文链接：https://blog.csdn.net/qq_73261465/article/details/134198482

版权

算法整理专栏收录该内容

11 篇文章 0 订阅

订阅专栏

简单

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

取模。

输入格式

第一行包含整数 n

。

接下来 n

行，每行包含一个整数 ai

。

输出格式

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

取模。

数据范围

1≤n≤100

,
1≤ai≤2×1e9

输入样例：

输出样例：

思路：

通过唯一分解定理就能得到质因子

一个整数N可以表示为:N=(p1^a1)*(p2^a2)*.......*(pk^ak)

那么此时约数的个数和就能表示为(a1+1)*(a2+1)*...*(ak+1)

具体举个例子

36=2^2 * 3^2
那么此时他的约数可以表示为
2^0 * 3^0 = 1
2^1 * 3^0 = 2
2^0 * 3^1 = 3
2^2 * 3^0 = 4
2^1 * 3^1 = 6
2^0 * 3^2 = 9
2^2 * 3^1 = 12
2^1 * 3^2 = 18
2^2 * 3^2 = 36

而这题让求的是约数之和所以要在以上基础上推到一下约数和的公示就是：

(p1^0+p1^1+....+p1^a1)*(p2^0+p2^1+...+p2^a2)*......*(pk^0+pk^1+...+pk^ak)

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<stack>
#include<string>
#include<algorithm>
#include<unordered_map>
#include<map>
#include<bitset>
#include<cstring>
#include <unordered_set>
//#include<priority_queue>
#include<queue>
#include<deque>
#include<set>
#include<stdlib.h>
#define dbug cout<<"*****hear*****"<<endl;
#define rep(a,b,c) for(ll a=b;a<=c;a++)
#define per(a,b,c) for(ll a=b;a>=c;a--)
#define no cout<<"NO"<<endl;
#define yes cout<<"YES"<<endl;
#define endl "\n"//交互题一定要关！！！！！！！！！
#define lowbit(x) (x&-x)
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
//priority_queue<int,vector<int>,greater<int> >q;
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll, ll> PII;
typedef pair<long double,long double> PDD;
 ll  INF = 0x3f3f3f3f;
//const ll LINF=LLONG_MAX;
// int get_len(int x1,int y1,int x2,int y2)
// {
//   return (x2-x1)*(x2-x1) + (y2-y1)*(y2-y1);
// }
const ll N = 2e5+ 10;
 const ll mod1 =998244353;
 const ll mod2 =1e9+7;
// const ll hash_num = 3e9+9;
ll n,m,ca;
ll arr[N],brr[N],crr[N],drr[N];
//ll h[N],ne[N],e[N],w[N],book[N],idx;
//ll idx;

// void add(ll a, ll b , ll c)
// {
//   e[idx] = b, w[idx] = c,ne[idx] = h[a], h[a] =idx ++ ; 
// }






void solve()
{
  cin >> n;
  unordered_map<ll ,ll>mp;
  rep(i,1,n)
  {
    ll x;
    cin >> x;
    rep(j,2,x/j)
    {
      while(x%j==0)
      {
        x/=j;
        mp[j]++;
      }
    }
    if(x>1)mp[x]++;
  }
  ll res=1;
  for(auto it : mp)
  {
    ll a=it.first,b=it.second;
    ll ant=1;
    while(b--)ant=(a*ant+1)%mod2;
    res= res*ant%mod2;
  }
  cout << res<<endl;
}



int main()
{
   IOS;
   ll _;
    _=1;
    //scanf("%lld",&_);
   // cin>>_;
  
    ca=1;
    while(_--)
    {
      solve(); 
      ca++;
    }    
    return 0;
}

普通

假设现在有两个自然数 A 和 B，S 是 AB

的所有约数之和。

请你求出 Smod9901

的值是多少。

输入格式

在一行中输入用空格隔开的两个整数 A

和 B

。

输出格式

输出一个整数，代表 Smod9901

的值。

数据范围

0≤A,B≤5×1e7

输入样例：

2 3

输出样例：

注意: A

和 B 不会同时为 0。

思路：

通过简单推导出的约数之和的公式为

(p1^0+p1^1+....+p1^a1)*(p2^0+p2^1+...+p2^a2)*......*(pk^0+pk^1+...+pk^ak)

这个时候可以通过分治的方法来解决p^0+p^1+...+p^k-1

在这里写一个sum函数：sum(p, k)表示p^0+p^1+...+p^k-1

当k为偶数的时候

将

p^0+p^1+...+p^k-1

变化为

p^0+p^1+…+p^k/2−1+p^k/2+p^k/2+1+…+p^k−1

再变为

p^0+p^1+…+p^k/2−1+p^k/2（p^0+p^1+…+p^k/2−1）

再变

（1+p^k/2)(p^0+p^1+…+p^k/2−1)

也就是
(p^k/2+1)∗sum(p,k/2)

如果k为奇数的话，那么直接把最后一项拿出来单算，剩下的当偶数算

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<stack>
#include<string>
#include<algorithm>
#include<unordered_map>
#include<map>
#include<bitset>
#include<cstring>
#include <unordered_set>
//#include<priority_queue>
#include<queue>
#include<deque>
#include<set>
#include<stdlib.h>
#define dbug cout<<"*****hear*****"<<endl;
#define rep(a,b,c) for(ll a=b;a<=c;a++)
#define per(a,b,c) for(ll a=b;a>=c;a--)
#define no cout<<"NO"<<endl;
#define yes cout<<"YES"<<endl;
#define endl "\n"//交互题一定要关！！！！！！！！！
#define lowbit(x) (x&-x)
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
//priority_queue<int,vector<int>,greater<int> >q;
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll, ll> PII;
typedef pair<long double,long double> PDD;
 ll  INF = 0x3f3f3f3f;
//const ll LINF=LLONG_MAX;
// int get_len(int x1,int y1,int x2,int y2)
// {
//   return (x2-x1)*(x2-x1) + (y2-y1)*(y2-y1);
// }
const ll N = 2e5+ 10;
 const ll mod1 =998244353;
 const ll mod2 =1e9+7;
// const ll hash_num = 3e9+9;
ll n,m,ca;
ll arr[N],brr[N],crr[N],drr[N];
//ll h[N],ne[N],e[N],w[N],book[N],idx;
//ll idx;

// void add(ll a, ll b , ll c)
// {
//   e[idx] = b, w[idx] = c,ne[idx] = h[a], h[a] =idx ++ ; 
// }
ll mod=9901;
unordered_map<ll,ll>prime;

ll fast_power(ll a,ll b)//快速幂
{
  ll res=1;
  while(b)
  {
    if(b&1)res=res*a%mod;
    b >>= 1;
    a=a*a%mod;
  }
  return res;
}

void get(ll x)//获得质因数
{
  for(ll i=2;i<=x/i;i++)
  {
    while(x%i==0)
    {
      x/=i;
      prime[i]++;
    }
  }
  if(x>1)prime[x]++;
}

ll sum(ll p,ll k)//sum函数
{
  if(k==1)return 1;
  if(k%2==0)
  {
    return (1+fast_power(p,k/2))*sum(p,k/2)%mod;
  }else
  {
    return (fast_power(p, k - 1) + sum(p, k - 1)) % mod;
  }
}


void solve()
{
  cin >> n >> m;
  get(n);
  ll ans=1;
  for(auto it : prime)
  {
    ll a = it.first, b = it.second * m;
    ans=ans*sum(a,b+1)%mod;
  }
  if(!n)ans=0;
  cout << ans;
}



int main()
{
   IOS;
   ll _;
    _=1;
    //scanf("%lld",&_);
   // cin>>_;
  
    ca=1;
    while(_--)
    {
      solve(); 
      ca++;
    }    
    return 0;
}