|
MATHEMATICA
|
FiveAlmostPrimePi[n_] := Sum[ PrimePi[n/(Prime@i*Prime@j*Prime@k*Prime@l)] - l + 1, {i, PrimePi[n^(1/5)]}, {j, i, PrimePi[(n/Prime@i)^(1/4)]}, {k, j, PrimePi[(n/(Prime@i*Prime@j)^(1/3))]}, {l, k, PrimePi@Sqrt[(n/(Prime@i*Prime@j*Prime@k))]}];
FiveAlmostPrime[n_] := Block[{e = Floor[Log[2, n] +4], a, b}, a = 2^e; Do[b = 2^p; While[FourAlmostPrimePi[a] < n, a = a + b]; a = a - b/2, {p, e, 0, -1}]; a + b/2]; Do[ Print@FiveAlmostPrime[10^n], {n, 0, 13}]
|