-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathDivisorAnalysis.cpp
More file actions
67 lines (54 loc) · 1.34 KB
/
DivisorAnalysis.cpp
File metadata and controls
67 lines (54 loc) · 1.34 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
//Divisor Analysis - https://cses.fi/problemset/task/2182
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD = 1e9 + 7;
const int MOD2 = MOD - 1;
const ll INF = 1e18;
ll expo(ll a, ll b, ll m) {
ll res = 1;
while (b > 0) {
if (b & 1) {
res *= a;
res %= m;
}
a *= a;
a %= m;
b >>= 1;
}
return res;
}
ll inverse(ll a, ll mod) {
return expo(a , mod - 2, mod);
}
void solve() {
int n;
cin >> n;
ll num_divisors = 1;
ll sum_of_divisors = 1;
ll product_of_divisors = 1;
ll num_divisors2 = 1;
for (int i = 0; i < n; i++) {
ll p, k;
cin >> p >> k;
num_divisors *= (k + 1);
num_divisors %= MOD;
sum_of_divisors *= ((expo(p, k + 1, MOD) - 1) * inverse(p - 1, MOD) % MOD);
sum_of_divisors %= MOD;
product_of_divisors = expo(product_of_divisors, k + 1, MOD) * expo(expo(p, k*(k+1)/2, MOD), num_divisors2, MOD);
product_of_divisors %= MOD;
num_divisors2 *= (k + 1);
num_divisors2 %= MOD2;
}
cout << num_divisors << " " << sum_of_divisors << " " << product_of_divisors << endl;
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
int T = 1;
// cin >> T;
while (T--) {
solve();
}
return 0;
}