OFFSET
0,4
COMMENTS
This sequence is the real 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.
We can also build the curve by successively applying the following substitution to an initial vector (1, 0):
.--->.
^ |
| v
.--->/<---/--->.
| ^
v |
.--->.
REFERENCES
Benoit B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman and Co., 1983, section 7, "Harnessing the Peano Monster Curves", page 62 description and plate 63 bottom right drawn with chamfered corners.
LINKS
Rémy Sigrist, Table of n, a(n) for n = 0..6561
Joerg Arndt, Plane-filling curves on all uniform grids, arXiv:1607.02433 [math.CO], 2016, 2018. Curve R9-1 drawn in figure 4.1-O (top row forms, vertical mirror image).
Donald E. Knuth, Selected Papers on Fun and Games, CSLI Lecture Notes Number 192, CSLI Publications, 2010, ISBN 978-1-57586-585-0, page 611 folding product DUUUDDDU drawn at 45 degrees in a labyrinth.
Walter Wunderlich, Über Peano-Kurven, Elemente der Mathematik, volume 28, number 1, 1973, pages 1-10. See section 4 serpentine type 010 101 010 as illustrated in figure 3, the coordinates here being diagonal steps across the unit squares there.
FORMULA
a(9^k) = 3^k 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, 86, print1 (real(z) ", "); z += I^vecsum(apply(d -> p[1+d], digits(n, #p)))) }
CROSSREFS
KEYWORD
AUTHOR
Rémy Sigrist, Feb 10 2020
STATUS
approved