This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A181545 G.f.: A(x) = Sum_{n>=0} (3n)!/(n!)^3 * x^(3n)/(1-x-x^2)^(3n+1). 9
 1, 1, 2, 9, 29, 92, 343, 1281, 4720, 17899, 68933, 266364, 1037423, 4072439, 16065148, 63658521, 253356763, 1012049086, 4055596343, 16299779331, 65683233938, 265310551667, 1073968967929, 4355988107100, 17699727361051 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,3 COMMENTS Limit_{n->oo} a(n+1)/a(n) = (Fibonacci(3)*sqrt(5) + Lucas(3))/2 = sqrt(5) + 2. REFERENCES C. Banderier, P. Hitczenko, Enumeration and asymptotics of restricted compositions having the same number of parts, Disc. Appl. Math. 160 (2012) 2542-2554 doi:10.1016/j.dam.2011.12.011, Proposition 3.1 LINKS Iain Fox, Table of n, a(n) for n = 0..1600 Steffen Eger, On the Number of Many-to-Many Alignments of N Sequences, arXiv:1511.00622 [math.CO], 2015. Edyta Hetmaniok, Barbara Smoleń, Roman Wituła, The Stirling triangles, Proceedings of the Symposium for Young Scientists in Technology, Engineering and Mathematics (SYSTEM 2017), Kaunas, Lithuania, April 28, 2017, p. 35-41. FORMULA a(n) = Sum_{k=0..floor(n/2)} binomial(n-k,k)^3. G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k)^3 * x^k. G.f.: A(x) = G( x^3/(1-x-x^2)^3 )/(1-x-x^2) where G(x) satisfies: * G(x^3) = G( x*(1+3*x+9*x^2)/(1+6*x)^3 )/(1+6*x) and G(x) is the g.f. of A006480. Recurrence: (n-3)*n^2*a(n) = (n-3)*(3*n^2 - 3*n + 1)*a(n-1) - (n-1)*a(n-2) + 2*(n-2)*(11*n^2 - 44*n + 34)*a(n-3) + (n-3)*a(n-4) + (n-1)*(3*n^2 - 21*n + 37)*a(n-5) + (n-4)^2*(n-1)*a(n-6). - Vaclav Kotesovec, Jul 31 2014 a(n) ~ sqrt((9+4*sqrt(5))/12) * (2+sqrt(5))^n / (Pi*n). - Vaclav Kotesovec, Jul 31 2014 EXAMPLE G.f. A(x) = 1 + x + 2*x^2 + 9*x^3 + 29*x^4 + 92*x^5 + 343*x^6 + ... which equals the series: A(x) = 1/(1-x-x^2) + 3!/1!^3*x^3/(1-x-x^2)^4 + 6!/2!^3*x^6/(1-x-x^2)^7 + 9!/3!^3*x^9/(1-x-x^2)^10 + 12!/4!^3*x^12/(1-x-x^2)^13 + ... The g.f. also equals the series: A(x) = 1 + x*(1 + x) + x^2*(1 + 2^3*x + x^2) + x^3*(1 + 3^3*x + 3^3*x^2 + x^3) + x^4*(1 + 4^3*x + 6^3*x^2 + 4^3*x^3 + x^4) + x^5*(1 + 5^3*x + 10^3*x^2 + 10^3*x^3 + 5^3*x^4 + x^5) + ... The terms begin: a(0) = a(1) = 1^3; a(2) = 1^3 + 1^3 = 2; a(3) = 1^3 + 2^3 = 9; a(4) = 1^3 + 3^3 + 1^3 = 29; a(5) = 1^3 + 4^3 + 3^3 = 92; a(6) = 1^3 + 5^3 + 6^3 + 1^3 = 343; a(7) = 1^3 + 6^3 + 10^3 + 4^3 = 1281; ... MATHEMATICA Table[Sum[Binomial[n-k, k]^3, {k, 0, Floor[n/2]}], {n, 0, 20}] (* Vaclav Kotesovec, Jul 31 2014 *) PROG (PARI) {a(n)=sum(k=0, n\2, binomial(n-k, k)^3)} (PARI) {a(n)=polcoeff(sum(m=0, n+1, x^m*sum(k=0, m, binomial(m, k)^3*x^k + x*O(x^n))), n)} (PARI) {a(n)=polcoeff(sum(m=0, n, x^(3*m)/(1-x-x^2+x*O(x^n))^(3*m+1)*(3*m)!/(m!)^3), n)} CROSSREFS Cf. A006480, A217421, A181546, A051286. Sequence in context: A062452 A265440 A152891 * A069006 A241774 A183383 Adjacent sequences:  A181542 A181543 A181544 * A181546 A181547 A181548 KEYWORD nonn AUTHOR Paul D. Hanna, Oct 29 2010 STATUS approved

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

Last modified October 23 08:31 EDT 2018. Contains 316522 sequences. (Running on oeis4.)