-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path10179.cpp
66 lines (57 loc) · 1.35 KB
/
10179.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
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
vector<int> prime({2, 3, 5, 7});
int seen = 9;
void explore(int x)
{
int sq = sqrt(seen);
int pw = sq * sq;
while (seen <= x) {
if (pw == seen) {
sq++;
pw = sq * sq;
}
else {
if (pw < seen) {
sq++;
pw = sq * sq;
}
bool pass = true;
for (int i = 1; pass && prime[i] < sq; i++) {
pass = seen % prime[i] != 0;
}
if (pass) prime.push_back(seen);
}
seen += 2;
}
}
int main(void)
{
int n, cnt, sq, frac;
explore(31623);
while (1) {
scanf("%d", &n);
if (!n) break;
sq = sqrt(n);
cnt = n;
frac = n;
for (int i = 0; i < int(prime.size()) && prime[i] <= sq; ++i)
if (n % prime[i] == 0) {
cnt = (cnt / prime[i]) * (prime[i] - 1);
while (frac > 1 && frac % prime[i] == 0)
frac /= prime[i];
}
if (frac > 1) cnt = cnt / frac * (frac - 1);
if (n == 1)
printf("1\n");
else if (cnt == n)
printf("%d\n", cnt - 1);
else
printf("%d\n", cnt);
}
return 0;
}