约数之和 (简单+普通)

目录

简单

输入格式

输出格式

数据范围

输入样例:

输出样例:

普通

输入格式

输出格式

数据范围

输入样例:

输出样例:


简单

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

取模。

输入格式

第一行包含整数 n

接下来 n

行,每行包含一个整数 ai

输出格式

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

取模。

数据范围

1≤n≤100

,
1≤ai≤2×1e9

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

思路:

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

        一个整数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
输出样例:
15

注意: 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;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值