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-10 of 11 results. Next

A332438 Decimal expansion of (2*cos(Pi/9))^2 = A332437^2.

Original entry on oeis.org

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

Views

Author

Wolfdieter Lang, Mar 31 2020

Keywords

Comments

This algebraic number rho(9)^2 of degree 3 is a root of its minimal polynomial x^3 - 6*x^2 + 9*x - 1.
The other two roots are x2 = (2*cos(5*Pi/9))^2 = (2*cos(4*Pi/9))^2 = (R(4,rho(9)))^2 = 2 - rho(9) = 0.120614758..., and x3 = (2*cos(7*Pi/9))^2 = (2*cos(7*Pi/9))^2 = (R(7,rho(9)))^2 = 4 + rho(9) - rho(9)^2 = 2.347296355... = A130880 + 2, with rho(9) = 2*cos(Pi/9) = A332437, the monic Chebyshev polynomials R (see A127672), and the computation is done modulo the minimal polynomial of rho(9) which is x^3 - 3*x - 1 (see A187360).
This gives the representation of these roots in the power basis of the simple field extension Q(rho(9)). See the linked W. Lang paper in A187360, sect. 4.
This number rho(9)^2 appears as limit of the quotient of consecutive numbers af various sequences, e.g., A094256 and A094829.
The algebraic number rho(9)^2 - 2 = 1.532088898... of Q(rho(9)) has minimal polynomial x^3 - 3*x + 1 over Q. The other roots are -rho(9) = -A332437 and 2 + rho(9) - rho(9)^2 = A130880. - Wolfdieter Lang, Sep 20 2022

Examples

			3.5320888862379560704047853011108333478716649...
		

Crossrefs

2 + 2*cos(2*Pi/n): A104457 (n = 5), A116425 (n = 7), A296184 (n = 10), A019973 (n = 12).

Programs

  • Mathematica
    RealDigits[(2*Cos[Pi/9])^2, 10, 100][[1]] (* Amiram Eldar, Mar 31 2020 *)
  • PARI
    (2*cos(Pi/9))^2 \\ Michel Marcus, Sep 23 2022

Formula

Equals (2*cos(Pi/9))^2 = rho(9)^2 = A332437^2.
Equals 2 + i^(4/9) - i^(14/9). - Peter Luschny, Apr 04 2020
Equals 2 + w1^(1/3) + w2^(1/3), where w1 = (-1 + sqrt(3)*i)/2 = exp(2*Pi*i/3) and w2 = (-1 - sqrt(3)*i)/2 are the complex roots of x^3 - 1. - Wolfdieter Lang, Sep 20 2022
Constant c = 2 + 2*cos(2*Pi/9). The linear fractional transformation z -> c - c/z has order 9, that is, z = c - c/(c - c/(c - c/(c - c/(c - c/(c - c/(c - c/(c - c/(c - c/(z))))))))). - Peter Bala, May 09 2024
From Amiram Eldar, Nov 22 2024: (Start)
Equals 3 + sec(Pi/9)/2 = 3 + 1/(2*A019879).
Equals 3 + Product_{k>=3} (1 + (-1)^k/A063289(k)). (End)

A094256 Expansion of x / ( (x-1)*(x^3 - 9*x^2 + 6*x - 1) ).

Original entry on oeis.org

1, 7, 34, 143, 560, 2108, 7752, 28101, 100947, 360526, 1282735, 4552624, 16131656, 57099056, 201962057, 714012495, 2523515514, 8916942687, 31504028992, 111295205284, 393151913464, 1388758662221, 4905479957435, 17327203698086, 61202661233823, 216176614077600
Offset: 1

Views

Author

Gary W. Adamson, Apr 25 2004

Keywords

Comments

Previous name was: Let M = the 4 X 4 matrix [0 1 0 0 / 0 0 1 0 / 0 0 0 1 / -1 10 -15 7]. Perform M^n * [1 0 0 0] = [p q r s]. Then a(n-3), a(n-2), a(n-1), a(n) = -p, -q, -r, -s respectively.
a(n)/a(n-1) tends to 3.53208888624... = 4*cos^2(Pi/9), which is an eigenvalue of the matrix and a root of the polynomial x^4 - 6x^3 + 15x^2 -10x + 1 = 0 (having roots 4*cos^2(r*Pi/9), with r = 1,2,3,4).
Number of (s(0), s(1), ..., s(2n+4)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+4, s(0) = 1, s(2n+4) = 7. - Herbert Kociemba, Jun 13 2004
From Wolfdieter Lang, Mar 27 2020: (Start)
This sequence, with offset -5, starting with -85, -10, -1, 0, 0, 0, 1, 7, ... appears in the formula for the n-th power of the 4 X 4 tridiagonal matrix given in A332602 as M_4 = matrix([1,1,0,0], [1,2,1,0], [0,1,2,1], [0,0,1,2]): (M_4)^n = a(n-2)*(M_4)^3 + b(n)*(M_4)^2 + c(n)*M_4 - a(n-3)*1_4, for n >= 0, with the 4 X 4 unit Matrix 1_4, b(n) = -15*a(n-3) + 10*a(n-4) - a(n-5), and c(n) = 10*a(n-3) - a(n-4). Proof from the characteristc polynomial of M_4 (see a comment in A332602) and the Cayley-Hamilton theorem.
From the proof that A094829(n+3)/A094829(n+2) -> rho(9)^2 = A332438 for n-> infinitiy, with rho(9) = 2*cos(Pi/9) = A332437 (see a comment in A094829), and a formula given below the same limit is obtained for a(n+1)/a(n) for n -> infinity, as stated in a comment above. (End)

Examples

			a(2), a(3), a(4), a(5) = 7, 34, 143, 560, since M^5 * [1 0 0 0] = [ -7 -34 -143 -560].
Cayley-Hamilton: (M_4)^5 = a(3)*(M_4)^3 + b(5)*(M_4)^2 + c(5)*M_4 - a(2)*1_4 = 34*(M_4)^3 - 95*(M_4)^2 + 69*M_4 - 7*1_4. - _Wolfdieter Lang_, Mar 27 2020
		

References

  • C. V. Durell and A. Robson, "Advanced Trigonometry", Dover 2003, p. 216.

Crossrefs

a(n) = A005023(n-1), n > 1. - R. J. Mathar, Sep 05 2008

Programs

  • Magma
    I:=[1,7,34,143]; [n le 4 select I[n] else 7*Self(n-1) - 15*Self(n-2) + 10*Self(n-3) - Self(n-4): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015
    
  • Mathematica
    Table[ (MatrixPower[{{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {-1, 10, -15, 7}}, n].{-1, 0, 0, 0})[[4]], {n, 24}] (* Robert G. Wilson v, Apr 28 2004 *)
    LinearRecurrence[{7, -15, 10, -1}, {1, 7, 34, 143}, 40] (* Vincenzo Librandi, Jul 25 2015 *)
  • PARI
    Vec(x / ( (x-1)*(x^3-9*x^2+6*x-1) ) + O(x^30)) \\ Michel Marcus, Jul 25 2015

Formula

From Herbert Kociemba, Jun 13 2004: (Start)
a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/9)*sin(7*r*Pi/9)*(2*cos(r*Pi/9))^(2n+4).
a(n) = 7*a(n-1) - 15*a(n-2) + 10*a(n-3) - a(n-4).
G.f.: x / ( (x-1)*(x^3 - 9*x^2 + 6*x - 1) ). (End)
3*a(n) = 1 - A094829(n+2) + 8*A094829(n+1) - A094829(n). - R. J. Mathar, Jun 29 2012 [offset corrected, and A094829(1) = 0. - Wolfdieter Lang, Mar 27 2020]
a(n) = (1/3)*(1 + 2*A094829(n+1) + 8*A094829(n) - A094829(n-1)), for n >= 1, with A094829(1) and A094829(0) = 0. - Wolfdieter Lang, Mar 27 2020

Extensions

More terms from Robert G. Wilson v, Apr 28 2004
a(25)-a(26) from Vincenzo Librandi, Jul 25 2015
New name (using g.f. from Herbert Kociemba) from Joerg Arndt, Jul 25 2015

A080938 Number of Catalan paths (nonnegative, starting and ending at 0, step +-1) of 2*n steps with all values less than or equal to 7.

Original entry on oeis.org

1, 1, 2, 5, 14, 42, 132, 429, 1429, 4846, 16645, 57686, 201158, 704420, 2473785, 8704089, 30664890, 108126325, 381478030, 1346396146, 4753200932, 16783118309, 59266297613, 209302921830, 739203970773, 2610763825782, 9221050139566, 32568630376132
Offset: 0

Views

Author

Henry Bottomley, Feb 25 2003

Keywords

Comments

From Wolfdieter Lang, Mar 27 2020: (Start)
a(n) also gives the upper left entry of the n-th power of the 4 X 4 tridiagonal matrix M_4, given in A332602: M_4 = Matrix([1,1,0,0], [1,2,1,0], [0,1,2,1], [0,0,1,2]): a(n) = (M_4)^n[1,1]. Proof from the formula for (M_4)^n, given in a comment in A094256, derived from the Cayley-Hamilton theorem, which leads to the recurrence. The formula for a(n) in terms of A094256 is given below.
For A094256(n+1)/A094256(n), like for A094829(n+1)/A094829(n), the limit for n -> infinity is rho(9)^2 = A332438 = 3.53208888..., with rho(9) = 2*cos(Pi/9) = A332437. Therefore the formula of a(n) in terms of A094256 shows that the same limit is reached for a(n+1)/a(n). See this conjecture by Gary W. Adamson in A332602.
(End)

Examples

			1 + x + 2*x^2 + 5*x^3 + 14*x^4 + 42*x^5 + 132*x^6 + 429*x^7 + ...
		

Crossrefs

Cf. A000007, A000012, A011782, A001519, A007051, A080937, A024175, A033191 which essentially provide the same sequence for different limits and tend to A000108.
Cf. A211216, A094826 (first differences), A094829, A094829, A332602, A332437, A332438.

Programs

  • Magma
    I:=[1,1,2,5]; [n le 4 select I[n] else 7*Self(n-1)-15*Self(n-2)+10*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Nov 30 2018
  • Mathematica
    CoefficientList[Series[(1 - 2 x) (2 x^2 - 4 x + 1) / ((x - 1) (x^3 - 9 x^2 + 6 x - 1)), {x, 0, 40}], x] (* Vincenzo Librandi, Nov 30 2018 *)
    LinearRecurrence[{7, -15, 10, -1}, {1, 1, 2, 5}, 30] (* Jean-François Alcover, Jan 07 2019 *)
  • PARI
    {a(n) = local(A); A = 1; for( i=1, 7, A = 1 / (1 - x*A)); polcoeff( A + x * O(x^n), n)} /* Michael Somos, May 12 2012 */
    

Formula

a(n) = A080934(n,7).
G.f.: -(2*x - 1)*(2*x^2 - 4*x + 1) / ( (x - 1)*(x^3 - 9*x^2 + 6*x - 1) ). - Ralf Stephan, May 13 2003
a(n) = 7*a(n-1) - 15*a(n-2) + 10*a(n-3) - a(n-4). - Herbert Kociemba, Jun 13 2004
G.f.: 1 / (1 - x / (1 - x / (1 - x / (1 - x / (1 - x / (1 - x / (1 - x))))))). - Michael Somos, May 12 2012
a(n) = 5*b(n-2) - 21*b(n-3) + 19*b(n-4) - 2*b(n-5), for n >= 0, with b(n) = A094256(n), for n >= -5. See a comment in A094256 for this offset, and the above comment. - Wolfdieter Lang, Mar 28 2020

A094833 Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 3, s(2n) = 5.

Original entry on oeis.org

1, 4, 15, 55, 199, 714, 2548, 9061, 32148, 113887, 403051, 1425471, 5039254, 17809336, 62928201, 222324436, 785402143, 2774421135, 9800231959, 34617003682, 122274355596, 431893332397, 1525507797700, 5388281150223
Offset: 1

Views

Author

Herbert Kociemba, Jun 13 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2n) counts (s(0), s(1), ..., s(2n)) such that 0 < s(i) < m and |s(i)-s(i-1)| = 1 for i = 1,2,...,2n, s(0) = j, s(2n) = k.

Programs

  • Mathematica
    Rest@ CoefficientList[Series[(-x + 2 x^2)/(-1 + 6 x - 9 x^2 + x^3), {x, 0, 24}], x] (* Michael De Vlieger, Jul 02 2021 *)

Formula

a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/3)*sin(5*r*Pi/9)*(2*cos(r*Pi/9))^(2n).
a(n) = 6a(n-1) - 9a(n-2) + a(n-3).
G.f.: (-x+2x^2)/(-1 + 6x - 9x^2 + x^3).
a(n+1) = 3*a(n) + A094832(n-1). - Philippe Deléham, Mar 20 2007
a(n) = A094829(n+1) - 2*A094829(n). - R. J. Mathar, Nov 14 2019

A005023 Number of walks of length 2n+7 in the path graph P_8 from one end to the other.

Original entry on oeis.org

7, 34, 143, 560, 2108, 7752, 28101, 100947, 360526, 1282735, 4552624, 16131656, 57099056, 201962057, 714012495, 2523515514, 8916942687, 31504028992, 111295205284, 393151913464, 1388758662221, 4905479957435, 17327203698086, 61202661233823, 216176614077600
Offset: 1

Views

Author

Keywords

References

  • W. Feller, An Introduction to Probability Theory and its Applications, 3rd ed, Wiley, New York, 1968, p. 96.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A094829 (first differences), A094256 (essentially the same).

Programs

  • Magma
    I:=[7, 34, 143, 560]; [n le 4 select I[n] else 7*Self(n-1)-15*Self(n-2)+10*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Jun 08 2013
  • Maple
    a:=k->sum(binomial(7+2*k,9*j+k-2),j=ceil((2-k)/9)..floor((9+k)/9))-sum(binomial(7+2*k,9*j+k-1),j=ceil((1-k)/9)..floor((8+k)/9)): seq(a(k),k=1..28);
    A005023:=-(-7+15*z-10*z**2+z**3)/(z-1)/(z**3-9*z**2+6*z-1); # Conjectured by Simon Plouffe in his 1992 dissertation.
  • Mathematica
    CoefficientList[Series[(-z^3 + 10 z^2 - 15 z + 7)/(z^4 - 10 z^3 + 15 z^2 - 7 z + 1), {z, 0, 100}], z] (* Vladimir Joseph Stephan Orlovsky, Jun 27 2011 *)
    LinearRecurrence[{7,-15,10,-1},{7,34,143,560},40] (* Harvey P. Dale, May 26 2013 *)
    CoefficientList[Series[(1 / x) (1 / (1 - 7 x + 15 x^2 - 10 x^3 + x^4) - 1), {x, 0, 50}], x] (* Vincenzo Librandi, Jun 08 2013 *)

Formula

G.f.: 1/(1-7x+15x^2-10x^3+x^4) - 1. a(n)=7a(n-1)-15a(n-2)+10a(n-3)-a(n-4). - Emeric Deutsch, Apr 02 2004
a(k) = sum(binomial(7+2k, 9j+k-2)-binomial(7+2k, 9j+k-1), j=-infinity..infinity) (a finite sum).

Extensions

Better definition from Emeric Deutsch, Apr 02 2004

A094826 Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 1, s(2n) = 3.

Original entry on oeis.org

1, 3, 9, 28, 90, 297, 1000, 3417, 11799, 41041, 143472, 503262, 1769365, 6230304, 21960801, 77461435, 273351705, 964918116, 3406804786, 12029917377, 42483179304, 150036624217, 529901048943, 1871559855009, 6610286313784
Offset: 1

Views

Author

Herbert Kociemba, Jun 13 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2n) counts (s(0), s(1), ..., s(2n)) such that 0 < s(i) < m and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = j, s(2n) = k.

Programs

  • Mathematica
    Rest@ CoefficientList[Series[x (-1 + 3 x)/(-1 + 6 x - 9 x^2 + x^3), {x, 0, 25}], x] (* Michael De Vlieger, Aug 05 2021 *)
    LinearRecurrence[{6,-9,1},{1,3,9},30] (* Harvey P. Dale, Dec 29 2021 *)

Formula

a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/9)*sin(3*r*Pi/9)*(2*cos(r*Pi/9))^(2n).
a(n) = 6*a(n-1) - 9*a(n-2) + a(n-3) = 7*a(n-1) - 15*a(n-2) + 10*a(n-3) - a(n-4).
G.f.: x(-1+3x)/(-1+6x-9x^2+x^3).
a(n) = A094829(n+1) - 3*A094829(n). - R. J. Mathar, Nov 14 2019

A094827 Number of (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = 1, s(2n+1) = 4.

Original entry on oeis.org

1, 4, 14, 48, 165, 571, 1988, 6953, 24396, 85786, 302104, 1064945, 3756519, 13256712, 46796545, 165225380, 583440086, 2060408640, 7276716445, 25700060995, 90770326604, 320598127113, 1132355884236, 3999522488002
Offset: 1

Views

Author

Herbert Kociemba, Jun 13 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2n+1) counts (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < m and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = j, s(2n+1) = k.

Programs

  • Mathematica
    LinearRecurrence[{7,-15,10,-1},{1,4,14,48},30] (* Harvey P. Dale, Jul 09 2020 *)

Formula

a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/9)*sin(4*r*Pi/9)*(2*cos(r*Pi/9))^(2*n+1).
a(n) = 7*a(n-1) - 15*a(n-2) + 10*a(n-3) - a(n-4).
G.f.: x*(1-3*x+x^2) / ( (x-1)*(x^3-9*x^2+6*x-1) ).
3*a(n) = A094829(n+2) -2*A094829(n+1) -2*A094829(n)-1. - R. J. Mathar, Nov 14 2019

A094828 Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 1, s(2n) = 5.

Original entry on oeis.org

1, 5, 20, 75, 274, 988, 3536, 12597, 44745, 158632, 561683, 1987154, 7026408, 24835744, 87763945, 310088381, 1095490524, 3869911659, 13670143618, 48287147300, 170561502896, 602454835293, 2127962632993, 7516243783216
Offset: 2

Views

Author

Herbert Kociemba, Jun 13 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2*n) counts (s(0), s(1), ..., s(2n)) such that 0 < s(i) < m and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = j, s(2n) = k.

Programs

  • Mathematica
    LinearRecurrence[{7,-15,10,-1},{1,5,20,75},30] (* Harvey P. Dale, Apr 27 2020 *)

Formula

a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/9)*sin(5*r*Pi/9)*(2*cos(r*Pi/9))^(2*n).
a(n) = 7*a(n-1) - 15*a(n-2) + 10*a(n-3) - a(n-4).
G.f.: -x^2*(-1+2*x) / ( (x-1)*(x^3-9*x^2+6*x-1) ).
a(n) = A094256(n-1) - 2*A094256(n-2). - R. J. Mathar, Nov 14 2019
3*a(n) = A094829(n+2) -5*A094829(n+1)+7*A094829(n)-1. - R. J. Mathar, Nov 14 2019

A094834 Number of (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < 9 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = 3, s(2n+1) = 6.

Original entry on oeis.org

1, 5, 21, 82, 308, 1131, 4096, 14705, 52497, 186733, 662630, 2347680, 8309143, 29388368, 103895601, 367187437, 1297452581, 4583924154, 16193659132, 57204089987, 202065531888, 713750040577, 2521114546457, 8905002445437
Offset: 1

Views

Author

Herbert Kociemba, Jun 13 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2n+1) counts (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < m and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = j, s(2n+1) = k.
This sequence is the odd bisection of A188048. - John Blythe Dobson, Jun 20 2015

Programs

  • Magma
    I:=[1,5,21]; [n le 3 select I[n] else 6*Self(n-1)-9*Self(n-2)+Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jun 21 2015
  • Mathematica
    CoefficientList[Series[(x - 1)/(- 1 + 6 x - 9 x^2 + x^3), {x, 0, 20}], x] (* Vincenzo Librandi, Jun 21 2015 *)
    LinearRecurrence[{6,-9,1},{1,5,21},30] (* Harvey P. Dale, Dec 27 2019 *)

Formula

a(n) = (2/9)*Sum_{r=1..8} sin(r*Pi/3)*sin(2*r*Pi/3)*(2*cos(r*Pi/9))^(2n+1).
a(n) = 6*a(n-1) - 9*a(n-2) + a(n-3).
G.f.: x(-1+x)/(-1 + 6x - 9x^2 + x^3).
a(n) = A094829(n+1) - A094829(n). - R. J. Mathar, Nov 15 2019

A217315 Square array T, read by antidiagonals: T(n,k) = 0 if n-k >= 1 or if k-n >= 8, T(0,k)= 1 if 0<=k<=7, T(n,k) = T(n-1,k) + T(n,k-1).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 0, 0, 1, 3, 2, 0, 0, 1, 4, 5, 0, 0, 0, 1, 5, 9, 5, 0, 0, 0, 1, 6, 14, 14, 0, 0, 0, 0, 0, 7, 20, 28, 14, 0, 0, 0, 0, 0, 7, 27, 48, 42, 0, 0, 0, 0, 0, 0, 0, 34, 75, 90, 42, 0, 0, 0, 0, 0, 0, 0, 34, 109, 165, 132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 143, 274, 297, 132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 143, 417, 571, 429, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Mar 17 2013

Keywords

Comments

A hexagon arithmetic of E. Lucas.

Examples

			Square array begins:
1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, ... row n=0
0, 1, 2, 3, 4, 5, 6, 7, 7, 0, 0, 0, 0, 0, 0, ... row n=1
0, 0, 2, 5, 9, 14, 20, 27, 34, 34, 0, 0, 0, ... row n=2
0, 0, 0, 5, 14, 28, 48, 75, 109, 143, 143, 0, 0, ... row n=3
0, 0, 0, 0, 14, 42, 90, 165, 274, 417, 560, 560, 0, ... row n=4
0, 0, 0, 0, 0, 42, 132, 297, 571, 988, 1548, 2108, 2108, 0, ... row n=5
...
		

Crossrefs

Cf. Similar sequence: A216230, A216228, A216226, A216238, A216054, A217257.

Programs

  • Mathematica
    t[0, k_ /; k <= 7] = 1; t[n_, k_] /; k < n || k > n+7 = 0; t[n_, k_] := t[n, k] = t[n-1, k] + t[n, k-1]; Table[t[n-k, k], {n, 0, 13}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Mar 18 2013 *)

Formula

T(n,n) = A080938(n).
T(n,n+1) = A080938(n+1).
T(n,n+2) = A094826(n+1).
T(n,n+3) = A094827(n+1).
T(n,n+4) = A094828(n+2).
T(n,n+5) = A094829(n+2).
T(n,n+6) = T(n,n+7) = A094256(n+1).
Sum_{k, 0<=k<=n} T(n-k,k) = A061551(n).
Showing 1-10 of 11 results. Next