OFFSET
0,13
COMMENTS
This sequence is the imaginary part of {f(n)} defined as:
- f(0) = 0,
- f(n+1) = f(n) + i^t(n)
where t(n) is the number of 1's and 7's minus the number of 3's and 5's
plus twice the number of 4's in the base 9 representation of n
and i denotes the imaginary unit.
LINKS
FORMULA
a(9^k) = 0 for any k >= 0.
PROG
(PARI) { [R, U, L, D]=[0..3]; p = [R, U, R, D, L, D, R, U, R]; z=0; for (n=0, 76, print1 (imag(z) ", "); z += I^vecsum(apply(d -> p[1+d], digits(n, #p)))) }
CROSSREFS
KEYWORD
sign,base
AUTHOR
Rémy Sigrist, Feb 10 2020
STATUS
approved