约数个数
给定 n 个正整数 ai,请你输出这些数的乘积的约数个数,答案对 109+7 取模。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一个整数 ai。
输出格式
输出一个整数,表示所给正整数的乘积的约数个数,答案需对 109+7 取模。
数据范围
1≤n≤100,
1≤ai≤2×109
输入样例:
3
2
6
8
输出样例:
12
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 110, mod = 1e9 + 7;
int n;
int main(){
cin >> n;
unordered_map<int, int> primes;
while(n -- ){
int x;
cin >> x;
for(int i = 2; i <= x / i; i ++ ){
while(x % i == 0){
x /= i;
primes[i] ++ ;
}
}
if(x > 1) primes[x] ++ ; // 最后 x > 1 是是特判大于 sqrt(x) 的那个唯一的质因子
}
LL res = 1;
for(auto p : primes)
res = res * (p.second + 1) % mod;
cout << res << endl;
return 0;
}
约数之和(基础版)
给定 n 个正整数 ai,请你输出这些数的乘积的约数之和,答案对 109+7 取模。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一个整数 ai。
输出格式
输出一个整数,表示所给正整数的乘积的约数之和,答案需对 109+7 取模。
数据范围
1≤n≤100,
1≤ai≤2×109
输入样例:
3
2
6
8
输出样例:
252
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 110, mod = 1e9 + 7;
int main(){
int n;
cin >> n;
unordered_map<int, int> primes;
while(n -- ){
int x;
cin >> x;
for(int i = 2; i <= x / i; i ++ ){
while(x % i == 0){
x /= i;
primes[i] ++ ;
}
}
if(x > 1) primes[x] ++ ;
}
LL res = 1;
for(auto p : primes){
LL a = p.first, b = p.second;
LL t = 1;
while(b -- ) t = (t * a + 1) % mod;
res = res * t % mod;
}
cout << res << endl;
return 0;
}
约数之和(提高版)
假设现在有两个自然数 A 和 B,S 是 A^B 的所有约数之和。
请你求出 Smod9901 的值是多少。
输入格式
在一行中输入用空格隔开的两个整数 A 和 B。
输出格式
输出一个整数,代表 Smod9901 的值。
数据范围
0≤A,B≤5×107
输入样例:
2 3
输出样例:
15
注意: A 和 B 不会同时为 0。
再次预习一下哈哈哈哈
#include<bits/stdc++.h>
#define x first
#define y second
using namespace std;
typedef long long LL;
const int mod = 9901;
int a, b;
unordered_map<int, int> primes;
void divide(int n){
for(int i = 2; i <= n / i; i ++ ){
if(n % i == 0){
while(n % i == 0){
primes[i] ++ ;
n /= i;
}
}
}
if(n > 1) primes[n] ++ ;
}
int qmi(int a, int b){
int res = 1;
while(b){
if(b & 1) res = (LL) res * a % mod;
b >>= 1;
a = (LL) a * a % mod;
}
return res;
}
int sum(int p, int k){
if(k == 1) return 1;
if(k % 2 == 0){
return (LL) (qmi(p, k / 2) + 1) * sum(p, k / 2) % mod;
}
return (LL) (qmi(p, k - 1) + sum(p, k - 1)) % mod;
}
int main(){
cin >> a >>b;
divide(a);
int res = 1;
for(auto it : primes){
int p = it.x, k = it.y * b;
res = (LL) res * sum(p, k + 1) % mod;
}
if(!a) res = 0;
cout << res << endl;
return 0;
}