login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A318479
For any n >= 0 with binary expansion Sum_{k=0..w} b_k * 2^k, let h(n) = Sum_{k=0..w} b_k * (i-1)^k (where i denotes the imaginary unit); a(n) is the square of the modulus of h(n).
3
0, 1, 2, 1, 4, 5, 2, 1, 8, 13, 10, 13, 4, 9, 2, 5, 16, 9, 26, 17, 20, 13, 26, 17, 8, 5, 18, 13, 4, 1, 10, 5, 32, 41, 18, 25, 52, 61, 34, 41, 40, 53, 26, 37, 52, 65, 34, 45, 16, 17, 10, 9, 36, 37, 26, 25, 8, 13, 2, 5, 20, 25, 10, 13, 64, 65, 82, 81, 36, 37, 50
OFFSET
0,3
COMMENTS
See A318438 for the real part of h and additional comments.
FORMULA
a(n) = A318438(n)^2 + A318439(n)^2.
a(2^k) = 2^k for any k >= 0.
a(3 * 2^k) = 2^k for any l >= 0.
PROG
(PARI) a(n) = my (d=Vecrev(digits(n, 2))); norm(sum(i=1, #d, d[i]*(I-1)^(i-1)))
CROSSREFS
Sequence in context: A229763 A163509 A161399 * A053985 A091564 A321301
KEYWORD
nonn,look,base
AUTHOR
Rémy Sigrist, Aug 27 2018
STATUS
approved