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.

A134861 Wythoff BAA numbers.

Original entry on oeis.org

2, 10, 15, 23, 31, 36, 44, 49, 57, 65, 70, 78, 86, 91, 99, 104, 112, 120, 125, 133, 138, 146, 154, 159, 167, 175, 180, 188, 193, 201, 209, 214, 222, 230, 235, 243, 248, 256, 264, 269, 277, 282, 290, 298, 303, 311, 319, 324, 332, 337, 345, 353, 358, 366, 371
Offset: 1

Views

Author

Clark Kimberling, Nov 14 2007

Keywords

Comments

The lower and upper Wythoff sequences, A and B, satisfy the complementary equation BAA = A+2B-3.
Also numbers with suffix string 0010, when written in Zeckendorf representation (with leading zero's for the first term). - A.H.M. Smeets, Mar 20 2024
The asymptotic density of this sequence is 1/phi^4 = A094214^4 = 0.145898... . - Amiram Eldar, Mar 24 2025

Crossrefs

Let A = A000201, B = A001950. Then AA = A003622, AB = A003623, BA = A035336, BB = A101864. The eight triples AAA, AAB, ..., BBB are A134859, A134860, A035337, A134862, A134861, A134863, A035338, A134864, resp.

Programs

  • Mathematica
    A[n_] := Floor[n * GoldenRatio]; B[n_] := Floor[n * GoldenRatio^2]; a[n_] := B[A[A[n]]]; Array[a, 100] (* Amiram Eldar, Mar 24 2025 *)
  • Python
    from sympy import floor
    from mpmath import phi
    def A(n): return floor(n*phi)
    def B(n): return floor(n*phi**2)
    def a(n): return B(A(A(n))) # Indranil Ghosh, Jun 10 2017
    
  • Python
    from math import isqrt
    def A134861(n): return 3*((n+isqrt(5*n**2)>>1)-1)+(n<<1) # Chai Wah Wu, Aug 10 2022

Formula

a(n) = B(A(A(n))), n>=1, with A=A000201, the lower Wythoff sequence and B=A001950, the upper Wythoff sequence.