约数之和 (简单+普通)

目录

简单

输入格式

输出格式

数据范围

输入样例:

输出样例:

普通

输入格式

输出格式

数据范围

输入样例:

输出样例:


简单

给定 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
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include 
//#include
#include
#include
#include
#include
#define dbug cout<<"*****hear*****"<=c;a--)
#define no cout<<"NO"<,greater >q;
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair PII;
typedef pair 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_mapmp;
  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<>_;
  
    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
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include 
//#include
#include
#include
#include
#include
#define dbug cout<<"*****hear*****"<=c;a--)
#define no cout<<"NO"<,greater >q;
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair PII;
typedef pair 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_mapprime;

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;
}

你可能感兴趣的:(算法整理,算法,数据结构,1024程序员节,c++)