OFFSET
-1,1
COMMENTS
If f(0)=A201942, f(n) will generate the odd numbers.
From Jon Maiga, May 02 2019: (Start)
The image shows f(n), note the symmetry around -LambertW(1) (~ -0.57). For example, all values in the lowest positive horizontal 'segment' (f(n) ~ 0.2) consists of the lesser twin primes (A001359). The corresponding segment at f(n)~-1.2 are n such that n+1 and n+3 are primes (A144834).
With 256-digit precision for f(0), f(n) will correctly decode all n up to 1069, which is the 180th prime number.
(End)
LINKS
FORMULA
Equals lim_{n->oo} sf(0)*exp(sf(1)*exp(...sf(n-1)*exp(sf(n)))) where sf(n)=2*(A010051(n)-1).
EXAMPLE
-0.02172268983561208427430192882333538885414410223300617362296942014938449977851...
MATHEMATICA
f0[n_] := (v = 0; For[i = n, i >= 0, i--, v = If[PrimeQ[i], 1, -1]*Exp[v]]; v);
RealDigits[f0[1000], 10, 100][[1]]
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Jon Maiga, Jul 30 2018
STATUS
approved