-
Notifications
You must be signed in to change notification settings - Fork 5
Expand file tree
/
Copy pathDistinctPowers.java
More file actions
95 lines (73 loc) · 2.68 KB
/
DistinctPowers.java
File metadata and controls
95 lines (73 loc) · 2.68 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
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
import java.io.*;
import java.util.*;
import java.lang.*;
public class Solution {
public static void main(String[] args) throws IOException{
BufferedReader f = new BufferedReader(new InputStreamReader(System.in));
long n = Long.parseLong(f.readLine());
Long output = (n - 1) * (n - 1);
ArrayList<Integer> lists = new ArrayList<>();
boolean[] checker = new boolean[200001];
Arrays.fill(checker, true);
checker[0] = false;
checker[1] = false;
for(int i = 0; i < checker.length; i++){
if(checker[i]){
lists.add(i);
for(int j = i * i; j < checker.length && j > 0; j*=i){
checker[j] = false;
}
}
}
int[] prime = new int[lists.size()];
for(int i = 0; i < prime.length; i++) prime[i] = lists.get(i);
//System.out.println(Arrays.toString(prime));
long num = 0;
long base = 1;
long largest = prime[0];
long times = 1;
while(largest*prime[0] <= n){
largest *= prime[0];
times++;
}
boolean[][] check = new boolean[(int)times + 1][(int)n + 1];
for(int k = 1; k < check.length; k++){
for(int j = 1; j < check[0].length; j++){
if(j != 0 && j != 1) check[k][j] = true;
}
}
while(base*prime[0] <= n){
base *= prime[0];
num++;
long cden = num + 1;
long cbase = base;
while(cbase*prime[0] <= n){
cbase *= prime[0];
//System.out.println(": " + cbase + " " + (int)exp);
for(int j = 2; cden*j <= num*n; j++){
if(cden * j % num == 0)check[(int)cden][j] = false;
}
cden++;
}
}
long[] sol = new long[check.length];
for(int i = 1; i < sol.length; i++){
int cur = 0;
for(int k = 1; k < check[0].length; k++){
if(k != 0 && k != 1 && !check[i][k]) sol[i]++;
}
}
for(int i = 0; prime[i] <= Math.sqrt(n) + 1; i++){
long clargest = prime[i];
int ctimes = 1;
while(clargest*prime[i] <= n){
clargest *= prime[i];
ctimes++;
}
for(int k = 1; k <= ctimes; k++){
output -= sol[k];
}
}
System.out.println(output);
}
}