The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A005425 a(n) = 2*a(n-1) + (n-1)*a(n-2).
(Formerly M1461)
33

%I M1461 #178 Jan 25 2024 08:36:39

%S 1,2,5,14,43,142,499,1850,7193,29186,123109,538078,2430355,11317646,

%T 54229907,266906858,1347262321,6965034370,36833528197,199037675054,

%U 1097912385851,6176578272782,35409316648435,206703355298074,1227820993510153,7416522514174082

%N a(n) = 2*a(n-1) + (n-1)*a(n-2).

%C Switchboard problem with n subscribers, where a subscriber who is not talking can be of either of two sexes. Subscribers who are talking cannot be distinguished by sex. See also A000085. - _Karol A. Penson_, Apr 15 2004

%C _John W. Layman_ observes that computationally this agrees with the binomial transform of A000085.

%C Number of self-inverse partial permutations.

%C Number of '12-3 and 214-3'-avoiding permutations.

%C Number of matchings of the corona K'(n) of the complete graph K(n) and the complete graph K(1); in other words, K'(n) is the graph constructed from K(n) by adding for each vertex v a new vertex v' and the edge vv'. Example: a(3)=14 because in the graph with vertex set {A,B,C,a,b,c} and edge set {AB,AC,BC,Aa,Bb,Cc} we have the following matchings: (i) the empty set (1); (ii) the edges as singletons (6); (iii) {Aa,BC},{Bb,AC},{Cc,AB},{Aa,Bb},{Aa,Cc}, {Bb,Cc} (6); (iv) {Aa,Bb,Cc} (1). Row sums of A100862. - _Emeric Deutsch_, Jan 10 2005

%C Consider finite sequences of positive integers <b(m)> of length n with b(1)=1 and with the constraint that b(m) <= max_{0<k<m} b(k)+k-m+2. The question is how many such sequences there are. (Note that when we consider only the term k=m-1, this becomes b(m) <= b(m-1)+1 and it is well known that the number of sequences under this constraint is the Catalan numbers.) This sequence starts (from n = 1): 1,2,5,14,43,142,499,1850,7193. This appears to be the present sequence. But I do not see any way to prove it. The number T(n,m) of sequences of length n which will limit the continuation to size n+1 to a maximum value of m+1 appears to be given by A111062. - _Franklin T. Adams-Watters_, Dec 21 2005, corrected Dec 31 2014

%C Number of n X n symmetric binary matrices with no row sum greater than 1. - _R. H. Hardin_, Jun 13 2008

%C Polynomials in A099174 evaluated at x=2 (see also formula by Deutsch below). - _Johannes W. Meijer_, Feb 04 2010

%C Equals eigensequence of triangle A128227. Example: a(5) = 142 = (1, 1, 2, 5, 14, 43) dot (1, 2, 3, 4, 5, 1) = (1 + 2 + 6 + 20 + 70 + 43); where (1, 2, 3, 4, 5, 1) = row 5 of triangle A128227. - _Gary W. Adamson_, Aug 27 2010

%C Number of words [d(1), d(2), ..., d(n)] where d(k) is either =0, or =k (a fixed point), or the only value repeating a previous fixed point, see example. - _Joerg Arndt_, Apr 18 2014

%C From _Robert A. Russell_, Apr 28 2018: (Start)

%C Stirling transform of this sequence is A002872;

%C Stirling transform of A005425(n-1) is A080337. (End)

%C Number of congruence orbits of upper-triangular n X n matrices on symmetric matrices, or the number of Borel orbits in largest sect of the type CI symmetric space Sp_{2n}(C)/GL_n(C). - _Aram Bingham_, Oct 10 2019

%C For a refined enumeration of the switchboard scenario presented by Penson above and in Donaghey and its relation to perfect matchings of simplices and an operator calculus, see A344678. - _Tom Copeland_, May 29 2021

%C Write [n] for {1, ..., n} and [n]^(k) for k-tuples without repeated entries. Then C [n]^(k) is naturally a complex S_n-representation, whose length is a(k) provided that n >= 2k. a(k) also gives the length of the countable dimensional Sym(N)-representation C N^(k), as remarked by Sam and Snowden (see link). - _Jingjie Yang_, Dec 28 2023

%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

%H T. D. Noe, <a href="/A005425/b005425.txt">Table of n, a(n) for n = 0..200</a>

%H E. Bagno and Y. Cherniavsky, <a href="https://doi.org/10.1016/j.disc.2011.12.018">Congruence B-orbits and the Bruhat poset of involutions of the symmetric group</a>, Discrete Math., 312 (1) (2012), pp. 1289-1299.

%H C. Banderier, M. Bousquet-Mélou, A. Denise, P. Flajolet, D. Gardy and D. Gouyou-Beauchamps, <a href="http://dx.doi.org/10.1016/S0012-365X(01)00250-3">Generating Functions for Generating Trees</a>, Discrete Mathematics 246(1-3), March 2002, pp. 29-55.

%H A. Bingham and O. Ugurlu, <a href="https://arxiv.org/abs/1903.07229">Sects and lattice paths over the Lagrangian Grassmannian</a>, arXiv preprint arXiv:1903.07229 [math.CO], 2019.

%H P. Blasiak, K. A. Penson, A. I. Solomon, A. Horzela and G. E. H. Duchamp, <a href="https://arXiv.org/abs/quant-ph/0405103">Combinatorial field theories via boson normal ordering</a>, arXiv:quant-ph/0405103, 2004-2006.

%H P. Blasiak, K. A. Penson, A. I. Solomon, A. Horzela and G. E. H. Duchamp, <a href="http://dx.doi.org/10.1063/1.1904161">Some useful combinatorial formulas for bosonic operators</a>, J. Math. Phys. 46, 052110 (2005) (6 pages).

%H R. Donaghey, <a href="http://dx.doi.org/10.1016/0097-3165(76)90060-1">Binomial self-inverse sequences and tangent coefficients</a>, J. Combin. Theory, Series A, 21 (1976), 155-163.

%H Adam M. Goyt and Lara K. Pudwell, <a href="https://arxiv.org/abs/1203.3786">Avoiding colored partitions of two elements in the pattern sense</a>, arXiv preprint arXiv:1203.3786 [math.CO], 2012. - From _N. J. A. Sloane_, Sep 17 2012

%H T. Halverson and M. Reeks, <a href="https://arxiv.org/abs/1302.6150">Gelfand Models for Diagram Algebras</a>, arXiv preprint arXiv:1302.6150 [math.RT], 2013-2014.

%H J. W. Layman, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL4/LAYMAN/hankel.html">The Hankel Transform and Some of its Properties</a>, J. Integer Sequences, 4 (2001), #01.1.5.

%H Mikhail Khovanov, Victor Ostrik and Yakov Kononov, <a href="https://arxiv.org/abs/2011.14758">Two-dimensional topological theories, rational functions and their tensor envelopes</a>, arXiv:2011.14758 [math.QA], 2020.

%H T. Mansour, <a href="https://arxiv.org/abs/math/0202219">Restricted permutations by patterns of type 2-1</a>, arXiv:math/0202219 [math.CO], 2002.

%H D. Panyushev, <a href="https://arxiv.org/abs/1407.6857">On the orbits of a Borel subgroup in abelian ideals</a>, arXiv preprint arXiv:1407.6857 [math.AG], 2014.

%H Ville H. Pettersson, <a href="https://doi.org/10.37236/4510">Enumerating Hamiltonian Cycles</a>, The Electronic Journal of Combinatorics, Volume 21, Issue 4, 2014.

%H K. Piejko, <a href="http://dx.doi.org/10.1155/2015/463650">Extremal trees with respect to number of (A, B, 2C)-edge colourings</a>, Journal of Applied Mathematics, Hindawi Publishing Corporation, Volume 2015, Article ID 463650, 5 pages.

%H Steven V Sam and Andrew Snowden, <a href="https://arxiv.org/abs/1302.5859">Stability patterns in representation theory</a>, arXiv:1302.5859 [math.RT], 2013-2015. See (1.3.4) p. 8.

%H R. P. Stanley, <a href="http://dx.doi.org/10.1016/S0196-8858(02)00537-7">On the enumeration of skew Young tableaux</a>, Advances in Applied Mathematics 30 (2003) 283-294, (see corollary 2.4).

%F E.g.f.: exp( 2*x + x^2/2 ).

%F a(n) = A027412(n+1)/2. - _N. J. A. Sloane_, Sep 13 2003

%F a(n) = Sum_{k=0..n} binomial(n, k)*A000085(n-k). - _Philippe Deléham_, Mar 07 2004

%F a(n) = (-i/sqrt(2))^n*H(n, i*sqrt(2)), where H(n, x) is the n-th Hermite polynomial and i = sqrt(-1). - _Emeric Deutsch_, Nov 24 2004

%F a(n) = Sum_{k=0..floor(n/2)} 2^{n-3*k}*n!/((n-2*k)!*k!). - Huajun Huang (hua_jun(AT)hotmail.com), Oct 10 2005

%F For all n, a(n) = [M_n]_1,1 = [M_n]_2,1, where M_n = A_n * A_n-1 * ... * A_1, being A_k the matrix A_k = [1, k;1, 1]. - _Simone Severini_, Apr 25 2007

%F a(n) = (1/sqrt(2*Pi))*Integral_{x=-infinity..infinity} exp(-x^2/2)*(x+2)^n. - _Groux Roland_, Mar 14 2011

%F G.f.: 1/(1-2x-x^2/(1-2x-2x^2/(1-2x-3x^2/(1-2x-4x^2/(1-... (continued fraction).

%F E.g.f.: G(0) where G(k) = 1 + x*(4+x)/(4*k + 2 - x*(4+x)*(4*k+2)/(x*(4+x) + 4*(k+1)/G(k+1))); (continued fraction). - _Sergei N. Gladkovskii_, Dec 17 2011

%F a(n) ~ exp(2*sqrt(n) - n/2 - 1)*n^(n/2)/sqrt(2). - _Vaclav Kotesovec_, Oct 08 2012

%F a(n) = 2^(n/2)*exp(-i*n*Pi/2)*KummerU(-n/2, 1/2, -2). - _Peter Luschny_, May 30 2021

%e From _Joerg Arndt_, Apr 18 2014: (Start)

%e The a(3) = 14 words [d(1), d(2), d(3)] where d(k) is either =0, or =k (a fixed point), or the only value repeating a previous fixed point are (dots for zeros):

%e # : word partial involution

%e 01: [ . . . ] ()

%e 02: [ . . 3 ] (3)

%e 03: [ . 2 . ] (2)

%e 04: [ . 2 2 ] (2 3)

%e 05: [ . 2 3 ] (2) (3)

%e 06: [ 1 . . ] (1)

%e 07: [ 1 . 1 ] (1 3)

%e 08: [ 1 . 3 ] (1) (3)

%e 09: [ 1 1 . ] (1 2)

%e 10: [ 1 1 3 ] (1 2) (3)

%e 11: [ 1 2 . ] (1) (2)

%e 12: [ 1 2 1 ] (1 3) (2)

%e 13: [ 1 2 2 ] (1) (2 3)

%e 14: [ 1 2 3 ] (1) (2) (3)

%e (End)

%p with(orthopoly): seq((-I/sqrt(2))^n*H(n,I*sqrt(2)),n=0..25);

%t a[0] = 1; a[1] = 2; a[n_]:= 2a[n-1] + (n-1)*a[n-2]; Table[ a[n], {n, 0, 25}] (* or *)

%t Range[0, 25]!CoefficientList[Series[Exp[2 x + x^2/2], {x, 0, 25}], x] (* or *)

%t f[n_] := Sum[2^(n - 3k)n!/((n - 2k)!k!), {k, 0, n}]; Table[ f[n], {n, 0, 25}] (* or *)

%t Table[(-I/Sqrt[2])^n*HermiteH[n, I*Sqrt[2]], {n, 0, 25}] (* _Robert G. Wilson v_, Nov 04 2005 *)

%t RecurrenceTable[{a[0]==1,a[1]==2,a[n]==2a[n-1]+(n-1)a[n-2]},a,{n,30}] (* _Harvey P. Dale_, Sep 30 2015 *)

%t a[n_] := 2^(n/2) Exp[- I n Pi/2] HypergeometricU[-n/2, 1/2, -2];

%t Table[a[n], {n, 0, 22}] (* _Peter Luschny_, May 30 2021 *)

%o (Haskell)

%o a005425 n = a005425_list !! n

%o a005425_list = 1 : 2 : zipWith (+)

%o (map (* 2) (tail a005425_list)) (zipWith (*) [1..] a005425_list)

%o -- _Reinhard Zumkeller_, Dec 18 2011

%o (PARI) A005425(n)=sum(k=0,n\2,2^(n-3*k)*n!/(n-2*k)!/k!) \\ _M. F. Hasler_, Jan 13 2012

%o (Maxima) makelist((%i/sqrt(2))^n*hermite(n,-%i*sqrt(2)),n,0,12); /* _Emanuele Munarini_, Mar 02 2016 */

%o (Magma) a:=[2,5];[1] cat [n le 2 select a[n] else 2*Self(n-1) + (n-1)*Self(n-2):n in [1..30]]; // _Marius A. Burtea_, Oct 10 2019

%o (SageMath) [(-i/sqrt(2))^n*hermite(n, i*sqrt(2)) for n in range(41)] # _G. C. Greubel_, Nov 19 2022

%Y Cf. A002872, A080337, A085483, A100862, A111062, A128227.

%Y Bisections give A093620, A100510.

%Y Row sums of A344678.

%K nonn,easy,nice

%O 0,2

%A _Simon Plouffe_

%E Recurrence and formula corrected by _Olivier Gérard_, Oct 15 1997

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 13 03:04 EDT 2024. Contains 372497 sequences. (Running on oeis4.)