OFFSET
0,2
COMMENTS
A related identity is Sum_{n>=0} (p + q^n)^n * r^n/n! = Sum_{n>=0} exp(p*q^n*r) * q^(n^2) * r^n/n!. Here, p = 2*Pi*i, q = 1/2, and r = 1.
FORMULA
Equals Im(Sum_{n>=0} (2*Pi*i + 1/2^n)^n / n!) (see A393771 for the real part).
Equals Sum_{n>=0} sin(2*Pi/2^n) / (n! * 2^(n^2)).
Equals Sum_{n>=0} sqrt((1 - cos(4*Pi/2^n))/2) / (n! * 2^(n^2)).
EXAMPLE
y = 0.03148042134069452309162958222422342516963974616558...
Constant y = 0 + 0/2 + 1/(2!*2^4) + sqrt(1/2)/(3!*2^9) + sqrt((1-sqrt(1/2))/2)/(4!*2^16) + sqrt((1-sqrt((1+sqrt(1/2))/2))/2)/(5!*2^25) + ... + sin(2*Pi/2^n)/(n!*2^(n^2)) + ...
equals the imaginary part of: Sum_{n>=0} (2*Pi*i + 1/2^n)^n/n! = (0.5002307656... + i*0.03148042134...).
PROG
(PARI) \\ y = imaginary( Sum_{n>=0} (2*Pi*i + 1/2^n)^n/n! ).
\p200 \\ set desired precision
{y = imag( suminf(n=0, (2*Pi*I + 1/2^n)^n/n! ) )}
for(n=1, 120, print1(floor(y*10^n)%10, ", "))
(PARI) \\ y = Sum_{n>=0} sin(2*Pi/2^n) / (n! * 2^(n^2)).
\p200 \\ set desired precision
{y = suminf(n=0, sin(2*Pi/2^n)/(n!*2^(n^2)) )}
for(n=1, 120, print1(floor(y*10^n)%10, ", "))
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Paul D. Hanna, Mar 12 2026
STATUS
approved
