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

A162007 Third left hand column of the EG1 triangle A162005.

Original entry on oeis.org

1, 270, 36096, 4766048, 704357760, 120536980224, 24060789342208, 5590122715250688, 1503080384197754880, 464520829174515630080, 163839204411117787938816, 65500849343294249018327040
Offset: 3

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Third left hand column of the EG1 triangle A162005.
Other left hand columns are A000182 and A162006.
Related to A094665, A083061 and A156919.
A000079, A036289 and A100381 appear in the a(n, 3) formula.
A001789, A003472, A054849, A002409, A054851, A140325 and A140354 (scaled by 2^(m-1)) appear one by one in the a(n, m) formulas for m= 4 and higher .

Programs

  • Maple
    nmax := 14; mmax := nmax: imax := nmax: T1(0, x) := 1: T1(0, x+1) := 1: for i from 1 to imax do T1(i, x) := expand((2*x+1)*(x+1)*T1(i-1, x+1) - 2*x^2*T1(i-1, x)): dx := degree(T1(i, x)): for k from 0 to dx do c(k) := coeff(T1(i, x), x, k) od: T1(i, x+1) := sum(c(j1)*(x+1)^(j1), j1 = 0..dx): od: for i from 0 to imax do for j from 0 to i do A083061(i, j) := coeff(T1(i, x), x, j) od: od: for n from 0 to nmax do for k from 0 to n do A094665(n+1, k+1) := A083061(n, k) od: od: A094665(0, 0) := 1: for n from 1 to nmax do A094665(n, 0) := 0 od: for m from 1 to mmax do A156919(0, m) := 0 end do: for n from 0 to nmax do A156919(n, 0) := 2^n end do: for n from 1 to nmax do for m from 1 to mmax do A156919(n, m) := (2*m+2)*A156919(n-1, m) + (2*n-2*m+1) * A156919(n-1, m-1) end do end do: m:=3; for n from m to nmax do a(n, m) := sum((-1)^(m-p1-1)*sum(2^(n-q-1)*binomial(n-q-1, m-p1-1) * A094665(n-1, q) * A156919(q, p1), q=1..n-m+p1), p1=0..m-1) od: seq(a(n, m), n = m..nmax);
    # Maple program edited by Johannes W. Meijer, Sep 25 2012

Formula

a(n) = sum((-1)^(m-p-1)*sum(2^(n-q-1)*binomial(n-q-1,m-p-1)*A094665(n-1,q)* A156919(q,p),q=1..n-m+p), p=0..m-1) with m = 3.

A162008 Third right hand column of the EG1 triangle A162005.

Original entry on oeis.org

16, 1032, 36096, 1035088, 27426960, 702812568, 17753262208, 445736371872, 11162877175440, 279268061007400, 6983654996144256, 174610650719469552, 4365455001524490256, 109138210900706764728, 2728473030627812279040
Offset: 1

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Third right hand column of the EG1 triangle A162005.
Other right hand columns are A004004 (2x), A162009 and A162010.
Cf. A162011 and A162014.

Formula

a(n) = (64*n^2+112*n+30-432*n*9^n-405*9^n+375*25^n)/128.
a(n) = 46*a(n-1)-663*a(n-2)+3748*a(n-3)-7711*a(n-4)+6606*a(n-5)-2025*a(n-6).
G.f.: (16+296*z-768*z^2-1080*z^3)/((1-z)^3*(1-9*z)^2*(1-25*z)).

A162009 Fourth right hand column of the EG1 triangle A162005.

Original entry on oeis.org

272, 52736, 4766048, 319830400, 18598875760, 1002968825344, 51882638754240, 2621627565515520, 130715075544000720, 6468157990602644480, 318685706549526508832, 15663443488266952501376, 768809642314801857986608
Offset: 1

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Fourth right hand column of the EG1 triangle A162005.
Other right hand columns are A004004 (2x), A162008 and A162010.
Cf. A162011 and A162014.

Formula

a(n) = (252105*49^n-525+76545*9^n-328125*25^n-225000*25^n*n-2320*n+134136*9^n*n-2112*n^2+46656*9^n*n^2-512*n^3)/3072
a(n) = 130*a(n-1)-6501*a(n-2)+163160*a(n-3)-2236466*a(n-4)+17123340*a(n-5)-71497186*a(n-6)+154127320*a(n-7)-174334221*a(n-8)+98986050*a(n-9)-22325625*a(n-10)
GF(z) = (272+17376*z-321360*z^2-1298624*z^3+8914800*z^4-11262240*z^5-10206000*z^6)/((1-z)^4*(1-9*z)^3*(1-25*z)^2*(1-49*z))

A162010 Fifth right hand column of the EG1 triangle A162005.

Original entry on oeis.org

7936, 3646208, 704357760, 93989648000, 10324483102720, 1013356176688128, 92857038223998720, 8148225153293502720, 695389790665420312320, 58282750219059501633280, 4827428305286309709508736
Offset: 1

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Fifth right hand column of the EG1 triangle A162005.
Other right hand columns are A004004 (2x), A162008 and A162009.
Cf. A162011 and A162014.

Formula

a(n) = (13230+75008*n^2+65184*n-778248135*49^n+30720*n^3+ 502211745*81^n+ 295312500*25^n-26034048*9^n*n^2-43600032*9^n*n-395300640*49^n*n-19289340*9^n+ 4096*n^4+ 352500000*25^n*n+90000000*25^n*n^2-4478976*9^n*n^3)/98304
a(n)= 295*a(n-1)-36729*a(n-2)+2549775*a(n-3)-109746165*a(n-4)+3080128275*a(n-5)-57713313405*a(n-6)+727045264875*a(n-7)-6122436806115*a(n-8)+33837597147925*a(n-9)-119061300168619*a(n-10)+257794693911405*a(n-11)-339251103039591*a(n-12)+264193039731825*a(n-13)-112000136889375*a(n-14)+19937341265625*a(n-15)
GF(z) = (7936+1305088*z-79792256*z^2-109331968*z^3+41828672000*z^4-460917924352*z^5+238697445120*z^6+5066784271872*z^7-14723693948160*z^8+12172737024000*z^9+8101522800000*z^10)/(1-z)^5/(1-9*z)^4/(1-25*z)^3/(1-49*z)^2/(1-81*z)

A162011 A sequence related to the recurrence relations of the right hand columns of the EG1 triangle A162005.

Original entry on oeis.org

1, -1, 1, -11, 19, -9, 1, -46, 663, -3748, 7711, -6606, 2025, 1, -130, 6501, -163160, 2236466, -17123340, 71497186, -154127320, 174334221, -98986050, 22325625, 1, -295, 36729, -2549775, 109746165, -3080128275, 57713313405, -727045264875
Offset: 1

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Comments

The recurrence relation RR(n) = 0 of the n-th right hand column can be found with RR(n) = expand( product((1-(2*k-1)^2*z)^(n-k+1),k=1..n),z) = 0 and replacing z^p by a(n-p).
The polynomials in the numerators of the generating functions GF(z) of the coefficients that precede the a(n), a(n-1), a(n-2) and a(n-3) sequences, see A000012, A006324, A162012 and A162013, are symmetrical. This phenomenon leads to the sequence [1, 1, 6, 1, 19, 492, 1218, 492, 19 , 9, 3631, 115138, 718465, 1282314, 718465, 115138, 3631, 9].

Examples

			The recurrence relations for the first few right hand columns:
n = 1: a(n) = 1*a(n-1)
n = 2: a(n) = 11*a(n-1)-19*a(n-2)+9*a(n-3)
n = 3: a(n) = 46*a(n-1)-663*a(n-2)+3748*a(n-3)-7711*a(n-4)+6606*a(n-5)-2025*a(n-6)
n = 4: a(n) = 130*a(n-1)-6501*a(n-2)+163160*a(n-3)-2236466*a(n-4)+17123340*a(n-5)-71497186*a(n-6)+154127320*a(n-7)-174334221*a(n-8)+98986050*a(n-9)-22325625*a(n-10)
		

Crossrefs

A000012, A004004 (2x), A162008, A162009 and A162010 are the first five right hand columns of EG1 triangle A162005.
A000124 (the Lazy Caterer's sequence) gives the number of terms of the RR(n).
A006324, A162012 and A162013 equal the absolute values of the coefficients that precede the a(n-1), a(n-2) and a(n-3) factors of the RR(n).

Programs

  • Maple
    nmax:=5; for n from 1 to nmax do RR(n) := expand(product((1-(2*k-1)^2*z)^(n-k+1), k=1..n), z) od: T:=1: for n from 1 to nmax do for m from 0 to(n)*(n+1)/2 do a(T):= coeff(RR(n), z, m): T:=T+1 od: od: seq(a(k), k=1..T-1);

Formula

RR(n) = expand( product((1-(2*k-1)^2*z)^(n-k+1),k=1..n),z) with n = 1, 2, 3, .. . The coefficients of these polynomials lead to the sequence given above.

A162006 Second left hand column of the EG1 triangle A162005.

Original entry on oeis.org

1, 28, 1032, 52736, 3646208, 330545664, 38188155904, 5488365862912, 961530104709120, 201865242068910080, 50052995352723193856, 14476381898608390176768, 4831399425299156001882112
Offset: 2

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Second left hand column of the EG1 triangle A162005.
Other left hand columns are A000182 and A162007.
Related to A094665, A083061 and A156919.
A000079 and A036289 appear in the Maple program.

Programs

  • Maple
    nmax := 14; mmax := nmax: imax := nmax: T1(0, x) := 1: T1(0, x+1) := 1: for i from 1 to imax do T1(i, x) := expand((2*x+1)*(x+1)*T1(i-1, x+1) - 2*x^2*T1(i-1, x)): dx := degree(T1(i, x)): for k from 0 to dx do c(k) := coeff(T1(i, x), x, k) od: T1(i, x+1) := sum(c(j1)*(x+1)^(j1), j1 = 0..dx): od: for i from 0 to imax do for j from 0 to i do A083061(i, j) := coeff(T1(i, x), x, j) od: od: for n from 0 to nmax do for k from 0 to n do A094665(n+1, k+1) := A083061(n, k) od: od: A094665(0, 0) := 1: for n from 1 to nmax do A094665(n, 0) := 0 od: for m from 1 to mmax do A156919(0, m) := 0 end do: for n from 0 to nmax do A156919(n, 0) := 2^n end do: for n from 1 to nmax do for m from 1 to mmax do A156919(n, m) := (2*m+2)*A156919(n-1, m) + (2*n-2*m+1) * A156919(n-1, m-1) end do end do: m:=2; for n from m to nmax do a(n, m) := sum((-1)^(m-p1-1)*sum(2^(n-q-1)*binomial(n-q-1, m-p1-1) * A094665(n-1, q) * A156919(q, p1), q=1..n-m+p1), p1=0..m-1) od: seq(a(n, m), n = m..nmax);
    # Maple program edited by Johannes W. Meijer, Sep 25 2012

Formula

a(n) = sum((-1)^(m-p-1)*sum(2^(n-q-1)*binomial(n-q-1,m-p-1)*A094665(n-1,q)*A156919(q,p),q=1..n-m+p), p=0..m-1) with m = 2.

A162014 Sequence related to the o.g.f.s. of the right hand columns of the EG1 triangle A162005.

Original entry on oeis.org

1, 8, -1536, -14155776, 10436770529280, 923378661099307008000, -13724698564186788948502118400000, -45695540009113634492156662349750599680000000
Offset: 1

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Comments

The a(n) are the sums of the coefficients of the polynomials that appear in the numerators of the o.g.f.s. of the right hand columns of the EG1 triangle A162005, see the examples.

Examples

			The polynomials in the numerators of the first few o.g.f.s are:
numer(GF(1)) = 1
numer(GF(2)) = 2+6*z
numer(GF(3)) = 16+296*z-768*z^2-1080*z^3
numer(GF(4)) = 272+17376*z-321360*z^2-1298624*z^3+8914800*z^4-11262240*z^5-10206000*z^6
numer(GF(5)) = 7936 + 1305088*z - 79792256*z^2 - 109331968*z^3 + 41828672000*z^4-460917924352*z^5 + 238697445120*z^6 + 5066784271872*z^7 - 14723693948160*z^8+ 12172737024000*z^9 + 8101522800000*z^10
		

Crossrefs

A000012, A004004 (2x), A162008, A162009 and A162010 are the first five right hand columns of the EG1 triangle A162005.
Cf. A055209 and A059332.

Formula

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

A000182 Tangent (or "Zag") numbers: e.g.f. tan(x), also (up to signs) e.g.f. tanh(x).

Original entry on oeis.org

1, 2, 16, 272, 7936, 353792, 22368256, 1903757312, 209865342976, 29088885112832, 4951498053124096, 1015423886506852352, 246921480190207983616, 70251601603943959887872, 23119184187809597841473536, 8713962757125169296170811392, 3729407703720529571097509625856
Offset: 1

Views

Author

Keywords

Comments

Number of Joyce trees with 2n-1 nodes. Number of tremolo permutations of {0,1,...,2n}. - Ralf Stephan, Mar 28 2003
The Hankel transform of this sequence is A000178(n) for n odd = 1, 12, 34560, ...; example: det([1, 2, 16; 2, 16, 272, 16, 272, 7936]) = 34560. - Philippe Deléham, Mar 07 2004
a(n) is the number of increasing labeled full binary trees with 2n-1 vertices. Full binary means every non-leaf vertex has two children, distinguished as left and right; labeled means the vertices are labeled 1,2,...,2n-1; increasing means every child has a label greater than its parent. - David Callan, Nov 29 2007
From Micha Hofri (hofri(AT)wpi.edu), May 27 2009: (Start)
a(n) was found to be the number of permutations of [2n] which when inserted in order, to form a binary search tree, yield the maximally full possible tree (with only one single-child node).
The e.g.f. is sec^2(x)=1+tan^2(x), and the same coefficients can be manufactured from the tan(x) itself, which is the e.g.f. for the number of trees as above for odd number of nodes. (End)
a(n) is the number of increasing strict binary trees with 2n-1 nodes. For more information about increasing strict binary trees with an associated permutation, see A245894. - Manda Riehl, Aug 07 2014
For relations to alternating permutations, Euler and Bernoulli polynomials, zigzag numbers, trigonometric functions, Fourier transform of a square wave, quantum algebras, and integrals over and in n-dimensional hypercubes and over Green functions, see Hodges and Sukumar. For further discussion on the quantum algebra, see the later Hodges and Sukumar reference and the paper by Hetyei presenting connections to the general combinatorial theory of Viennot on orthogonal polynomials, inverse polynomials, tridiagonal matrices, and lattice paths (thereby related to continued fractions and cumulants). - Tom Copeland, Nov 30 2014
The Zigzag Hankel transform is A000178. That is, A000178(2*n - k) = det( [a(i+j - k)]{i,j = 1..n} ) for n>0 and k=0,1. - _Michael Somos, Mar 12 2015
a(n) is the number of standard Young tableaux of skew shape (n,n,n-1,n-2,...,3,2)/(n-1,n-2,n-3,...,2,1). - Ran Pan, Apr 10 2015
For relations to the Sheffer Appell operator calculus and a Riccati differential equation for generating the Meixner-Pollaczek and Krawtchouk orthogonal polynomials, see page 45 of the Feinsilver link and Rzadkowski. - Tom Copeland, Sep 28 2015
For relations to an elliptic curve, a Weierstrass elliptic function, the Lorentz formal group law, a Lie infinitesimal generator, and the Eulerian numbers A008292, see A155585. - Tom Copeland, Sep 30 2015
Absolute values of the alternating sums of the odd-numbered rows (where the single 1 at the apex of the triangle is counted as row #1) of the Eulerian triangle, A008292. The actual alternating sums alternate in sign, e.g., 1, -2, 16, -272, etc. (Even-numbered rows have alternating sums always 0.) - Gregory Gerard Wojnar, Sep 28 2018
The sequence is periodic modulo any odd prime p. The minimal period is (p-1)/2 if p == 1 mod 4 and p-1 if p == 3 mod 4 [Knuth & Buckholtz, 1967, Theorem 1]. - Allen Stenger, Aug 03 2020
From Peter Bala, Dec 24 2021: (Start)
Conjectures:
1) The sequence taken modulo any integer k eventually becomes periodic with period dividing phi(k).
2) The Gauss congruences a(n*p^k) == a(n*p^(k-1)) ( mod p^k ) hold for all prime p and positive integers n and k, except when p = 2, n = 1 and k = 1 or 2.
3) For i >= 1 define a_i(n) = a(n+i). The Gauss congruences a_i(n*p^k) == a_i(n*p^(k-1)) ( mod p^k ) hold for all prime p and positive integers n and k. If true, then for each i >= 1 the expansion of exp(Sum_{n >= 1} a_i(n)*x^n/n) has integer coefficients. For an example, see A262145.(End)

Examples

			tan(x) = x + 2*x^3/3! + 16*x^5/5! + 272*x^7/7! + ... = x + 1/3*x^3 + 2/15*x^5 + 17/315*x^7 + 62/2835*x^9 + O(x^11).
tanh(x) = x - 1/3*x^3 + 2/15*x^5 - 17/315*x^7 + 62/2835*x^9 - 1382/155925*x^11 + ...
(sec x)^2 = 1 + x^2 + 2/3*x^4 + 17/45*x^6 + ...
a(3)=16 because we have: {1, 3, 2, 5, 4}, {1, 4, 2, 5, 3}, {1, 4, 3, 5, 2},
  {1, 5, 2, 4, 3}, {1, 5, 3, 4, 2}, {2, 3, 1, 5, 4}, {2, 4, 1, 5, 3},
  {2, 4, 3, 5, 1}, {2, 5, 1, 4, 3}, {2, 5, 3, 4, 1}, {3, 4, 1, 5, 2},
  {3, 4, 2, 5, 1}, {3, 5, 1, 4, 2}, {3, 5, 2, 4, 1}, {4, 5, 1, 3, 2},
  {4, 5, 2, 3, 1}. - _Geoffrey Critzer_, May 19 2013
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 932.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 88.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 111.
  • H. Doerrie, 100 Great Problems of Elementary Mathematics, Dover, NY, 1965, p. 69.
  • L. M. Milne-Thompson, Calculus of Finite Differences, 1951, p. 148 (the numbers |C^{2n-1}|).
  • J. W. Milnor and J. D. Stasheff, Characteristic Classes, Princeton, 1974, p. 282.
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 444.
  • H. Rademacher, Topics in Analytic Number Theory, Springer, 1973, Chap. 1, p. 20.
  • L. Seidel, Über eine einfache Entstehungsweise der Bernoullischen Zahlen und einiger verwandten Reihen, Sitzungsberichte der mathematisch-physikalischen Classe der königlich bayerischen Akademie der Wissenschaften zu München, volume 7 (1877), 157-187.
  • 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).
  • E. van Fossen Conrad, Some continued fraction expansions of elliptic functions, PhD thesis, The Ohio State University, 2002, p. 28.
  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, pp. 267-268.

Crossrefs

A350972 is essentially the same sequence.
a(n)=2^(n-1)*A002105(n). Apart from signs, 2^(2n-2)*A001469(n) = n*a(n).
Cf. A001469, A002430, A036279, A000364 (secant numbers), A000111 (secant-tangent numbers), A024283, A009764. First diagonal of A059419 and of A064190.
Equals A002425(n) * 2^A101921(n).
Equals leftmost column of A162005. - Johannes W. Meijer, Jun 27 2009

Programs

  • Maple
    series(tan(x),x,40);
    with(numtheory): a := n-> abs(2^(2*n)*(2^(2*n)-1)*bernoulli(2*n)/(2*n));
    A000182_list := proc(n) local T,k,j; T[1] := 1;
    for k from 2 to n do T[k] := (k-1)*T[k-1] od;
       for k from 2 to n do
           for j from k to n do
               T[j] := (j-k)*T[j-1]+(j-k+2)*T[j] od od;
    seq(T[j], j=1..n)  end:
    A000182_list(15);  # Peter Luschny, Apr 02 2012
  • Mathematica
    Table[ Sum[2^(2*n + 1 - k)*(-1)^(n + k + 1)*k!*StirlingS2[2*n + 1, k], {k, 1, 2*n + 1}], {n, 0, 7}] (* Victor Adamchik, Oct 05 2005 *)
    v[1] = 2; v[n_] /; n >= 2 := v[n] = Sum[ Binomial[2 n - 3, 2 k - 2] v[k] v[n - k], {k, n - 1}]; Table[ v[n]/2, {n, 15}] (* Zerinvary Lajos, Jul 08 2009 *)
    Rest@ Union[ Range[0, 29]! CoefficientList[ Series[ Tan[x], {x, 0, 30}], x]] (* Harvey P. Dale, Oct 19 2011; modified by Robert G. Wilson v, Apr 02 2012 *)
    t[1, 1] = 1; t[1, 0] = 0; t[n_ /; n > 1, m_] := t[n, m] = m*(m+1)*Sum[t[n-1, k], {k, m-1, n-1}]; a[n_] := t[n, 1]; Table[a[n], {n, 1, 15}]  (* Jean-François Alcover, Jan 02 2013, after A064190 *)
    a[ n_] := If[ n < 1, 0, With[{m = 2 n - 1}, m! SeriesCoefficient[ Tan[x], {x, 0, m}]]]; (* Michael Somos, Mar 12 2015 *)
    a[ n_] := If[ n < 1, 0, ((-16)^n - (-4)^n) Zeta[1 - 2 n]]; (* Michael Somos, Mar 12 2015 *)
    Table[2 PolyGamma[2n - 1, 1/2]/Pi^(2n), {n, 1, 10}] (* Vladimir Reshetnikov, Oct 18 2015 *)
    a[ n_] := a[n] = If[ n < 2, Boole[n == 1], Sum[Binomial[2 n - 2, 2 k - 1] a[k] a[n - k], {k, n - 1}]]; (* Michael Somos, Aug 02 2018 *)
    a[n_] := (2^(2*n)*(2^(2*n) - 1)*Abs[BernoulliB[2*n]])/(2*n); a /@  Range[20]  (* Stan Wagon, Nov 21 2022 *)
  • Maxima
    a(n):=sum(sum(binomial(k,r)*sum(sum(binomial(l,j)/2^(j-1)*sum((-1)^(n)*binomial(j,i)*(j-2*i)^(2*n),i,0,floor((j-1)/2))*(-1)^(l-j),j,1,l)*(-1)^l*binomial(r+l-1,r-1),l,1,2*n)*(-1)^(1-r),r,1,k)/k,k,1,2*n); /* Vladimir Kruchinin, Aug 23 2010 */
    
  • Maxima
    a[n]:=if n=1 then 1 else 2*sum(sum(binomial(2*j,j+k)*(-4*k^2)^(n-1)*(-1)^k/(4^j),k,1,j),j,1,n-1);
    makelist(a[n],n,1,30); /* Tani Akinari, Sep 20 2023 */
    
  • PARI
    {a(n) = if( n<1, 0, ((-4)^n - (-16)^n) * bernfrac(2*n) / (2*n))};
    
  • PARI
    {a(n) = my(an); if( n<2, n==1, an = vector(n, m, 1); for( m=2, n, an[m] = sum( k=1, m-1, binomial(2*m - 2, 2*k - 1) * an[k] * an[m-k])); an[n])}; /* Michael Somos */
    
  • PARI
    {a(n) = if( n<1, 0, (2*n - 1)! * polcoeff( tan(x + O(x^(2*n + 2))), 2*n - 1))}; /* Michael Somos */
    
  • PARI
    {a(n) = my(X=x+x*O(x^n),Egf); Egf = x*sum(m=0,n, prod(k=1,m, tanh(2*k*X))); (n-1)!*polcoeff(Egf,n)} /* Paul D. Hanna, May 11 2010 */
    
  • PARI
    /* Continued Fraction for the e.g.f. tan(x), from Paul D. Hanna: */
    {a(n)=local(CF=1+O(x)); for(i=1, n, CF=1/(2*(n-i+1)-1-x^2*CF)); (2*n-1)!*polcoeff(x*CF, 2*n-1)}
    
  • PARI
    /* O.g.f. Sum_{n>=1} a(n)*x^n, from Paul D. Hanna Feb 05 2013: */
    {a(n)=polcoeff( x+2*x*sum(m=1, n, x^m*prod(k=1, m, (2*k-1)^2/(1+(2*k-1)^2*x +x*O(x^n))) ), n)}
    
  • Python
    # The objective of this implementation is efficiency.
    # n -> [0, a(1), a(2), ..., a(n)] for n > 0.
    def A000182_list(n):
        T = [0 for i in range(1, n+2)]
        T[1] = 1
        for k in range(2, n+1):
            T[k] = (k-1)*T[k-1]
        for k in range(2, n+1):
            for j in range(k, n+1):
                T[j] = (j-k)*T[j-1]+(j-k+2)*T[j]
        return T
    print(A000182_list(100)) # Peter Luschny, Aug 07 2011
    
  • Python
    from sympy import bernoulli
    def A000182(n): return abs(((2-(2<<(m:=n<<1)))*bernoulli(m)<Chai Wah Wu, Apr 14 2023
    
  • Sage
    # Algorithm of L. Seidel (1877)
    # n -> [a(1), ..., a(n)] for n >= 1.
    def A000182_list(len) :
        R = []; A = {-1:0, 0:1}; k = 0; e = 1
        for i in (0..2*len-1) :
            Am = 0; A[k + e] = 0; e = -e
            for j in (0..i) : Am += A[k]; A[k] = Am; k += e
            if e > 0 : R.append(A[i//2])
        return R
    A000182_list(15) # Peter Luschny, Mar 31 2012

Formula

E.g.f.: log(sec x) = Sum_{n > 0} a(n)*x^(2*n)/(2*n)!.
E.g.f.: tan x = Sum_{n >= 0} a(n+1)*x^(2*n+1)/(2*n+1)!.
E.g.f.: (sec x)^2 = Sum_{n >= 0} a(n+1)*x^(2*n)/(2*n)!.
2/(exp(2x)+1) = 1 + Sum_{n>=1} (-1)^(n+1) a(n) x^(2n-1)/(2n-1)! = 1 - x + x^3/3 - 2*x^5/15 + 17*x^7/315 - 62*x^9/2835 + ...
a(n) = 2^(2*n) (2^(2*n) - 1) |B_(2*n)| / (2*n) where B_n are the Bernoulli numbers (A000367/A002445 or A027641/A027642).
Asymptotics: a(n) ~ 2^(2*n+1)*(2*n-1)!/Pi^(2*n).
Sum[2^(2*n + 1 - k)*(-1)^(n + k + 1)*k!*StirlingS2[2*n + 1, k], {k, 1, 2*n + 1}]. - Victor Adamchik, Oct 05 2005
a(n) = abs[c(2*n-1)] where c(n)= 2^(n+1) * (1-2^(n+1)) * Ber(n+1)/(n+1) = 2^(n+1) * (1-2^(n+1)) * (-1)^n * Zeta(-n) = [ -(1+EN(.))]^n = 2^n * GN(n+1)/(n+1) = 2^n * EP(n,0) = (-1)^n * E(n,-1) = (-2)^n * n! * Lag[n,-P(.,-1)/2] umbrally = (-2)^n * n! * C{T[.,P(.,-1)/2] + n, n} umbrally for the signed Euler numbers EN(n), the Bernoulli numbers Ber(n), the Genocchi numbers GN(n), the Euler polynomials EP(n,t), the Eulerian polynomials E(n,t), the Touchard / Bell polynomials T(n,t), the binomial function C(x,y) = x!/[(x-y)!*y! ] and the polynomials P(j,t) of A131758. - Tom Copeland, Oct 05 2007
a(1) = A094665(0,0)*A156919(0,0) and a(n) = Sum_{k=1..n-1} 2^(n-k-1)*A094665(n-1, k)*A156919(k,0) for n = 2, 3, .., see A162005. - Johannes W. Meijer, Jun 27 2009
G.f.: 1/(1-1*2*x/(1-2*3*x/(1-3*4*x/(1-4*5*x/(1-5*6*x/(1-... (continued fraction). - Paul Barry, Feb 24 2010
From Paul Barry, Mar 29 2010: (Start)
G.f.: 1/(1-2x-12x^2/(1-18x-240x^2/(1-50x-1260x^2/(1-98x-4032x^2/(1-162x-9900x^2/(1-... (continued fraction);
coefficient sequences given by 4*(n+1)^2*(2n+1)*(2n+3) and 2(2n+1)^2 (see Van Fossen Conrad reference). (End)
E.g.f.: x*Sum_{n>=0} Product_{k=1..n} tanh(2*k*x) = Sum_{n>=1} a(n)*x^n/(n-1)!. - Paul D. Hanna, May 11 2010 [corrected by Paul D. Hanna, Sep 28 2023]
a(n) = (-1)^(n+1)*Sum_{j=1..2*n+1} j!*Stirling2(2*n+1,j)*2^(2*n+1-j)*(-1)^j for n >= 0. Vladimir Kruchinin, Aug 23 2010: (Start)
If n is odd such that 2*n-1 is prime, then a(n) == 1 (mod (2*n-1)); if n is even such that 2*n-1 is prime, then a(n) == -1 (mod (2*n-1)). - Vladimir Shevelev, Sep 01 2010
Recursion: a(n) = (-1)^(n-1) + Sum_{i=1..n-1} (-1)^(n-i+1)*C(2*n-1,2*i-1)* a(i). - Vladimir Shevelev, Aug 08 2011
E.g.f.: tan(x) = Sum_{n>=1} a(n)*x^(2*n-1)/(2*n-1)! = x/(1 - x^2/(3 - x^2/(5 - x^2/(7 - x^2/(9 - x^2/(11 - x^2/(13 -...))))))) (continued fraction from J. H. Lambert - 1761). - Paul D. Hanna, Sep 21 2011
From Sergei N. Gladkovskii, Oct 31 2011 to Oct 09 2013: (Start)
Continued fractions:
E.g.f.: (sec(x))^2 = 1+x^2/(x^2+U(0)) where U(k) = (k+1)*(2k+1) - 2x^2 + 2x^2*(k+1)*(2k+1)/U(k+1).
E.g.f.: tan(x) = x*T(0) where T(k) = 1-x^2/(x^2-(2k+1)*(2k+3)/T(k+1)).
E.g.f.: tan(x) = x/(G(0)+x) where G(k) = 2*k+1 - 2*x + x/(1 + x/G(k+1)).
E.g.f.: tanh(x) = x/(G(0)-x) where G(k) = k+1 + 2*x - 2*x*(k+1)/G(k+1).
E.g.f.: tan(x) = 2*x - x/W(0) where W(k) = 1 + x^2*(4*k+5)/((4*k+1)*(4*k+3)*(4*k+5) - 4*x^2*(4*k+3) + x^2*(4*k+1)/W(k+1)).
E.g.f.: tan(x) = x/T(0) where T(k) = 1 - 4*k^2 + x^2*(1 - 4*k^2)/T(k+1).
E.g.f.: tan(x) = -3*x/(T(0)+3*x^2) where T(k)= 64*k^3 + 48*k^2 - 4*k*(2*x^2 + 1) - 2*x^2 - 3 - x^4*(4*k -1)*(4*k+7)/T(k+1).
G.f.: 1/G(0) where G(k) = 1 - 2*x*(2*k+1)^2 - x^2*(2*k+1)*(2*k+2)^2*(2*k+3)/G(k+1).
G.f.: 2*Q(0) - 1 where Q(k) = 1 + x^2*(4*k + 1)^2/(x + x^2*(4*k + 1)^2 - x^2*(4*k + 3)^2*(x + x^2*(4*k + 1)^2)/(x^2*(4*k + 3)^2 + (x + x^2*(4*k + 3)^2)/Q(k+1) )).
G.f.: (1 - 1/G(0))*sqrt(-x), where G(k) = 1 + sqrt(-x) - x*(k+1)^2/G(k+1).
G.f.: Q(0), where Q(k) = 1 - x*(k+1)*(k+2)/( x*(k+1)*(k+2) - 1/Q(k+1)). (End)
O.g.f.: x + 2*x*Sum_{n>=1} x^n * Product_{k=1..n} (2*k-1)^2 / (1 + (2*k-1)^2*x). - Paul D. Hanna, Feb 05 2013
a(n) = (-4)^n*Li_{1-2*n}(-1). - Peter Luschny, Jun 28 2012
a(n) = (-4)^n*(4^n-1)*Zeta(1-2*n). - Jean-François Alcover, Dec 05 2013
Asymptotic expansion: 4*((2*(2*n-1))/(Pi*e))^(2*n-1/2)*exp(1/2+1/(12*(2*n-1))-1/(360*(2*n-1)^3)+1/(1260*(2*n-1)^5)-...). (See Luschny link.) - Peter Luschny, Jul 14 2015
From Peter Bala, Sep 11 2015: (Start)
The e.g.f. A(x) = tan(x) satisfies the differential equation A''(x) = 2*A(x)*A'(x) with A(0) = 0 and A'(0) = 1, leading to the recurrence a(0) = 0, a(1) = 1, else a(n) = 2*Sum_{i = 0..n-2} binomial(n-2,i)*a(i)*a(n-1-i) for the aerated sequence [0, 1, 0, 2, 0, 16, 0, 272, ...].
Note, the same recurrence, but with the initial conditions a(0) = 1 and a(1) = 1, produces the sequence n! and with a(0) = 1/2 and a(1) = 1 produces A080635. Cf. A002105, A234797. (End)
a(n) = 2*polygamma(2*n-1, 1/2)/Pi^(2*n). - Vladimir Reshetnikov, Oct 18 2015
a(n) = 2^(2n-2)*|p(2n-1,-1/2)|, where p_n(x) are the shifted row polynomials of A019538. E.g., a(2) = 2 = 2^2 * |1 + 6(-1/2) + 6(-1/2)^2|. - Tom Copeland, Oct 19 2016
From Peter Bala, May 05 2017: (Start)
With offset 0, the o.g.f. A(x) = 1 + 2*x + 16*x^2 + 272*x^3 + ... has the property that its 4th binomial transform 1/(1 - 4*x) A(x/(1 - 4*x)) has the S-fraction representation 1/(1 - 6*x/(1 - 2*x/(1 - 20*x/(1 - 12*x/(1 - 42*x/(1 - 30*x/(1 - ...))))))), where the coefficients [6, 2, 20, 12, 42, 30, ...] in the partial numerators of the continued fraction are obtained from the sequence [2, 6, 12, 20, ..., n*(n + 1), ...] by swapping adjacent terms. Compare with the S-fraction associated with A(x) given above by Paul Barry.
A(x) = 1/(1 + x - 3*x/(1 - 4*x/(1 + x - 15*x/(1 - 16*x/(1 + x - 35*x/(1 - 36*x/(1 + x - ...))))))), where the unsigned coefficients in the partial numerators [3, 4, 15, 16, 35, 36,...] come in pairs of the form 4*n^2 - 1, 4*n^2 for n = 1,2,.... (End)
a(n) = Sum_{k = 1..n-1} binomial(2*n-2, 2*k-1) * a(k) * a(n-k), with a(1) = 1. - Michael Somos, Aug 02 2018
a(n) = 2^(2*n-1) * |Euler(2*n-1, 0)|, where Euler(n,x) are the Euler polynomials. - Daniel Suteu, Nov 21 2018 (restatement of one of Copeland's 2007 formulas.)
x - Sum_{n >= 1} (-1)^n*a(n)*x^(2*n)/(2*n)! = x - log(cosh(x)). The series reversion of x - log(cosh(x)) is (1/2)*x - (1/2)*log(2 - exp(x)) = Sum_{n >= 0} A000670(n)*x^(n+1)/(n+1)!. - Peter Bala, Jul 11 2022
For n > 1, a(n) = 2*Sum_{j=1..n-1} Sum_{k=1..j} binomial(2*j,j+k)*(-4*k^2)^(n-1)*(-1)^k/(4^j). - Tani Akinari, Sep 20 2023
a(n) = A110501(n) * 4^(n-1) / n (Han and Liu, 2018). - Amiram Eldar, May 17 2024

A079484 a(n) = (2n-1)!! * (2n+1)!!, where the double factorial is A006882.

Original entry on oeis.org

1, 3, 45, 1575, 99225, 9823275, 1404728325, 273922023375, 69850115960625, 22561587455281875, 9002073394657468125, 4348001449619557104375, 2500100833531245335015625, 1687568062633590601135546875, 1321365793042101440689133203125
Offset: 0

Views

Author

Benoit Cloitre, Jan 17 2003

Keywords

Comments

a(n) is the determinant of M(2n+1) where M(k) is the k X k matrix with m(i,j)=j if i+j=k m(i,j)=i otherwise. - Adapted to offset 0, Rainer Rosenthal, Jun 19 2024
In the following two comments on the calculation of the terms using permanents, offset 1 is assumed. In the corresponding PARI code, this is implemented with offset 0. - Hugo Pfoertner, Jun 23 2024
(-1)^n*a(n)/2^(2n-1) is the permanent of the (m X m) matrix {1/(x_i-y_j), 1<=i<=m, 1<=j<=m}, where x_1,x_2,...,x_m are the zeros of x^m-1 and y_1,y_2,...,y_m the zeros of y^m+1 and m=2n-1.
In 1881, R. F. Scott posed a conjecture that the absolute value of permanent of square matrix with elements a(i,j)= (x_i - y_j)^(-1), where x_1,...,x_n are roots of x^n=1, while y_1,...,y_n are roots of y^n=-1, equals a((n-1)/2)/2^n, if n>=1 is odd, and 0, if n>=2 is even. After a century (in 1979), the conjecture was proved by H. Minc. - Vladimir Shevelev, Dec 01 2013
a(n) is the number of permutations in S_{2n+1} in which all cycles have odd length. - José H. Nieto S., Jan 09 2012
Number of 3-bundled increasing bilabeled trees with 2n labels. - Markus Kuba, Nov 18 2014
a(n) is the number of rooted, binary, leaf-labeled topologies with 2n+2 leaves that have n+1 cherry nodes. - Noah A Rosenberg, Feb 12 2019

Examples

			G.f. = 1 + 3*x + 45*x^2 + 1575*x^3 + 99225*x^4 + 9823275*x^5 + ...
M(5) =
  [1, 2, 3, 1, 5]
  [1, 2, 2, 4, 5]
  [1, 3, 3, 4, 5]
  [4, 2, 3, 4, 5]
  [1, 2, 3, 4, 5].
Integral_{x=0..oo} x^3*BesselK(1, sqrt(x)) = 1575*Pi. - _Olivier Gérard_, May 20 2009
		

References

  • Miklós Bóna, A walk through combinatorics, World Scientific, 2006.

Crossrefs

Bisection of A000246, A053195, |A013069|, |A046126|. Cf. A000909.
Cf. A001044, A010791, |A129464|, A114779, are also values of similar moments.
Equals the row sums of A162005.
Cf. A316728.
Diagonal elements of A306364 in even-numbered rows.

Programs

  • Magma
    I:=[1, 3]; [n le 2 select I[n] else (4*n^2-8*n+3)*Self(n-1): n in [1..20]]; // Vincenzo Librandi, Nov 18 2014
    
  • Maple
    a:= n-> (d-> d(2*n-1)*d(2*n+1))(doublefactorial):
    seq(a(n), n=0..15);  # Alois P. Heinz, Jan 30 2013
    # second Maple program:
    A079484 := n-> LinearAlgebra[Determinant](Matrix(2*n+1, (i, j)-> `if`(i+j=2*n+1, j, i))): seq(A079484(n), n=0..14); # Rainer Rosenthal, Jun 18 2024
  • Mathematica
    a[n_] := (2n - 1)!!*(2n + 1)!!; Table[a[n], {n, 0, 13}] (* Jean-François Alcover, Jan 30 2013 *)
  • PARI
    /* Formula using the zeta function and a log integral:*/
    L(n)= intnum(t=0, 1, log(1-1/t)^n);
    Zetai(n)= -I*I^n*(2*Pi)^(n-1)/(n-1)*L(1-n);
    a(m)={my(n=m+1);round(real(-I*2^(2*n-1)*Zetai(1/2-n)*L(-1/2+n)/(Zetai(-1/2+n)*L(1/2-n))))};
    /* Gerry Martens, Mar 07 2011, adapted to offset 0 by Hugo Pfoertner, Jun 19 2024 */
    
  • PARI
    {a(n) = if( n<0, -1 / self()(-1-n), (2*n + 1)! * (2*n)! / (n! * 2^n)^2 )}; /* Michael Somos, May 04 2017 */
    
  • PARI
    {a(n) = if( n<0, -1 / self()(-1-n), my(m = 2*n + 1); m! * polcoeff( x / sqrt( 1 - x^2 + x * O(x^m) ), m))}; /* Michael Somos, May 04 2017 */
    
  • PARI
    \\ using the Pochhammer symbol
    a(n) = {my(P(x,k)=gamma(x+k)/gamma(x)); 4^n*round(P(1/2,n)*P(3/2,n))} \\ Hugo Pfoertner, Jun 20 2024
    
  • PARI
    \\ Scott's (1881) method
    a(n) = {my(m=2*n+1, X = polroots(x^m-1), Y = polroots(x^m+1), M = matrix(m, m, i, j, 1/(X[i]-Y[j]))); (-1)^n * round(2^m * real(matpermanent(M)))}; \\ Hugo Pfoertner, Jun 23 2024

Formula

D-finite with recurrence a(n) = (4*n^2 - 1) * a(n-1) for all n in Z.
a(n) = A001147(n)*A001147(n+1).
E.g.f.: 1/(1-x^2)^(3/2) (with interpolated zeros). - Paul Barry, May 26 2003
a(n) = (2n+1)! * C(2n, n) / 2^(2n). - Ralf Stephan, Mar 22 2004.
Alternatingly signed values have e.g.f. sqrt(1+x^2).
a(n) is the value of the n-th moment of (1/Pi)*BesselK(1, sqrt(x)) on the positive part of the real line. - Olivier Gérard, May 20 2009
a(n) = -2^(2*n-1)*exp(i*n*Pi)*gamma(1/2+n)/gamma(3/2-n). - Gerry Martens, Mar 07 2011
E.g.f. (odd powers) tan(arcsin(x)) = Sum_{n>=0} (2n-1)!!*(2n+1)!!*x^(2*n+1)/(2*n+1)!. - Vladimir Kruchinin, Apr 22 2011
G.f.: 1 + x*(G(0) - 1)/(x-1) where G(k) = 1 - ((2*k+2)^2-1)/(1-x/(x - 1/G(k+1))); ( continued fraction ). - Sergei N. Gladkovskii, Jan 15 2013
a(n) = (2^(2*n+3)*Gamma(n+3/2)*Gamma(n+5/2))/Pi. - Jean-François Alcover, Jul 20 2015
Limit_{n->oo} 4^n*(n!)^2/a(n) = Pi/2. - Daniel Suteu, Feb 05 2017
From Michael Somos, May 04 2017: (Start)
a(n) = (2*n + 1) * A001818(n).
E.g.f.: Sum_{n>=0} a(n) * x^(2*n+1) / (2*n+1)! = x / sqrt(1 - x^2) = tan(arcsin(x)).
Given e.g.f. A(x) = y, then x * y' = y + y^3.
a(n) = -1 / a(-1-n) for all n in Z.
0 = +a(n)*(+288*a(n+2) -60*a(n+3) +a(n+4)) +a(n+1)*(-36*a(n+2) -4*a(n+3)) +a(n+2)*(+3*a(n+2)) for all n in Z. (End)
a(n) = Sum_{k=0..2n} (k+1) * A316728(n,k). - Alois P. Heinz, Jul 12 2018
From Amiram Eldar, Mar 18 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 + L_1(1)*Pi/2, where L is the modified Struve function.
Sum_{n>=0} (-1)^n/a(n) = 1 - H_1(1)*Pi/2, where H is the Struve function. (End)

Extensions

Simpler description from Daniel Flath (deflath(AT)yahoo.com), Mar 05 2004

A156919 Table of coefficients of polynomials related to the Dirichlet eta function.

Original entry on oeis.org

1, 2, 1, 4, 10, 1, 8, 60, 36, 1, 16, 296, 516, 116, 1, 32, 1328, 5168, 3508, 358, 1, 64, 5664, 42960, 64240, 21120, 1086, 1, 128, 23488, 320064, 900560, 660880, 118632, 3272, 1, 256, 95872, 2225728
Offset: 0

Views

Author

Johannes W. Meijer, Feb 20 2009, Jun 24 2009

Keywords

Comments

Essentially the same as A185411. Row reverse of A185410. - Peter Bala, Jul 24 2012
The SF(z; n) formulas, see below, were discovered while studying certain properties of the Dirichlet eta function.
From Peter Bala, Apr 03 2011: (Start)
Let D be the differential operator 2*x*d/dx. The row polynomials of this table come from repeated application of the operator D to the function g(x) = 1/sqrt(1 - x). For example,
D(g) = x*g^3
D^2(g) = x*(2 + x)*g^5
D^3(g) = x*(4 + 10*x + x^2)*g^7
D^4(g) = x*(8 + 60*x + 36*x^2 + x^3)*g^9.
Thus this triangle is analogous to the triangle of Eulerian numbers A008292, whose row polynomials come from the repeated application of the operator x*d/dx to the function 1/(1 - x). (End)

Examples

			The first few rows of the triangle are:
  [1]
  [2, 1]
  [4, 10, 1]
  [8, 60, 36, 1]
  [16, 296, 516, 116, 1]
The first few P(z;n) are:
  P(z; n=0) = 1
  P(z; n=1) = 2 + z
  P(z; n=2) = 4 + 10*z + z^2
  P(z; n=3) = 8 + 60*z + 36*z^2 + z^3
The first few SF(z;n) are:
  SF(z; n=0) = (1/2)*(1)/(1-z)^(3/2);
  SF(z; n=1) = (1/4)*(2+z)/(1-z)^(5/2);
  SF(z; n=2) = (1/8)*(4+10*z+z^2)/(1-z)^(7/2);
  SF(z; n=3) = (1/16)*(8+60*z+36*z^2+z^3)/(1-z)^(9/2);
In the Savage-Viswanathan paper, the coefficients appear as
  1;
  1,    2;
  1,   10,     4;
  1,   36,    60,     8;
  1,  116,   516,   296,    16;
  1,  358,  3508,  5168,  1328,   32;
  1, 1086, 21120, 64240, 42960, 5664, 64;
  ...
		

Crossrefs

A142963 and this sequence can be mapped onto the A156920 triangle.
FP1 sequences A000340, A156922, A156923, A156924.
FP2 sequences A050488, A142965, A142966, A142968.
Appears in A162005, A000182, A162006 and A162007.
Cf. A185410 (row reverse), A185411.

Programs

  • Maple
    A156919 := proc(n,m) if n=m then 1; elif m=0 then 2^n ; elif m<0 or m>n then 0; else 2*(m+1)*procname(n-1,m)+(2*n-2*m+1)*procname(n-1,m-1) ; end if; end proc: seq(seq(A156919(n,m), m=0..n), n=0..7); # R. J. Mathar, Feb 03 2011
  • Mathematica
    g[0] = 1/Sqrt[1-x]; g[n_] := g[n] = 2x*D[g[n-1], x]; p[n_] := g[n] / g[0]^(2n+1) // Cancel; row[n_] := CoefficientList[p[n], x] // Rest; Table[row[n], {n, 0, 9}] // Flatten (* Jean-François Alcover, Aug 09 2012, after Peter Bala *)
    Flatten[Table[Rest[CoefficientList[Nest[2 x D[#, x] &, (1 - x)^(-1/2), k] (1 - x)^(k + 1/2), x]], {k, 10}]] (* Jan Mangaldan, Mar 15 2013 *)

Formula

SF(z; n) = Sum_{m >= 1} m^(n-1)*4^(-m)*z^(m-1)*Gamma(2*m+1)/(Gamma(m)^2) = P(z;n) / (2^(n+1)*(1-z)^((2*n+3)/2)) for n >= 0. The polynomials P(z;n) = Sum_{k = 0..n} a(k)*z^k generate the a(n) sequence.
If we write the sequence as a triangle the following relation holds: T(n,m) = (2*m+2)*T(n-1,m) + (2*n-2*m+1)*T(n-1,m-1) with T(n,m=0) = 2^n and T(n,n) = 1, n >= 0 and 0 <= m <= n.
G.f.: 1/(1-xy-2x/(1-3xy/(1-4x/(1-5xy/(1-6x/(1-7xy/(1-8x/(1-... (continued fraction). - Paul Barry, Jan 26 2011
From Peter Bala, Apr 03 2011 (Start)
E.g.f.: exp(z*(x + 2)) * (1 - x)/(exp(2*x*z) - x*exp(2*z))^(3/2) = Sum_{n >= 0} P(x,n)*z^n/n! = 1 + (2 + x)*z + (4 + 10*x + x^2)*z^2/2! + (8 + 60*x + 36*x^2 + x^3)*z^3/3! + ... .
Explicit formula for the row polynomials:
P(x,n-1) = Sum_{k = 1..n} 2^(n-2*k)*binomial(2k,k)*k!*Stirling2(n,k)*x^(k-1)*(1 - x)^(n-k).
The polynomials x*(1 + x)^n * P(x/(x + 1),n) are the row polynomials of A187075.
The polynomials x^(n+1) * P((x + 1)/x,n) are the row polynomials of A186695.
Row sums are A001147(n+1). (End)
Sum_{k = 0..n} (-1)^k*T(n,k) = (-1)^binomial(n,2)*A012259(n+1). - Johannes W. Meijer, Sep 27 2011

Extensions

Minor edits from Johannes W. Meijer, Sep 27 2011
Showing 1-10 of 18 results. Next