-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathG.cpp
120 lines (115 loc) · 3.26 KB
/
G.cpp
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
#include<bits/stdc++.h>
using namespace std;
using LL = long long;
const int maxn = 2e6 + 5;
bool isPrime[maxn];
int primes[maxn], minp[maxn], id[maxn], cnt;
int mx1[maxn], mx2[maxn], c[maxn];
void init(){
isPrime[1] = 1;
for(int i = 2; i < maxn; i++){
if (!isPrime[i]){
id[i] = cnt;
primes[cnt++] = i;
minp[i] = i;
}
for(int j = 0; i * primes[j] < maxn; j++){
isPrime[i * primes[j]] = 1;
minp[i * primes[j]] = primes[j];
if (i % primes[j] == 0) break;
}
}
}
int main(){
#ifdef LOCAL
freopen("data.in", "r", stdin);
freopen("data.out", "w", stdout);
#endif
cin.tie(0);
cout.tie(0);
ios::sync_with_stdio(0);
const int mod = 1e9 + 7;
init();
int n;
cin >> n;
for(int i = 0; i < n; i++){
int x;
cin >> x;
c[id[x]] += 1;
}
int ans = 1;
int add = 0;
for(int i = cnt - 1; i >= 0; i--){
if (mx1[i] != 0){
for(int j = 0; j < mx1[i]; j++){
ans = 1LL * ans * primes[i] % mod;
}
if (c[i] >= 1){
int t = primes[i] - 1;
while(t != 1){
int p = minp[t], c = 0;
while(t % p == 0) t /= p, c += 1;
if (c >= mx1[id[p]]){
mx2[id[p]] = mx1[id[p]];
mx1[id[p]] = c;
}
else if (c > mx2[id[p]]){
mx2[id[p]] = c;
}
}
}
}
else{
if (c[i] >= 1){
ans = 1LL * ans * primes[i] % mod;
}
if (c[i] >= 2){
int t = primes[i] - 1;
while(t != 1){
int p = minp[t], c = 0;
while(t % p == 0) t /= p, c += 1;
if (c >= mx1[id[p]]){
mx2[id[p]] = mx1[id[p]];
mx1[id[p]] = c;
}
else if (c > mx2[id[p]]){
mx2[id[p]] = c;
}
}
}
}
}
for(int i = cnt - 1; i >= 0; i--){
if (mx1[i] != 0){
if (c[i] >= 1){
bool ok = true;
int t = primes[i] - 1;
while(t != 1){
int p = minp[t], c = 0;
while(t % p == 0) t /= p, c += 1;
if (c > mx2[id[p]]){
ok = false;
}
}
if (ok) add = 1;
}
if (c[i] >= 2) add = 1;
}
else{
if (c[i] >= 2){
bool ok = true;
int t = primes[i] - 1;
while(t != 1){
int p = minp[t], c = 0;
while(t % p == 0) t /= p, c += 1;
if (c > mx2[id[p]]){
ok = false;
}
}
if (ok) add = 1;
}
if (c[i] >= 3) add = 1;
}
}
cout << ans + add << '\n';
}