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 16 results. Next

A167552 A triangle related to the a(n) formulas of the rows of the ED1 array A167546.

Original entry on oeis.org

1, 3, -2, 5, -5, 2, 7, -7, 14, -8, 9, -6, 63, -66, 24, 11, 0, 209, -264, 308, -144, 13, 13, 559, -689, 2236, -2132, 720, 15, 35, 1281, -1255, 11640, -14980, 14064, -5760, 17, 68, 2618, -1360, 47753, -68068, 145452, -126480, 40320
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009, Nov 23 2009

Keywords

Comments

The a(n) formulas given below correspond to the first ten rows of the ED1 array A167546.
The recurrence relations of the a(n) formulas for the left hand triangle columns, see the cross-references below, lead to the sequences A003148 and A007318.

Examples

			Row 1: a(n) = 1.
Row 2: a(n) = 3*n - 2.
Row 3: a(n) = 5*n^2 - 5*n + 2.
Row 4: a(n) = 7*n^3 - 7*n^2 + 14*n - 8.
Row 5: a(n) = 9*n^4 - 6*n^3 + 63*n^2 - 66*n + 24.
Row 6: a(n) = 11*n^5 + 0*n^4 + 209*n^3 - 264*n^2 + 308*n - 144.
Row 7: a(n) = 13*n^6 +13*n^5 +559*n^4 -689*n^3 +2236*n^2 -2132*n +720.
Row 8: a(n) = 15*n^7 + 35*n^6 + 1281*n^5 - 1255*n^4 + 11640*n^3 - 14980*n^2 + 14064*n - 5760.
Row 9: a(n) = 17*n^8 + 68*n^7 + 2618*n^6 - 1360*n^5 + 47753*n^4 - 68068*n^3 + 145452*n^2 - 126480*n + 40320.
Row 10: a(n) = 19*n^9 + 114*n^8 + 4902*n^7 + 684*n^6 + 163419*n^5 - 224694*n^4 + 1048268*n^3 - 1308264*n^2 + 1081632*n - 403200.
		

Crossrefs

A167546 is the ED1 array.
A000012, A016777, 2*A005891, A167547, A167548 and A167549 are the first sixth ED1 array rows.
A098557 and A167553 equal the first two right hand columns of this triangle.
A005408, A167554 and A167555, A168302 and A168303 equal the first five left hand columns of this triangle.
A000142 equals the row sums.
Cf. A003148 and A007318.

A167556 A triangle related to the GF(z) formulas of the rows of the ED1 array A167546.

Original entry on oeis.org

1, 1, 2, 2, 6, 2, 6, 24, 4, 8, 24, 120, 0, 48, 24, 120, 720, -120, 384, 72, 144, 720, 5040, -1680, 3696, -432, 1296, 720, 5040, 40320, -20160, 40320, -15840, 17280, 2880, 5760, 40320, 362880, -241920, 483840, -311040, 288000, -46080, 69120, 40320
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

The GF(z) formulas given below correspond to the first ten rows of the ED1 array A167546. The polynomials in their numerators lead to the triangle given above.

Examples

			Row 1: GF(z) = 1/(1-z).
Row 2: GF(z) = (1 + 2*z)/(1-z)^2.
Row 3: GF(z) = (2 + 6*z + 2*z^2)/(1-z)^3.
Row 4: GF(z) = (6 + 24*z + 4*z^2 + 8*z^3)/(1-z)^4.
Row 5: GF(z) = (24 + 120*z + 0*z^2 + 48*z^3 + 24*z^4)/(1-z)^5.
Row 6: GF(z) = (120 + 720*z - 120*z^2 + 384*z^3 + 72*z^4 + 144*z^5)/ (1-z)^6.
Row 7: GF(z) = (720 + 5040*z - 1680*z^2 + 3696*z^3 - 432*z^4 + 1296*z^5 + 720*z^6)/(1-z)^7.
Row 8: GF(z) = (5040 + 40320*z - 20160*z^2 + 40320*z^3 - 15840*z^4 + 17280*z^5 + 2880*z^6 + 5760*z^7)/(1-z)^8.
Row 9: GF(z) = (40320 +362880*z -241920*z^2 + 483840*z^3 - 311040*z^4 + 288000*z^5 - 46080*z^6 + 69120*z^7 + 40320*z^8)/(1-z)^9.
Row 10: GF(z) = (362880 +3628800*z -3024000*z^2 +6289920*z^3 -5495040*z^4 + 5276160*z^5 - 2131200*z^6 + 1382400*z^7 + 201600*z^8 + 403200*z^9)/(1-z)^10;
		

Crossrefs

A167546 is the ED1 array.
A000142, A000142 (n=>2) and 120*A062148 (with three extra terms at the beginning of the sequence) equal the first three left hand triangle columns.
A098557(n) and A098557(n)*A064455(n) equal the first two right hand triangle columns.
A007680 equals the row sums.

A167557 The lower left triangle of the ED1 array A167546.

Original entry on oeis.org

1, 1, 4, 2, 12, 32, 6, 48, 160, 384, 24, 240, 960, 2688, 6144, 120, 1440, 6720, 21504, 55296, 122880, 720, 10080, 53760, 193536, 552960, 1351680, 2949120, 5040, 80640, 483840, 1935360, 6082560, 16220160, 38338560, 82575360
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

We discovered that the numbers that appear in the lower left triangle of the ED1 array A167546 (m <= n) behave in a regular way, see the formula below. This rather simple regularity doesn't show up in the upper right triangle of the ED1 array (m > n).

Examples

			The first few triangle rows are:
[1]
[1, 4]
[2, 12, 32]
[6, 48, 160, 384]
[24, 240, 960, 2688, 6144]
[120, 1440, 6720, 21504, 55296, 122880]
		

Crossrefs

A167546 is the ED1 array.
A047053 and A167558 are the first two right hand triangle columns.
A000142, 4*A001710 (n>=2), 32*A001720, 384*A001730, 6144*A049389, 122880*A051431 are the first six left hand triangle columns.
A167559 equals the row sums.

Programs

  • Maple
    a := proc(n, m): 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)! end: seq(seq(a(n, m), m=1..n), n=1..8);  # Johannes W. Meijer, revised Nov 23 2012
  • Mathematica
    Flatten[Table[(4^(m-1) (m-1)!(n+m-2)!)/(2m-2)!,{n,10},{m,n}]] (* Harvey P. Dale, Sep 29 2013 *)

Formula

a(n,m) = 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)!.

A167547 The fourth row of the ED1 array A167546.

Original entry on oeis.org

6, 48, 160, 384, 762, 1336, 2148, 3240, 4654, 6432, 8616, 11248, 14370, 18024, 22252, 27096, 32598, 38800, 45744, 53472, 62026, 71448, 81780, 93064, 105342, 118656, 133048, 148560, 165234, 183112, 202236, 222648, 244390, 267504, 292032, 318016
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Crossrefs

Cf. A167546.

Programs

  • Mathematica
    Table[7*n^3 - 7*n^2 + 14*n - 8, {n, 1, 100}] (* or *) LinearRecurrence[{4,-6,4,-1},{6, 48, 160, 384}, 100] (* G. C. Greubel, Jun 15 2016 *)

Formula

a(n) = 7*n^3 - 7*n^2 + 14*n - 8.
G.f.: (6+24*z+4*z^2+8*z^3)/(1-z)^4.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - G. C. Greubel, Jun 15 2016

A167548 The fifth row of the ED1 array A167546.

Original entry on oeis.org

24, 240, 960, 2688, 6144, 12264, 22200, 37320, 59208, 89664, 130704, 184560, 253680, 340728, 448584, 580344, 739320, 929040, 1153248, 1415904, 1721184, 2073480, 2477400, 2937768, 3459624, 4048224, 4709040, 5447760, 6270288
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Crossrefs

Cf. A167546.

Programs

  • Mathematica
    Table[9*n^4 - 6*n^3 + 63*n^2 - 66*n + 24, {n, 1, 100}] (* or *) LinearRecurrence[{5,-10,10,-5,1}, {24, 240, 960, 2688, 6144}, 100] (* G. C. Greubel, Jun 15 2016 *)

Formula

a(n) = 9*n^4 - 6*n^3 + 63*n^2 - 66*n + 24.
G.f.: (24+120*z+0*z^2+48*z^3+24*z^4)/(1-z)^5.
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - G. C. Greubel, Jun 15 2016

A167549 The sixth row of the ED1 array A167546.

Original entry on oeis.org

120, 1440, 6720, 21504, 55296, 122880, 245640, 452880, 783144, 1285536, 2021040, 3063840, 4502640, 6441984, 9003576, 12327600, 16574040, 21924000, 28581024, 36772416, 46750560, 58794240, 73209960, 90333264, 110530056
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Crossrefs

Cf. A167546.

Programs

  • Mathematica
    Table[11*n^5 + 209*n^3 - 264*n^2 + 308*n - 144, {n, 1, 100}] (* or *) LinearRecurrence[{6,-15,20,-15,6,-1},{120, 1440, 6720, 21504, 55296, 122880}, 100] (* G. C. Greubel, Jun 15 2016 *)

Formula

a(n) = 11*n^5 + 0*n^4 + 209*n^3 - 264*n^2 + 308*n - 144.
G.f.: (120+720*z-120*z^2+384*z^3+72*z^4+144*z^5)/(1-z)^6.
a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6). - G. C. Greubel, Jun 15 2016

A167551 The row sums of the ED1 array A167546 read by antidiagonals.

Original entry on oeis.org

1, 2, 7, 26, 115, 596, 3623, 25462, 203939, 1836032, 18361831, 201984098, 2423820227, 31509694204, 441135811847, 6617037458606, 105872600222083, 1799834206605272, 32397015728269319, 615543298868687962
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Comments

The first six terms coincide with the known first six terms of A006373.

Crossrefs

Cf. A006373.
A167546 is the ED1 array.

A167550 The a(n,n+1) diagonal of the ED1 array A167546.

Original entry on oeis.org

1, 7, 62, 762, 12264, 245640, 5897520, 165145680, 5284782720, 190253266560, 7610141548800, 334846347897600, 16072626136089600, 835776577757721600, 46803488615967283200, 2808209320881060096000, 179725396599156215808000
Offset: 1

Views

Author

Johannes W. Meijer, Nov 10 2009

Keywords

Crossrefs

Equals the a(n, n+1) diagonal of the ED1 array A167546.

Formula

a(n,n+1) = 4*(n-1)*a(n-1,n)+3*(n-1)! with a(1,2) = 1.

A005891 Centered pentagonal numbers: (5n^2+5n+2)/2; crystal ball sequence for 3.3.3.4.4. planar net.

Original entry on oeis.org

1, 6, 16, 31, 51, 76, 106, 141, 181, 226, 276, 331, 391, 456, 526, 601, 681, 766, 856, 951, 1051, 1156, 1266, 1381, 1501, 1626, 1756, 1891, 2031, 2176, 2326, 2481, 2641, 2806, 2976, 3151, 3331, 3516, 3706, 3901, 4101, 4306, 4516, 4731, 4951, 5176, 5406
Offset: 0

Views

Author

Keywords

Comments

Equals the triangular numbers convolved with [1, 3, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
From Ant King, Jun 15 2012: (Start)
a(n) == 1 (mod 5) for all n.
The digital roots of the a(n) form a purely periodic palindromic 9-cycle 1, 6, 7, 4, 6, 4, 7, 6, 1.
The units' digits of the a(n) form a purely periodic palindromic 4-cycle 1, 6, 6, 1.
(End)
Binomial transform of (1, 5, 5, 0, 0, 0, ...) and second partial sum of (1, 4, 5, 5, 5, ...). - Gary W. Adamson, Sep 09 2015
a(n) = a(-1-n) for all n in Z. - Michael Somos, Jan 25 2019
On the plane start with a single regular pentagon, and repeat the following procedure, "For each edge of any pentagon not already connected to an existing pentagon create a mirror image such that the mirror image does not overlap with an existing pentagon." a(n) is the number of pentagons occupying the plane after n repetitions. - Torlach Rush, Sep 14 2022

Examples

			a(2)= 5*T(2) + 1 = 5*3 + 1 = 16, a(4) = 5*T(4) + 1 = 5*10 + 1 = 51. - _Thomas M. Green_, Nov 16 2009
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.

Crossrefs

Cf. A028895, A001844, A003215, A004068 (partial sums), A006322, A001263.
Partial sums of A008706.
Equals second row of A167546 divided by 2.

Programs

Formula

G.f.: (1 + 3*x + x^2)/(1 - x)^3. Simon Plouffe in his 1992 dissertation
Narayana transform (A001263) of [1, 5, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=1, a(1)=6, a(2)=16. - Jaume Oliver Lafont, Dec 02 2008
a(n) = 5*A000217(n) + 1 = 5*T(n) + 1, for n = 0, 1, 2, 3, ... and where T(n) = n*(n+1)/2 = n-th triangular number. - Thomas M. Green, Nov 25 2009
a(n) = a(n-1) + 5*n, with a(0)=1. - Vincenzo Librandi, Nov 18 2010
a(n) = A028895(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = 2*a(n-1) - a(n-2) + 5. - Ant King, Jun 12 2012
Sum_{n>=0} 1/a(n) = 2*Pi /sqrt(15) *tanh(Pi/2*sqrt(3/5)) = 1.360613169863... - Ant King, Jun 15 2012
a(n) = A101321(5,n). - R. J. Mathar, Jul 28 2016
E.g.f.: (2 + 10*x + 5*x^2)*exp(x)/2. - Ilya Gutkovskiy, Jul 28 2016
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 17*e/2.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 3/(2*e). (End)

Extensions

Formula corrected and relocated by Johannes W. Meijer, Nov 07 2009

A007680 a(n) = (2n+1)*n!.

Original entry on oeis.org

1, 3, 10, 42, 216, 1320, 9360, 75600, 685440, 6894720, 76204800, 918086400, 11975040000, 168129561600, 2528170444800, 40537905408000, 690452066304000, 12449059983360000, 236887827111936000, 4744158915944448000, 99748982335242240000
Offset: 0

Views

Author

Keywords

Comments

Denominators in series for sqrt(Pi/4)*erf(x): sqrt(Pi/4)*erf(x)= x/1 - x^3/3 + x^5/10 - x^7/42 + x^9/216 -+ ...
Appears to be the BinomialMean transform of A000354 (after truncating the first term of A000354). (See A075271 for the definition of BinomialMean.) - John W. Layman, Apr 16 2003
Number of permutations p of {1,2,...,n+2} such that max|p(i)-i|=n+1. Example: a(1)=3 since only the permutations 312,231 and 321 of {1,2,3} satisfy the given condition. - Emeric Deutsch, Jun 04 2003
Stirling transform of A000670(n+1) = [3, 13, 75, 541, ...] is a(n) = [3, 10, 42, 216, ...]. - Michael Somos, Mar 04 2004
Stirling transform of a(n) = [2, 10, 42, 216, ...] is A052875(n+1) = [2, 12, 74, ...]. - Michael Somos, Mar 04 2004
A related sequence also arises in evaluating indefinite integrals of sec(x)^(2k+1), k=0, 1, 2, ... Letting u = sec(x) and d = sqrt(u^2-1), one obtains a(0) = log(u+d) 2*k*a(k) = (2*k-1)*u^(2*k-1)*d + a(k-1). Viewing these as polynomials in u gives 2^k*k!*a(k) = a(0) + d*Sum(i=0..k-1){ (2*i+1)*i!*2^i*u^(2*i+1) }, which is easily proved by induction. Apart from the power of 2, which could be incorporated into the definition of u (or by looking at erf(ix/2)/ i (i=sqrt(-1)), the sum's coefficients form our series and are the reciprocals of the power series terms for -sqrt(-Pi/4)*erf(ix/2)). This yields a direct but somewhat mysterious relationship between the power series of erf(x) and integrals involving sec(x). - William A. Huber (whuber(AT)quantdec.com), Mar 14 2002
When written in factoradic ("factorial base"), this sequence from a(1) onwards gives the smallest number containing two adjacent digits, increasing when read from left to right, whose difference is n-1. - Christian Perfect, May 03 2016
a(n-1)^2 is the number of permutations p of [1..2n] such that Sum_{i=1..2n} abs(p(i)-i) = 2n^2-2. - Fang Lixing, Dec 07 2018
A standard series for the calculation of coordinates on a clothoid (also called cornuspiral):
x = s*(a(0) - (tau^2/a(2)) + (tau^4/a(4)) - (tau^6/a(6)) + ...)
y = s*((tau/a(1)) + (tau^3/a(3)) - (tau^5/a(5)) + ...).
s is the arclength from the clothoids origin to the desired point p(x,y). The tangent at the clothoids origin intersects with the tangent at the point p(x,y) with an angle of tau. - Thomas Scheuerle, Oct 13 2021
a(n) = P_n(1) where P_n(x) is the Pidduck polynomials. - Michael Somos, May 27 2023

Examples

			G.f. = 1 + 3*x + 10*x^2 + 42*x^3 + 216*x^4 + 1320*x^5 + 9360*x^6 + ... - _Michael Somos_, Jan 01 2019
		

References

  • H. W. Gould, A class of binomial sums and a series transform, Utilitas Math., 45 (1994), 71-83.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • N. Wirth, Systematisches Programmieren, 1975, exercise 9.3

Crossrefs

From Johannes W. Meijer, Nov 12 2009: (Start)
Appears in A167546.
Equals the rows sums of A167556.
(End)

Programs

  • GAP
    a:=List([0..20],n->(2*n+1)*Factorial(n));; Print(a); # Muniru A Asiru, Jan 01 2019
  • Magma
    [(2*n+1)*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Aug 20 2011
    
  • Maple
    [(2*n+1)*factorial(n)$n=0..20]; # Muniru A Asiru, Jan 01 2019
  • Mathematica
    Table[(2n + 1)*n!, {n, 0, 20}] (* Stefan Steinerberger, Apr 08 2006 *)
  • PARI
    {a(n) = if( n<0, 0, (2*n+1) * n!)}; /* Michael Somos, Mar 04 2004 */
    

Formula

E.g.f.: (1+x)/(1-x)^2.
This is the binomial mean transform of A000354 (after truncating the first term). See Spivey and Steil (2006). - Michael Z. Spivey (mspivey(AT)ups.edu), Feb 26 2006
E.g.f.: (of aerated sequence) 1+x^2/2+sqrt(pi)*(x+x^3/4)*exp(x^2/4)*ERF(x/2). - Paul Barry, Apr 11 2010
G.f.: 1 + x*G(0), where G(k)= 1 + x*(k+1)/(1 - (k+2)/(k+2 + (k+1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 08 2013
a(n-2) = (A208528(n)+A208529(n))/2, for n>=2. - Luis Manuel Rivera Martínez, Mar 05 2014
D-finite with recurrence: (-2*n+1)*a(n) +n*(2*n+1)*a(n-1)=0. - R. J. Mathar, Jan 27 2020
Sum_{n>=0} 1/a(n) = sqrt(Pi)*erfi(1)/2 = A019704 * A099288 = A347910. - Amiram Eldar, Oct 07 2020
Sum_{n>=0} (-1)^n/a(n) = A347909 . - R. J. Mathar, Sep 30 2021
Showing 1-10 of 16 results. Next