OFFSET
1,1
COMMENTS
The constant c is given in the article [Plouffe, 2018] with 2600 digits of precision.
LINKS
Simon Plouffe, A set of formulas for primes, arXiv:1901.01849 [math.NT], 2019.
Simon Plouffe, The calculation of p(n) and pi(n), arXiv:2002.12137 [math.NT], 2020.
FORMULA
a(n) = round(c^((5/4)^n)), where c is a real constant starting 43.80468771580293481859664562569089495081037087137495184074061328752670419506...
EXAMPLE
a(1) = round(c^((5/4)^1)) = round(112.69...) = 113,
a(2) = round(c^((5/4)^2)) = round(367.17...) = 367,
a(3) = round(c^((5/4)^3)) = round(1607.2...) = 1607, etc..
MAPLE
# Computes the values according to the formula, v = 43.804..., e = 5/4, m the
# number of terms. Returns the real and the rounded values (primes).
val := proc(s, e, m)
local ll, v, n, kk;
v := s;
ll := [];
for n to m do
v := v^e; ll := [op(ll), v]
end do;
return [ll, map(round, ll)]
end;
CROSSREFS
KEYWORD
nonn
AUTHOR
Simon Plouffe, Jan 05 2019
STATUS
approved