cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A332246 a(n) is the X-coordinate of the n-th point of the Minkowski sausage (or Minkowski curve). Sequence A332247 gives Y-coordinates.

Original entry on oeis.org

0, 1, 1, 2, 2, 2, 3, 3, 4, 4, 3, 3, 4, 5, 5, 4, 4, 5, 5, 6, 6, 6, 7, 7, 8, 8, 9, 9, 8, 7, 7, 8, 8, 8, 9, 9, 8, 7, 7, 8, 8, 9, 9, 10, 10, 10, 11, 11, 12, 12, 11, 11, 12, 13, 13, 12, 12, 13, 13, 14, 14, 14, 15, 15, 16, 16, 15, 15, 16, 17, 17, 16, 16, 15, 15, 14
Offset: 0

Views

Author

Rémy Sigrist, Feb 08 2020

Keywords

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 6's minus the number of 3's and 4's
in the base 8 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 |
.--->.

Crossrefs

See A163528, A323258 and A332204 for similar sequences.
Cf. A332247 (Y-coordinates).

Programs

  • PARI
    { dd = [0,1,0,-1,-1,0,1,0]; z=0; for (n=0, 75, print1 (real(z)", "); z += I^vecsum(apply(d -> dd[1+d], digits(n, #dd)))) }

Formula

a(8^k-m) + a(m) = 4^k for any k >= 0 and m = 0..8^k.

A332205 a(n) is the imaginary part of f(n) defined by f(0) = 0, and f(n+1) = f(n) + g((1+i)^(A065359(n) mod 8)) (where g(z) = z/gcd(Re(z), Im(z)) and i denotes the imaginary unit).

Original entry on oeis.org

0, 0, 1, 0, 0, 1, 2, 2, 3, 2, 2, 1, 0, 0, 1, 0, 0, 1, 2, 2, 3, 4, 5, 6, 7, 7, 8, 7, 7, 8, 9, 9, 10, 9, 9, 8, 7, 7, 8, 7, 7, 6, 5, 4, 3, 2, 2, 1, 0, 0, 1, 0, 0, 1, 2, 2, 3, 2, 2, 1, 0, 0, 1, 0, 0, 1, 2, 2, 3, 4, 5, 6, 7, 7, 8, 7, 7, 8, 9, 9, 10, 11, 12, 13, 14
Offset: 0

Views

Author

Rémy Sigrist, Feb 07 2020

Keywords

Comments

Looks much like A005536, in particular in respect of its symmetries of scale (compare the scatterplots). - Peter Munn, Jun 21 2021

Crossrefs

Cf. A005536, A007052, A065359, A332204 (real part and additional comments), A332206 (positions of 0's, cf. A001196).

Programs

  • Mathematica
    A065359[0] = 0;
    A065359[n_] := -Total[(-1)^PositionIndex[Reverse[IntegerDigits[n, 2]]][1]];
    g[z_] := z/GCD[Re[z], Im[z]];
    Module[{n = 0}, Im[NestList[# + g[(1+I)^A065359[n++]] &, 0, 100]]] (* Paolo Xausa, Aug 28 2024 *)
  • PARI
    \\ See Links section.

Formula

a(2^(2*k-1)) = A007052(k) for any k >= 0.
a(4^k-m) = a(m) for any k >= 0 and m = 0..4^k.
Showing 1-2 of 2 results.