6 lines
245 B
Text
6 lines
245 B
Text
ClearAll[Phi,pi]
|
|
$RecursionLimit = 10^6;
|
|
Phi[x_, 0] := x
|
|
Phi[x_, a_] := Phi[x, a] = Phi[x, a - 1] - Phi[Floor[x/Prime[a]], a - 1]
|
|
pi[n_] := Module[{a}, If[n < 2, 0, a = pi[Floor[Sqrt[n]]]; Phi[n, a] + a - 1]]
|
|
Scan[Print[pi[10^#]] &, Range[0,9]]
|