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-8 of 8 results.

A035338 4th column of Wythoff array.

Original entry on oeis.org

5, 18, 26, 39, 52, 60, 73, 81, 94, 107, 115, 128, 141, 149, 162, 170, 183, 196, 204, 217, 225, 238, 251, 259, 272, 285, 293, 306, 314, 327, 340, 348, 361, 374, 382, 395, 403, 416, 429, 437, 450, 458, 471, 484, 492, 505, 518, 526, 539, 547, 560, 573, 581, 594
Offset: 0

Views

Author

Keywords

Comments

The asymptotic density of this sequence is 1/phi^5 = phi^5 - 11 = A244593 - 4 = 0.0901699... . - 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

  • Maple
    t := (1+sqrt(5))/2 ; [ seq(5*floor((n+1)*t)+3*n,n=0..80) ];
  • Mathematica
    f[n_] := 5 Floor[(n + 1) GoldenRatio] + 3n; Array[f, 54, 0] (* Robert G. Wilson v, Dec 11 2017 *)
  • Python
    from math import isqrt
    def A035338(n): return 5*(n+1+isqrt(5*(n+1)**2)>>1)+3*n # Chai Wah Wu, Aug 11 2022

A134862 Wythoff ABB numbers.

Original entry on oeis.org

8, 21, 29, 42, 55, 63, 76, 84, 97, 110, 118, 131, 144, 152, 165, 173, 186, 199, 207, 220, 228, 241, 254, 262, 275, 288, 296, 309, 317, 330, 343, 351, 364, 377, 385, 398, 406, 419, 432, 440, 453, 461, 474, 487, 495, 508, 521, 529, 542, 550, 563, 576, 584, 597
Offset: 1

Views

Author

Clark Kimberling, Nov 14 2007

Keywords

Comments

The lower and upper Wythoff sequences, A and B, satisfy the complementary equation ABB = 2A+3B.
The asymptotic density of this sequence is 1/phi^5 = phi^5 - 11 = A244593 - 4 = 0.0901699... . - 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_] := A[B[B[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 A(B(B(n))) # Indranil Ghosh, Jun 10 2017
    
  • Python
    from math import isqrt
    def A134862(n): return 5*(n+isqrt(5*n**2)>>1)+3*n # Chai Wah Wu, Aug 10 2022

Formula

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

A134863 Wythoff BAB numbers.

Original entry on oeis.org

7, 20, 28, 41, 54, 62, 75, 83, 96, 109, 117, 130, 143, 151, 164, 172, 185, 198, 206, 219, 227, 240, 253, 261, 274, 287, 295, 308, 316, 329, 342, 350, 363, 376, 384, 397, 405, 418, 431, 439, 452, 460, 473, 486, 494, 507, 520, 528, 541, 549, 562, 575, 583, 596
Offset: 1

Views

Author

Clark Kimberling, Nov 14 2007

Keywords

Comments

The lower and upper Wythoff sequences, A and B, satisfy the complementary equation BAB = 2A+3B-1.
Also numbers with suffix string 1010, when written in Zeckendorf representation. - A.H.M. Smeets, Mar 24 2024
The asymptotic density of this sequence is 1/phi^5 = phi^5 - 11 = A244593 - 4 = 0.0901699... . - 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[B[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(B(n))) # Indranil Ghosh, Jun 10 2017
    
  • Python
    from math import isqrt
    def A134863(n): return 5*(n+isqrt(5*n**2)>>1)+3*n-1 # Chai Wah Wu, Aug 11 2022

Formula

a(n) = B(A(B(n))), n>=1, with A=A000201, the lower Wythoff sequence and B=A001950, the upper Wythoff sequence.
From A.H.M. Smeets, Mar 24 2024: (Start)
a(n) = 2*A(n) + 3*B(n) - 1 (see Clark Kimberling 2008), with A=A000201, B=A001950, the lower and upper Wythoff sequences, respectively.
Equals {A035336}\{A134861} (= Wythoff BA \ Wythoff BAA). (End)

A244847 Decimal expansion of rho_c = (5-sqrt(5))/10, the asymptotic critical density for the hard hexagon model.

Original entry on oeis.org

2, 7, 6, 3, 9, 3, 2, 0, 2, 2, 5, 0, 0, 2, 1, 0, 3, 0, 3, 5, 9, 0, 8, 2, 6, 3, 3, 1, 2, 6, 8, 7, 2, 3, 7, 6, 4, 5, 5, 9, 3, 8, 1, 6, 4, 0, 3, 8, 8, 4, 7, 4, 2, 7, 5, 7, 2, 9, 1, 0, 2, 7, 5, 4, 5, 8, 9, 4, 7, 9, 0, 7, 4, 3, 6, 2, 1, 9, 5, 1, 0, 0, 5, 8, 5, 5, 8, 5, 5, 9, 1, 6, 2, 1, 2, 1, 7, 7, 2, 5, 0, 3
Offset: 0

Views

Author

Jean-François Alcover, Nov 12 2014

Keywords

Comments

The vertical distance between the accumulation point and the outermost point of a golden spiral inscribed inside a golden rectangle with dimensions phi and 1 along the x and y axes, respectively (the horizontal distance is A176015). - Amiram Eldar, May 18 2021
The limiting frequency of the digit 1 in the base phi representation of real numbers in the range [0,1], where phi is the golden ratio (A001622) (Rényi, 1957). - Amiram Eldar, Mar 18 2025

Examples

			0.2763932022500210303590826331268723764559381640388474275729102754589479...
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 1.2 The Golden Mean, phi, p. 7.
  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 5.12.1 Phase transitions in Lattice Gas Models, p. 347.

Crossrefs

Essentially the same sequence of digits as A229760 and A187799.

Programs

  • Mathematica
    RealDigits[(5 - Sqrt[5])/10, 10, 102] // First

Formula

Equals 1/(sqrt(5)*phi), where phi = (1+sqrt(5))/2 = A001622. - Vaclav Kotesovec, Nov 13 2014
Equals lim_{n -> infinity} A000045(n)/A000032(n+1). - Bruno Berselli, Jan 22 2018
Equals Sum_{n>=1} A000045(3^(n-1))/A000032(3^n) = Sum_{n>=1} A045529(n-1)/A006267(n). - Amiram Eldar, Dec 20 2018
Equals 1 - A242671. - Amiram Eldar, Mar 18 2025

A224776 Number of lattice paths from (0,0) to (n,n) that do not go below the x-axis or above the diagonal x=y and consist of steps D=(1,-1), H=(1,0) and S=(0,1).

Original entry on oeis.org

1, 1, 3, 14, 83, 568, 4271, 34296, 288946, 2524676, 22695611, 208713400, 1955285936, 18601484936, 179267898087, 1746795785272, 17183086302528, 170427862676296, 1702621483524154, 17118538010217472, 173092651634957516, 1759113081143064184, 17959329720442879275
Offset: 0

Views

Author

Alois P. Heinz, Apr 25 2013

Keywords

Examples

			a(0) = 1: the empty path.
a(1) = 1: HS.
a(2) = 3: HSHS, HHSS, HSDSS.
a(3) = 14: HSHSHS, HHSSHS, HSDSSHS, HSHHSS, HHSHSS, HSDSHSS, HHHSSS, HSDHSSS, HSHDSSS, HHSDSSS, HSDSDSSS, HSHSDSS, HHSSDSS, HSDSSDSS.
		

Crossrefs

Cf. A000108 (without D-steps), A114296 (paths to (n,0)), A225042 (with additional U-steps), A244593, A286764.

Programs

  • Maple
    b:= proc(x, y) option remember; `if`(y>x, 0, `if`(x=0, 1,
           b(x-1, y)+`if`(y>0, b(x, y-1), 0)+b(x-1, y+1)))
        end:
    a:= n-> b(n, n):
    seq(a(n), n=0..25);
  • Mathematica
    b[x_, y_] := b[x, y] = If[y > x, 0, If[x == 0, 1, b[x - 1, y] + If[y > 0, b[x, y - 1], 0] + b[x - 1, y + 1]]];
    a[n_] := b[n, n];
    a /@ Range[0, 25] (* Jean-François Alcover, Dec 18 2020, after Alois P. Heinz *)

Formula

a(n) ~ c * ((11+5*sqrt(5))/2)^n / n^(3/2), where c = 0.01403940208697420741365874329992235342402687... . - Vaclav Kotesovec, Sep 07 2014

A108554 Main diagonal of square array A108553, in which row n equals the crystal ball sequence for D_n lattice.

Original entry on oeis.org

1, 3, 13, 147, 1681, 18733, 204205, 2195399, 23429185, 249258777, 2650436845, 28209301211, 300745906449, 3212758605189, 34392151363021, 368908132170511, 3964622452867329, 42682164804034609, 460242445806795853
Offset: 0

Views

Author

Paul D. Hanna, Jun 10 2005

Keywords

Comments

Limit a(n+1)/a(n) ~ 11.08...
Limit is equal to A001622^5 = A244593 = 11.090169943749474241... - Vaclav Kotesovec, Aug 31 2025

Crossrefs

Programs

  • Mathematica
    Join[{1, 3}, Table[Sum[Binomial[2*n - j, n - j]*(Binomial[2*n, 2*j] - 2*j*(n - j)*Binomial[n, j]/(n - 1)), {j, 0, n}], {n, 2, 20}]] (* Vaclav Kotesovec, Aug 31 2025 *)
  • PARI
    {a(n)=if(n<0,0,if(n==0,1,if(n==1,3,sum(j=0,n,binomial(2*n-j,n-j)* (binomial(2*n,2*j)-2*j*(n-j)*binomial(n,j)/(n-1))))))}

Formula

a(n) ~ phi^(5*n + 1) / (2*5^(1/4)*sqrt(Pi*n)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Aug 31 2025

A374149 Decimal expansion of the minimum volume of an axis-aligned bounding box which includes the shortest minimum-link polygonal chain joining all the vertices of the cube {0,1}^3.

Original entry on oeis.org

5, 5, 4, 5, 0, 8, 4, 9, 7, 1, 8, 7, 4, 7, 3, 7, 1, 2, 0, 5, 1, 1, 4, 6, 7, 0, 8, 5, 9, 1, 4, 0, 9, 5, 2, 9, 4, 3, 0, 0, 7, 7, 2, 9, 4, 9, 5, 1, 4, 4, 0, 7, 1, 5, 5, 3, 3, 8, 6, 2, 1, 5, 5, 6, 7, 6, 3, 1, 5, 1, 1, 5, 7, 0, 4, 7, 2, 5, 6, 1, 2, 4, 2, 6, 8, 0, 1
Offset: 1

Views

Author

Marco Ripà, Jun 29 2024

Keywords

Comments

It has been proved that it is not possible to join the 8 vertices of a cube with a polygonal chain that has fewer than 6 edges (see Links, Optimal cycles enclosing all the nodes of a k-dimensional hypercube, Theorem 2.2).
Here we are considering the additional constraint that asks to minimize the volume of the Axis-Aligned Bounding Box (AABB) including the above-mentioned optimal polygonal chain consisting of only 6 connected line segments and that joins all the vertices of the cube [0,1] X [0,1] X [0,1].
Given phi = (1+sqrt(5))/2, the well-known golden ratio (see A001622), a valid polygonal chain is (0, 1, 0)-(0, 0, 0)-((1+phi)/2, 0, (1+phi)/2)-(1/2, 1+phi, 1/2)-((1-phi)/2, 0, (1+phi)/2)-(1, 0, 0)-(1, 1, 0) (see Links, p. 164), and consequently the minimum volume AABB is [(1-phi)/2, (1+phi)/2] X [0, 1+phi] X [0, (1+phi)/2].
As noted by Hugo Pfoertner, the present sequence is also given by phi^5/2 (i.e., A244593/2).

Examples

			5.5450849718747371205114670859140952943...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[(11+5*Sqrt[5])/4, 10, 100][[1]]

Formula

Equals phi*(1+phi)*((1+phi)/2), where phi := (1+sqrt(5))/2 is the golden ratio.
Equals (11+5*sqrt(5))/4.
Equals phi^5/2.
Equals 10*A134944 + 3/2.

A351898 Decimal expansion of metallic ratio for N = 14.

Original entry on oeis.org

1, 4, 0, 7, 1, 0, 6, 7, 8, 1, 1, 8, 6, 5, 4, 7, 5, 2, 4, 4, 0, 0, 8, 4, 4, 3, 6, 2, 1, 0, 4, 8, 4, 9, 0, 3, 9, 2, 8, 4, 8, 3, 5, 9, 3, 7, 6, 8, 8, 4, 7, 4, 0, 3, 6, 5, 8, 8, 3, 3, 9, 8, 6, 8, 9, 9, 5, 3, 6, 6, 2, 3, 9, 2, 3, 1, 0, 5, 3, 5, 1, 9, 4, 2, 5, 1, 9
Offset: 2

Views

Author

A.H.M. Smeets, Feb 24 2022

Keywords

Comments

Decimal expansion of continued fraction [14; 14, 14, 14, ...].
Also largest solution of x^2 - 14 x - 1 = 0.
Essentially the same digit sequence as A010503, A157214, A174968 and A268683.
The metallic ratio's for N = A077444(n) are equal to powers of the silver ratio, i.e., A014166^(2n-1); this constant represents the special case for N = A077444(2).

Examples

			14.0710678118654752440084436210484903928483593...
		

Crossrefs

Metallic ratios: A001622 (N=1), A014176 (N=2), A098316 (N=3), A098317 (N=4), A098318 (N=5), A176398 (N=6), A176439 (N=7), A176458 (N=8), A176522 (N=9), A176537 (N=10), A244593 (N=11).

Programs

  • Mathematica
    RealDigits[7 + 5*Sqrt[2], 10, 100][[1]] (* Amiram Eldar, Feb 24 2022 *)
  • PARI
    (1+sqrt(2))^3

Formula

Equals 2 + 5*A014176.
Equals A014176^3.
Equals exp(arcsinh(7)). - Amiram Eldar, Jul 04 2023
Showing 1-8 of 8 results.