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.

Previous Showing 11-20 of 141 results. Next

A026935 a(n) = Sum_{k=0..n-2} T(n,k) * T(n,k+2), with T given by A008288.

Original entry on oeis.org

1, 10, 75, 500, 3149, 19214, 115031, 680424, 3992921, 23305234, 135514019, 785892316, 4549048229, 26295995926, 151857925039, 876366840784, 5055045581745, 29148894792730, 168045778127355, 968679251764676, 5583525654107645, 32183666525389086, 185514611981021959
Offset: 2

Views

Author

Keywords

Crossrefs

Cf. A008288.

Programs

Extensions

More terms from Sean A. Irvine, Oct 17 2019

A026936 a(n) = Sum_{k=0..n-3} T(n,k) * T(n,k+3), with T given by A008288.

Original entry on oeis.org

1, 14, 131, 1028, 7333, 49362, 319943, 2020552, 12526153, 76603094, 463676491, 2784550604, 16619355501, 98706690714, 583943891087, 3443573355152, 20254044891793, 118871039217822, 696401860980499, 4073713439171732, 23799679604123189, 138895001298962786, 809854583626750039
Offset: 3

Views

Author

Keywords

Crossrefs

Cf. A008288.

Programs

Extensions

More terms from Sean A. Irvine, Oct 17 2019

A204135 Array: row n shows the coefficients of the characteristic polynomial of the n-th principal submatrix of the Delannoy matrix (A008288).

Original entry on oeis.org

1, -1, 2, -4, 1, 8, -28, 17, -1, 64, -384, 424, -80, 1, 1024, -10624, 19400, -7700, 401, -1, 32768, -598016, 1748224, -1225536, 161618, -2084, 1, 2097152, -68550656, 319410176, -363159040, 95891872
Offset: 1

Views

Author

Clark Kimberling, Jan 12 2012

Keywords

Comments

Let p(n)=p(n,x) be the characteristic polynomial of the n-th principal submatrix. The zeros of p(n) are real, and they interlace the zeros of p(n+1). See A202605 and A204016 for guides to related sequences.

Examples

			Top of the array:
1....-1
2....-4.....1
8....-28....17....-1
64...-384...424...-80...1
The interlacing of zeros is illustrated by these zeros (truncated):
p(1):  1
p(2): .58, 3.41
p(3): .36, 1.44, 15.19
p(4): .21, .87, 4.53, 74.3
p(5): .12, .59, 2.14, 17.22, 380.91
		

References

  • (For references regarding interlacing roots, see A202605.)

Crossrefs

Programs

  • Mathematica
    f[i_, 1] := 1; f[1, j_] := 1;
    f[i_, j_] := f[i, j - 1] + f[i - 1, j - 1] + f[i - 1, j]
    m[n_] := Table[f[i, j], {i, 1, n}, {j, 1, n}]
    TableForm[m[8]] (* 8x8 principal submatrix *)
    Flatten[Table[f[i, n + 1 - i],
      {n, 1, 15}, {i, 1, n}]]  (*  Delannoy, A008288 *)
    p[n_] := CharacteristicPolynomial[m[n], x];
    c[n_] := CoefficientList[p[n], x]
    TableForm[Flatten[Table[p[n], {n, 1, 10}]]]
    Table[c[n], {n, 1, 8}]
    Flatten[%]                 (* 204135 *)
    TableForm[Table[c[n], {n, 1, 6}]]

A000129 Pell numbers: a(0) = 0, a(1) = 1; for n > 1, a(n) = 2*a(n-1) + a(n-2).

Original entry on oeis.org

0, 1, 2, 5, 12, 29, 70, 169, 408, 985, 2378, 5741, 13860, 33461, 80782, 195025, 470832, 1136689, 2744210, 6625109, 15994428, 38613965, 93222358, 225058681, 543339720, 1311738121, 3166815962, 7645370045, 18457556052, 44560482149, 107578520350, 259717522849
Offset: 0

Views

Author

Keywords

Comments

Sometimes also called lambda numbers.
Also denominators of continued fraction convergents to sqrt(2): 1, 3/2, 7/5, 17/12, 41/29, 99/70, 239/169, 577/408, 1393/985, 3363/2378, 8119/5741, 19601/13860, 47321/33461, 114243/80782, ... = A001333/A000129.
Number of lattice paths from (0,0) to the line x=n-1 consisting of U=(1,1), D=(1,-1) and H=(2,0) steps (i.e., left factors of Grand Schroeder paths); for example, a(3)=5, counting the paths H, UD, UU, DU and DD. - Emeric Deutsch, Oct 27 2002
a(2*n) with b(2*n) := A001333(2*n), n >= 1, give all (positive integer) solutions to Pell equation b^2 - 2*a^2 = +1 (see Emerson reference). a(2*n+1) with b(2*n+1) := A001333(2*n+1), n >= 0, give all (positive integer) solutions to Pell equation b^2 - 2*a^2 = -1.
Bisection: a(2*n+1) = T(2*n+1, sqrt(2))/sqrt(2) = A001653(n), n >= 0 and a(2*n) = 2*S(n-1,6) = 2*A001109(n), n >= 0, with T(n,x), resp. S(n,x), Chebyshev's polynomials of the first, resp. second kind. S(-1,x)=0. See A053120, resp. A049310. - Wolfdieter Lang, Jan 10 2003
Consider the mapping f(a/b) = (a + 2b)/(a + b). Taking a = b = 1 to start with and carrying out this mapping repeatedly on each new (reduced) rational number gives the following sequence 1/1, 3/2, 7/5, 17/12, 41/29, ... converging to 2^(1/2). Sequence contains the denominators. - Amarnath Murthy, Mar 22 2003
This is also the Horadam sequence (0,1,1,2). Limit_{n->oo} a(n)/a(n-1) = sqrt(2) + 1 = A014176. - Ross La Haye, Aug 18 2003
Number of 132-avoiding two-stack sortable permutations.
From Herbert Kociemba, Jun 02 2004: (Start)
For n > 0, the number of (s(0), s(1), ..., s(n)) such that 0 < s(i) < 4 and |s(i) - s(i-1)| <= 1 for i = 1,2,...,n, s(0) = 2, s(n) = 3.
Number of (s(0), s(1), ..., s(n)) such that 0 < s(i) < 4 and |s(i) - s(i-1)| <= 1 for i = 1,2,...,n, s(0) = 1, s(n) = 2. (End)
Counts walks of length n from a vertex of a triangle to another vertex to which a loop has been added. - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
Apart from initial terms, Pisot sequence P(2,5). See A008776 for definition of Pisot sequences. - David W. Wilson
Sums of antidiagonals of A038207 [Pascal's triangle squared]. - Ross La Haye, Oct 28 2004
The Pell primality test is "If N is an odd prime, then P(N)-Kronecker(2,N) is divisible by N". "Most" composite numbers fail this test, so it makes a useful pseudoprimality test. The odd composite numbers which are Pell pseudoprimes (i.e., that pass the above test) are in A099011. - Jack Brennen, Nov 13 2004
a(n) = sum of n-th row of triangle in A008288 = A094706(n) + A000079(n). - Reinhard Zumkeller, Dec 03 2004
Pell trapezoids (cf. A084158); for n > 0, A001109(n) = (a(n-1) + a(n+1))*a(n)/2; e.g., 1189 = (12+70)*29/2. - Charlie Marion, Apr 01 2006
(0!a(1), 1!a(2), 2!a(3), 3!a(4), ...) and (1,-2,-2,0,0,0,...) form a reciprocal pair under the list partition transform and associated operations described in A133314. - Tom Copeland, Oct 29 2007
Let C = (sqrt(2)+1) = 2.414213562..., then for n > 1, C^n = a(n)*(1/C) + a(n+1). Example: C^3 = 14.0710678... = 5*(0.414213562...) + 12. Let X = the 2 X 2 matrix [0, 1; 1, 2]; then X^n * [1, 0] = [a(n-1), a(n); a(n), a(n+1)]. a(n) = numerator of n-th convergent to (sqrt(2)-1) = 0.414213562... = [2, 2, 2, ...], the convergents being [1/2, 2/5, 5/12, ...]. - Gary W. Adamson, Dec 21 2007
A = sqrt(2) = 2/2 + 2/5 + 2/(5*29) + 2/(29*169) + 2/(169*985) + ...; B = ((5/2) - sqrt(2)) = 2/2 + 2/(2*12) + 2/(12*70) + 2/(70*408) + 2/(408*2378) + ...; A+B = 5/2. C = 1/2 = 2/(1*5) + 2/(2*12) + 2/(5*29) + 2/(12*70) + 2/(29*169) + ... - Gary W. Adamson, Mar 16 2008
From Clark Kimberling, Aug 27 2008: (Start)
Related convergents (numerator/denominator):
lower principal convergents: A002315/A001653
upper principal convergents: A001541/A001542
intermediate convergents: A052542/A001333
lower intermediate convergents: A005319/A001541
upper intermediate convergents: A075870/A002315
principal and intermediate convergents: A143607/A002965
lower principal and intermediate convergents: A143608/A079496
upper principal and intermediate convergents: A143609/A084068. (End)
Equals row sums of triangle A143808 starting with offset 1. - Gary W. Adamson, Sep 01 2008
Binomial transform of the sequence:= 0,1,0,2,0,4,0,8,0,16,..., powers of 2 alternating with zeros. - Philippe Deléham, Oct 28 2008
a(n) is also the sum of the n-th row of the triangle formed by starting with the top two rows of Pascal's triangle and then each next row has a 1 at both ends and the interior values are the sum of the three numbers in the triangle above that position. - Patrick Costello (pat.costello(AT)eku.edu), Dec 07 2008
Starting with offset 1 = eigensequence of triangle A135387 (an infinite lower triangular matrix with (2,2,2,...) in the main diagonal and (1,1,1,...) in the subdiagonal). - Gary W. Adamson, Dec 29 2008
Starting with offset 1 = row sums of triangle A153345. - Gary W. Adamson, Dec 24 2008
From Charlie Marion, Jan 07 2009: (Start)
In general, denominators, a(k,n) and numerators, b(k,n), of continued fraction convergents to sqrt((k+1)/k) may be found as follows:
let a(k,0) = 1, a(k,1) = 2k; for n > 0, a(k,2n) = 2*a(k,2n-1) + a(k,2n-2)
and a(k,2n+1) = (2k)*a(k,2n) + a(k,2n-1);
let b(k,0) = 1, b(k,1) = 2k+1; for n > 0, b(k,2n) = 2*b(k,2n-1) + b(k,2n-2)
and b(k,2n+1) = (2k)*b(k,2n) + b(k,2n-1).
For example, the convergents to sqrt(2/1) start 1/1, 3/2, 7/5, 17/12, 41/29.
In general, if a(k,n) and b(k,n) are the denominators and numerators, respectively, of continued fraction convergents to sqrt((k+1)/k) as defined above, then
k*a(k,2n)^2 - a(k,2n-1)*a(k,2n+1) = k = k*a(k,2n-2)*a(k,2n) - a(k,2n-1)^2 and
b(k,2n-1)*b(k,2n+1) - k*b(k,2n)^2 = k+1 = b(k,2n-1)^2 - k*b(k,2n-2)*b(k,2n);
for example, if k=1 and n=3, then a(1,n) = a(n+1) and
1*a(1,6)^2 - a(1,5)*a(1,7) = 1*169^2 - 70*408 = 1;
1*a(1,4)*a(1,6) - a(1,5)^2 = 1*29*169 - 70^2 = 1;
b(1,5)*b(1,7) - 1*b(1,6)^2 = 99*577 - 1*239^2 = 2;
b(1,5)^2 - 1*b(1,4)*b(1,6) = 99^2 - 1*41*239 = 2.
(End)
Starting with offset 1 = row sums of triangle A155002, equivalent to the statement that the Fibonacci sequence convolved with the Pell sequence prefaced with a "1": (1, 1, 2, 5, 12, 29, ...) = (1, 2, 5, 12, 29, ...). - Gary W. Adamson, Jan 18 2009
It appears that P(p) == 8^((p-1)/2) (mod p), p = prime; analogous to [Schroeder, p. 90]: Fp == 5^((p-1)/2) (mod p). Example: Given P(11) = 5741, == 8^5 (mod 11). Given P(17) = 11336689, == 8^8 (mod 17) since 17 divides (8^8 - P(17)). - Gary W. Adamson, Feb 21 2009
Equals eigensequence of triangle A154325. - Gary W. Adamson, Feb 12 2009
Another combinatorial interpretation of a(n-1) arises from a simple tiling scenario. Namely, a(n-1) gives the number of ways of tiling a 1 X n rectangle with indistinguishable 1 X 2 rectangles and 1 X 1 squares that come in two varieties, say, A and B. For example, with C representing the 1 X 2 rectangle, we obtain a(4)=12 from AAA, AAB, ABA, BAA, ABB, BAB, BBA, BBB, AC, BC, CA and CB. - Martin Griffiths, Apr 25 2009
a(n+1) = 2*a(n) + a(n-1), a(1)=1, a(2)=2 was used by Theon from Smyrna. - Sture Sjöstedt, May 29 2009
The n-th Pell number counts the perfect matchings of the edge-labeled graph C_2 x P_(n-1), or equivalently, the number of domino tilings of a 2 X (n-1) cylindrical grid. - Sarah-Marie Belcastro, Jul 04 2009
As a fraction: 1/79 = 0.0126582278481... or 1/9799 = 0.000102051229...(1/119 and 1/10199 for sequence in reverse). - Mark Dols, May 18 2010
Limit_{n->oo} (a(n)/a(n-1) - a(n-1)/a(n)) tends to 2.0. Example: a(7)/a(6) - a(6)/a(7) = 169/70 - 70/169 = 2.0000845... - Gary W. Adamson, Jul 16 2010
Numbers k such that 2*k^2 +- 1 is a square. - Vincenzo Librandi, Jul 18 2010
Starting (1, 2, 5, ...) = INVERTi transform of A006190: (1, 3, 10, 33, 109, ...). - Gary W. Adamson, Aug 06 2010
[u,v] = [a(n), a(n-1)] generates all Pythagorean triples [u^2-v^2, 2uv, u^2+v^2] whose legs differ by 1. - James R. Buddenhagen, Aug 14 2010
An elephant sequence, see A175654. For the corner squares six A[5] vectors, with decimal values between 21 and 336, lead to this sequence (without the leading 0). For the central square these vectors lead to the companion sequence A078057. - Johannes W. Meijer, Aug 15 2010
Let the 2 X 2 square matrix A=[2, 1; 1, 0] then a(n) = the (1,1) element of A^(n-1). - Carmine Suriano, Jan 14 2011
Define a t-circle to be a first-quadrant circle tangent to the x- and y-axes. Such a circle has coordinates equal to its radius. Let C(0) be the t-circle with radius 1. Then for n > 0, define C(n) to be the next larger t-circle which is tangent to C(n - 1). C(n) has radius A001333(2n) + a(2n)*sqrt(2) and each of the coordinates of its point of intersection with C(n + 1) is a(2n + 1) + (A001333(2n + 1)*sqrt(2))/2. See similar Comments for A001109 and A001653, Sep 14 2005. - Charlie Marion, Jan 18 2012
A001333 and A000129 give the diagonal numbers described by Theon from Smyrna. - Sture Sjöstedt, Oct 20 2012
Pell numbers could also be called "silver Fibonacci numbers", since, for n >= 1, F(n+1) = ceiling(phi*F(n)), if n is even and F(n+1) = floor(phi*F(n)), if n is odd, where phi is the golden ratio, while a(n+1) = ceiling(delta*a(n)), if n is even and a(n+1) = floor(delta*a(n)), if n is odd, where delta = delta_S = 1+sqrt(2) is the silver ratio. - Vladimir Shevelev, Feb 22 2013
a(n) is the number of compositions (ordered partitions) of n-1 into two sorts of 1's and one sort of 2's. Example: the a(3)=5 compositions of 3-1=2 are 1+1, 1+1', 1'+1, 1'+1', and 2. - Bob Selcoe, Jun 21 2013
Between every two consecutive squares of a 1 X n array there is a flap that can be folded over one of the two squares. Two flaps can be lowered over the same square in 2 ways, depending on which one is on top. The n-th Pell number counts the ways n-1 flaps can be lowered. For example, a sideway representation for the case n = 3 squares and 2 flaps is \\., .//, \./, ./., .\., where . is an empty square. - Jean M. Morales, Sep 18 2013
Define a(-n) to be a(n) for n odd and -a(n) for n even. Then a(n) = A005319(k)*(a(n-2k+1) - a(n-2k)) + a(n-4k) = A075870(k)*(a(n-2k+2) - a(n-2k+1)) - a(n-4k+2). - Charlie Marion, Nov 26 2013
An alternative formulation of the combinatorial tiling interpretation listed above: Except for n=0, a(n-1) is the number of ways of partial tiling a 1 X n board with 1 X 1 squares and 1 X 2 dominoes. - Matthew Lehman, Dec 25 2013
Define a(-n) to be a(n) for n odd and -a(n) for n even. Then a(n) = A077444(k)*a(n-2k+1) + a(n-4k+2). This formula generalizes the formula used to define this sequence. - Charlie Marion, Jan 30 2014
a(n-1) is the top left entry of the n-th power of any of the 3 X 3 matrices [0, 1, 1; 1, 1, 1; 0, 1, 1], [0, 1, 1; 0, 1, 1; 1, 1, 1], [0, 1, 0; 1, 1, 1; 1, 1, 1] or [0, 0, 1; 1, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
a(n+1) counts closed walks on K2 containing two loops on the other vertex. Equivalently the (1,1) entry of A^(n+1) where the adjacency matrix of digraph is A=(0,1;1,2). - David Neil McGrath, Oct 28 2014
For n >= 1, a(n) equals the number of ternary words of length n-1 avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
This is a divisibility sequence (i.e., if n|m then a(n)|a(m)). - Tom Edgar, Jan 28 2015
A strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n, m)) for all positive integers n and m. - Michael Somos, Jan 03 2017
a(n) is the number of compositions (ordered partitions) of n-1 into two kinds of parts, n and n', when the order of the 1 does not matter, or equivalently, when the order of the 1' does not matter. Example: When the order of the 1 does not matter, the a(3)=5 compositions of 3-1=2 are 1+1, 1+1'=1+1, 1'+1', 2 and 2'. (Contrast with entry from Bob Selcoe dated Jun 21 2013.) - Gregory L. Simay, Sep 07 2017
Number of weak orderings R on {1,...,n} that are weakly single-peaked w.r.t. the total ordering 1 < ... < n and for which {1,...,n} has exactly one minimal element for the weak ordering R. - J. Devillet, Sep 28 2017
Also the number of matchings in the (n-1)-centipede graph. - Eric W. Weisstein, Sep 30 2017
Let A(r,n) be the total number of ordered arrangements of an n+r tiling of r red squares and white tiles of total length n, where the individual tile lengths can range from 1 to n. A(r,0) corresponds to a tiling of r red squares only, and so A(r,0)=1. Let A_1(r,n) = Sum_{j=0..n} A(r,j) and let A_s(r,n) = Sum_{j=0..n} A_(s-1)(r,j). Then A_0(1,n) + A_2(3,n-4) + A_4(5,n-8) + ... + A_(2j) (2j+1, n-4j) = a(n) without the initial 0. - Gregory L. Simay, May 25 2018
(1, 2, 5, 12, 29, ...) is the fourth INVERT transform of (1, -2, 5, -12, 29, ...), as shown in A073133. - Gary W. Adamson, Jul 17 2019
Number of 2-compositions of n restricted to odd parts (and allowed zeros); see Hopkins & Ouvry reference. - Brian Hopkins, Aug 17 2020
Also called the 2-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence. - Michael A. Allen, Jan 23 2023
Named by Lucas (1878) after the English mathematician John Pell (1611-1685). - Amiram Eldar, Oct 02 2023
a(n) is the number of compositions of n when there are F(i) parts of size i, with i,n >= 1, F(n) the Fibonacci numbers, A000045(n) (see example below). - Enrique Navarrete, Dec 15 2023

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 12*x^4 + 29*x^5 + 70*x^6 + 169*x^7 + 408*x^8 + 985*x^9 + ...
From _Enrique Navarrete_, Dec 15 2023: (Start)
From the comment on compositions with Fibonacci number of parts, F(n), there are F(1)=1 type of 1, F(2)=1 type of 2, F(3)=2 types of 3, F(4)=3 types of 4, F(5)=5 types of 5 and F(6)=8 types of 6.
The following table gives the number of compositions of n=6 with Fibonacci number of parts:
Composition, number of such compositions, number of compositions of this type:
6,           1,     8;
5+1,         2,    10;
4+2,         2,     6;
3+3,         1,     4;
4+1+1,       3,     9;
3+2+1,       6,    12;
2+2+2,       1,     1;
3+1+1+1,     4,     8;
2+2+1+1,     6,     6;
2+1+1+1+1,   5,     5;
1+1+1+1+1+1, 1,     1;
for a total of a(6)=70 compositions of n=6. (End).
		

References

  • J. Austin and L. Schneider, Generalized Fibonacci sequences in Pythagorean triple preserving sequences, Fib. Q., 58:1 (2020), 340-350.
  • P. Bachmann, Niedere Zahlentheorie (1902, 1910), reprinted Chelsea, NY, 1968, vol. 2, p. 76.
  • A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, pp. 122-125, 1964.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 941.
  • J. M. Borwein, D. H. Bailey, and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 53.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 204.
  • John Derbyshire, Prime Obsession, Joseph Henry Press, 2004, see p. 16.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, Section 1.1.
  • Shaun Giberson and Thomas J. Osler, Extending Theon's Ladder to Any Square Root, Problem 3858, Elementa, No. 4 1996.
  • R. P. Grimaldi, Ternary strings with no consecutive 0's and no consecutive 1's, Congressus Numerantium, 205 (2011), 129-149.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.5 The Fibonacci and Related Sequences, p. 288.
  • Thomas Koshy, Pell and Pell-Lucas Numbers with Applications, Springer, New York, 2014.
  • Serge Lang, Introduction to Diophantine Approximations, Addison-Wesley, New York, 1966.
  • Paulo Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 43.
  • Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 3.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 46, 61.
  • J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 224.
  • Manfred R. Schroeder, "Number Theory in Science and Communication", 5th ed., Springer-Verlag, 2009, p. 90.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 34.
  • D. B. West, Combinatorial Mathematics, Cambridge, 2021, p. 62.

Crossrefs

Partial sums of A001333.
2nd row of A172236.
a(n) = A054456(n-1, 0), n>=1 (first column of triangle).
Cf. A175181 (Pisano periods), A214028 (Entry points), A214027 (number of zeros in a fundamental period).
A077985 is a signed version.
INVERT transform of Fibonacci numbers (A000045).
Cf. A038207.
The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.
Cf. A048739.
Cf. A073133.
Cf. A041085.
Sequences with g.f. 1/(1-k*x-x^2) or x/(1-k*x-x^2): A000045 (k=1), this sequence (k=2), A006190 (k=3), A001076 (k=4), A052918 (k=5), A005668 (k=6), A054413 (k=7), A041025 (k=8), A099371 (k=9), A041041 (k=10), A049666 (k=11), A041061 (k=12), A140455 (k=13), A041085 (k=14), A154597 (k=15), A041113 (k=16), A178765 (k=17), A041145 (k=18), A243399 (k=19), A041181 (k=20).

Programs

  • GAP
    a := [0,1];; for n in [3..10^3] do a[n] := 2 * a[n-1] + a[n-2]; od; A000129 := a; # Muniru A Asiru, Oct 16 2017
    
  • Haskell
    a000129 n = a000129_list !! n
    a000129_list = 0 : 1 : zipWith (+) a000129_list (map (2 *) $ tail a000129_list)
    -- Reinhard Zumkeller, Jan 05 2012, Feb 05 2011
    
  • Magma
    [0] cat [n le 2 select n else 2*Self(n-1) + Self(n-2): n in [1..35]]; // Vincenzo Librandi, Aug 08 2015
    
  • Maple
    A000129 := proc(n) option remember; if n <=1 then n; else 2*procname(n-1)+procname(n-2); fi; end;
    a:= n-> (<<2|1>, <1|0>>^n)[1, 2]: seq(a(n), n=0..40); # Alois P. Heinz, Aug 01 2008
    A000129 := n -> `if`(n<2, n, 2^(n-1)*hypergeom([1-n/2, (1-n)/2], [1-n], -1)):
    seq(simplify(A000129(n)), n=0..31); # Peter Luschny, Dec 17 2015
  • Mathematica
    CoefficientList[Series[x/(1 - 2*x - x^2), {x, 0, 60}], x] (* Stefan Steinerberger, Apr 08 2006 *)
    Expand[Table[((1 + Sqrt[2])^n - (1 - Sqrt[2])^n)/(2Sqrt[2]), {n, 0, 30}]] (* Artur Jasinski, Dec 10 2006 *)
    LinearRecurrence[{2, 1}, {0, 1}, 60] (* Harvey P. Dale, Jan 04 2012 *)
    a[ n_] := With[ {s = Sqrt@2}, ((1 + s)^n - (1 - s)^n) / (2 s)] // Simplify; (* Michael Somos, Jun 01 2013 *)
    Table[Fibonacci[n, 2], {n, 0, 20}] (* Vladimir Reshetnikov, May 08 2016 *)
    Fibonacci[Range[0, 20], 2] (* Eric W. Weisstein, Sep 30 2017 *)
    a[ n_] := ChebyshevU[n - 1, I] / I^(n - 1); (* Michael Somos, Oct 30 2021 *)
  • Maxima
    a[0]:0$
    a[1]:1$
    a[n]:=2*a[n-1]+a[n-2]$
    A000129(n):=a[n]$
    makelist(A000129(n),n,0,30); /* Martin Ettl, Nov 03 2012 */
    
  • Maxima
    makelist((%i)^(n-1)*ultraspherical(n-1,1,-%i),n,0,24),expand; /* Emanuele Munarini, Mar 07 2018 */
    
  • PARI
    for (n=0, 4000, a=contfracpnqn(vector(n, i, 1+(i>1)))[2, 1]; if (a > 10^(10^3 - 6), break); write("b000129.txt", n, " ", a)); \\ Harry J. Smith, Jun 12 2009
    
  • PARI
    {a(n) = imag( (1 + quadgen( 8))^n )}; /* Michael Somos, Jun 01 2013 */
    
  • PARI
    {a(n) = if( n<0, -(-1)^n, 1) * contfracpnqn( vector( abs(n), i, 1 + (i>1))) [2, 1]}; /* Michael Somos, Jun 01 2013 */
    
  • PARI
    a(n)=([2, 1; 1, 0]^n)[2,1] \\ Charles R Greathouse IV, Mar 04 2014
    
  • PARI
    {a(n) = polchebyshev(n-1, 2, I) / I^(n-1)}; /* Michael Somos, Oct 30 2021 */
    
  • Python
    from itertools import islice
    def A000129_gen(): # generator of terms
        a, b = 0, 1
        yield from [a,b]
        while True:
            a, b = b, a+2*b
            yield b
    A000129_list = list(islice(A000129_gen(),20)) # Chai Wah Wu, Jan 11 2022
  • Sage
    [lucas_number1(n, 2, -1) for n in range(30)]  # Zerinvary Lajos, Apr 22 2009
    

Formula

G.f.: x/(1 - 2*x - x^2). - Simon Plouffe in his 1992 dissertation.
a(2n+1)=A001653(n). a(2n)=A001542(n). - Ira Gessel, Sep 27 2002
G.f.: Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (2*k + x)/(1 + 2*k*x) ) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (x + 1 + k)/(1 + k*x) ) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (x + 3 - k)/(1 - k*x) ) may all be proved using telescoping series. - Peter Bala, Jan 04 2015
a(n) = 2*a(n-1) + a(n-2), a(0)=0, a(1)=1.
a(n) = ((1 + sqrt(2))^n - (1 - sqrt(2))^n)/(2*sqrt(2)).
For initial values a(0) and a(1), a(n) = ((a(0)*sqrt(2)+a(1)-a(0))*(1+sqrt(2))^n + (a(0)*sqrt(2)-a(1)+a(0))*(1-sqrt(2))^n)/(2*sqrt(2)). - Shahreer Al Hossain, Aug 18 2019
a(n) = integer nearest a(n-1)/(sqrt(2) - 1), where a(0) = 1. - Clark Kimberling
a(n) = Sum_{i, j, k >= 0: i+j+2k = n} (i+j+k)!/(i!*j!*k!).
a(n)^2 + a(n+1)^2 = a(2n+1) (1999 Putnam examination).
a(2n) = 2*a(n)*A001333(n). - John McNamara, Oct 30 2002
a(n) = ((-i)^(n-1))*S(n-1, 2*i), with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. See A049310. S(-1, x)=0, S(-2, x)= -1.
Binomial transform of expansion of sinh(sqrt(2)x)/sqrt(2). E.g.f.: exp(x)sinh(sqrt(2)x)/sqrt(2). - Paul Barry, May 09 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2k+1)*2^k. - Paul Barry, May 13 2003
a(n-2) + a(n) = (1 + sqrt(2))^(n-1) + (1 - sqrt(2))^(n-1) = A002203(n-1). (A002203(n))^2 - 8(a(n))^2 = 4(-1)^n. - Gary W. Adamson, Jun 15 2003
Unreduced g.f.: x(1+x)/(1 - x - 3x^2 - x^3); a(n) = a(n-1) + 3*a(n-2) + a(n-2). - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
a(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k, k)*2^(n-2k). - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
Apart from initial terms, inverse binomial transform of A052955. - Paul Barry, May 23 2004
a(n)^2 + a(n+2k+1)^2 = A001653(k)*A001653(n+k); e.g., 5^2 + 70^2 = 5*985. - Charlie Marion Aug 03 2005
a(n+1) = Sum_{k=0..n} binomial((n+k)/2, (n-k)/2)*(1+(-1)^(n-k))*2^k/2. - Paul Barry, Aug 28 2005
a(n) = a(n-1) + A001333(n-1) = A001333(n) - a(n-1) = A001109(n)/A001333(n) = sqrt(A001110(n)/A001333(n)^2) = ceiling(sqrt(A001108(n)/2)). - Henry Bottomley, Apr 18 2000
a(n) = F(n, 2), the n-th Fibonacci polynomial evaluated at x=2. - T. D. Noe, Jan 19 2006
Define c(2n) = -A001108(n), c(2n+1) = -A001108(n+1) and d(2n) = d(2n+1) = A001652(n); then ((-1)^n)*(c(n) + d(n)) = a(n). [Proof given by Max Alekseyev.] - Creighton Dement, Jul 21 2005
a(r+s) = a(r)*a(s+1) + a(r-1)*a(s). - Lekraj Beedassy, Sep 03 2006
a(n) = (b(n+1) + b(n-1))/n where {b(n)} is the sequence A006645. - Sergio Falcon, Nov 22 2006
From Miklos Kristof, Mar 19 2007: (Start)
Let F(n) = a(n) = Pell numbers, L(n) = A002203 = companion Pell numbers (A002203):
For a >= b and odd b, F(a+b) + F(a-b) = L(a)*F(b).
For a >= b and even b, F(a+b) + F(a-b) = F(a)*L(b).
For a >= b and odd b, F(a+b) - F(a-b) = F(a)*L(b).
For a >= b and even b, F(a+b) - F(a-b) = L(a)*F(b).
F(n+m) + (-1)^m*F(n-m) = F(n)*L(m).
F(n+m) - (-1)^m*F(n-m) = L(n)*F(m).
F(n+m+k) + (-1)^k*F(n+m-k) + (-1)^m*(F(n-m+k) + (-1)^k*F(n-m-k)) = F(n)*L(m)*L(k).
F(n+m+k) - (-1)^k*F(n+m-k) + (-1)^m*(F(n-m+k) - (-1)^k*F(n-m-k)) = L(n)*L(m)*F(k).
F(n+m+k) + (-1)^k*F(n+m-k) - (-1)^m*(F(n-m+k) + (-1)^k*F(n-m-k)) = L(n)*F(m)*L(k).
F(n+m+k) - (-1)^k*F(n+m-k) - (-1)^m*(F(n-m+k) - (-1)^k*F(n-m-k)) = 8*F(n)*F(m)*F(k). (End)
a(n+1)*a(n) = 2*Sum_{k=0..n} a(k)^2 (a similar relation holds for A001333). - Creighton Dement, Aug 28 2007
a(n+1) = Sum_{k=0..n} binomial(n+1,2k+1) * 2^k = Sum_{k=0..n} A034867(n,k) * 2^k = (1/n!) * Sum_{k=0..n} A131980(n,k) * 2^k. - Tom Copeland, Nov 30 2007
Equals row sums of unsigned triangle A133156. - Gary W. Adamson, Apr 21 2008
a(n) (n >= 3) is the determinant of the (n-1) X (n-1) tridiagonal matrix with diagonal entries 2, superdiagonal entries 1 and subdiagonal entries -1. - Emeric Deutsch, Aug 29 2008
a(n) = A000045(n) + Sum_{k=1..n-1} A000045(k)*a(n-k). - Roger L. Bagula and Gary W. Adamson, Sep 07 2008
From Hieronymus Fischer, Jan 02 2009: (Start)
fract((1+sqrt(2))^n) = (1/2)*(1 + (-1)^n) - (-1)^n*(1+sqrt(2))^(-n) = (1/2)*(1 + (-1)^n) - (1-sqrt(2))^n.
See A001622 for a general formula concerning the fractional parts of powers of numbers x > 1, which satisfy x - x^(-1) = floor(x).
a(n) = round((1+sqrt(2))^n/(2*sqrt(2))) for n > 0. (End) [last formula corrected by Josh Inman, Mar 05 2024]
a(n) = ((4+sqrt(18))*(1+sqrt(2))^n + (4-sqrt(18))*(1-sqrt(2))^n)/4 offset 0. - Al Hakanson (hawkuu(AT)gmail.com), Aug 08 2009
If p[i] = Fibonacci(i) and if A is the Hessenberg matrix of order n defined by A[i,j] = p[j-i+1] when i<=j, A[i,j]=-1 when i=j+1, and A[i,j]=0 otherwise, then, for n >= 1, a(n) = det A. - Milan Janjic, May 08 2010
a(n) = 3*a(n-1) - a(n-2) - a(n-3), n > 2. - Gary Detlefs, Sep 09 2010
From Charlie Marion, Apr 13 2011: (Start)
a(n) = 2*(a(2k-1) + a(2k))*a(n-2k) - a(n-4k).
a(n) = 2*(a(2k) + a(2k+1))*a(n-2k-1) + a(n-4k-2). (End)
G.f.: x/(1 - 2*x - x^2) = sqrt(2)*G(0)/4; G(k) = ((-1)^k) - 1/(((sqrt(2) + 1)^(2*k)) - x*((sqrt(2) + 1)^(2*k))/(x + ((sqrt(2) - 1)^(2*k + 1))/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 02 2011
In general, for n > k, a(n) = a(k+1)*a(n-k) + a(k)*a(n-k-1). See definition of Pell numbers and the formula for Sep 04 2008. - Charlie Marion, Jan 17 2012
Sum{n>=1} (-1)^(n-1)/(a(n)*a(n+1)) = sqrt(2) - 1. - Vladimir Shevelev, Feb 22 2013
From Vladimir Shevelev, Feb 24 2013: (Start)
(1) Expression a(n+1) via a(n): a(n+1) = a(n) + sqrt(2*a^2(n) + (-1)^n);
(2) a(n+1)^2 - a(n)*a(n+2) = (-1)^n;
(3) Sum_{k=1..n} (-1)^(k-1)/(a(k)*a(k+1)) = a(n)/a(n+1);
(4) a(n)/a(n+1) = sqrt(2) - 1 + r(n), where |r(n)| < 1/(a(n+1)*a(n+2)). (End)
a(-n) = -(-1)^n * a(n). - Michael Somos, Jun 01 2013
G.f.: G(0)/(2+2*x) - 1/(1+x), where G(k) = 1 + 1/(1 - x*(2*k-1)/(x*(2*k+1) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Aug 10 2013
G.f.: Q(0)*x/2, where Q(k) = 1 + 1/(1 - x*(4*k+2 + x)/( x*(4*k+4 + x) + 1/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 30 2013
a(n) = Sum_{r=0..n-1} Sum_{k=0..n-r-1} binomial(r+k,k)*binomial(k,n-k-r-1). - Peter Luschny, Nov 16 2013
a(n) = Sum_{k=1,3,5,...<=n} C(n,k)*2^((k-1)/2). - Vladimir Shevelev, Feb 06 2014
a(2n) = 2*a(n)*(a(n-1) + a(n)). - John Blythe Dobson, Mar 08 2014
a(k*n) = a(k)*a(k*n-k+1) + a(k-1)*a(k*n-k). - Charlie Marion, Mar 27 2014
a(k*n) = 2*a(k)*(a(k*n-k)+a(k*n-k-1)) + (-1)^k*a(k*n-2k). - Charlie Marion, Mar 30 2014
a(n+1) = (1+sqrt(2))*a(n) + (1-sqrt(2))^n. - Art DuPre, Apr 04 2014
a(n+1) = (1-sqrt(2))*a(n) + (1+sqrt(2))^n. - Art DuPre, Apr 04 2014
a(n) = F(n) + Sum_{k=1..n} F(k)*a(n-k), n >= 0 where F(n) the Fibonacci numbers A000045. - Ralf Stephan, May 23 2014
a(n) = round(sqrt(a(2n) + a(2n-1)))/2. - Richard R. Forberg, Jun 22 2014
a(n) = Product_{k divides n} A008555(k). - Tom Edgar, Jan 28 2015
a(n+k)^2 - A002203(k)*a(n)*a(n+k) + (-1)^k*a(n)^2 = (-1)^n*a(k)^2. - Alexander Samokrutov, Aug 06 2015
a(n) = 2^(n-1)*hypergeom([1-n/2, (1-n)/2], [1-n], -1) for n >= 2. - Peter Luschny, Dec 17 2015
a(n+1) = Sum_{k=0..n} binomial(n,k)*2^floor(k/2). - Tony Foster III, May 07 2017
a(n) = exp((i*Pi*n)/2)*sinh(n*arccosh(-i))/sqrt(2). - Peter Luschny, Mar 07 2018
From Rogério Serôdio, Mar 30 2018: (Start)
Some properties:
(1) a(n)^2 - a(n-2)^2 = 2*a(n-1)*(a(n) + a(n-2)) (see A005319);
(2) a(n-k)*a(n+k) = a(n)^2 + (-1)^(n+k+1)*a(k)^2;
(3) Sum_{k=2..n+1} a(k)*a(k-1) = a(n+1)^2 if n is odd, else a(n+1)^2 - 1 if n is even;
(4) a(n) - a(n-2*k+1) = (A077444(k) - 1)*a(n-2*k+1) + a(n-4*k+2);
(5) Sum_{k=n..n+9} a(k) = 41*A001333(n+5). (End)
From Kai Wang, Dec 30 2019: (Start)
a(m+r)*a(n+s) - a(m+s)*a(n+r) = -(-1)^(n+s)*a(m-n)*a(r-s).
a(m+r)*a(n+s) + a(m+s)*a(n+r) = (2*A002203(m+n+r+s) - (-1)^(n+s)*A002203(m-n)*A002203(r-s))/8.
A002203(m+r)*A002203(n+s) - A002203(m+s)*A002203(n+r) = (-1)^(n+s)*8*a(m-n)*a(r-s).
A002203(m+r)*A002203(n+s) - 8*a(m+s)*a(n+r) = (-1)^(n+s)*A002203(m-n)*A002203(r-s).
A002203(m+r)*A002203(n+s) + 8*a(m+s)*a(n+r) = 2*A002203(m+n+r+s)+ (-1)^(n+s)*8*a(m-n)*a(r-s). (End)
From Kai Wang, Jan 12 2020: (Start)
a(n)^2 - a(n+1)*a(n-1) = (-1)^(n-1).
a(n)^2 - a(n+r)*a(n-r) = (-1)^(n-r)*a(r)^2.
a(m)*a(n+1) - a(m+1)*a(n) = (-1)^n*a(m-n).
a(m-n) = (-1)^n (a(m)*A002203(n) - A002203(m)*a(n))/2.
a(m+n) = (a(m)*A002203(n) + A002203(m)*a(n))/2.
A002203(n)^2 - A002203(n+r)*A002203(n-r) = (-1)^(n-r-1)*8*a(r)^2.
A002203(m)*A002203(n+1) - A002203(m+1)*A002203(n) = (-1)^(n-1)*8*a(m-n).
A002203(m-n) = (-1)^(n)*(A002203(m)*A002203(n) - 8*a(m)*a(n) )/2.
A002203(m+n) = (A002203(m)*A002203(n) + 8*a(m)*a(n) )/2. (End)
From Kai Wang, Mar 03 2020: (Start)
Sum_{m>=1} arctan(2/a(2*m+1)) = arctan(1/2).
Sum_{m>=2} arctan(2/a(2*m+1)) = arctan(1/12).
In general, for n > 0,
Sum_{m>=n} arctan(2/a(2*m+1)) = arctan(1/a(2*n)). (End)
a(n) = (A001333(n+3*k) + (-1)^(k-1)*A001333(n-3*k)) / (20*A041085(k-1)) for any k>=1. - Paul Curtz, Jun 23 2021
Sum_{i=0..n} a(i)*J(n-i) = (a(n+1) + a(n) - J(n+2))/2 for J(n) = A001045(n). - Greg Dresden, Jan 05 2022
From Peter Bala, Aug 20 2022: (Start)
Sum_{n >= 1} 1/(a(2*n) + 1/a(2*n)) = 1/2.
Sum_{n >= 1} 1/(a(2*n+1) - 1/a(2*n+1)) = 1/4. Both series telescope - see A075870 and A005319.
Product_{n >= 1} ( 1 + 2/a(2*n) ) = 1 + sqrt(2).
Product_{n >= 2} ( 1 - 2/a(2*n) ) = (1/3)*(1 + sqrt(2)). (End)
G.f. = 1/(1 - Sum_{k>=1} Fibonacci(k)*x^k). - Enrique Navarrete, Dec 17 2023
Sum_{n >=1} 1/a(n) = 1.84220304982752858079237158327980838... - R. J. Mathar, Feb 05 2024
a(n) = ((3^(n+1) + 1)^(n-1) mod (9^(n+1) - 2)) mod (3^(n+1) - 1). - Joseph M. Shunia, Jun 06 2024

A000073 Tribonacci numbers: a(n) = a(n-1) + a(n-2) + a(n-3) for n >= 3 with a(0) = a(1) = 0 and a(2) = 1.

Original entry on oeis.org

0, 0, 1, 1, 2, 4, 7, 13, 24, 44, 81, 149, 274, 504, 927, 1705, 3136, 5768, 10609, 19513, 35890, 66012, 121415, 223317, 410744, 755476, 1389537, 2555757, 4700770, 8646064, 15902591, 29249425, 53798080, 98950096, 181997601, 334745777, 615693474, 1132436852
Offset: 0

Views

Author

Keywords

Comments

The name "tribonacci number" is less well-defined than "Fibonacci number". The sequence A000073 (which begins 0, 0, 1) is probably the most important version, but the name has also been applied to A000213, A001590, and A081172. - N. J. A. Sloane, Jul 25 2024
Also (for n > 1) number of ordered trees with n+1 edges and having all leaves at level three. Example: a(4)=2 because we have two ordered trees with 5 edges and having all leaves at level three: (i) one edge emanating from the root, at the end of which two paths of length two are hanging and (ii) one path of length two emanating from the root, at the end of which three edges are hanging. - Emeric Deutsch, Jan 03 2004
a(n) is the number of compositions of n-2 with no part greater than 3. Example: a(5)=4 because we have 1+1+1 = 1+2 = 2+1 = 3. - Emeric Deutsch, Mar 10 2004
Let A denote the 3 X 3 matrix [0,0,1;1,1,1;0,1,0]. a(n) corresponds to both the (1,2) and (3,1) entries in A^n. - Paul Barry, Oct 15 2004
Number of permutations satisfying -k <= p(i)-i <= r, i=1..n-2, with k=1, r=2. - Vladimir Baltic, Jan 17 2005
Number of binary sequences of length n-3 that have no three consecutive 0's. Example: a(7)=13 because among the 16 binary sequences of length 4 only 0000, 0001 and 1000 have 3 consecutive 0's. - Emeric Deutsch, Apr 27 2006
Therefore, the complementary sequence to A050231 (n coin tosses with a run of three heads). a(n) = 2^(n-3) - A050231(n-3) - Toby Gottfried, Nov 21 2010
Convolved with the Padovan sequence = row sums of triangle A153462. - Gary W. Adamson, Dec 27 2008
For n > 1: row sums of the triangle in A157897. - Reinhard Zumkeller, Jun 25 2009
a(n+2) is the top left entry of the n-th power of any of the 3 X 3 matrices [1, 1, 1; 0, 0, 1; 1, 0, 0] or [1, 1, 0; 1, 0, 1; 1, 0, 0] or [1, 1, 1; 1, 0, 0; 0, 1, 0] or [1, 0, 1; 1, 0, 0; 1, 1, 0]. - R. J. Mathar, Feb 03 2014
a(n-1) is the top left entry of the n-th power of any of the 3 X 3 matrices [0, 0, 1; 1, 1, 1; 0, 1, 0], [0, 1, 0; 0, 1, 1; 1, 1, 0], [0, 0, 1; 1, 0, 1; 0, 1, 1] or [0, 1, 0; 0, 0, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
Also row sums of A082601 and of A082870. - Reinhard Zumkeller, Apr 13 2014
Least significant bits are given in A021913 (a(n) mod 2 = A021913(n)). - Andres Cicuttin, Apr 04 2016
The nonnegative powers of the tribonacci constant t = A058265 are t^n = a(n)*t^2 + (a(n-1) + a(n-2))*t + a(n-1)*1, for n >= 0, with a(-1) = 1 and a(-2) = -1. This follows from the recurrences derived from t^3 = t^2 + t + 1. See the example in A058265 for the first nonnegative powers. For the negative powers see A319200. - Wolfdieter Lang, Oct 23 2018
The term "tribonacci number" was coined by Mark Feinberg (1963), a 14-year-old student in the 9th grade of the Susquehanna Township Junior High School in Pennsylvania. He died in 1967 in a motorcycle accident. - Amiram Eldar, Apr 16 2021
Andrews, Just, and Simay (2021, 2022) remark that it has been suggested that this sequence is mentioned in Charles Darwin's Origin of Species as bearing the same relation to elephant populations as the Fibonacci numbers do to rabbit populations. - N. J. A. Sloane, Jul 12 2022

Examples

			G.f. = x^2 + x^3 + 2*x^4 + 4*x^5 + 7*x^6 + 13*x^7 + 24*x^8 + 44*x^9 + 81*x^10 + ...
		

References

  • M. Agronomof, Sur une suite récurrente, Mathesis (Series 4), Vol. 4 (1914), pp. 125-126.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 47, ex. 4.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, Section 1.2.2.
  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • J. Riordan, An Introduction to Combinatorial Analysis, Princeton University Press, Princeton, NJ, 1978.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000045, A000078, A000213, A000931, A001590 (first differences, also a(n)+a(n+1)), A001644, A008288 (tribonacci triangle), A008937 (partial sums), A021913, A027024, A027083, A027084, A046738 (Pisano periods), A050231, A054668, A062544, A063401, A077902, A081172, A089068, A118390, A145027, A153462, A230216.
A057597 is this sequence run backwards: A057597(n) = a(1-n).
Row 3 of arrays A048887 and A092921 (k-generalized Fibonacci numbers).
Partitions: A240844 and A117546.
Cf. also A092836 (subsequence of primes), A299399 = A092835 + 1 (indices of primes).

Programs

  • GAP
    a:=[0,0,1];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # Muniru A Asiru, Oct 24 2018
  • Haskell
    a000073 n = a000073_list !! n
    a000073_list = 0 : 0 : 1 : zipWith (+) a000073_list (tail
                              (zipWith (+) a000073_list $ tail a000073_list))
    -- Reinhard Zumkeller, Dec 12 2011
    
  • Magma
    [n le 3 select Floor(n/3) else Self(n-1)+Self(n-2)+Self(n-3): n in [1..70]]; // Vincenzo Librandi, Jan 29 2016
    
  • Maple
    a:= n-> (<<0|1|0>, <0|0|1>, <1|1|1>>^n)[1,3]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Dec 19 2016
    # second Maple program:
    A000073:=proc(n) option remember; if n <= 1 then 0 elif n=2 then 1 else procname(n-1)+procname(n-2)+procname(n-3); fi; end; # N. J. A. Sloane, Aug 06 2018
  • Mathematica
    CoefficientList[Series[x^2/(1 - x - x^2 - x^3), {x, 0, 50}], x]
    a[0] = a[1] = 0; a[2] = 1; a[n_] := a[n] = a[n - 1] + a[n - 2] + a[n - 3]; Array[a, 36, 0] (* Robert G. Wilson v, Nov 07 2010 *)
    LinearRecurrence[{1, 1, 1}, {0, 0, 1}, 60] (* Vladimir Joseph Stephan Orlovsky, May 24 2011 *)
    a[n_] := SeriesCoefficient[If[ n < 0, x/(1 + x + x^2 - x^3), x^2/(1 - x - x^2 - x^3)], {x, 0, Abs @ n}] (* Michael Somos, Jun 01 2013 *)
    Table[-RootSum[-1 - # - #^2 + #^3 &, -#^n - 9 #^(n + 1) + 4 #^(n + 2) &]/22, {n, 0, 20}] (* Eric W. Weisstein, Nov 09 2017 *)
  • Maxima
    A000073[0]:0$
    A000073[1]:0$
    A000073[2]:1$
    A000073[n]:=A000073[n-1]+A000073[n-2]+A000073[n-3]$
      makelist(A000073[n], n, 0, 40);  /* Emanuele Munarini, Mar 01 2011 */
    
  • PARI
    {a(n) = polcoeff( if( n<0, x / ( 1 + x + x^2 - x^3), x^2 / ( 1 - x - x^2 - x^3) ) + x * O(x^abs(n)), abs(n))}; /* Michael Somos, Sep 03 2007 */
    
  • PARI
    my(x='x+O('x^99)); concat([0, 0], Vec(x^2/(1-x-x^2-x^3))) \\ Altug Alkan, Apr 04 2016
    
  • PARI
    a(n)=([0,1,0;0,0,1;1,1,1]^n)[1,3] \\ Charles R Greathouse IV, Apr 18 2016, simplified by M. F. Hasler, Apr 18 2018
    
  • Python
    def a(n, adict={0:0, 1:0, 2:1}):
        if n in adict:
            return adict[n]
        adict[n]=a(n-1)+a(n-2)+a(n-3)
        return adict[n] # David Nacin, Mar 07 2012
    from functools import cache
    @cache
    def A000073(n: int) -> int:
        if n <= 1: return 0
        if n == 2: return 1
        return A000073(n-1) + A000073(n-2) + A000073(n-3) # Peter Luschny, Nov 21 2022
    

Formula

G.f.: x^2/(1 - x - x^2 - x^3).
G.f.: x^2 / (1 - x / (1 - x / (1 + x^2 / (1 + x)))). - Michael Somos, May 12 2012
G.f.: Sum_{n >= 0} x^(n+2) *[ Product_{k = 1..n} (k + k*x + x^2)/(1 + k*x + k*x^2) ] = x^2 + x^3 + 2*x^4 + 4*x^5 + 7*x^6 + 13*x^7 + ... may be proved by the method of telescoping sums. - Peter Bala, Jan 04 2015
a(n+1)/a(n) -> A058265. a(n-1)/a(n) -> A192918.
a(n) = central term in M^n * [1 0 0] where M = the 3 X 3 matrix [0 1 0 / 0 0 1 / 1 1 1]. (M^n * [1 0 0] = [a(n-1) a(n) a(n+1)].) a(n)/a(n-1) tends to the tribonacci constant, 1.839286755... = A058265, an eigenvalue of M and a root of x^3 - x^2 - x - 1 = 0. - Gary W. Adamson, Dec 17 2004
a(n+2) = Sum_{k=0..n} T(n-k, k), where T(n, k) = trinomial coefficients (A027907). - Paul Barry, Feb 15 2005
A001590(n) = a(n+1) - a(n); A001590(n) = a(n-1) + a(n-2) for n > 1; a(n) = (A000213(n+1) - A000213(n))/2; A000213(n-1) = a(n+2) - a(n) for n > 0. - Reinhard Zumkeller, May 22 2006
Let C = the tribonacci constant, 1.83928675...; then C^n = a(n)*(1/C) + a(n+1)*(1/C + 1/C^2) + a(n+2)*(1/C + 1/C^2 + 1/C^3). Example: C^4 = 11.444...= 2*(1/C) + 4*(1/C + 1/C^2) + 7*(1/C + 1/C^2 + 1/C^3). - Gary W. Adamson, Nov 05 2006
a(n) = j*C^n + k*r1^n + L*r2^n where C is the tribonacci constant (C = 1.8392867552...), real root of x^3-x^2-x-1=0, and r1 and r2 are the two other roots (which are complex), r1 = m+p*i and r2 = m-p*i, where i = sqrt(-1), m = (1-C)/2 (m = -0.4196433776...) and p = ((3*C-5)*(C+1)/4)^(1/2) = 0.6062907292..., and where j = 1/((C-m)^2 + p^2) = 0.1828035330..., k = a+b*i, and L = a-b*i, where a = -j/2 = -0.0914017665... and b = (C-m)/(2*p*((C-m)^2 + p^2)) = 0.3405465308... . - Philippe LALLOUET (philip.lallouet(AT)wanadoo.fr), Jun 23 2007
a(n+1) = 3*c*((1/3)*(a+b+1))^n/(c^2-2*c+4) where a=(19+3*sqrt(33))^(1/3), b=(19-3*sqrt(33))^(1/3), c=(586+102*sqrt(33))^(1/3). Round to the nearest integer. - Al Hakanson (hawkuu(AT)gmail.com), Feb 02 2009
a(n) = round(3*((a+b+1)/3)^n/(a^2+b^2+4)) where a=(19+3*sqrt(33))^(1/3), b=(19-3*sqrt(33))^(1/3).. - Anton Nikonov
Another form of the g.f.: f(z) = (z^2-z^3)/(1-2*z+z^4). Then we obtain a(n) as a sum: a(n) = Sum_{i=0..floor((n-2)/4)} ((-1)^i*binomial(n-2-3*i,i)*2^(n-2-4*i)) - Sum_{i=0..floor((n-3)/4)} ((-1)^i*binomial(n-3-3*i,i)*2^(n-3-4*i)) with natural convention: Sum_{i=m..n} alpha(i) = 0 for m > n. - Richard Choulet, Feb 22 2010
a(n+2) = Sum_{k=0..n} Sum_{i=k..n, mod(4*k-i,3)=0} binomial(k,(4*k-i)/3)*(-1)^((i-k)/3)*binomial(n-i+k-1,k-1). - Vladimir Kruchinin, Aug 18 2010
a(n) = 2*a(n-2) + 2*a(n-3) + a(n-4). - Gary Detlefs, Sep 13 2010
Sum_{k=0..2*n} a(k+b)*A027907(n,k) = a(3*n+b), b >= 0 (see A099464, A074581).
a(n) = 2*a(n-1) - a(n-4), with a(0)=a(1)=0, a(2)=a(3)=1. - Vincenzo Librandi, Dec 20 2010
Starting (1, 2, 4, 7, ...) is the INVERT transform of (1, 1, 1, 0, 0, 0, ...). - Gary W. Adamson, May 13 2013
G.f.: Q(0)*x^2/2, where Q(k) = 1 + 1/(1 - x*(4*k+1 + x + x^2)/( x*(4*k+3 + x + x^2) + 1/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Sep 09 2013
a(n+2) = Sum_{j=0..floor(n/2)} Sum_{k=0..j} binomial(n-2*j,k)*binomial(j,k)*2^k. - Tony Foster III, Sep 08 2017
Sum_{k=0..n} (n-k)*a(k) = (a(n+2) + a(n+1) - n - 1)/2. See A062544. - Yichen Wang, Aug 20 2020
a(n) = A008937(n-1) - A008937(n-2) for n >= 2. - Peter Luschny, Aug 20 2020
From Yichen Wang, Aug 27 2020: (Start)
Sum_{k=0..n} a(k) = (a(n+2) + a(n) - 1)/2. See A008937.
Sum_{k=0..n} k*a(k) = ((n-1)*a(n+2) - a(n+1) + n*a(n) + 1)/2. See A337282. (End)
For n > 1, a(n) = b(n) where b(1) = 1 and then b(n) = Sum_{k=1..n-1} b(n-k)*A000931(k+2). - J. Conrad, Nov 24 2022
Conjecture: the congruence a(n*p^(k+1)) + a(n*p^k) + a(n*p^(k-1)) == 0 (mod p^k) holds for positive integers k and n and for all the primes p listed in A106282. - Peter Bala, Dec 28 2022
Sum_{k=0..n} k^2*a(k) = ((n^2-4*n+6)*a(n+1) - (2*n^2-2*n+5)*a(n) + (n^2-2*n+3)*a(n-1) - 3)/2. - Prabha Sivaramannair, Feb 10 2024
a(n) = Sum_{r root of x^3-x^2-x-1} r^n/(3*r^2-2*r-1). - Fabian Pereyra, Nov 23 2024

Extensions

Minor edits by M. F. Hasler, Apr 18 2018
Deleted certain dangerous or potentially dangerous links. - N. J. A. Sloane, Jan 30 2021

A208510 Triangle of coefficients of polynomials u(n,x) jointly generated with A029653; see the Formula section.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 4, 1, 1, 7, 9, 5, 1, 1, 9, 16, 14, 6, 1, 1, 11, 25, 30, 20, 7, 1, 1, 13, 36, 55, 50, 27, 8, 1, 1, 15, 49, 91, 105, 77, 35, 9, 1, 1, 17, 64, 140, 196, 182, 112, 44, 10, 1, 1, 19, 81, 204, 336, 378, 294, 156, 54, 11, 1, 1, 21, 100, 285, 540, 714, 672, 450, 210, 65, 12, 1
Offset: 1

Views

Author

Clark Kimberling, Feb 28 2012

Keywords

Comments

Row sums: A083329
Alternating row sums: 1,0,-1,-1,-1,-1,-1,-1,-1,-1,...
Antidiagonal sums: A000071 (-1+Fibonacci numbers)
col 1: A000012
col 2: A005408
col 3: A000290
col 4: A000330
col 5: A002415
col 6: A005585
col 7: A040977
col 8: A050486
col 9: A053347
col 10: A054333
col 11: A054334
col 12: A057788
col 2n-1 of A208510 is column n of A208508
col 2n of A208510 is column n of A208509.
...
GENERAL DISCUSSION:
A208510 typifies arrays generated by paired recurrence equations of the following form:
u(n,x)=a(n,x)*u(n-1,x)+b(n,x)*v(n-1,x)+c(n,x)
v(n,x)=d(n,x)*u(n-1,x)+e(n,x)*v(n-1,x)+f(n,x).
...
These first-order recurrences imply separate second-order recurrences. In order to show them, the six functions a(n,x),...,f(n,x) are abbreviated as a,b,c,d,e,f.
Then, starting with initial values u(1,x)=1 and u(2,x)=a+b+c: u(n,x) = (a+e)u(n-1,x) + (bd-ae)u(n-2,x) + bf-ce+c.
With initial values v(1,x)=1 and v(2,x)=d+e+f: v(n,x) = (a+e)v(n-1,x) + (bd-ae)v(n-2,x) + cd-af+f.
...
In the guide below, the last column codes certain sequences that occur in one of these ways: row, column, edge, row sum, alternating row sum. Coding:
A: 1,-1,1,-1,1,-1,1.... A033999
B: 1,2,4,8,16,32,64,... powers of 2
C: 1,1,1,1,1,1,1,1,.... A000012
D: 2,2,2,2,2,2,2,2,.... A007395
E: 2,4,6,8,10,12,14,... even numbers
F: 1,1,2,3,5,8,13,21,.. Fibonacci numbers
N: 1,2,3,4,5,6,7,8,.... A000027
O: 1,3,5,7,9,11,13,.... odd numbers
P: 1,3,9,27,81,243,.... powers of 3
S: 1,4,9,16,25,36,49,.. squares
T: 1,3,6,10,15,21,38,.. triangular numbers
Z: 1,0,0,0,0,0,0,0,0,.. A000007
*: (eventually) periodic alternating row sums
^: has a limiting row; i.e., the polynomials "approach" a power series
This coding includes indirect and repeated occurrences; e.g. F occurs thrice at A094441: in column 1 directly as Fibonacci numbers, in row sums as odd-indexed Fibonacci numbers, and in alternating row sums as signed Fibonacci numbers.
......... a....b....c....d....e....f....code
A034839 u 1....1....0....1....x....0....CCOT
A034867 v 1....1....0....1....x....0....CEN
A210221 u 1....1....0....1....2x...0....BBFF
A210596 v 1....1....0....1....2x...0....BBFF
A105070 v 1....2x...0....1....1....0....BN
A207605 u 1....1....0....1....x+1..0....BCFFN
A106195 v 1....1....0....1....x+1..0....BCFFN
A207606 u 1....1....0....x....x+1..0....DNT
A207607 v 1....1....0....x....x+1..0....DNT
A207608 u 1....1....0....2x...x+1..0....N
A207609 v 1....1....0....2x...x+1..0....C
A207610 u 1....1....0....1....x....1....CF
A207611 v 1....1....0....1....x....1....BCF
A207612 u 1....1....0....1....2x...1....BF
A207613 v 1....1....0....1....2x...1....BF
A207614 u 1....1....0....1....x+1..1....CN
A207615 v 1....1....0....1....x+1..1....CFN
A207616 u 1....1....0....x....1....1....CE
A207617 v 1....1....0....x....1....1....CNO
A029638 u 1....1....0....x....x....1....CDNO
A029635 v 1....1....0....x....x....1....CDNOZ
A207618 u 1....1....0....x....2x...1....N
A207619 v 1....1....0....x....2x...1....CFN
A207620 u 1....1....0....x....x+1..1....DET
A207621 v 1....1....0....x....x+1..1....DNO
A207622 u 1....1....0....2x...1....1....BT
A207623 v 1....1....0....2x...1....1....BN
A207624 u 1....1....0....2x...x....1....N
A102662 v 1....1....0....2x...x....1....CO
A207625 u 1....1....0....2x...x+1..1....T
A207626 v 1....1....0....2x...x+1..1....N
A207627 u 1....1....0....2x...2x...1....BN
A207628 v 1....1....0....2x...2x...1....BCE
A207629 u 1....1....0....x+1..1....1....CET
A207630 v 1....1....0....x+1..1....1....CO
A207631 u 1....1....0....x+1..x....1....DF
A207632 v 1....1....0....x+1..x....1....DEF
A207633 u 1....1....0....x+1..2x...1....F
A207634 v 1....1....0....x+1..2x...1....F
A207635 u 1....1....0....x+1..x+1..1....DN
A207636 v 1....1....0....x+1..x+1..1....CD
A160232 u 1....x....0....1....2x...0....BCFN
A208341 v 1....x....0....1....2x...0....BCFFN
A085478 u 1....x....0....1....x+1..0....CCOFT*
A078812 v 1....x....0....1....x+1..0....CEFN*
A208342 u 1....x....0....x....x....0....CCFNO
A208343 v 1....x....0....x....x....0....BBCDFZ
A208344 u 1....x....0....x....2x...0....CCFN
A208345 v 1....x....0....x....2x...0....CFZ
A094436 u 1....x....0....x....x+1..0....CFFN
A094437 v 1....x....0....x....x+1..0....CEFF
A117919 u 1....x....0....2x...1....0....BCNT
A135837 v 1....x....0....2x...1....0....BCET
A208328 u 1....x....0....2x...x....0....CCOP
A208329 v 1....x....0....2x...x....0....DPZ
A208330 u 1....x....0....2x...x+1..0....CNPT
A208331 v 1....x....0....2x...x+1..0....CN
A208332 u 1....x....0....2x...2x...0....CCE
A208333 v 1....x....0....2x...2x...0....DZ
A208334 u 1....x....0....x+1..1....0....CCNT
A208335 v 1....x....0....x+1..1....0....CCN*
A208336 u 1....x....0....x+1..x....0....CFNT*
A208337 v 1....x....0....x+1..x....0....ACFN*
A208338 u 1....x....0....x+1..2x...0....CNP
A208339 v 1....x....0....x+1..2x...0....BCNP
A202390 u 1....x....0....x+1..x+1..0....CFPTZ*
A208340 v 1....x....0....x+1..x+1..0....FNPZ*
A208508 u 1....x....0....1....1....1....CCES
A208509 v 1....x....0....1....1....1....BCO
A208510 u 1....x....0....1....x....1....CCCNOS*
A029653 v 1....x....0....1....x....1....BCDOSZ*
A208511 u 1....x....0....1....2x...1....BCFO
A208512 v 1....x....0....1....2x...1....BDFO
A208513 u 1....x....0....1....x+1..1....CCES*
A111125 v 1....x....0....1....x+1..1....COO*
A133567 u 1....x....0....x....1....1....CCOTT
A133084 v 1....x....0....x....1....1....BBCEN
A208514 u 1....x....0....x....x....1....CEFN
A208515 v 1....x....0....x....x....1....BCDFN
A208516 u 1....x....0....x....2x...1....CNN
A208517 v 1....x....0....x....2x...1....CCN
A208518 u 1....x....0....x....x+1..1....CFNT
A208519 v 1....x....0....x....x+1..1....NFFT
A208520 u 1....x....0....2x...1....1....BCTT
A208521 v 1....x....0....2x...1....1....BEN
A208522 u 1....x....0....2x...x....1....CCN
A208523 v 1....x....0....2x...x....1....CCO
A208524 u 1....x....0....2x...x+1..1....CT*
A208525 v 1....x....0....2x...x+1..1....ACNP*
A208526 u 1....x....0....2x...2x...1....CEN
A208527 v 1....x....0....2x...2x...1....CCE
A208606 u 1....x....0....x+1..1....1....CCS
A208607 v 1....x....0....x+1..1....1....CNO
A208608 u 1....x....0....x+1..x....1....CFOT
A208609 v 1....x....0....x+1..x....1....DEN*
A208610 u 1....x....0....x+1..2x...1....CO
A208611 v 1....x....0....x+1..2x...1....DE
A208612 u 1....x....0....x+1..x+1..1....CFNS
A208613 v 1....x....0....x+1..x+1..1....CFN*
A105070 u 1....2x...0....1....1....0....BN
A207536 u 1....2x...0....1....1....0....BCT
A208751 u 1....2x...0....1....x+1..0....CDPT
A208752 v 1....2x...0....1....x+1..0....CNP
A135837 u 1....2x...0....x....1....0....BCNT
A117919 v 1....2x...0....x....1....0....BCNT
A208755 u 1....2x...0....x....x....0....BCDEP
A208756 v 1....2x...0....x....x....0....BCCOZ
A208757 u 1....2x...0....x....2x...0....CDEP
A208758 v 1....2x...0....x....2x...0....CCEPZ
A208763 u 1....2x...0....2x...x....0....CDOP
A208764 v 1....2x...0....2x...x....0....CCCP
A208765 u 1....2x...0....2x...x+1..0....CE
A208766 v 1....2x...0....2x...x+1..0....CC
A208747 u 1....2x...0....2x...2x...0....CDE
A208748 v 1....2x...0....2x...2x...0....CCZ
A208749 u 1....2x...0....x+1..1....0....BCOPT
A208750 v 1....2x...0....x+1..1....0....BCNP*
A208759 u 1....2x...0....x+1..2x....0...CE
A208760 v 1....2x...0....x+1..2x....0...BCO
A208761 u 1....2x...0....x+1..x+1...0...BCCT*
A208762 v 1....2x...0....x+1..x+1...0...BNZ*
A208753 u 1....2x...0....1....1.....1...BCS
A208754 v 1....2x...0....1....1.....1...BO
A105045 u 1....2x...0....1....2x....1...BCCOS*
A208659 v 1....2x...0....1....2x....1...BDOSZ*
A208660 u 1....2x...0....1....x+1...1...CDS
A208904 v 1....2x...0....1....x+1...1...CNO
A208905 u 1....2x...0....x....1.....1...BCT
A208906 v 1....2x...0....x....1.....1...BNN
A208907 u 1....2x...0....x....x.....1...BCN
A208756 v 1....2x...0....x....x.....1...BCCE
A208755 u 1....2x...0....x....2x....1...CEN
A208910 v 1....2x...0....x....2x....1...CCE
A208911 u 1....2x...0....x....x+1...1...BCT
A208912 v 1....2x...0....x....x+1...1...BNT
A208913 u 1....2x...0....2x...1.....1...BCT
A208914 v 1....2x...0....2x...1.....1...BEN
A208915 u 1....2x...0....2x...x.....1...CE
A208916 v 1....2x...0....2x...x.....1...CCO
A208919 u 1....2x...0....2x...x+1...1...CT
A208920 v 1....2x...0....2x...x+1...1...N
A208917 u 1....2x...0....2x...2x....1...CEN
A208918 v 1....2x...0....2x...2x....1...CCNP
A208921 u 1....2x...0....x+1..1.....1...BC
A208922 v 1....2x...0....x+1..1.....1...BON
A208923 u 1....2x...0....x+1..x.....1...BCNO
A208908 v 1....2x...0....x+1..x.....1...BDN*
A208909 u 1....2x...0....x+1..2x....1...BN
A208930 v 1....2x...0....x+1..2x....1...DN
A208931 u 1....2x...0....x+1..x+1...1...BCOS
A208932 v 1....2x...0....x+1..x+1...1...BCO*
A207537 u 1....x+1..0....1....1.....0...BCO
A207538 v 1....x+1..0....1....1.....0...BCE
A122075 u 1....x+1..0....1....x.....0...CCFN*
A037027 v 1....x+1..0....1....x.....0...CCFN*
A209125 u 1....x+1..0....1....2x....0...BCFN*
A164975 v 1....x+1..0....1....2x....0...BF
A209126 u 1....x+1..0....x....x.....0...CDFO*
A209127 v 1....x+1..0....x....x.....0...DFOZ*
A209128 u 1....x+1..0....x....2x....0...CDE*
A209129 v 1....x+1..0....x....2x....0...DEZ
A102756 u 1....x+1..0....x....x+1...0...CFNP*
A209130 v 1....x+1..0....x....x+1...0...CCFNP*
A209131 u 1....x+1..0....2x...x.....0...CDEP*
A209132 v 1....x+1..0....2x...x.....0...CNPZ*
A209133 u 1....x+1..0....2x...2x....0...CDN
A209134 v 1....x+1..0....2x...2x....0...CCN*
A209135 u 1....x+1..0....2x...x+1...0...CN*
A209136 v 1....x+1..0....2x...x+1...0...CCS*
A209137 u 1....x+1..0....x+1..x.....0...CFFP*
A209138 v 1....x+1..0....x+1..x.....0...AFFP*
A209139 u 1....x+1..0....x+1..2x....0...CF*
A209140 v 1....x+1..0....x+1..2x....0...BF
A209141 u 1....x+1..0....x+1..x+1...0...BCF*
A209142 v 1....x+1..0....x+1..x+1...0...BFZ*
A209143 u 1....x+1..0....1....1.....1...CCE*
A209144 v 1....x+1..0....1....1.....1...COO*
A209145 u 1....x+1..0....1....x.....1...CCFN*
A122075 v 1....x+1..0....1....x.....1...CCFN*
A209146 u 1....x+1..0....1....2x....1...BCF*
A209147 v 1....x+1..0....1....2x....1...BF
A209148 u 1....x+1..0....1....x+1...1...CCO*
A209149 v 1....x+1..0....1....x+1...1...CDO*
A209150 u 1....x+1..0....x....1.....1...CCNT*
A208335 v 1....x+1..0....x....1.....1...CDNN*
A209151 u 1....x+1..0....x....x.....1...CFN*
A208337 v 1....x+1..0....x....x.....1...ACFN*
A209152 u 1....x+1..0....x....2x....1...CN*
A208339 v 1....x+1..0....x....x.....1...BCN
A209153 u 1....x+1..0....x....x+1...1...CFT*
A208340 v 1....x+1..0....x....x.....1...FNZ*
A209154 u 1....x+1..0....2x...1.....1...BCT*
A209157 v 1....x+1..0....2x...1.....1...BNN
A209158 u 1....x+1..0....2x...x.....1...CN*
A209159 v 1....x+1..0....2x...x.....1...CO*
A209160 u 1....x+1..0....2x...2x....1...CN*
A209161 v 1....x+1..0....2x...2x....1...CE
A209162 u 1....x+1..0....2x...x+1...1...CT*
A209163 v 1....x+1..0....2x...x+1...1...CO*
A209164 u 1....x+1..0....x+1..1.....1...CC*
A209165 v 1....x+1..0....x+1..1.....1...CCN
A209166 u 1....x+1..0....x+1..x.....1...CFF*
A209167 v 1....x+1..0....x+1..x.....1...FF*
A209168 u 1....x+1..0....x+1..2x....1...CF*
A209169 v 1....x+1..0....x+1..2x....1...CF
A209170 u 1....x+1..0....x+1..x+1...1...CF*
A209171 v 1....x+1..0....x+1..x+1...1...CF*
A053538 u x....1....0....1....1.....0...BBCCFN
A076791 v x....1....0....1....1.....0...BBCDF
A209172 u x....1....0....1....2x....0...BCCFF
A209413 v x....1....0....1....2x....0...BCCFF
A094441 u x....1....0....1....x+1...0...CFFFN
A094442 v x....1....0....1....x+1...0...CEFFF
A054142 u x....1....0....x....x+1...0...CCFOT*
A172431 v x....1....0....x....x+1...0...CEFN*
A008288 u x....1....0....2x...1.....0...CCOO*
A035607 v x....1....0....2x...1.....0...ACDE*
A209414 u x....1....0....2x...x+1...0...CCS
A112351 v x....1....0....2x...x+1...0...CON
A209415 u x....1....0....x+1..x.....0...CCTN
A209416 v x....1....0....x+1..x.....0...ACN*
A209417 u x....1....0....x+1..2x....0...CC
A209418 v x....1....0....x+1..2x....0...BBC
A209419 u x....1....0....x+1..x+1...0...CFTZ*
A209420 v x....1....0....x+1..x+1...0...FNZ*
A209421 u x....1....0....1....1.....1...CCN
A209422 v x....1....0....1....1.....1...CD
A209555 u x....1....0....1....x.....1...CNN
A209556 v x....1....0....1....x.....1...CNN
A209557 u x....1....0....1....2x....1...BCN
A209558 v x....1....0....1....2x....1...BN
A209559 u x....1....0....1....x+1...1...CN
A209560 v x....1....0....1....x+1...1...CN
A209561 u x....1....0....x....1.....1...CCNNT*
A209562 v x....1....0....x....1.....1...CDNNT*
A209563 u x....1....0....x....x.....1...CCFT^
A209564 v x....1....0....x....x.....1...CFN^
A209565 u x....1....0....x....2x....1...CC^
A209566 v x....1....0....x....2x....1...BC^
A209567 u x....1....0....x....x+1...1...CNT*
A209568 v x....1....0....x....x+1...1...NNS*
A209569 u x....1....0....2x...1.....1...CNO*
A209570 v x....1....0....2x...1.....1...DNN*
A209571 u x....1....0....2x...x.....1...CCS^
A209572 v x....1....0....2x...x.....1...CN^
A209573 u x....1....0....2x...x+1...1...CNS
A209574 v x....1....0....2x...x+1...1...NO
A209575 u x....1....0....2x...2x....1...CC
A209576 v x....1....0....2x...2x....1...C
A209577 u x....1....0....x+1..1.....1...CNNT
A209578 v x....1....0....x+1..1.....1...CNN
A209579 u x....1....0....x+1..x.....1...CNNT
A209580 v x....1....0....x+1..x.....1...NN*
A209581 u x....1....0....x+1..2x....1...CN
A209582 v x....1....0....x+1..2x....1...BN
A209583 u x....1....0....x+1..x+1...1...CT*
A209584 v x....1....0....x+1..x+1...1...CN*
A121462 u x....x....0....x....x+1...0...BCFFNZ
A208341 v x....x....0....x....x+1...0...BCFFN
A209687 u x....x....0....2x...x+1...0...BCNZ
A208339 v x....x....0....2x...x+1...0...BCN
A115241 u x....x....0....1....1.....1...CDNZ*
A209688 v x....x....0....1....1.....1...DDN*
A209689 u x....x....0....1....x.....1...FNZ^
A209690 v x....x....0....1....x.....1...FN^
A209691 u x....x....0....1....2x....1...BCZ^
A209692 v x....x....0....1....2x....1...BCC^
A209693 u x....x....0....1....x+1...1...NNZ*
A209694 v x....x....0....1....x+1...1...CN*
A209697 u x....x....0....x....x+1...1...BNZ
A209698 v x....x....0....x....x+1...1...BNT
A209699 u x....x....0....2x...1.....1...BNNZ
A209700 v x....x....0....2x...1.....1...BDN
A209701 u x....x....0....2x...x+1...1...NZ
A209702 v x....x....0....2x...x+1...1...N
A209703 u x....x....0....x+1..1.....1...FNTZ
A209704 v x....x....0....x+1..1.....1...FNNT
A209705 u x....x....0....x+1..x+1...1...BNZ*
A209706 v x....x....0....x+1..x+1...1...BCN*
A209695 u x....x+1..0....2x...x+1...0...ACN*
A209696 v x....x+1..0....2x...x+1...0...CDN*
A209830 u x....x+1..0....x+1..2x....0...ACF
A209831 v x....x+1..0....x+1..2x....0...BCF*
A209745 u x....x+1..0....x+1..x+1...0...ABF*
A209746 v x....x+1..0....x+1..x+1...0...BFZ*
A209747 u x....x+1..0....1....1.....1...ADE*
A209748 v x....x+1..0....1....1.....1...DEO
A209749 u x....x+1..0....1....x.....1...ANN*
A209750 v x....x+1..0....1....x.....1...CNO
A209751 u x....x+1..0....1....2x....1...ABN*
A209752 v x....x+1..0....1....2x....1...BN
A209753 u x....x+1..0....1....x+1...1...AN*
A209754 v x....x+1..0....1....x+1...1...NT*
A209755 u x....x+1..0....x....1.....1...AFN
A209756 v x....x+1..0....x....1.....1...FNO*
A209759 u x....x+1..0....x....2x....1...ACF^
A209760 v x....x+1..0....x....2x....1...CF^*
A209761 u x....x+1..0....x.....x+1..1...ABNS*
A209762 v x....x+1..0....x.....x+1..1...BNS*
A209763 u x....x+1..0....2x....1....1...ABN*
A209764 v x....x+1..0....2x....1....1...BNN
A209765 u x....x+1..0....2x....x....1...ACF^*
A209766 v x....x+1..0....2x....x....1...CF^
A209767 u x....x+1..0....2x....x+1..1...AN*
A209768 v x....x+1..0....2x....x+1..1...N*
A209769 u x....x+1..0....x+1...1....1...AF*
A209770 v x....x+1..0....x+1...1....1...FN
A209771 u x....x+1..0....x+1...x....1...ABN*
A209772 v x....x+1..0....x+1...x....1...BN*
A209773 u x....x+1..0....x+1...2x...1...AF
A209774 v x....x+1..0....x+1...2x...1...FN*
A209775 u x....x+1..0....x+1...x+1..1...AB*
A209776 v x....x+1..0....x+1...x+1..1...BC*
A210033 u 1....1....1....1.....x....1...BCN
A210034 v 1....1....1....1.....x....1...BCDFN
A210035 u 1....1....1....1.....2x...1...BBF
A210036 v 1....1....1....1.....2x...1...BBFF
A210037 u 1....1....1....1.....x+1..1...BCFFN
A210038 v 1....1....1....1.....x+1..1...BCFFN
A210039 u 1....1....1....x.....1....1...BCOT
A210040 v 1....1....1....x.....1....1...BCEN
A210042 u 1....1....1....x.....x....1...BCDEOT*
A124927 v 1....1....1....x.....x....1...BCDET*
A210041 u 1....1....1....x.....2x...1...BFO
A209758 v 1....1....1....x.....2x...1...BCFO
A210187 u 1....1....1....x.....x+1..1...DTF*
A210188 v 1....1....1....x.....x+1..1...DNF*
A210189 u 1....1....1....2x....1....1...BT
A210190 v 1....1....1....2x....1....1...BN
A210191 u 1....1....1....2x....x....1...CO*
A210192 v 1....1....1....2x....x....1...CCO*
A210193 u 1....1....1....2x....x+1..1...CPT
A210194 v 1....1....1....2x....x+1..1...CN
A210195 u 1....1....1....2x....2x...1...BOPT*
A210196 v 1....1....1....2x....2x...1...BCC*
A210197 u 1....1....1....x+1...1....1...BCOT
A210198 v 1....1....1....x+1...1....1...BCEN
A210199 u 1....1....1....x+1...x....1...DFT
A210200 v 1....1....1....x+1...x....1...DFO*
A210201 u 1....1....1....x+1...2x...1...BFP
A210202 v 1....1....1....x+1...2x...1...BF
A210203 u 1....1....1....x+1...x+1..1...BDOP
A210204 v 1....1....1....x+1...x+1..1...BCDN*
A210211 u x....1....1....1.....2x...1...BCFN
A210212 v x....1....1....1.....2x...1...BFN
A210213 u x....1....1....1.....x+1..1...CFFN
A210214 v x....1....1....1.....x+1..1...CFFO
A210215 u x....1....1....x.....x....1...BCDFT^
A210216 v x....1....1....x.....x....1...BCFO^
A210217 u x....1....1....x.....2x...1...CDF^
A210218 v x....1....1....x.....2x...1...BCF^
A210219 u x....1....1....x.....x+1..1...CNSTF*
A210220 v x....1....1....x.....x+1..1...FNNT*
A104698 u x....1....1....2x......1..1...CENS*
A210220 v x....1....1....2x....x+1..1...DNNT*
A210223 u x....1....1....2x....x....1...CD^
A210224 v x....1....1....2x....x....1...CO^
A210225 u x....1....1....2x....x+1..1...CNP
A210226 v x....1....1....2x....x+1..1...NOT
A210227 u x....1....1....2x....2x...1...CDP^
A210228 v x....1....1....2x....2x...1...C^
A210229 u x....1....1....x+1...1....1...CFNN
A210230 v x....1....1....x+1...1....1...CCN
A210231 u x....1....1....x+1...x....1...CNT
A210232 v x....1....1....x+1...x....1...NN*
A210233 u x....1....1....x+1...2x...1...CNP
A210234 v x....1....1....x+1...2x...1...BN
A210235 u x....1....1....x+1...x+1..1...CCFPT*
A210236 v x....1....1....x+1...x+1..1...CFN*
A124927 u x....x....1....1.....1....1...BCDEET*
A210042 v x....1....1....x+1...x+1..1...BDEOT*
A210216 u x....x....1....1.....x....1...BCFO^
A210215 v x....x....1....1.....x....1...BCDFT^
A210549 u x....x....1....1.....2x...1...BCF^
A210550 v x....x....1....1.....2x...1...BDF^
A172431 u x....x....1....1.....x+1..1...CEFN*
A210551 v x....x....1....1.....x+1..1...CFOT*
A210552 u x....x....1....x.....1....1...BBCFNO
A210553 v x....x....1....x.....1....1...BNNFB
A208341 u x....x....1....x.....x+1..1...BCFFN
A210554 v x....x....1....x.....x+1..1...BNFFT
A210555 u x....x....1....2x....1....1...BCNN
A210556 v x....x....1....2x....1....1...BENP
A210557 u x....x....1....2x....x+1..1...CNP
A210558 v x....x....1....2x....x+1..1...N
A210559 u x....x....1....x+1...1....1...CEF
A210560 v x....x....1....x+1...1....1...OFNS
A210561 u x....x....1....x+1...x....1...BCNP^
A210562 v x....x....1....x+1...x....1...BDP*^
A210563 u x....x....1....x+1...2x...1...CFP^
A210564 v x....x....1....x+1...2x...1...DF^
A013609 u x....x....1....x+1...x+1..1...BCEPT*
A209757 v x....x....1....x+1...x+1..1...BCOS*
A209819 u x....2x...1....x+1...x....1...CFN^
A209820 v x....2x...1....x+1...x....1...DF^
A209996 u x....2x...1....x+1...2x...1...CP^
A209998 v x....2x...1....x+1...2x...1...DP^
A209999 u x....x+1..1....1.....x+1..1...FN*
A210287 v x....x+1..1....1.....x+1..1...CFT*
A210565 u x....x+1..1....x.....1....1...FNT*
A210595 v x....x+1..1....x.....1....1...FNNT
A210598 u x....x+1..1....x+1...2x...1...FN*
A210599 v x....x+1..1....x+1...2x...1...FN
A210600 u x....x+1..1....x+1...x+1..1...BF*
A210601 v x....x+1..1....x+1...x+1..1...BF*
A210597 u 2x...1....1....x+1...1....1...BF
A210601 v 2x...1....1....x+1...1....1...BFN*
A210603 u 2x...1....1....x+1...x+1..1...BF
A210738 v 2x...1....1....x+1...x+1..1...CBF*
A210739 u 2x...x....1....x+1...x....1...CF^
A210740 v 2x...x....1....x+1...x....1...DF*^
A210741 u 2x...x....1....x+1...x+1..1...BCFO
A210742 v 2x...x....1....x+1...x+1..1...CFO*
A210743 u 2x...x+1..1....x+1...1....1...F
A210744 v 2x...x+1..1....x+1...1....1...FN
A210747 u 2x...x+1..1....x+1...x+1..1...FF
A210748 v 2x...x+1..1....x+1...x+1..1...CFF*
A210749 u x+1..1....1....x+1...2x...1...BCF
A210750 v x+1..1....1....x+1...2x...1...BF
A210751 u x+1..x....1....x+1...2x...1...FNT
A210752 v x+1..x....1....x+1...2x...1...FN
A210753 u x+1..x....1....x+1...x+1..1...BNZ*
A210754 v x+1..x....1....x+1...x+1..1...BCT*
A210755 u x+1..2x...1....x+1...x+1..1...N*
A210756 v x+1..2x...1....x+1...x+1..1...CT*
A210789 u 1....x....0....x+2...x-1..0...CFFN
A210790 v 1....x....0....x+2...x-1..0...CEFF
A210791 u 1....x....0....x-1...x+2..0...CFNP
A210792 v 1....x....0....x-1...x+2..0...CF
A210793 u 1....x+1..0....x+2...x-1..0...CFNP
A210794 v 1....x+1..0....x+2...x-1..0...FPP
A210795 u 1....x....1....x+2...x-1..0...FN
A210796 v 1....x....1....x+2...x-1..0...FO
A210797 u 1....x....0....x+2...x-1..1...CF
A210798 v 1....x....0....x+2...x-1..1...F
A210799 u 1....x+1..1....x+2...x-1..0...FN
A210800 v 1....x+1..1....x+2...x-1..0...F
A210801 u 1....x+1..1....x+2...x-1..1...FN
A210802 v 1....x+1..1....x+2...x-1..1...F
A210803 u 1....x....0....x-1...x+3..0...F*
A210804 v 1....x....0....x-1...x+3..0...F*
A210805 u 1....x....0....x+2...x-1.-1...CFFN
A210806 v 1....x....0....x+2...x-1.-1...FF
A210858 u 1....x....0....x+n...x....0...CFT*
A210859 v 1....x....0....x+n...x....0...FN*
A210860 u 1....x+1..0....x+n...x....0...F
A210861 v 1....x+1..0....x+n...x....0...F*
A210862 u 1....x....1....x+n-1.x....0...FN
A210863 v 1....x....1....x+n-1.x....0...FS
A210864 u 1....x....1....x+n...x....0...FN
A210865 v 1....x....1....x+n...x....0...FT
A210866 u 1....x....0....x+n...x...-x...CFT
A210867 v 1....x....0....x+n...x...-x...FN
A210868 u 1....x....0....x+1...x-1..0...BCFN
A210869 v 1....x....0....x+1...x-1..0...BBCFNZ
A210870 u 1....x....0....x+1...x-1..1...CFFN
A210871 v 1....x....0....x+1...x-1..1...CFF
A210872 u x....1...-1....x.....x....1...BDFZ^
A210873 v x....1...-1....x.....x....1...BCFN^
A210876 u x....1....1....x.....x....x...BCCF^
A210877 v x....1....1....x.....x....x...BDFNZ^
A210878 u x....2x...0....x+1...x....1...DFZ^
A210879 v x....2x...0....x+1...x....1...FC*^
Some of these triangles have irregular row lengths, making it difficult to retrieve individual rows/columns/diagonals without actually computing the recurrence. - Georg Fischer, Sep 04 2021

Examples

			First five rows:
1
1...1
1...3...1
1...5...4...1
1...7...9...5...1
First five polynomials u(n,x):
1
1 + x
1 + 3x + x^2
1 + 5x + 4x^2 + x^3
1 + 7x + 9x^2 + 5x^3 + x^4
		

Crossrefs

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x];
    v[n_, x_] := u[n - 1, x] + x*v[n - 1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]   (* A208510 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A029653 *)
  • Python
    from sympy import Poly
    from sympy.abc import x
    def u(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x)
    def v(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x) + 1
    def a(n): return Poly(u(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 27 2017

Formula

u(n,x)=u(n-1,x)+x*v(n-1,x),
v(n,x)=u(n-1,x)+x*v(n-1,x)+1,
where u(1,x)=1, v(1,x)=1.
Also, u(n,x)=(x+1)*u(n-1,x)+x for n>2, with u(n,2)=x+1.

Extensions

Corrected by Philippe Deléham, Apr 10 2012
Corrections and additions by Clark Kimberling, May 09 2012
Corrections in the overview by Georg Fischer, Sep 04 2021

A001844 Centered square numbers: a(n) = 2*n*(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.

Original entry on oeis.org

1, 5, 13, 25, 41, 61, 85, 113, 145, 181, 221, 265, 313, 365, 421, 481, 545, 613, 685, 761, 841, 925, 1013, 1105, 1201, 1301, 1405, 1513, 1625, 1741, 1861, 1985, 2113, 2245, 2381, 2521, 2665, 2813, 2965, 3121, 3281, 3445, 3613, 3785, 3961, 4141, 4325, 4513
Offset: 0

Views

Author

Keywords

Comments

These are Hogben's central polygonal numbers denoted by
...2...
....P..
...4.n.
Numbers of the form (k^2+1)/2 for k odd.
(y(2x+1))^2 + (y(2x^2+2x))^2 = (y(2x^2+2x+1))^2. E.g., let y = 2, x = 1; (2(2+1))^2 + (2(2+2))^2 = (2(2+2+1))^2, (2(3))^2 + (2(4))^2 = (2(5))^2, 6^2 + 8^2 = 10^2, 36 + 64 = 100. - Glenn B. Cox (igloos_r_us(AT)canada.com), Apr 08 2002
a(n) is also the number of 3 X 3 magic squares with sum 3(n+1). - Sharon Sela (sharonsela(AT)hotmail.com), May 11 2002
For n > 0, a(n) is the smallest k such that zeta(2) - Sum_{i=1..k} 1/i^2 <= zeta(3) - Sum_{i=1..n} 1/i^3. - Benoit Cloitre, May 17 2002
Number of convex polyominoes with a 2 X (n+1) minimal bounding rectangle.
The prime terms are given by A027862. - Lekraj Beedassy, Jul 09 2004
First difference of a(n) is 4n = A008586(n). Any entry k of the sequence is followed by k + 2*(1 + sqrt(2k - 1)). - Lekraj Beedassy, Jun 04 2006
Integers of the form 1 + x + x^2/2 (generating polynomial is Schur's polynomial as in A127876). - Artur Jasinski, Feb 04 2007
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n-4) is equal to the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Row sums of triangle A132778. - Gary W. Adamson, Sep 02 2007
Binomial transform of [1, 4, 4, 0, 0, 0, ...]; = inverse binomial transform of A001788: (1, 6, 24, 80, 240, ...). - Gary W. Adamson, Sep 02 2007
Narayana transform (A001263) of [1, 4, 0, 0, 0, ...]. Equals A128064 (unsigned) * [1, 2, 3, ...]. - Gary W. Adamson, Dec 29 2007
k such that the Diophantine equation x^3 - y^3 = x*y + k has a solution with y = x-1. If that solution is (x,y) = (m+1,m) then m^2 + (m+1)^2 = k. Note that this Diophantine equation is an elliptic curve and (m+1,m) is an integer point on it. - James R. Buddenhagen, Aug 12 2008
Numbers k such that (k, k, 2*k-2) are the sides of an isosceles triangle with integer area. Also, k such that 2*k-1 is a square. - James R. Buddenhagen, Oct 17 2008
a(n) is also the least weight of self-conjugate partitions having n+1 different odd parts. - Augustine O. Munagi, Dec 18 2008
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) = A153869 * (1, 2, 3, ...). - Gary W. Adamson, Jan 03 2009
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) where a(n) = 2n*(n-1)+1, all tuples of square numbers (X-Y, X, X+Y) are produced by ((m*(a(n)-2n))^2, (m*a(n))^2, (m*(a(n)+2n-2))^2) where m is a whole number. - Doug Bell, Feb 27 2009
Equals (1, 2, 3, ...) convolved with (1, 3, 4, 4, 4, ...). E.g., a(3) = 25 = (1, 2, 3, 4) dot (4, 4, 3, 1) = (4 + 8 + 9 + 4). - Gary W. Adamson, May 01 2009
The running sum of squares taken two at a time. - Al Hakanson (hawkuu(AT)gmail.com), May 18 2009
Equals the odd integers convolved with (1, 2, 2, 2, ...). - Gary W. Adamson, May 25 2009
Equals the triangular numbers convolved with [1, 2, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
When the positive integers are written in a square array by diagonals as in A038722, a(n) gives the numbers appearing on the main diagonal. - Joshua Zucker, Jul 07 2009
The finite continued fraction [n,1,1,n] = (2n+1)/(2n^2 + 2n + 1) = (2n+1)/a(n); and the squares of the first two denominators of the convergents = a(n). E.g., the convergents and value of [4,1,1,4] = 1/4, 1/5, 2/9, 9/41 where 4^2 + 5^2 = 41. - Gary W. Adamson, Jul 15 2010
From Keith Tyler, Aug 10 2010: (Start)
Running sum of A008574.
Square open pyramidal number; that is, the number of elements in a square pyramid of height (n) with only surface and no bottom nodes. (End)
For k>0, x^4 + x^2 + k factors over the integers iff sqrt(k) is in this sequence. - James R. Buddenhagen, Aug 15 2010
Create the simple continued fraction from Pythagorean triples to get [2n + 1; 2n^2 + 2n, 2n^2 + 2n + 1]; its value equals the rational number 2n + 1 + a(n) / (4n^4 + 8n^3 + 6n^2 + 2n + 1). - J. M. Bergot, Sep 30 2011
a(n), n >= 1, has in its prime number factorization only primes of the form 4*k+1, i.e., congruent to 1 (mod 4) (see A002144). This follows from the fact that a(n) is a primitive sum of two squares and odd. See Theorem 3.20, p. 164, in the given Niven-Zuckerman-Montgomery reference. E.g., a(3) = 25 = 5^2, a(6) = 85 = 5*17. - Wolfdieter Lang, Mar 08 2012
From Ant King, Jun 15 2012: (Start)
a(n) is congruent to 1 (mod 4) for all n.
The digital roots of the a(n) form a purely periodic palindromic 9-cycle 1, 5, 4, 7, 5, 7, 4, 5, 1.
The units' digits of the a(n) form a purely periodic palindromic 5-cycle 1, 5, 3, 5, 1.
(End)
Number of integer solutions (x,y) of |x| + |y| <= n. Geometrically: number of lattice points inside a square with vertices (n,0), (0,-n), (-n,0), (0,n). - César Eliud Lozada, Sep 18 2012
(a(n)-1)/a(n) = 2*x / (1+x^2) where x = n/(n+1). Note that in this form, this is the velocity-addition formula according to the special theory of relativity (two objects traveling at 1/(n+1) slower than c relative to each other appear to travel at 1/a(n) less than c to a stationary observer). - Christian N. K. Anderson, May 20 2013 [Corrected by Rémi Guillaume, May 22 2025]
A geometric curiosity: the envelope of the circles x^2 + (y-a(n)/2)^2 = ((2n+1)/2)^2 is the parabola y = x^2, the n=0 circle being the osculating circle at the parabola vertex. - Jean-François Alcover, Dec 02 2013
Draw n ellipses in the plane (n>0), any 2 meeting in 4 points; a(n-1) gives the number of internal regions into which the plane is divided (cf. A051890, A046092); a(n-1) = A051890(n) - 1 = A046092(n-1) + 1. - Jaroslav Krizek, Dec 27 2013
a(n) is also, of course, the scalar product of the 2-vector (n, n+1) (or (n+1, n)) with itself. The unique inverse of (n, n+1) as vector in the Clifford algebra over the Euclidean 2-space is (1/a(n))(0, n, n+1, 0) (similarly for the other vector). In general the unique inverse of such a nonzero vector v (odd element in Cl_2) is v^(-1) = (1/|v|^2) v. Note that the inverse with respect to the scalar product is not unique for any nonzero vector. See the P. Lounesto reference, sects. 1.7 - 1.12, pp. 7-14. See also the Oct 15 2014 comment in A147973. - Wolfdieter Lang, Nov 06 2014
Subsequence of A004431, for n >= 1. - Bob Selcoe, Mar 23 2016
Numbers k such that 2k - 1 is a perfect square. - Juri-Stepan Gerasimov, Apr 06 2016
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 574", based on the 5-celled von Neumann neighborhood. - Robert Price, May 13 2016
a(n) is the first integer in a sum of (2*n + 1)^2 consecutive integers that equals (2*n + 1)^4. - Patrick J. McNab, Dec 24 2016
Central elements of odd-length rows of the triangular array of positive integers. a(n) is the mean of the numbers in the (2*n + 1)-th row of this triangle. - David James Sycamore, Aug 01 2018
Intersection of A000982 and A080827. - David James Sycamore, Aug 07 2018
An off-diagonal of the array of Delannoy numbers, A008288, (or a row/column when the array is shown as a square). As such, this is one of the crystal ball sequences. - Jack W Grahl, Feb 15 2021 and Shel Kaphan, Jan 18 2023
a(n) appears as a solution to a "Riddler Express" puzzle on the FiveThirtyEight website. The Jan 21 2022 issue (problem) and the Jan 28 2022 issue (solution) present the following puzzle and include a proof. - Fold a square piece of paper in half, obtaining a rectangle. Fold again to obtain a square with 1/4 the size of the original square. Then make n cuts through the folded paper. a(n) is the greatest number of pieces of the unfolded paper after the cutting. - Manfred Boergens, Feb 22 2022
a(n) is (1/6) times the number of 2 X 2 triangles in the n-th order hexagram with 12*n^2 cells. - Donghwi Park, Feb 06 2024
If k is a centered square number, its index in this sequence is n = (sqrt(2k-1)-1)/2. - Rémi Guillaume, Mar 30 2025.
Row sums of the symmetric triangle of odd numbers [1]; [1, 3, 1]; [1, 3, 5, 3, 1]; [1, 3, 5, 7, 5, 3, 1]; .... - Marco Zárate, Jun 15 2025

Examples

			G.f.: 1 + 5*x + 13*x^2 + 25*x^3 + 41*x^4 + 61*x^5 + 85*x^6 + 113*x^7 + 145*x^8 + ...
The first few triples are (1,0,1), (3,4,5), (5,12,13), (7,24,25), ...
The first four such partitions, corresponding to n = 0,1,2,3, i.e., to a(n) = 1,5,13,25, are 1, 3+1+1, 5+3+3+1+1, 7+5+5+3+3+1+1. - _Augustine O. Munagi_, Dec 18 2008
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 3.
  • A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, p. 125, 1964.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 50.
  • Pertti Lounesto, Clifford Algebras and Spinors, second edition, Cambridge University Press, 2001.
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 483.
  • Ivan Niven, Herbert S. Zuckerman and Hugh L. Montgomery, An Introduction to the Theory Of Numbers, Fifth Edition, John Wiley and Sons, Inc., NY 1991.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Travers et al., The Mysterious Lost Proof, Using Advanced Algebra, (1976), pp. 27.

Crossrefs

X values are A005408; Y values are A046092.
Cf. A008586 (first differences), A005900 (partial sums), A254373 (digital roots).
Subsequence of A004431.
Right edge of A055096; main diagonal of A069480, A078475, A129312.
Row n=2 (or column k=2) of A008288.
Cf. A016754.

Programs

  • Haskell
    a001844 n = 2 * n * (n + 1) + 1
    a001844_list = zipWith (+) a000290_list $ tail a000290_list
    -- Reinhard Zumkeller, Dec 04 2012
    
  • Magma
    [2*n^2 + 2*n + 1: n in [0..50]]; // Vincenzo Librandi, Jan 19 2013
    
  • Magma
    [n: n in [0..4400] | IsSquare(2*n-1)]; // Juri-Stepan Gerasimov, Apr 06 2016
    
  • Maple
    A001844:=-(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation
  • Mathematica
    Table[2n(n + 1) + 1, {n, 0, 50}]
    FoldList[#1 + #2 &, 1, 4 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
    maxn := 47; Flatten[Table[SeriesCoefficient[Series[(n + (n - 1)*x)/(1 - x)^2, {x, 0, maxn}], k], {n, maxn}, {k, n - 1, n - 1}]] (* L. Edson Jeffery, Aug 24 2014 *)
    CoefficientList[ Series[-(x^2 + 2x + 1)/(x - 1)^3, {x, 0, 48}], x] (* or *)
    LinearRecurrence[{3, -3, 1}, {1, 5, 13}, 48] (* Robert G. Wilson v, Aug 01 2018 *)
    Total/@Partition[Range[0,50]^2,2,1] (* Harvey P. Dale, Dec 05 2020 *)
    Table[ j! Coefficient[Series[Exp[x]*(1 + 4*x + 2*x^2), {x, 0, 20}], x,
    j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 07 2023 *)
  • PARI
    {a(n) = 2*n*(n+1) + 1};
    
  • PARI
    x='x+O('x^200); Vec((1+x)^2/(1-x)^3) \\ Altug Alkan, Mar 23 2016
    
  • Python
    print([2*n*(n+1)+1 for n in range(48)]) # Michael S. Branicky, Jan 05 2021
  • Sage
    [i**2 + (i + 1)**2 for i in range(46)] # Zerinvary Lajos, Jun 27 2008
    

Formula

a(n) = 2*n^2 + 2*n + 1 = n^2 + (n+1)^2.
a(n) = 1 + 3 + 5 + ... + 2*n-1 + 2*n+1 + 2*n-1 + ... + 3 + 1. - Amarnath Murthy, May 28 2001
a(n) = 1/real(z(n+1)) where z(1)=i, (i^2=-1), z(k+1) = 1/(z(k)+2i). - Benoit Cloitre, Aug 06 2002
Nearest integer to 1/Sum_{k>n} 1/k^3. - Benoit Cloitre, Jun 12 2003
G.f.: (1+x)^2/(1-x)^3.
E.g.f.: exp(x)*(1+4x+2x^2).
a(n) = a(n-1) + 4n.
a(-n) = a(n-1).
a(n) = A064094(n+3, n) (fourth diagonal).
a(n) = 1 + Sum_{j=0..n} 4*j. - Xavier Acloque, Oct 08 2003
a(n) = A046092(n)+1 = (A016754(n)+1)/2. - Lekraj Beedassy, May 25 2004
a(n) = Sum_{k=0..n+1} (-1)^k*binomial(n, k)*Sum_{j=0..n-k+1} binomial(n-k+1, j)*j^2. - Paul Barry, Dec 22 2004
a(n) = ceiling((2n+1)^2/2). - Paul Barry, Jul 16 2006
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=1, a(1)=5, a(2)=13. - Jaume Oliver Lafont, Dec 02 2008
a(n)*a(n-1) = 4*n^4 + 1 for n > 0. - Reinhard Zumkeller, Feb 12 2009
Prefaced with a "1" (1, 1, 5, 13, 25, 41, ...): a(n) = 2*n*(n-1)+1. - Doug Bell, Feb 27 2009
a(n) = sqrt((A056220(n)^2 + A056220(n+1)^2) / 2). - Doug Bell, Mar 08 2009
a(n) = floor(2*(n+1)^3/(n+2)). - Gary Detlefs, May 20 2010
a(n) = A000330(n) - A000330(n-2). - Keith Tyler, Aug 10 2010
a(n) = A069894(n)/2. - J. M. Bergot, Jun 11 2012
a(n) = 2*a(n-1) - a(n-2) + 4. - Ant King, Jun 12 2012
Sum_{n>=0} 1/a(n) = (Pi/2)*tanh(Pi/2) = 1.4406595199775... = A228048. - Ant King, Jun 15 2012
a(n) = A209297(2*n+1,n+1). - Reinhard Zumkeller, Jan 19 2013
a(n)^3 = A048395(n)^2 + A048395(-n-1)^2. - Vincenzo Librandi, Jan 19 2013
a(n) = A000217(2n+1) - n. - Ivan N. Ianakiev, Nov 08 2013
a(n) = A251599(3*n+1). - Reinhard Zumkeller, Dec 13 2014
a(n) = A101321(4,n). - R. J. Mathar, Jul 28 2016
From Ilya Gutkovskiy, Jul 30 2016: (Start)
a(n) = Sum_{k=0..n} A008574(k).
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = exp(-1) = A068985. (End)
a(n) = 4 * A000217(n) + 1. - Bruce J. Nicholson, Jul 10 2017
a(n) = A002522(n) + A005563(n) = A002522(n+1) + A005563(n-1). - Bruce J. Nicholson, Aug 05 2017
Sum_{n>=0} a(n)/n! = 7*e. Sum_{n>=0} 1/a(n) = A228048. - Amiram Eldar, Jun 20 2020
a(n) = A000326(n+1) + A000217(n-1). - Charlie Marion, Nov 16 2020
a(n) = Integral_{x=0..2n+2} |1-x| dx. - Pedro Caceres, Dec 29 2020
From Amiram Eldar, Feb 17 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)*sech(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi*csch(Pi)*sinh(Pi/2). (End)
a(n) = A001651(n+1) + 1 - A028242(n). - Charlie Marion, Apr 05 2022
a(n) = A016754(n) - A046092(n). - Leo Tavares, Sep 16 2022
For n>0, a(n) = A101096(n+2) / 30. - Andy Nicol, Feb 06 2025
From Rémi Guillaume, Apr 21 2025: (Start)
a(n) = (2*A003215(n)+1)/3.
a(n) = (4*A005448(n+1)-1)/3.
a(n) + a(n-1) = A001845(n) - A001845(n-1), for n >= 1.
a(n) = (A005917(n+1))/(2n+1). (End)

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

A002162 Decimal expansion of the natural logarithm of 2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Newton calculated the first 16 terms of this sequence.
Area bounded by y = tan x, y = cot x, y = 0. - Clark Kimberling, Jun 26 2020
Choose four values independently and uniformly at random from the unit interval [0,1]. Sort them, and label them a,b,c,d from least to greatest (so that a b^2+c^2. - Akiva Weinberger, Dec 02 2024
Define the trihyperboloid to be the intersection of the three solid hyperboloids x^2+y^2-z^2<1, x^2-y^2+z^2<1, and -x^2+y^2+z^2<1. This fits perfectly within the cube [-1,1]^3. Then this is the ratio of the volume of the trihyperboloid to its bounding cube. - Akiva Weinberger, Dec 02 2024

Examples

			0.693147180559945309417232121458176568075500134360255254120680009493393...
		

References

  • G. Boros and V. H. Moll, Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals, Cambridge University Press, 2004.
  • Calvin C. Clawson, Mathematical Mysteries: The Beauty and Magic of Numbers, Springer, 2013. See p. 227.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 24, 250.
  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Sections 1.3.3, 2.21, 6.2, and 7.2.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 25 and appendix A, equations 25:14:3 and A:7:3 at pages 232, 670.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 29.

Crossrefs

Cf. A016730 (continued fraction), A002939, A008288, A142979, A142992.

Programs

  • Mathematica
    RealDigits[N[Log[2],200]][[1]] (* Vladimir Joseph Stephan Orlovsky, Feb 21 2011 *)
    RealDigits[Log[2],10,120][[1]] (* Harvey P. Dale, Jan 25 2024 *)
  • PARI
    { default(realprecision, 20080); x=10*log(2); for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b002162.txt", n, " ", d)); } \\ Harry J. Smith, Apr 21 2009

Formula

log(2) = Sum_{k>=1} 1/(k*2^k) = Sum_{j>=1} (-1)^(j+1)/j.
log(2) = Integral_{t=0..1} dt/(1+t).
log(2) = (2/3) * (1 + Sum_{k>=1} 2/((4*k)^3-4*k)) (Ramanujan).
log(2) = 4*Sum_{k>=0} (3-2*sqrt(2))^(2*k+1)/(2*k+1) (Y. Luke). - R. J. Mathar, Jul 13 2006
log(2) = 1 - (1/2)*Sum_{k>=1} 1/(k*(2*k+1)). - Jaume Oliver Lafont, Jan 06 2009, Jan 08 2009
log(2) = 4*Sum_{k>=0} 1/((4*k+1)*(4*k+2)*(4*k+3)). - Jaume Oliver Lafont, Jan 08 2009
log(2) = 7/12 + 24*Sum_{k>=1} 1/(A052787(k+4)*A000079(k)). - R. J. Mathar, Jan 23 2009
From Alexander R. Povolotsky, Jul 04 2009: (Start)
log(2) = (1/4)*(3 - Sum_{n>=1} 1/(n*(n+1)*(2*n+1))).
log(2) = (230166911/9240 - Sum_{k>=1} (1/2)^k*(11/k + 10/(k+1) + 9/(k+2) + 8/(k+3) + 7/(k+4) + 6/(k+5) - 6/(k+7) - 7/(k+8) - 8/(k+9) - 9/(k+10) - 10/(k+11)))/35917. (End)
log(2) = A052882/A000670. - Mats Granvik, Aug 10 2009
From log(1-x-x^2) at x=1/2, log(2) = (1/2)*Sum_{k>=1} L(k)/(k*2^k), where L(n) is the n-th Lucas number (A000032). - Jaume Oliver Lafont, Oct 24 2009
log(2) = Sum_{k>=1} 1/(cos(k*Pi/3)*k*2^k) (cf. A176900). - Jaume Oliver Lafont, Apr 29 2010
log(2) = (Sum_{n>=1} 1/(n^2*(n+1)^2*(2*n+1)) + 11)/16. - Alexander R. Povolotsky, Jan 13 2011
log(2) = ((Sum_{n>=1} (2*n+1)/(Sum_{k=1..n} k^2)^2)+396)/576. - Alexander R. Povolotsky, Jan 14 2011
From Alexander R. Povolotsky, Dec 16 2008: (Start)
log(2) = 105*(319/44100 - Sum_{n>=1} 1/(2*n*(2*n+1)*(2*n+3)*(2*n+5)*(2*n+7))).
log(2) = 319/420 - (3/2)*Sum_{n>=1} 1/(6*n^2+39*n+63). (End)
log(2) = Sum_{k>=1} A191907(2,k)/k. - Mats Granvik, Jun 19 2011
log(2) = Integral_{x=0..oo} 1/(1 + e^x) dx. - Jean-François Alcover, Mar 21 2013
log(2) = lim_{s->1} zeta(s)*(1-1/2^(s-1)). - Mats Granvik, Jun 18 2013
From Peter Bala, Dec 10 2013: (Start)
log(2) = 2*Sum_{n>=1} 1/( n*A008288(n-1,n-1)*A008288(n,n) ), a result due to Burnside.
log(2) = (1/3)*Sum_{n >= 0} (5*n+4)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(1/2)^n = (1/12)*Sum_{n >= 0} (28*n+17)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(-1/4)^n.
log(2) = (3/16)*Sum_{n >= 0} (14*n+11)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(1/4)^n = (1/12)*Sum_{n >= 0} (34*n+25)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(-1/18)^n. For more series of this type see the Bala link.
See A142979 for series acceleration formulas for log(2) obtained from the Mercator series log(2) = Sum_{n >= 1} (-1)^(n+1)/n. See A142992 for series for log(2) related to the root lattice C_n. (End)
log(2) = lim_{n->oo} Sum_{k=2^n..2^(n+1)-1} 1/k. - Richard R. Forberg, Aug 16 2014
From Peter Bala, Feb 03 2015: (Start)
log(2) = (2/3)*Sum_{k >= 0} 1/((2*k + 1)*9^k).
Define a pair of integer sequences A(n) = 9^n*(2*n + 1)!/n! and B(n) = A(n)*Sum_{k = 0..n} 1/((2*k + 1)*9^k). Both satisfy the same second-order recurrence equation u(n) = (40*n + 16)*u(n-1) - 36*(2*n - 1)^2*u(n-2). From this observation we obtain the continued fraction expansion log(2) = (2/3)*(1 + 2/(54 - 36*3^2/(96 - 36*5^2/(136 - ... - 36*(2*n - 1)^2/((40*n + 16) - ... ))))). Cf. A002391, A073000 and A105531 for similar expansions. (End)
log(2) = Sum_{n>=1} (Zeta(2*n)-1)/n. - Vaclav Kotesovec, Dec 11 2015
From Peter Bala, Oct 30 2016: (Start)
Asymptotic expansions:
for N even, log(2) - Sum_{k = 1..N/2} (-1)^(k-1)/k ~ (-1)^(N/2)*(1/N - 1/N^2 + 2/N^4 - 16/N^6 + 272/N^8 - ...), where the sequence of unsigned coefficients [1, 1, 2, 16, 272, ...] is A000182 with an extra initial term of 1. See Borwein et al., Theorem 1 (b);
for N odd, log(2) - Sum_{k = 1..(N-1)/2} (-1)^(k-1)/k ~ (-1)^((N-1)/2)*(1/N - 1/N^3 + 5/N^5 - 61/N^7 + 1385/N^9 - ...), by Borwein et al., Lemma 2 with f(x) := 1/(x + 1/2), h := 1/2 and then set x = (N - 1)/2, where the sequence of unsigned coefficients [1, 1, 5, 61, 1385, ...] is A000364. (End)
log(2) = lim_{n->oo} Sum_{k=1..n} sin(1/(n+k)). See Mathematical Reflections link. - Michel Marcus, Jan 07 2017
log(2) = Sum_{n>=1} A006519(n) / ((1 + 2^A006519(n)) * A000265(n) * (1 + A000265(n))). - Nicolas Nagel, Mar 19 2018
From Amiram Eldar, Jul 02 2020: (Start)
Equals Sum_{k>=2} zeta(k)/2^k.
Equals -Sum_{k>=2} log(1 - 1/k^2).
Equals Sum_{k>=1} 1/A002939(k).
Equals Integral_{x=0..Pi/3} tan(x) dx. (End)
log(2) = Integral_{x=0..Pi/2} (sec(x) - tan(x)) dx. - Clark Kimberling, Jul 08 2020
From Peter Bala, Nov 14 2020: (Start)
log(2) = Integral_{x = 0..1} (x - 1)/log(x) dx (Boros and Moll, p. 97).
log(2) = (1/2)*Integral_{x = 0..1} (x + 2)*(x - 1)^2/log(x)^2 dx.
log(2) = (1/4)*Integral_{x = 0..1} (x^2 + 3*x + 4)*(x - 1)^3/log(x)^3 dx. (End)
log(2) = 2*arcsinh(sqrt(2)/4) = 2*sqrt(2)*Sum_{n >= 0} (-1)^n*C(2*n,n)/ ((8*n+4)*32^n) = 3*Sum_{n >= 0} (-1)^n/((8*n+4)*(2^n)*C(2*n,n)). - Peter Bala, Jan 14 2022
log(2) = Integral_{x=0..oo} ( e^(-x) * (1-e^(-2x)) * (1-e^(-4x)) * (1-e^(-6x)) ) / ( x * (1-e^(-14x)) ) dx (see Crux Mathematicorum link). - Bernard Schott, Jul 11 2022
From Peter Bala, Oct 22 2023: (Start)
log(2) = 23/32 + 2!^3/16 * Sum_{n >= 1} (-1)^n * (n + 1)/(n*(n + 1)*(n + 2))^2 = 707/1024 - 4!^3/(16^2 * 2!^2) * Sum_{n >= 1} (-1)^n * (n + 2)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4))^2 = 42611/61440 + 6!^3/(16^3 * 3!^2) * Sum_{n >= 1} (-1)^n * (n + 3)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*(n + 6))^2.
More generally, it appears that for k >= 0, log(2) = c(k) + (2*k)!^3/(16^k * k!^2) * Sum_{n >= 1} (-1)^(n+k+1) * (n + k)/(n*(n + 1)*...*(n + 2*k))^2 , where c(k) is a rational approximation to log(2). The first few values of c(k) are [0, 23/32, 707/1024, 42611/61440, 38154331/55050240, 76317139/110100480, 26863086823/38755368960, ...].
Let P(n,k) = n*(n + 1)*...*(n + k).
Conjecture: for k >= 0 and r even with r - 1 <= k, the series Sum_{n >= 1} (-1)^n * (d/dn)^r (P(n,k)) / (P(n,k)^2 = A(r,k)*log(2) + B(r,k), where A(r,k) and B(r,k) are both rational numbers. (End)
From Peter Bala, Nov 13 2023: (Start)
log(2) = 5/8 + (1/8)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)^2 / ( k*(k + 1) )^4
= 257/384 + (3!^5/2^9)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)^2*(2*k + 5) / ( k*(k + 1)*(k + 2)*(k + 3) )^4
= 267515/393216 + (5!^5/2^19)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)*(2*k + 5)^2*(2*k + 7)*(2*k + 9) / ( k*(k + 1)*(k + 2)*(k + 3)*(k + 4)*(k + 5) )^4
log(2) = 3/4 - 1/128 * Sum_{k >= 0} (-1/16)^k * (10*k + 12)*binomial(2*k+2,k+1)/ ((k + 1)*(2*k + 3)). The terms of the series are O(1/(k^(3/2)*4^n)). (End)
log(2) = eta(1) is a period, where eta(x) is the Dirichlet eta function. - Andrea Pinos, Mar 19 2024
log(2) = K_{n>=0} (n^2 + [n=0])/1, where K is the Gauss notation for an infinite continued fraction. In the expanded form, log(2) = 1/(1 + 1/(1 + 4/(1 + 9/1 + 16/(1 + 25/(1 + ... (see Clawson at p. 227). - Stefano Spezia, Jul 01 2024
log(2) = lim_{n->oo} Sum_{k=1..n} 1/(n + k) = lim_{x->0} (2^x - 1)/x = lim_{x->0} (2^x - 2^(-x))/(2*x) (see Finch). - Stefano Spezia, Oct 19 2024
From Colin Linzer, Nov 08 2024: (Start)
log(2) = Integral_{t=0...oo} (1 - tanh(t)) dt.
log(2) = Integral_{t=0...1} arctanh(t) dt.
log(2) = (1/2) * Integral_{t=-1...1} |arctanh(t)| dt. (End)
log(2) = 1 + Sum_{n >= 1} (-1)^n/(n*(4*n^2 - 1)) = 1/2 + (1/2)*Sum_{n >= 1} 1/(n*(4*n^2 - 1)). - Peter Bala, Jan 07 2025
log(2) = Integral_{x=0..1} Integral_{y=0..1} 1/((1 - x*y)*(1 + x)*(1 + y)) dy dx. - Kritsada Moomuang, May 22 2025

A001850 Central Delannoy numbers: a(n) = Sum_{k=0..n} C(n,k)*C(n+k,k).

Original entry on oeis.org

1, 3, 13, 63, 321, 1683, 8989, 48639, 265729, 1462563, 8097453, 45046719, 251595969, 1409933619, 7923848253, 44642381823, 252055236609, 1425834724419, 8079317057869, 45849429914943, 260543813797441, 1482376214227923, 8443414161166173, 48141245001931263
Offset: 0

Views

Author

Keywords

Comments

Number of paths from (0,0) to (n,n) in an n X n grid using only steps north, northeast and east (i.e., steps (1,0), (1,1), and (0,1)).
Also the number of ways of aligning two sequences (e.g., of nucleotides or amino acids) of length n, with at most 2*n gaps (-) inserted, so that while unnecessary gappings: - -a a- - are forbidden, both b- and -b are allowed. (If only other of the latter is allowed, then the sequence A000984 gives the number of alignments.) There is an easy bijection from grid walks given by Dickau to such set of alignments (e.g., the straight diagonal corresponds to the perfect alignment with no gaps). - Antti Karttunen, Oct 10 2001
Also main diagonal of array A008288 defined by m(i,1) = m(1,j) = 1, m(i,j) = m(i-1,j-1) + m(i-1,j) + m(i,j-1). - Benoit Cloitre, May 03 2002
So, as a special case of Dmitry Zaitsev's Dec 10 2015 comment on A008288, a(n) is the number of points in Z^n that are L1 (Manhattan) distance <= n from any given point. These terms occur in the crystal ball sequences: a(n) here is the n-th term in the sequence for the n-dimensional cubic lattice. See A008288 for a list of crystal ball sequences (rows or columns of A008288). - Shel Kaphan, Dec 26 2022
a(n) is the number of n-matchings of a comb-like graph with 2*n teeth. Example: a(2) = 13 because the graph consisting of a horizontal path ABCD and the teeth Aa, Bb, Cc, Dd has 13 2-matchings: any of the six possible pairs of teeth and {Aa, BC}, {Aa, CD}, {Bb, CD}, {Cc, AB}, {Dd, AB}, {Dd, BC}, {AB, CD}. - Emeric Deutsch, Jul 02 2002
Number of ordered trees with 2*n+1 edges, having root of odd degree, nonroot nodes of outdegree at most 2 and branches of odd length. - Emeric Deutsch, Aug 02 2002
The sum of the first n coefficients of ((1 - x) / (1 - 2*x))^n is a(n-1). - Michael Somos, Sep 28 2003
Row sums of A063007 and A105870. - Paul Barry, Apr 23 2005
The Hankel transform (see A001906 for definition) of this sequence is A036442: 1, 4, 32, 512, 16384, ... . - Philippe Deléham, Jul 03 2005
Also number of paths from (0,0) to (n,0) using only steps U = (1,1), H = (1,0) and D =(1,-1), U can have 2 colors and H can have 3 colors. - N-E. Fahssi, Jan 27 2008
Equals row sums of triangle A152250 and INVERT transform of A109980: (1, 2, 8, 36, 172, 852, ...). - Gary W. Adamson, Nov 30 2008
Number of overpartitions in the n X n box (treat a walk of the type in the first comment as an overpartition, by interpreting a NE step as N, E with the part thus created being overlined). - William J. Keith, May 19 2017
Diagonal of rational functions 1/(1 - x - y - x*y), 1/(1 - x - y*z - x*y*z). - Gheorghe Coserea, Jul 03 2018
Dimensions of endomorphism algebras End(R^{(n)}) in the Delannoy category attached to the oligomorphic group of order preserving self-bijections of the real line. - Noah Snyder, Mar 22 2023
a(n) is the number of ways to tile a strip of length n with white squares, black squares, and red dominos, where we must have an equal number of white and black squares. - Greg Dresden and Leo Zhang, Jul 11 2025

Examples

			G.f. = 1 + 3*x + 13*x^2 + 63*x^3 + 321*x^4 + 1683*x^5 + 8989*x^6 + ...
		

References

  • Frits Beukers, Arithmetic properties of Picard-Fuchs equations, Séminaire de Théorie des nombres de Paris, 1982-83, Birkhäuser Boston, Inc.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 593.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
  • L. Moser and W. Zayachkowski, Lattice paths with diagonal steps, Scripta Math., 26 (1961), 223-229.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Wadsworth, Vol. 2, 1999; see Example 6.3.8 and Problem 6.49.
  • D. B. West, Combinatorial Mathematics, Cambridge, 2021, p. 28.

Crossrefs

Main diagonal of A064861.
Column k=2 of A262809 and A263159.

Programs

  • Maple
    seq(add(multinomial(n+k,n-k,k,k),k=0..n),n=0..20); # Zerinvary Lajos, Oct 18 2006
    seq(orthopoly[P](n,3), n=0..100); # Robert Israel, Nov 03 2015
  • Mathematica
    f[n_] := Sum[ Binomial[n, k] Binomial[n + k, k], {k, 0, n}]; Array[f, 21, 0] (* Or *)
    a[0] = 1; a[1] = 3; a[n_] := a[n] = (3(2 n - 1)a[n - 1] - (n - 1)a[n - 2])/n; Array[a, 21, 0] (* Or *)
    CoefficientList[ Series[1/Sqrt[1 - 6x + x^2], {x, 0, 20}], x] (* Robert G. Wilson v *)
    Table[LegendreP[n, 3], {n, 0, 22}] (* Jean-François Alcover, Jul 16 2012, from first formula *)
    a[n_] := Hypergeometric2F1[-n, n+1, 1, -1]; Table[a[n], {n, 0, 22}] (* Jean-François Alcover, Feb 26 2013 *)
    a[ n_] := With[ {m = If[n < 0, -1 - n, n]}, SeriesCoefficient[ (1 - 6 x + x^2)^(-1/2), {x, 0, m}]]; (* Michael Somos, Jun 10 2015 *)
  • Maxima
    a(n):=coeff(expand((1+3*x+2*x^2)^n),x,n);
    makelist(a(n),n,0,12); /* Emanuele Munarini, Mar 02 2011 */
    
  • PARI
    {a(n) = if( n<0, n = -1 - n); polcoeff( 1 / sqrt(1 - 6*x + x^2 + x * O(x^n)), n)}; /* Michael Somos, Sep 23 2006 */
    
  • PARI
    {a(n) = if( n<0, n = -1 - n); subst( pollegendre(n), x, 3)}; /* Michael Somos, Sep 23 2006 */
    
  • PARI
    {a(n) = if( n<0, n = -1 - n); n++; subst( Pol(((1 - x) / (1 - 2*x) + O(x^n))^n), x, 1);} /* Michael Somos, Sep 23 2006 */
    
  • PARI
    a(n)=if(n<0, 0, polcoeff((1+3*x+2*x^2)^n, n)) \\ Paul Barry, Aug 22 2007
    
  • PARI
    /* same as in A092566 but use */
    steps=[[1,0], [0,1], [1,1]]; /* Joerg Arndt, Jun 30 2011 */
    
  • PARI
    a(n)=sum(k=0,n,binomial(n,k)*binomial(n+k,k)); \\ Joerg Arndt, May 11 2013
    
  • PARI
    my(x='x+O('x^30)); Vec(1/sqrt(1 - 6*x + x^2)) \\ Altug Alkan, Oct 17 2015
    
  • Python
    # from Nick Hobson.
    def f(a, b):
        if a == 0 or b == 0:
            return 1
        return f(a, b - 1) + f(a - 1, b) + f(a - 1, b - 1)
    [f(n, n) for n in range(7)]
    
  • Python
    from gmpy2 import divexact
    A001850 = [1, 3]
    for n in range(2,10**3):
        A001850.append(divexact(A001850[-1]*(6*n-3)-(n-1)*A001850[-2],n))
    # Chai Wah Wu, Sep 01 2014
    
  • Sage
    a = lambda n: hypergeometric([-n, -n], [1], 2)
    [simplify(a(n)) for n in range(23)] # Peter Luschny, Nov 19 2014

Formula

a(n) = P_n(3), where P_n is n-th Legendre polynomial.
G.f.: 1 / sqrt(1 - 6*x + x^2).
a(n) = a(n-1) + 2*A002002(n) = Sum_{j} A063007(n, j). - Henry Bottomley, Jul 02 2001
Dominant term in asymptotic expansion is binomial(2*n, n)/2^(1/4)*((sqrt(2) + 1)/2)^(2*n + 1)*(1 + c_1/n + c_2/n^2 + ...). - Michael David Hirschhorn
a(n) = Sum_{i=0..n} (A000079(i)*A008459(n, i)) = Sum_{i=0..n} (2^i * C(n, i)^2). - Antti Karttunen, Oct 10 2001
a(n) = Sum_{k=0..n} C(n+k, n-k)*C(2*k, k). - Benoit Cloitre, Feb 13 2003
a(n) = Sum_{k=0..n} C(n, k)^2 * 2^k. - Michael Somos, Oct 08 2003
a(n - 1) = coefficient of x^n in A120588(x)^n if n>=0. - Michael Somos, Apr 11 2012
G.f. of a(n-1) = 1 / (1 - x / (1 - 2*x / (1 - 2*x / (1 - x / (1 - 2*x / (1 - x / ...)))))). - Michael Somos, May 11 2012
INVERT transform is A109980. BINOMIAL transform is A080609. BINOMIAL transform of A006139. PSUM transform is A089165. PSUMSIGN transform is A026933. First backward difference is A110170. - Michael Somos, May 11 2012
E.g.f.: exp(3*x)*BesselI(0, 2*sqrt(2)*x). - Vladeta Jovovic, Mar 21 2004
a(n) = Sum_{k=0..n} C(2*n-k, n)*C(n, k). - Paul Barry, Apr 23 2005
a(n) = Sum_{k>=n} binomial(k, n)^2/2^(k+1). - Vladeta Jovovic, Aug 25 2006
a(n) = a(-1 - n) for all n in Z. - Michael Somos, Sep 23 2006
D-finite with recurrence: a(-1) = a(0) = 1; n*a(n) = 3*(2*n-1)*a(n-1) - (n-1)*a(n-2). Eq (4) in T. D. Noe's article in JIS 9 (2006) #06.2.7.
Define general Delannoy numbers by (i,j > 0): d(i,0) = d(0,j) = 1 =: d(0,0) and d(i,j) = d(i-1,j-1) + d(i-2,j-1) + d(i-1,j). Then a(k) = Sum_{j >= 0} d(k,j)^2 + d(k-1,j)^2 = A026933(n)+A026933(n-1). This is a special case of the following formula for general Delannoy numbers: d(k,j) = Sum_{i >= 0, p=0..n} d(p, i) * d(n-p, j-i) + d(p-1, i) * d(n-p-1, j-i-1). - Peter E John, Oct 19 2006
Coefficient of x^n in (1 + 3*x + 2*x^2)^n. - N-E. Fahssi, Jan 11 2008
a(n) = A008288(A046092(n)). - Philippe Deléham, Apr 08 2009
G.f.: 1/(1 - x - 2*x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - ... (continued fraction). - Paul Barry, May 28 2009
G.f.: d/dx log(1/(1 - x*A001003(x))). - Vladimir Kruchinin, Apr 19 2011
G.f.: 1/(2*Q(0) + x - 1) where Q(k) = 1 + k*(1-x) - x - x*(k + 1)*(k + 2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Mar 14 2013
a(n) = Sum_{k=0..n} C(n,k) * C(n+k,k). - Joerg Arndt, May 11 2013
G.f.: G(0), where G(k) = 1 + x*(6 - x)*(4*k + 1)/(4*k + 2 - 2*x*(6-x)*(2*k + 1)*(4*k + 3)/(x*(6 - x)*(4*k + 3) + 4*(k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 22 2013
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - x*(6 - x)*(2*k - 1)/(x*(6 - x)*(2*k - 1) + 2*(k + 1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(6 - x)*(2*k + 1)/(x*(6 - x)*(2*k + 1) + 2*(k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jul 17 2013
a(n)^2 = Sum_{k=0..n} 2^k * C(2*k, k)^2 * C(n+k, n-k) = A243949(n). - Paul D. Hanna, Aug 17 2014
a(n) = hypergeom([-n, -n], [1], 2). - Peter Luschny, Nov 19 2014
a(n) = Sum_{k=0..n/2} C(n-k,k) * 3^(n-2*k) * 2^k * C(n,k). - Vladimir Kruchinin, Jun 29 2015
a(n) = A049600(n, n-1).
a(n) = Sum_{0 <= j, k <= n} (-1)^(n+j)*C(n,k)*C(n,j)*C(n+k,k)*C(n+k+j,k+j). Cf. A126086 and A274668. - Peter Bala, Jan 15 2020
a(n) ~ c * (3 + 2*sqrt(2))^n / sqrt(n), where c = 1/sqrt(4*Pi*(3*sqrt(2)-4)) = 0.572681... (Banderier and Schwer, 2005). - Amiram Eldar, Jun 07 2020
a(n+1) = 3*a(n) + 2*Sum_{l=1..n} A006318(l)*a(n-l). [Eq. (1.16) in Qi-Shi-Guo (2016)]
a(n) ~ (1 + sqrt(2))^(2*n+1) / (2^(5/4) * sqrt(Pi*n)). - Vaclav Kotesovec, Jan 09 2023
a(n-1) + a(n) = A241023(n) for n >= 1. - Peter Bala, Sep 18 2024
a(n) = Sum_{k=0..n} C(n+k, 2*k) * C(2*k, k). - Greg Dresden and Leo Zhang, Jul 11 2025

Extensions

New name and reference Sep 15 1995
Formula and more references from Don Knuth, May 15 1996

A180662 The Golden Triangle: T(n,k) = A001654(k) for n>=0 and 0<=k<=n.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 2, 6, 0, 1, 2, 6, 15, 0, 1, 2, 6, 15, 40, 0, 1, 2, 6, 15, 40, 104, 0, 1, 2, 6, 15, 40, 104, 273, 0, 1, 2, 6, 15, 40, 104, 273, 714, 0, 1, 2, 6, 15, 40, 104, 273, 714, 1870, 0, 1, 2, 6, 15, 40, 104, 273, 714, 1870, 4895
Offset: 0

Views

Author

Johannes W. Meijer, Sep 21 2010

Keywords

Comments

The terms in the n-th row of the Golden Triangle are the first (n+1) golden rectangle numbers. The golden rectangle numbers are A001654(n)=F(n)*F(n+1), with F(n) the Fibonacci numbers. The mirror image of the Golden Triangle is A180663.
We define below 24 mostly new triangle sums. The Row1 and Row2 sums are the ordinary and alternating row sums respectively and the Kn11 and Kn12 sums are commonly known as antidiagonal sums. Each of the names of these sums, except for the row sums, comes from a (fairy) chess piece that moves in its own peculiar way over a chessboard, see Hooper and Whyld. All pieces are leapers: knight (sqrt(5) or 1,2), fil (sqrt(8) or 2,2), camel (sqrt(10) or 3,1), giraffe (sqrt(17) or 4,1) and zebra (sqrt(13) or 3,2). Information about the origin of these chess sums can be found in "Famous numbers on a chessboard", see Meijer.
Each triangle or chess sum formula adds up numbers on a chessboard using the moves of its namesake. Converting a number triangle to a square array of numbers shows this most clearly (use the table button!). The formulas given below are for number triangles.
The chess sums of the Golden Triangle lead to six different sequences, see the crossrefs. As could be expected all these sums are related to the golden rectangle numbers.
Some triangles with complete sets of triangle sums are: A002260 (Natural Numbers), A007318 (Pascal), A008288 (Delannoy) A013609 (Pell-Jacobsthal), A036561 (Nicomachus), A104763 (Fibonacci(n)), A158405 (Odd Numbers) and of course A180662 (Golden Triangle).
#..Name....Type..Code....Definition of triangle sums.
1. Row......1....Row1.. a(n) = Sum_{k=0..n} T(n, k).
2. Row Alt..2....Row2.. a(n) = Sum_{k=0..n} (-1)^(n+k)*T(n, k).
3. Knight...1....Kn11.. a(n) = Sum_{k=0..floor(n/2)} T(n-k, k).
4. Knight...1....Kn12.. a(n) = Sum_{k=0..floor(n/2)} T(n-k+1, k+1).
5. Knight...1....Kn13.. a(n) = Sum_{k=0..floor(n/2)} T(n-k+2, k+2).
6. Knight...2....Kn21.. a(n) = Sum_{k=0..floor(n/2)} T(n-k, n-2*k).
7. Knight...2....Kn22.. a(n) = Sum_{k=0..floor(n/2)} T(n-k+1, n-2*k).
8. Knight...2....Kn23.. a(n) = Sum_{k=0..floor(n/2)} T(n-k+2, n-2*k).
9. Knight...3....Kn3... a(n) = Sum_{k=0..n} T(n+k, 2*k).
10. Knight...4....Kn4... a(n) = Sum_{k=0..n} T(n+k, n-k).
11. Fil......1....Fi1... a(n) = Sum_{k=0..floor(n/2)} T(n, 2*k).
12. Fil......2....Fi2... a(n) = Sum_{k=0..floor(n/2)} T(n, n-2*k).
13. Camel....1....Ca1... a(n) = Sum_{k=0..floor(n/3)} T(n-2*k, k).
14. Camel....2....Ca2... a(n) = Sum_{k=0..floor(n/3)} T(n-2*k, n-3*k).
15. Camel....3....Ca3... a(n) = Sum_{k=0..n} T(n+2*k, 3*k).
16. Camel....4....Ca4... a(n) = Sum_{k=0..n} T(n+2*k, n-k).
17. Giraffe..1....Gi1... a(n) = Sum_{k=0..floor(n/4)} T(n-3*k, k).
18. Giraffe..2....Gi2... a(n) = Sum_{k=0..floor(n/4)} T(n-3*k, n-4*k).
19. Giraffe..3....Gi3... a(n) = Sum_{k=0..n} T(n+3*k, 4*k).
20. Giraffe..4....Gi4... a(n) = Sum_{k=0..n} T(n+3*k, n-k).
21. Zebra....1....Ze1... a(n) = Sum_{k=0..floor(n/2)} T(n+k, 3*k).
22. Zebra....2....Ze2... a(n) = Sum_{k=0..floor(n/2)} T(n+k, n-2*k).
23. Zebra....3....Ze3... a(n) = Sum_{k=0..floor(n/3)} T(n-k, 2*k).
24. Zebra....4....Ze4... a(n) = Sum_{k=0..floor(n/3)} T(n-k, n-3*k).

Examples

			The first few rows of the Golden Triangle are:
  0;
  0, 1;
  0, 1, 2;
  0, 1, 2, 6;
  0, 1, 2, 6, 15;
  0, 1, 2, 6, 15, 40;
		

References

  • David Hooper and Kenneth Whyld, The Oxford Companion to Chess, p. 221, 1992.

Crossrefs

Cf. A180663 (Mirror), A001654 (Golden Rectangle), A000045 (F(n)).
Triangle sums: A064831 (Row1, Kn11, Kn12, Kn4, Ca1, Ca4, Gi1, Gi4), A077916 (Row2), A180664 (Kn13), A180665 (Kn21, Kn22, Kn23, Fi2, Ze2), A180665(2*n) (Kn3, Fi1, Ze3), A115730(n+1) (Ca2, Ze4), A115730(3*n+1) (Ca3, Ze1), A180666 (Gi2), A180666(4*n) (Gi3).

Programs

  • Haskell
    import Data.List (inits)
    a180662 n k = a180662_tabl !! n !! k
    a180662_row n = a180662_tabl !! n
    a180662_tabl = tail $ inits a001654_list
    -- Reinhard Zumkeller, Jun 08 2013
    
  • Magma
    [Fibonacci(k)*Fibonacci(k+1): k in [0..n], n in [0..12]]; // G. C. Greubel, May 25 2021
    
  • Maple
    F:= combinat[fibonacci]:
    T:= (n, k)-> F(k)*F(k+1):
    seq(seq(T(n, k), k=0..n), n=0..10); # revised Johannes W. Meijer, Sep 13 2012
  • Mathematica
    Table[Times @@ Fibonacci@ {k, k + 1}, {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Aug 18 2016 *)
    Module[{nn=20,f},f=Times@@@Partition[Fibonacci[Range[0,nn]],2,1];Table[Take[f,n],{n,nn}]]//Flatten (* Harvey P. Dale, Nov 26 2022 *)
  • PARI
    T(n,k)=fibonacci(k)*fibonacci(k+1) \\ Charles R Greathouse IV, Nov 07 2016
    
  • Sage
    flatten([[fibonacci(k)*fibonacci(k+1) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 25 2021

Formula

T(n, k) = F(k)*F(k+1) with F(n) = A000045(n), for n>=0 and 0<=k<=n.
From Johannes W. Meijer, Jun 22 2015: (Start)
Kn1p(n) = Sum_{k=0..floor(n/2)} T(n-k+p-1, k+p-1), p >= 1.
Kn1p(n) = Kn11(n+2*p-2) - Sum_{k=0..p-2} T(n-k+2*p-2, k), p >= 2.
Kn2p(n) = Sum_{k=0..floor(n/2)} T(n-k+p-1, n-2*k), p >= 1.
Kn2p(n) = Kn21(n+2*p-2) - Sum_{k=0..p-2} T(n+k+p, n+2*k+2), p >= 2. (End)
G.f. as triangle: xy/((1-x)(1+xy)(1-3xy+x^2 y^2)). - Robert Israel, Sep 06 2015
Previous Showing 11-20 of 141 results. Next