cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 15 results. Next

A000567 Octagonal numbers: n*(3*n-2). Also called star numbers.

Original entry on oeis.org

0, 1, 8, 21, 40, 65, 96, 133, 176, 225, 280, 341, 408, 481, 560, 645, 736, 833, 936, 1045, 1160, 1281, 1408, 1541, 1680, 1825, 1976, 2133, 2296, 2465, 2640, 2821, 3008, 3201, 3400, 3605, 3816, 4033, 4256, 4485, 4720, 4961, 5208, 5461
Offset: 0

Views

Author

Keywords

Comments

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,1,....
The spiral begins:
.
85--84--83--82--81--80
/ \
86 56--55--54--53--52 79
/ / \ \
87 57 33--32--31--30 51 78
/ / / \ \ \
88 58 34 16--15--14 29 50 77
/ / / / \ \ \ \
89 59 35 17 5---4 13 28 49 76
/ / / / / \ \ \ \ \
90 60 36 18 6 0 3 12 27 48 75
/ / / / / / / / / / /
91 61 37 19 7 1---2 11 26 47 74
\ \ \ \ \ . / / / /
92 62 38 20 8---9--10 25 46 73
\ \ \ \ . / / /
93 63 39 21--22--23--24 45 72
\ \ \ . / /
94 64 40--41--42--43--44 71
\ \ . /
95 65--66--67--68--69--70
\ .
96
.
(End)
From Lekraj Beedassy, Oct 02 2003: (Start)
Also the number of distinct three-cell blocks that may be removed out of A000217(n+1) square cells arranged in a stepping triangular array of side (n+1). A 5-layer triangular array of square cells, for instance, has vertices outlined thus:
x x
x x x
x x x x
x x x x x
x x x x x x
x x x x x x (End)
First derivative at n of A045991. - Ross La Haye, Oct 23 2004
Starting from n=1, the sequence corresponds to the Wiener index of K_{n,n} (the complete bipartite graph wherein each independent set has n vertices). - Kailasam Viswanathan Iyer, Mar 11 2009
Number of divisors of 24^(n-1) for n > 0 (cf A009968). - J. Lowell, Aug 30 2008
a(n) = A001399(6n-5), number of partitions of 6*n - 5 into parts < 4. For example a(2)=8 and partitions of 6*2 - 5 = 7 into parts < 4 are: [1,1,1,1,1,1,1], [1,1,1,1,1,2],[1,1,1,1,3], [1,1,1,2,2], [1,1,2,3], [1,2,2,2], [1,3,3], [2,2,3]. - Adi Dani, Jun 07 2011
Also, sequence found by reading the line from 0 in the direction 0, 8, ..., and the parallel line from 1 in the direction 1, 21, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. - Omar E. Pol, Sep 10 2011
Partial sums give A002414. - Omar E. Pol, Jan 12 2013
Generate a Pythagorean triple using Euclid's formula with (n, n-1) to give A,B,C. a(n) = B + (A + C)/2. - J. M. Bergot, Jul 13 2013
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n-4; {1, 2n-2, 3, 2n-2, 1, 18n-8}]. For n=1, this collapses to [5; {5, 10}]. - Magus K. Chu, Oct 10 2022
a(n)*a(n+1) + 1 = (3n^2 + n - 1)^2. In general, a(n)*a(n+k) + k^2 = (3n^2 + (3k-2)n - k)^2. - Charlie Marion, May 23 2023

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 38.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 1.
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 19-20.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 123.

Crossrefs

Cf. A014641, A014642, A014793, A014794, A001835, A016777, A045944, A093563 ((6, 1) Pascal, column m=2). A016921 (differences).
Cf. A005408 (the odd numbers).

Programs

  • GAP
    List([0..50], n -> n*(3*n-2)); # G. C. Greubel, Nov 15 2018
    
  • Haskell
    a000567 n = n * (3 * n - 2)  -- Reinhard Zumkeller, Dec 20 2012
    
  • Magma
    [n*(3*n-2) : n in [0..50]]; // Wesley Ivan Hurt, Oct 10 2021
  • Maple
    A000567 := proc(n)
        n*(3*n-2) ;
    end proc:
    seq(A000567(n), n=1..50) ;
  • Mathematica
    Table[n (3 n - 2), {n, 0, 50}] (* Harvey P. Dale, May 06 2012 *)
    Table[PolygonalNumber[RegularPolygon[8], n], {n, 0, 43}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    PolygonalNumber[8, Range[0, 20]] (* Eric W. Weisstein, Sep 07 2017 *)
    LinearRecurrence[{3, -3, 1}, {1, 8, 21}, {0, 20}] (* Eric W. Weisstein, Sep 07 2017 *)
  • PARI
    a(n)=n*(3*n-2) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    vector(50, n, n--; n*(3*n-2)) \\ G. C. Greubel, Nov 15 2018
    
  • Python
    # Intended to compute the initial segment of the sequence, not isolated terms.
    def aList():
         x, y = 1, 1
         yield 0
         while True:
             yield x
             x, y = x + y + 6, y + 6
    A000567 = aList()
    print([next(A000567) for i in range(49)]) # Peter Luschny, Aug 04 2019
    
  • Python
    [n*(3*n-2) for n in range(50)] # Gennady Eremin, Mar 10 2022
    
  • Sage
    [n*(3*n-2) for n in range(50)] # G. C. Greubel, Nov 15 2018
    

Formula

a(n) = n*(3*n-2).
a(n) = (3n-2)*(3n-1)*(3n)/((3n-1) + (3n-2) + (3n)), i.e., (the product of three consecutive numbers)/(their sum). a(1) = 1*2*3/(1+2+3), a(2) = 4*5*6/(4+5+6), etc. - Amarnath Murthy, Aug 29 2002
E.g.f.: exp(x)*(x+3*x^2). - Paul Barry, Jul 23 2003
G.f.: x*(1+5*x)/(1-x)^3. Simon Plouffe in his 1992 dissertation
a(n) = Sum_{k=1..n} (5*n - 4*k). - Paul Barry, Sep 06 2005
a(n) = n + 6*A000217(n-1). - Floor van Lamoen, Oct 14 2005
a(n) = C(n+1,2) + 5*C(n,2).
Starting (1, 8, 21, 40, 65, ...) = binomial transform of [1, 7, 6, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=0, a(1)=1, a(2)=8. - Jaume Oliver Lafont, Dec 02 2008
a(n) = A000578(n) - A007531(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = a(n-1) + 6*n - 5 (with a(0)=0). - Vincenzo Librandi, Nov 20 2010
a(n) = 2*a(n-1) - a(n-2) + 6. - Ant King, Sep 01 2011
a(n) = A000217(n) + 5*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = (A185212(n) - 1) / 4. - Reinhard Zumkeller, Dec 20 2012
a(n) = A174709(6n). - Philippe Deléham, Mar 26 2013
a(n) = (2*n-1)^2 - (n-1)^2. - Ivan N. Ianakiev, Apr 10 2013
a(6*a(n) + 16*n + 1) = a(6*a(n) + 16*n) + a(6*n + 1). - Vladimir Shevelev, Jan 24 2014
a(0) = 0, a(n) = Sum_{k=0..n-1} A005408(A051162(n-1,k)), n >= 1. - L. Edson Jeffery, Jul 28 2014
Sum_{n>=1} 1/a(n) = (sqrt(3)*Pi + 9*log(3))/12 = 1.2774090575596367311949534921... . - Vaclav Kotesovec, Apr 27 2016
From Ilya Gutkovskiy, Jul 29 2016: (Start)
Inverse binomial transform of A084857.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(2*sqrt(3)) = A093766. (End)
a(n) = n * A016777(n-1) = A053755(n) - A000290(n+1). - Bruce J. Nicholson, Aug 10 2017
Product_{n>=2} (1 - 1/a(n)) = 3/4. - Amiram Eldar, Jan 21 2021
P(4k+4,n) = ((k+1)*n - k)^2 - (k*n - k)^2 where P(m,n) is the n-th m-gonal number (a generalization of the Apr 10 2013 formula, a(n) = (2*n-1)^2 - (n-1)^2). - Charlie Marion, Oct 07 2021
From Leo Tavares, Oct 31 2021: (Start)
a(n) = A000290(n) + 4*A000217(n-1). See Square Rays illustration.
a(n) = A000290(n) + A046092(n-1)
a(n) = A000384(n) + 2*A000217(n-1). See Twin Rectangular Rays illustration.
a(n) = A000384(n) + A002378(n-1)
a(n) = A003154(n) - A045944(n-1). See Star Rows illustration. (End)

Extensions

Incorrect example removed by Joerg Arndt, Mar 11 2010

A016921 a(n) = 6*n + 1.

Original entry on oeis.org

1, 7, 13, 19, 25, 31, 37, 43, 49, 55, 61, 67, 73, 79, 85, 91, 97, 103, 109, 115, 121, 127, 133, 139, 145, 151, 157, 163, 169, 175, 181, 187, 193, 199, 205, 211, 217, 223, 229, 235, 241, 247, 253, 259, 265, 271, 277, 283, 289, 295, 301, 307, 313, 319, 325, 331
Offset: 0

Views

Author

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 22 ).
Also solutions to 2^x + 3^x == 5 (mod 7). - Cino Hilliard, May 10 2003
Except for 1, exponents n > 1 such that x^n - x^2 - 1 is reducible. - N. J. A. Sloane, Jul 19 2005
Let M(n) be the n X n matrix m(i,j) = min(i,j); then the trace of M(n)^(-2) is a(n-1) = 6*n - 5. - Benoit Cloitre, Feb 09 2006
If Y is a 3-subset of an (2n+1)-set X then, for n >= 3, a(n-1) is the number of 3-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 16 2007
All composite terms belong to A269345 as shown in there. - Waldemar Puszkarz, Apr 13 2016
First differences of the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
For b(n) = A103221(n) one has b(a(n)-1) = b(a(n)+1) = b(a(n)+2) = b(a(n)+3) = b(a(n)+4) = n+1 but b(a(n)) = n. So-called "dips" in A103221. See the Avner and Gross remark on p. 178. - Wolfdieter Lang, Sep 16 2016
A (n+1,n) pebbling move involves removing n + 1 pebbles from a vertex in a simple graph and placing n pebbles on an adjacent vertex. A two-player impartial (n+1,n) pebbling game involves two players alternating (n+1,n) pebbling moves. The first player unable to make a move loses. The sequence a(n) is also the minimum number of pebbles such that any assignment of those pebbles on a complete graph with 3 vertices is a next-player winning game in the two player impartial (k+1,k) pebbling game. These games are represented by A347637(3,n). - Joe Miller, Oct 18 2021
Interleaving of A017533 and A017605. - Leo Tavares, Nov 16 2021

Examples

			From _Ilya Gutkovskiy_, Apr 15 2016: (Start)
Illustration of initial terms:
                      o
                    o o o
              o     o o o
            o o o   o o o
      o     o o o   o o o
    o o o   o o o   o o o
o   o o o   o o o   o o o
n=0  n=1     n=2     n=3
(End)
		

References

  • Avner Ash and Robert Gross, Summing it up, Princeton University Press, 2016, p. 178.

Crossrefs

Cf. A093563 ((6, 1) Pascal, column m=1).
a(n) = A007310(2*(n+1)); complement of A016969 with respect to A007310.
Cf. A287326 (second column).

Programs

Formula

a(n) = 6*n + 1, n >= 0 (see the name).
G.f.: (1+5*x)/(1-x)^2.
A008615(a(n)) = n. - Reinhard Zumkeller, Feb 27 2008
A157176(a(n)) = A013730(n). - Reinhard Zumkeller, Feb 24 2009
a(n) = 4*(3*n-1) - a(n-1) (with a(0)=1). - Vincenzo Librandi, Nov 20 2010
E.g.f.: (1 + 6*x)*exp(x). - G. C. Greubel, Sep 18 2019
a(n) = A003215(n) - 6*A000217(n-1). See Hexagonal Lines illustration. - Leo Tavares, Sep 10 2021
From Leo Tavares, Oct 27 2021: (Start)
a(n) = 6*A001477(n-1) + 7
a(n) = A016813(n) + 2*A001477(n)
a(n) = A017605(n-1) + A008588(n-1)
a(n) = A016933(n) - 1
a(n) = A008588(n) + 1. (End)
Sum_{n>=0} (-1)^n/a(n) = Pi/6 + sqrt(3)*arccoth(sqrt(3))/3. - Amiram Eldar, Dec 10 2021

A029653 Numbers in (2,1)-Pascal triangle (by row).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Reverse of A029635. Row sums are A003945. Diagonal sums are Fibonacci(n+2) = Sum_{k=0..floor(n/2)} (2n-3k)*C(n-k,n-2k)/(n-k). - Paul Barry, Jan 30 2005
Riordan array ((1+x)/(1-x), x/(1-x)). The signed triangle (-1)^(n-k)T(n,k) or ((1-x)/(1+x), x/(1+x)) is the inverse of A055248. Row sums are A003945. Diagonal sums are F(n+2). - Paul Barry, Feb 03 2005
Row sums = A003945: (1, 3, 6, 12, 24, 48, 96, ...) = (1, 3, 7, 15, 31, 63, 127, ...) - (0, 0, 1, 3, 7, 15, 31, ...); where (1, 3, 7, 15, ...) = A000225. - Gary W. Adamson, Apr 22 2007
Triangle T(n,k), read by rows, given by (2,-1,0,0,0,0,0,0,0,...) DELTA (1,0,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 17 2011
A029653 is jointly generated with A208510 as an array of coefficients of polynomials v(n,x): initially, u(1,x)=v(1,x)=1; for n>1, u(n,x)=u(n-1,x)+x*v(n-1)x and v(n,x)=u(n-1,x)+x*v(n-1,x)+1. See the Mathematica section. - Clark Kimberling, Feb 28 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle, see A228196. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle, see A228576. - Boris Putievskiy, Sep 04 2013
The n-th row polynomial is (2 + x)*(1 + x)^(n-1) for n >= 1. More generally, the n-th row polynomial of the Riordan array ( (1-a*x)/(1-b*x), x/(1-b*x) ) is (b - a + x)*(b + x)^(n-1) for n >= 1. - Peter Bala, Feb 25 2018

Examples

			The triangle T(n,k) begins:
n\k 0  1  2   3   4   5   6   7  8  9 10 ...
0:  1
1:  2  1
2:  2  3  1
3:  2  5  4   1
4:  2  7  9   5   1
5:  2  9 16  14   6   1
6:  2 11 25  30  20   7   1
7:  2 13 36  55  50  27   8   1
8:  2 15 49  91 105  77  35   9  1
9:  2 17 64 140 196 182 112  44 10  1
10: 2 19 81 204 336 378 294 156 54 11  1
... Reformatted. - _Wolfdieter Lang_, Jan 09 2015
With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins
/1        \/1         \/1        \      /1        \
|2 1      ||0 1       ||0 1      |      |2 1      |
|2 1 1    ||0 2 1     ||0 0 1    |... = |2 3 1    |
|2 1 1 1  ||0 2 1 1   ||0 0 2 1  |      |2 5 4 1  |
|2 1 1 1 1||0 2 1 1 1 ||0 0 2 1 1|      |2 7 9 5 1|
|...      ||...       ||...      |      |...      |
- _Peter Bala_, Dec 27 2014
		

References

  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.

Crossrefs

(d, 1) Pascal triangles: A007318(d=1), A093560(3), A093561(4), A093562(5), A093563(6), A093564(7), A093565(8), A093644(9), A093645(10).

Programs

  • Haskell
    a029653 n k = a029653_tabl !! n !! k
    a029653_row n = a029653_tabl !! n
    a029653_tabl = [1] : iterate
                   (\xs -> zipWith (+) ([0] ++ xs) (xs ++ [0])) [2, 1]
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Maple
    A029653 :=  proc(n,k)
    if n = 0 then
      1;
    else
      binomial(n-1, k)+binomial(n, k)
    fi
    end proc: # R. J. Mathar, Jun 30 2013
  • 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 *)
    (* Clark Kimberling, Feb 28 2012 *)
  • 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(v(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 27 2017
    
  • Python
    from math import comb, isqrt
    def A029653(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*((r<<1)-a)//r if n else 1 # Chai Wah Wu, Nov 12 2024

Formula

T(n, k) = C(n-2, k-1) + C(n-2, k) + C(n-1, k-1) + C(n-1, k) except for n=0.
G.f.: (1 + x + y + xy)/(1 - y - xy). - Ralf Stephan, May 17 2004
T(n, k) = (2n-k)*binomial(n, n-k)/n, n, k > 0. - Paul Barry, Jan 30 2005
Sum_{k=0..n} T(n, k)*x^k gives A003945-A003954 for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10. - Philippe Deléham, Jul 10 2005
T(n, k) = C(n-1, k) + C(n, k). - Philippe Deléham, Jul 10 2005
Equals A097806 * A007318, i.e., the pairwise operator * Pascal's Triangle as infinite lower triangular matrices. - Gary W. Adamson, Apr 22 2007
From Peter Bala, Dec 27 2014: (Start)
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(2 + 5*x + 4*x^2/2! + x^3/3!) = 2 + 7*x + 16*x^2/2! + 30*x^3/3! + 50*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the lower unit triangular array with 1's on the main diagonal and 1's everywhere else below the main diagonal except for the first column which consists of the sequence [1,2,2,2,...]. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle equals the infinite product M(0)*M(1)*M(2)*... (which is clearly well-defined). See the Example section. (End)

Extensions

More terms from James Sellers

A002414 Octagonal pyramidal numbers: a(n) = n*(n+1)*(2*n-1)/2.

Original entry on oeis.org

1, 9, 30, 70, 135, 231, 364, 540, 765, 1045, 1386, 1794, 2275, 2835, 3480, 4216, 5049, 5985, 7030, 8190, 9471, 10879, 12420, 14100, 15925, 17901, 20034, 22330, 24795, 27435, 30256, 33264, 36465, 39865, 43470, 47286, 51319, 55575, 60060, 64780
Offset: 1

Views

Author

Keywords

Comments

Number of ways of covering a 2n X 2n lattice with 2n^2 dominoes with exactly 2 horizontal dominoes. - Yong Kong (ykong@curagen.com), May 06 2000
Equals binomial transform of [0, 1, 7, 6, 0, 0, 0, ...]. - Gary W. Adamson, Jun 14 2008, corrected Oct 25 2012
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF3 denominators of A156927. See A157704 for background information. - Johannes W. Meijer, Mar 07 2009
This sequence is related to A000326 by a(n) = n*A000326(n) - Sum_{i=0..n-1} A000326(i) and this is the case d=3 in the identity n*(n*(d*n-d+2)/2)-Sum_{k=0..n-1} k*(d*k-d+2)/2 = n*(n+1)*(2*d*n-2*d+3)/6. - Bruno Berselli, Apr 21 2010
2*a(n) gives the principal diagonal of the convolution array A213819. - Clark Kimberling, Jul 04 2012
Partial sums of the figurate octagonal numbers A000567. For each sequence with a linear recurrence with constant coefficients, the values reduced modulo some constant m generate a periodic sequence. For this sequence, these Pisano periods have length 1, 4, 3, 8, 5, 12, 7, 16, 9, 20, 11, 24, 13, 28, 15, 32, 17, ... for m >= 1. - Ant King, Oct 26 2012
Partial sums of the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
On a square grid of side length n+1, the number of embedded rectangles (where each side is greater than 1). For example, in a 2 X 2 square there is one rectangle, in a 3 X 3 square there are nine rectangles, etc. - Peter Woodward, Nov 26 2017
a(n) is the sum of the numbers in the n X n square array A204154(n). - Ali Sada, Jun 21 2019
Sum of all multiples of n and n+1 that are <= n^2. - Wesley Ivan Hurt, May 25 2023

Examples

			a(2) = 9 since there are 9 ways to cover a 4 X 4 lattice with 8 dominoes, 2 of which is horizontal and the other 6 are vertical. - Yong Kong (ykong@curagen.com), May 06 2000
G.f. = x + 9*x^2 + 30*x^3 + 70*x^4 + 135*x^5 + 231*x^6 + 364*x^7 + 540*x^8 + 765*x^9 + ...
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 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).

Crossrefs

Cf. A093563 (( 6, 1) Pascal, column m=3). A000567 (differences).
Cf. A156927, A157704. - Johannes W. Meijer, Mar 07 2009
Cf. A000326.
Cf. similar sequences listed in A237616.
Cf. A260234 (largest prime factor of a(n+1)).

Programs

  • Magma
    [n*(n+1)*(2*n-1)/2: n in [1..50]]; // Vincenzo Librandi, May 24 2016
  • Maple
    A002414 := n-> 1/2*n*(n+1)*(2*n-1): seq(A002414(n), n=1..100);
  • Mathematica
    LinearRecurrence[{4,-6,4,-1},{1,9,30,70},40] (* Harvey P. Dale, Apr 12 2013 *)
  • PARI
    {a(n) = (2*n - 1) * n * (n + 1) / 2} \\ Michael Somos, Mar 17 2011
    

Formula

a(n) = odd numbers * triangular numbers = (2*n-1)* binomial(n+1,2). - Xavier Acloque, Oct 27 2003
G.f.: x*(1+5*x)/(1-x)^4. [Conjectured by Simon Plouffe in his 1992 dissertation.]
a(n) = A000578(n) + A000217(n-1). - Kieren MacMillan, Mar 19 2007
a(-n) = -A160378(n). - Michael Somos, Mar 17 2011
From Ant King, Oct 26 2012: (Start)
a(n) = a(n-1) + n*(3*n-2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) = n*A000326(n) - A002411(n-1), see Berselli's comment.
a(n) = (n+1)*(2*A000567(n)+n)/6.
a(n) = A000292(n) + 5*A000292(n-1) = binomial(n+2,3)+5*binomial(n+1,3).
a(n) = A002413(n) + A000292(n-1).
a(n) = A000217(n) + 6*A000292(n-1).
Sum_{n>=1} 1/a(n) = 2*(4*log(2)-1)/3 = 1.1817258148265...
(End)
a(n) = Sum_{i=0..n-1} (n-i)*(6*i+1), with a(0)=0. - Bruno Berselli, Feb 10 2014
E.g.f.: x*(2 + 7*x + 2*x^2)*exp(x)/2. - Ilya Gutkovskiy, May 23 2016
a(n) = A080851(6,n-1). - R. J. Mathar, Jul 28 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(Pi + 1 - 4*log(2))/3. - Amiram Eldar, Jul 02 2020

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), May 09 2000
Incorrect formula deleted by Ant King, Oct 04 2012

A005009 a(n) = 7*2^n.

Original entry on oeis.org

7, 14, 28, 56, 112, 224, 448, 896, 1792, 3584, 7168, 14336, 28672, 57344, 114688, 229376, 458752, 917504, 1835008, 3670016, 7340032, 14680064, 29360128, 58720256, 117440512, 234881024, 469762048, 939524096, 1879048192, 3758096384
Offset: 0

Views

Author

Keywords

Comments

The first differences are the sequence itself. - Alexandre Wajnberg & Eric Angelini, Sep 07 2005

Crossrefs

Sequences of the form (2*m+1)*2^n: A000079 (m=0), A007283 (m=1), A020714 (m=2), this sequence (m=3), A005010 (m=4), A005015 (m=5), A005029 (m=6), A110286 (m=7), A110287 (m=8), A110288 (m=9), A175805 (m=10), A248646 (m=11), A164161 (m=12), A175806 (m=13), A257548 (m=15).
Row sums of (6, 1)-Pascal triangle A093563 and of (1, 6)-Pascal triangle A096956, n>=1.

Programs

Formula

G.f.: 7/(1-2*x).
a(n) = A118416(n+1,4) for n > 3. - Reinhard Zumkeller, Apr 27 2006
a(n) = 2*a(n-1), for n > 0, with a(0)=7 . - Philippe Deléham, Nov 23 2008
a(n) = 7 * A000079(n). - Omar E. Pol, Dec 16 2008
a(n) = A173787(n+3,n). - Reinhard Zumkeller, Feb 28 2010
Intersection of A014311 and A212191: all terms and their squares are the sum of exactly three distinct powers of 2, A000120(a(n)) = A000120(a(n)^2) = 3. - Reinhard Zumkeller, May 03 2012
G.f.: 2/x/G(0) - 1/x + 9, where G(k)= 1 + 1/(1 - x*(7*k+2)/(x*(7*k+9) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 03 2013
E.g.f.: 7*exp(2*x). - Stefano Spezia, May 15 2021

A228196 A triangle formed like Pascal's triangle, but with n^2 on the left border and 2^n on the right border instead of 1.

Original entry on oeis.org

0, 1, 2, 4, 3, 4, 9, 7, 7, 8, 16, 16, 14, 15, 16, 25, 32, 30, 29, 31, 32, 36, 57, 62, 59, 60, 63, 64, 49, 93, 119, 121, 119, 123, 127, 128, 64, 142, 212, 240, 240, 242, 250, 255, 256, 81, 206, 354, 452, 480, 482, 492, 505, 511, 512, 100, 287, 560, 806, 932, 962, 974, 997, 1016, 1023, 1024
Offset: 1

Views

Author

Boris Putievskiy, Aug 15 2013

Keywords

Comments

The third row is (n^4 - n^2 + 24*n + 24)/12.
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013

Examples

			The start of the sequence as a triangular array read by rows:
   0;
   1,  2;
   4,  3,  4;
   9,  7,  7,  8;
  16, 16, 14, 15, 16;
  25, 32, 30, 29, 31, 32;
  36, 57, 62, 59, 60, 63, 64;
		

Crossrefs

Cf. We denote Pascal-like triangle with L(n) on the left border and R(n) on the right border by (L(n),R(n)). A007318 (1,1), A008949 (1,2^n), A029600 (2,3), A029618 (3,2), A029635 (1,2), A029653 (2,1), A037027 (Fibonacci(n),1), A051601 (n,n) n>=0, A051597 (n,n) n>0, A051666 (n^2,n^2), A071919 (1,0), A074829 (Fibonacci(n), Fibonacci(n)), A074909 (1,n), A093560 (3,1), A093561 (4,1), A093562 (5,1), A093563 (6,1), A093564 (7,1), A093565 (8,1), A093644 (9,1), A093645 (10,1), A095660 (1,3), A095666 (1,4), A096940 (1,5), A096956 (1,6), A106516 (3^n,1), A108561(1,(-1)^n), A132200 (4,4), A134636 (2n+1,2n+1), A137688 (2^n,2^n), A160760 (3^(n-1),1), A164844(1,10^n), A164847 (100^n,1), A164855 (101*100^n,1), A164866 (101^n,1), A172171 (1,9), A172185 (9,11), A172283 (-9,11), A177954 (int(n/2),1), A193820 (1,2^n), A214292 (n,-n), A227074 (4^n,4^n), A227075 (3^n,3^n), A227076 (5^n,5^n), A227550 (n!,n!), A228053 ((-1)^n,(-1)^n), A228074 (Fibonacci(n), n).
Cf. A000290 (row 1), A153056 (row 2), A000079 (column 1), A000225 (column 2), A132753 (column 3), A118885 (row sums of triangle array + 1), A228576 (generalized Pascal's triangle).

Programs

  • GAP
    T:= function(n,k)
        if k=0 then return n^2;
        elif k=n then return 2^n;
        else return T(n-1,k-1) + T(n-1,k);
        fi;
      end;
    Flat(List([0..12], n-> List([0..n], k-> T(n,k) ))); # G. C. Greubel, Nov 12 2019
  • Maple
    T:= proc(n, k) option remember;
          if k=0 then n^2
        elif k=n then 2^k
        else T(n-1, k-1) + T(n-1, k)
          fi
        end:
    seq(seq(T(n, k), k=0..n), n=0..10); # G. C. Greubel, Nov 12 2019
  • Mathematica
    T[n_, k_]:= T[n, k] = If[k==0, n^2, If[k==n, 2^k, T[n-1, k-1] + T[n-1, k]]]; Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 12 2019 *)
    Flatten[Table[Sum[i^2 Binomial[n-1-i, n-k-i], {i,1,n-k}] + Sum[2^i Binomial[n-1-i, k-i], {i,1,k}], {n,0,10}, {k,0,n}]] (* Greg Dresden, Aug 06 2022 *)
  • PARI
    T(n,k) = if(k==0, n^2, if(k==n, 2^k, T(n-1, k-1) + T(n-1, k) )); \\ G. C. Greubel, Nov 12 2019
    
  • Python
    def funcL(n):
       q = n**2
       return q
    def funcR(n):
       q = 2**n
       return q
    for n in range (1,9871):
       t=int((math.sqrt(8*n-7) - 1)/ 2)
       i=n-t*(t+1)/2-1
       j=(t*t+3*t+4)/2-n-1
       sum1=0
       sum2=0
       for m1 in range (1,i+1):
          sum1=sum1+funcR(m1)*binomial(i+j-m1-1,i-m1)
       for m2 in range (1,j+1):
          sum2=sum2+funcL(m2)*binomial(i+j-m2-1,j-m2)
       sum=sum1+sum2
    
  • Sage
    @CachedFunction
    def T(n, k):
        if (k==0): return n^2
        elif (k==n): return 2^n
        else: return T(n-1, k-1) + T(n-1, k)
    [[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 12 2019
    

Formula

T(n,0) = n^2, n>0; T(0,k) = 2^k; T(n, k) = T(n-1, k-1) + T(n-1, k) for n,k > 0. [corrected by G. C. Greubel, Nov 12 2019]
Closed-form formula for general case. Let L(m) and R(m) be the left border and the right border of Pascal like triangle, respectively. We denote binomial(n,k) by C(n,k).
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} R(m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} L(m2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} R(m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} L(m2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2); n>0.
Some special cases. If L(m)={b,b,b...} b*A000012, then the second sum takes form b*C(n+k-1,j). If L(m) is {0,b,2b,...} b*A001477, then the second sum takes form b*C(n+k,n-1). Similarly for R(m) and the first sum.
For this sequence L(m)=m^2 and R(m)=2^m.
As table read by antidiagonals T(n,k) = Sum_{m1=1..n} (2^m1)*C(n+k-m1-1, n-m1) + Sum_{m2=1..k} (m2^2)*C(n+k-m2-1, k-m2); n,k >=0.
As linear sequence a(n) = Sum_{m1=1..i} (2^m1)*C(i+j-m1-1, i-m1) + Sum_{m2=1..j} (m2^2)*C(i+j-m2-1, j-m2), where i=n-t*(t+1)/2-1, j=(t*t+3*t+4)/2-n-1, t=floor((-1+sqrt(8*n-7))/2).
As a triangular array read by rows, T(n,k) = Sum_{i=1..n-k} i^2*C(n-1-i, n-k-i) + Sum_{i=1..k} 2^i*C(n-1-i, k-i); n,k >=0. - Greg Dresden, Aug 06 2022

Extensions

Cross-references corrected and extended by Philippe Deléham, Dec 27 2013

A002419 4-dimensional figurate numbers: a(n) = (6*n-2)*binomial(n+2,3)/4.

Original entry on oeis.org

1, 10, 40, 110, 245, 476, 840, 1380, 2145, 3190, 4576, 6370, 8645, 11480, 14960, 19176, 24225, 30210, 37240, 45430, 54901, 65780, 78200, 92300, 108225, 126126, 146160, 168490, 193285, 220720, 250976, 284240, 320705, 360570, 404040, 451326, 502645, 558220
Offset: 1

Views

Author

Keywords

Comments

a(n) is the n-th antidiagonal sum of the convolution array A213761. - Clark Kimberling, Jul 04 2012
Convolution of A000027 with A000567 (excluding 0). - Bruno Berselli, Dec 07 2012
a(n) = the sum of all the ways of adding the k-tuples of A016777(0) to A016777(n-1). For n=4, the terms are 1,4,7,10 giving (1)+(4)+(7)+(10)=22; (1+4)+(4+7)+(7+10)=33; (1+4+7)+(4+7+10)=33; (1+4+7+10)=22; adding 22+33+33+22=110. - J. M. Bergot, Jun 26 2017
Also the number of chordless cycles in the (n+2)-crown graph. - Eric W. Weisstein, Jan 02 2018

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • 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. A093563 ((6, 1) Pascal, column m=4).
Cf. A000027, A000567, A002414 (first differences), A016777, A080852, A213761.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.

Programs

  • GAP
    List([1..40], n-> n*(n+1)*(n+2)*(3*n-1)/12); # G. C. Greubel, Jul 03 2019
  • Magma
    /* A000027 convolved with A000567 (excluding 0): */ A000567:=func; [&+[(n-i+1)*A000567(i): i in [1..n]]: n in [1..40]]; // Bruno Berselli, Dec 07 2012
    
  • Mathematica
    CoefficientList[Series[(1+5*x)/(1-x)^5, {x,0,40}], x] (* Vincenzo Librandi, Jun 20 2013 *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {1, 10, 40, 110, 245}, 40] (* Harvey P. Dale, Nov 30 2014 *)
    Table[n(n+1)(n+2)(3n-1)/12, {n, 40}] (* Eric W. Weisstein, Jan 02 2018 *)
    Table[Sum[2 x + 3 x^2 - 2 y, {x, 0, g}, {y, x, g}], {g, 1, 20}] (* Horst H. Manninger, Jun 20 2025 *)
  • PARI
    a(n)=(3*n-1)*binomial(n+2,3)/2 \\ Charles R Greathouse IV, Sep 24 2015
    
  • Python
    A002419_list, m = [], [6, 1, 1, 1, 1]
    for _ in range(10**2):
        A002419_list.append(m[-1])
        for i in range(4):
            m[i+1] += m[i] # Chai Wah Wu, Jan 24 2016
    
  • Sage
    [n*(n+1)*(n+2)*(3*n-1)/12 for n in (1..40)] # G. C. Greubel, Jul 03 2019
    

Formula

a(n) = (3*n-1)*binomial(n+2, 3)/2.
G.f.: x*(1+5*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation.
Sum_{n>=1} 1/a(n) = (-24+81*log(3) -9*Pi*sqrt(3))/14 = 1.143929... - R. J. Mathar, Mar 29 2011
a(n) = (3*n^4 + 8*n^3 + 3*n^2 - 2*n)/12. - Chai Wah Wu, Jan 24 2016
a(n) = A080852(6,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: x*(12 + 48*x + 26*x^2 + 3*x^3)*exp(x)/12. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*(3*sqrt(3)*Pi - 32*log(2) + 8)/7. - Amiram Eldar, Feb 11 2022

A022096 Fibonacci sequence beginning 1, 6.

Original entry on oeis.org

1, 6, 7, 13, 20, 33, 53, 86, 139, 225, 364, 589, 953, 1542, 2495, 4037, 6532, 10569, 17101, 27670, 44771, 72441, 117212, 189653, 306865, 496518, 803383, 1299901, 2103284, 3403185, 5506469, 8909654, 14416123, 23325777, 37741900, 61067677, 98809577, 159877254
Offset: 0

Views

Author

Keywords

Comments

a(n-1) = Sum_{k=0..ceiling((n-1)/2)} P(6;n-1-k,k), n>=1, with a(-1)=5. These are the sums of the SW-NE diagonals in P(6;n,k), the (6,1) Pascal triangle A093563. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs. Also sums of SW-NE diagonals in (1,5)-Pascal triangle A096940.
Subsequence of primes: 7, 13, 53, 139, 953, 44771, 189653, 1494692464747, ... - R. J. Mathar, Aug 09 2012
a(n) is the sum of seven consecutive Fibonacci numbers. a(n) = F(n-4) + F(n-3) + F(n-2) + F(n-1) + F(n) + F(n+1) + F(n+2), where F(n)=A000045(n), extended so that F(-1)=1, F(-2)=-1, F(-3)=2, and F(-4)=-3. - Graeme McRae, Apr 24 2014

Crossrefs

Programs

Formula

a(n) = a(n-1) + a(n-2), n>=2, a(0)=1, a(1)=6.
G.f.: (1+5*x)/(1-x-x^2).
a(n) = A109754(5, n+1).
a(n) = 5*Fibonacci(n+2) - 4*Fibonacci(n+1). - Gary Detlefs, Dec 21 2010
a(n) = (2^(-1-n)*((1 - sqrt(5))^n*(-11 + sqrt(5)) + (1 + sqrt(5))^n*(11 + sqrt(5))))/sqrt(5). - Herbert Kociemba, Dec 18 2011
a(n) = Fibonacci(n+3) - Fibonacci(n-4). - Greg Dresden and Sam Neale, Mar 08 2022

Extensions

Spelling correction by Jason G. Wurtzel, Aug 22 2010

A093562 (5,1) Pascal triangle.

Original entry on oeis.org

1, 5, 1, 5, 6, 1, 5, 11, 7, 1, 5, 16, 18, 8, 1, 5, 21, 34, 26, 9, 1, 5, 26, 55, 60, 35, 10, 1, 5, 31, 81, 115, 95, 45, 11, 1, 5, 36, 112, 196, 210, 140, 56, 12, 1, 5, 41, 148, 308, 406, 350, 196, 68, 13, 1, 5, 46, 189, 456, 714, 756, 546, 264, 81, 14, 1, 5, 51, 235, 645, 1170
Offset: 0

Views

Author

Wolfdieter Lang, Apr 22 2004

Keywords

Comments

This is the fifth member, d=5, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-1, for d=1..4.
This is an example of a Riordan triangle (see A093560 for a comment and A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group). Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=(1+4*z)/(1-(1+x)*z).
The SW-NE diagonals give A022095(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 4. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
The array F(5;n,m) gives in the columns m >= 1 the figurate numbers based on A016861, including the heptagonal numbers A000566 (see the W. Lang link).
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
The n-th row polynomial is (4 + x)*(1 + x)^(n-1) for n >= 1. More generally, the n-th row polynomial of the Riordan array ( (1-a*x)/(1-b*x), x/(1-b*x) ) is (b - a + x)*(b + x)^(n-1) for n >= 1. - Peter Bala, Mar 02 2018

Examples

			Triangle begins
  [1];
  [5,  1];
  [5,  6,  1];
  [5, 11,  7,  1];
  ...
		

References

  • Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
  • Ivo Schneider, Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.

Crossrefs

Cf. Row sums: A007283(n-1), n>=1, 1 for n=0. A082505(n+1), alternating row sums are 1 for n=0, 4 for n=2 and 0 else.
Column sequences give for m=1..9: A016861, A000566 (heptagonal), A002413, A002418, A027800, A051946, A050484, A052255, A055844.

Programs

  • Haskell
    a093562 n k = a093562_tabl !! n !! k
    a093562_row n = a093562_tabl !! n
    a093562_tabl = [1] : iterate
                   (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [5, 1]
    -- Reinhard Zumkeller, Aug 31 2014
    
  • Python
    from math import comb, isqrt
    def A093562(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*(r+(r-a<<2))//r if n else 1 # Chai Wah Wu, Nov 12 2024

Formula

a(n, m) = F(5;n-m, m) for 0<= m <= n, otherwise 0, with F(5;0, 0)=1, F(5;n, 0)=5 if n>=1 and F(5;n, m):=(5*n+m)*binomial(n+m-1, m-1)/m if m>=1.
G.f. column m (without leading zeros): (1+4*x)/(1-x)^(m+1), m>=0.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=5 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
T(n, k) = C(n, k) + 4*C(n-1, k). - Philippe Deléham, Aug 28 2005
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(5 + 11*x + 7*x^2/2! + x^3/3!) = 5 + 16*x + 34*x^2/2! + 60*x^3/3! + 95*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014

A027810 a(n) = (n+1)*binomial(n+5, 5).

Original entry on oeis.org

1, 12, 63, 224, 630, 1512, 3234, 6336, 11583, 20020, 33033, 52416, 80444, 119952, 174420, 248064, 345933, 474012, 639331, 850080, 1115730, 1447160, 1856790, 2358720, 2968875, 3705156, 4587597, 5638528, 6882744, 8347680, 10063592, 12063744, 14384601, 17066028
Offset: 0

Views

Author

Thi Ngoc Dinh (via R. K. Guy)

Keywords

Comments

Number of 11-subsequences of [ 1, n ] with just 5 contiguous pairs.

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.
  • Herbert John Ryser, Combinatorial Mathematics, Carus Mathematical Monographs No. 14, John Wiley and Sons, 1963, pp. 1-8.

Crossrefs

Partial sums of A051843.
Cf. A093563 ((6, 1) Pascal, column m=6).

Programs

  • Haskell
    a027810 n = (n + 1) * a007318' (n + 5) 5
    -- Reinhard Zumkeller, Aug 31 2014
    
  • Magma
    [(n+1)*Binomial(n+5,5): n in [0..40]]; // Vincenzo Librandi, Jul 30 2014
    
  • Magma
    [n*(n-1)*(n-2)*(n-3)*(n-4)^2/120: n in [5..40]]; // Vincenzo Librandi, Jul 30 2014
    
  • Maple
    [seq(n*(n-1)*(n-2)*(n-3)*(n-4)^2/5!,n=5..33)]; # Zerinvary Lajos, Oct 19 2006
  • Mathematica
    Table[(n+1)Binomial[n+5,5],{n,0,30}] (* Harvey P. Dale, Jul 29 2014 *)
    CoefficientList[Series[(1 + 5 x)/(1 - x)^7, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
  • PARI
    a(n)=n*(n^5+16*n^4+100*n^3+310*n^2+499*n+394)/120+1 \\ Charles R Greathouse IV, Sep 28 2015

Formula

G.f.: (1+5*x)/(1-x)^7.
a(n) = A245334(n+5, 5)/A000142(5). - Reinhard Zumkeller, Aug 31 2014
From Amiram Eldar, Jan 28 2022: (Start)
Sum_{n>=0} 1/a(n) = 5*Pi^2/6 - 1025/144.
Sum_{n>=0} (-1)^n/a(n) = 5*Pi^2/12 - 160*log(2)/3 + 4865/144. (End)

Extensions

Two redundant formulas deleted by N. J. A. Sloane, Jul 30 2014
Showing 1-10 of 15 results. Next