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.

A274773 a(n) = floor(sqrt(2*n-1) + 1/2) - abs(2*(n-1) - (floor(sqrt(2*n-1) + 1/2))^2) + 1.

Original entry on oeis.org

1, 1, 3, 1, 3, 3, 1, 3, 5, 3, 1, 3, 5, 5, 3, 1, 3, 5, 7, 5, 3, 1, 3, 5, 7, 7, 5, 3, 1, 3, 5, 7, 9, 7, 5, 3, 1, 3, 5, 7, 9, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 11, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 13, 11, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 13, 13, 11, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 13, 15, 13, 11, 9, 7, 5, 3, 1, 3, 5, 7, 9, 11, 13, 15, 15, 13, 11, 9, 7, 5, 3
Offset: 1

Views

Author

Ilya Gutkovskiy, Aug 11 2016

Keywords

Comments

First bisection of A004738.
All terms are odd.

Examples

			Triangle begins:
  1;
  1, 3;
  1, 3, 3;
  1, 3, 5, 3;
  1, 3, 5, 5, 3;
  1, 3, 5, 7, 5, 3;
  1, 3, 5, 7, 7, 5, 3;
  1, 3, 5, 7, 9, 7, 5, 3;
  1, 3, 5, 7, 9, 9, 7, 5, 3;
  ...
Read like the Ulam spiral, starting with 1:
    x  7  x  5  x  3  x  1
    7  x  5  x  3  x  1  x
    x  5  x  3  x  1  x  3
    5  x  3  x  1  x  3  x
    x  3  x  1  x  3  x  5
    3  x  1  x  3  x  5  x
    x  1  x  3  x  5  x  7
    1  x  3  x  5  x  7  x
		

Crossrefs

Programs

  • Mathematica
    Table[Floor[Sqrt[2 n - 1] + 1/2] - Abs[2 (n - 1) - Floor[Sqrt[2 n - 1] + 1/2]^2] + 1, {n, 1, 120}]
  • Python
    from gmpy2 import isqrt_rem
    def A274773(n):
        i, j = isqrt_rem(2*n-1)
        return int(i+2 - abs(j-2*(i+1)) if 4*(i-j) + 1 <= 0 else i+1 - abs(j-1)) # Chai Wah Wu, Aug 15 2016