A163333 -id:A163333 - cp's OEIS frontend

A163336 Peano curve in an n X n grid, starting downwards from the top left corner, listed antidiagonally as A(0,0), A(0,1), A(1,0), A(0,2), A(1,1), A(2,0), ...

0, 5, 1, 6, 4, 2, 47, 7, 3, 15, 48, 46, 8, 14, 16, 53, 49, 45, 9, 13, 17, 54, 52, 50, 44, 10, 12, 18, 59, 55, 51, 39, 43, 11, 23, 19, 60, 58, 56, 38, 40, 42, 24, 22, 20, 425, 61, 57, 69, 37, 41, 29, 25, 21, 141, 426, 424, 62, 68, 70, 36, 30, 28, 26, 140, 142, 431, 427

Offset: 0

Antti Karttunen, Jul 29 2009

			The top left 9 X 9 corner of the array shows how this surjective self-avoiding walk begins (connect the terms in numerical order, 0-1-2-3-...):
   0  5  6 47 48 53 54 59 60
   1  4  7 46 49 52 55 58 61
   2  3  8 45 50 51 56 57 62
  15 14  9 44 39 38 69 68 63
  16 13 10 43 40 37 70 67 64
  17 12 11 42 41 36 71 66 65
  18 23 24 29 30 35 72 77 78
  19 22 25 28 31 34 73 76 79
  20 21 26 27 32 33 74 75 80

A. Karttunen, Table of n, a(n) for n = 0..3320
E. H. Moore, On Certain Crinkly Curves, Transactions of the American Mathematical Society, volume 1, number 1, 1900, pages 72-90. (And errata.) See section 7 (figure 3 with Y downwards is the table here).
Giuseppe Peano, Sur une courbe, qui remplit toute une aire plane, Mathematische Annalen, volume 36, number 1, 1890, pages 157-160. Also EUDML (link to GDZ).
Rémy Sigrist, Colored scatterplot of (x, y) such that 0 <= x, y < 3^6 (where the hue is function of T(x, y))
Eric Weisstein's World of Mathematics, Hilbert curve (this curve called "Hilbert II").
Wikipedia, Self-avoiding walk
Wikipedia, Space-filling curve
Index entries for sequences that are permutations of the natural numbers

Transpose: A163334. Inverse: A163337. a(n) = A163332(A163330(n)) = A163327(A163333(A163328(n))) = A163334(A061579(n)). One-based version: A163340. Row sums: A163342. Row 0: A163481. Column 0: A163480. Central diagonal: A163343.

See A163357 and A163359 for the Hilbert curve.

b[{n_, k_}, {m_}] := (A[n, k] = m - 1);
MapIndexed[b, List @@ PeanoCurve[4][[1]]];
Table[A[n - k, k], {n, 0, 12}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Mar 07 2021 *)

Name corrected by Kevin Ryde, Aug 28 2020

A163332 Self-inverse permutation of the integers for constructing the Peano curve in an N X N grid.

Original entry on oeis.org

0, 1, 2, 5, 4, 3, 6, 7, 8, 15, 16, 17, 14, 13, 12, 9, 10, 11, 18, 19, 20, 23, 22, 21, 24, 25, 26, 47, 46, 45, 48, 49, 50, 53, 52, 51, 44, 43, 42, 39, 40, 41, 38, 37, 36, 29, 28, 27, 30, 31, 32, 35, 34, 33, 54, 55, 56, 59, 58, 57, 60, 61, 62, 69, 70, 71, 68, 67, 66, 63, 64, 65

Offset: 0

Antti Karttunen, Jul 29 2009

nonn

The integers [0,(3^k)-1] are confined to range [0,(3^k)-1].
From Kevin Ryde, Sep 04 2020: (Start)
To calculate a(n), write n in ternary digits n[k],..,n[0], where n[0] is the least significant digit.  Then the ternary digits of a(n) are a[j] = k^{n[j+1]+n[j+3]+...}(n[j]) where Peano's complement operator k^{s}(d) = d if s even or 2-d if s odd.
A single complement is k(d) = 2-d and the "exponent" is repeats k(k(k(...))).  Sum s = n[j+1] + n[j+3] + ... is every second digit above j, so digit j flips 0 <-> 2 when an odd number of odd digits (1's) among these.  The complement does not change digit parity so a second transformation re-complements back to the original digits and so self-inverse a(a(n)) = n.
Peano's curve is formed by digits of a(n) put alternately to x and y coordinates, so a(n) maps between the Peano curve the ternary Z-order curve per the formulas in A163528, A163529.
(End)

Antti Karttunen, Table of n, a(n) for n = 0..59048
Giuseppe Peano, Sur une courbe, qui remplit toute une aire plane, Mathematische Annalen, volume 36, number 1, 1890, pages 157-160. Also EUDML (link to GDZ).
Index entries for sequences that are permutations of the natural numbers

Coordinates using this transform: A163528, A163529.
A163334 & A163336 give two variants of the Peano curve in an N X N grid.
Cf. A163355 (Hilbert curve).

a[n_] := FromDigits[With[{d = Reverse@IntegerDigits[n, 3]}, Reverse@Table[
  If[EvenQ@Total@d[[j+1 ;; ;; 2]], d[[j]], 2-d[[j]]], {j, Length@d}]], 3];
Array[a, 100] (* Andrey Zabolotskiy, Apr 08 2021, after Kevin Ryde *)

a(n) = my(v=digits(n,3)); for(start=2,3, my(s=0); forstep(i=start,#v,2, s+=v[i-1]; if(s%2,v[i]=2-v[i]))); fromdigits(v,3); \\ Kevin Ryde, Sep 04 2020

a(n) = A163327(A163333(A163327(n))).

Name corrected by Kevin Ryde, Aug 27 2020

cp's OEIS Frontend

A163336 Peano curve in an n X n grid, starting downwards from the top left corner, listed antidiagonally as A(0,0), A(0,1), A(1,0), A(0,2), A(1,1), A(2,0), ...

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Extensions