Displaying 1-10 of 64 results found.
G.f.: A(x) = exp( Sum_{n>=1} A005260(n)*x^n/n ) where A005260(n) = Sum_{k=0..n} C(n,k)^4.
+20
6
1, 2, 11, 74, 621, 5850, 60212, 659712, 7583514, 90494068, 1112755389, 14022849582, 180362150901, 2360201899690, 31344689243344, 421621652965160, 5734850816825046, 78773961705345324, 1091497852618784390
FORMULA
a(n) ~ c * 16^n / n^(5/2), where c = 0.30919827904959014083681667605470681109347914449671378054261267779... - Vaclav Kotesovec, Nov 27 2017
EXAMPLE
G.f.: A(x) = 1 + 2*x + 11*x^2 + 74*x^3 + 621*x^4 + 5850*x^5 + 60212*x^6 +...
log(A(x)) = 2*x + 18*x^2/2 + 164*x^3/3 + 1810*x^4/4 + 21252*x^5/5 + 263844*x^6/6 + 3395016*x^7/7 +...+ A005260(n)*x^n/n +...
MATHEMATICA
a[n_] := Sum[(Binomial[n, k])^4, {k, 0, n}]; f[x_] := Sum[a[n]*x^n/(n), {n, 1, 75}]; CoefficientList[Series[Exp[f[x]], {x, 0, 50}], x] (* G. C. Greubel, May 30 2016 *)
PROG
(PARI) {a(n)=polcoeff(exp(sum(m=1, n, sum(k=0, m, binomial(m, k)^4)*x^m/m)+x*O(x^n)), n)}
G.f.: A(x) = exp( Sum_{n>=1} A005260(n)*x^n/(2*n) ), where A005260(n) = Sum_{k=0..n} C(n,k)^4.
+20
5
1, 1, 5, 32, 266, 2499, 25765, 283084, 3264502, 39077898, 481942608, 6089941550, 78523226064, 1029859481949, 13704960309415, 184688556173542, 2516342539576510, 34617557176739174, 480336524752492608
FORMULA
a(n) ~ c * 16^n / n^(5/2), where c = 0.14011467789446087641913961305130549385534145578464604013551918158... - Vaclav Kotesovec, Nov 27 2017
EXAMPLE
G.f.: A(x) = 1 + x + 5*x^2 + 32*x^3 + 266*x^4 + 2499*x^5 + 25765*x^6 +...
log(A(x)) = x + 9*x^2/2 + 82*x^3/3 + 905*x^4/4 + 10626*x^5/5 + 131922*x^6/6 + 1697508*x^7/7 +...+ A005260(n)/2*x^n/n +...
MATHEMATICA
a[n_] := Sum[(Binomial[n, k])^4, {k, 0, n}]; f[x_] := Sum[a[n]*x^n/(2*n), {n, 1, 75}]; CoefficientList[Series[Exp[f[x]], {x, 0, 50}], x] (* G. C. Greubel, May 30 2016 *)
PROG
(PARI) {a(n)=polcoeff(exp(sum(m=1, n, sum(k=0, m, binomial(m, k)^4)/2*x^m/m)+x*O(x^n)), n)}
Least prime divisor of A005260(n) which does not divide any previous term A005260(k) with k < n, or 1 if such a primitive prime divisor of A005260(n) does not exist.
+20
4
2, 3, 41, 5, 7, 349, 61, 75617, 31, 13, 499, 643897693, 17, 19, 1729774061, 101, 2859112064587, 138407, 83, 167, 59, 29, 653, 257, 997540809461453561581, 347, 13679, 37, 160449179727717672892660463, 211, 151, 43, 97, 73, 47
COMMENTS
Conjecture: a(n) is prime for any n > 0. In general, for any r > 2, if n is large enough then f_r(n) = sum_{k=0..n}C(n,k)^r has a prime divisor which does not divide any previous terms f_r(k) with k < n.
MATHEMATICA
u[n_]:=Sum[Binomial[n, k]^4, {k, 0, n}]
f[n_]:=FactorInteger[u[n]]
p[n_]:=Table[Part[Part[f[n], k], 1], {k, 1, Length[f[n]]}]
Do[If[u[n]<2, Goto[cc]]; Do[Do[If[Mod[u[i], Part[p[n], k]]==0, Goto[aa]], {i, 1, n-1}]; Print[n, " ", Part[p[n], k]]; Goto[bb]; Label[aa]; Continue, {k, 1, Length[p[n]]}]; Label[cc]; Print[n, " ", 1]; Label[bb]; Continue, {n, 1, 35}]
G.f. A(x) satisfies: coefficient of x^n in A(x)^(2*n) equals A005260(n) = Sum_{k=0..n} C(n,k)^4.
+20
1
1, 1, 3, 9, 60, 417, 3430, 29927, 278316, 2693437, 26976407, 277394148, 2916106328, 31220964707, 339508802940, 3741551907530, 41714692453164, 469827584596185, 5339334757945439, 61165396353689573, 705720529604453193, 8195208178337460065, 95724512701573485819, 1124070800784913396731
FORMULA
G.f.: sqrt( x / Series_Reversion( x*exp( Sum_{n>=1} A005260(n)*x^n/n ) ) ), where A005260(n) = Sum_{k=0..n} C(n,k)^4. a(n) ~ c * d^n / n^(5/2), where d= 13.142352254618115022093263384837224..., c = 0.051491668112404252102416729094836... . - Vaclav Kotesovec, Jun 05 2014
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 9*x^3 + 60*x^4 + 417*x^5 + 3430*x^6 +...
Form a table of coefficients in A(x)^(2*n) for n>=0, which begins:
[1, 0, 0, 0, 0, 0, 0, 0, 0, ...];
[1, 2, 7, 24, 147, 1008, 8135, 70296, 648172, ...];
[1, 4, 18, 76, 439, 2940, 22936, 194300, 1761411, ...];
[1, 6, 33, 164, 960, 6378, 48526, 403440, 3598050, ...];
[1, 8, 52, 296, 1810, 12128, 90972, 744656, 6542519, ...];
[1, 10, 75, 480, 3105, 21252, 158845, 1286240, 11157705, ...];
[1, 12, 102, 724, 4977, 35100, 263844, 2125020, 18253680, ...];
[1, 14, 133, 1036, 7574, 55342, 421484, 3395016, 28975933, ...];
[1, 16, 168, 1424, 11060, 84000, 651848, 5277696, 44916498, ...]; ...
then the main diagonal forms A005260(n) = Sum_{k=0..n} C(n,k)^4.
PROG
(PARI) {a(n)=polcoeff(sqrt(x/serreverse(x*exp(sum(m=1, n+1, sum(k=0, m, binomial(m, k)^4)*x^m/m +x^2*O(x^n))))), n)}
for(n=0, 30, print1(a(n), ", "))
Central binomial coefficients: binomial(2*n,n) = (2*n)!/(n!)^2.
(Formerly M1645 N0643)
+10
1082
1, 2, 6, 20, 70, 252, 924, 3432, 12870, 48620, 184756, 705432, 2704156, 10400600, 40116600, 155117520, 601080390, 2333606220, 9075135300, 35345263800, 137846528820, 538257874440, 2104098963720, 8233430727600, 32247603683100, 126410606437752, 495918532948104, 1946939425648112
COMMENTS
Devadoss refers to these numbers as type B Catalan numbers (cf. A000108). Equal to the binomial coefficient sum Sum_{k=0..n} binomial(n,k)^2.
Number of possible interleavings of a program with n atomic instructions when executed by two processes. - Manuel Carro (mcarro(AT)fi.upm.es), Sep 22 2001
Convolving a(n) with itself yields A000302, the powers of 4. - T. D. Noe, Jun 11 2002 Number of ordered trees with 2n+1 edges, having root of odd degree and nonroot nodes of outdegree 0 or 2. - Emeric Deutsch, Aug 02 2002 Also number of directed, convex polyominoes having semiperimeter n+2.
Also number of diagonally symmetric, directed, convex polyominoes having semiperimeter 2n+2. - Emeric Deutsch, Aug 03 2002 The second inverse binomial transform of this sequence is this sequence with interpolated zeros. Its g.f. is (1 - 4*x^2)^(-1/2), with n-th term C(n,n/2)(1+(-1)^n)/2. - Paul Barry, Jul 01 2003 Number of possible values of a 2n-bit binary number for which half the bits are on and half are off. - Gavin Scott (gavin(AT)allegro.com), Aug 09 2003
Ordered partitions of n with zeros to n+1, e.g., for n=4 we consider the ordered partitions of 11110 (5), 11200 (30), 13000 (20), 40000 (5) and 22000 (10), total 70 and a(4)=70. See A001700 (esp. Mambetov Bektur's comment). - Jon Perry, Aug 10 2003 Number of nondecreasing sequences of n integers from 0 to n: a(n) = Sum_{i_1=0..n} Sum_{i_2=i_1..n}...Sum_{i_n=i_{n-1}..n}(1). - J. N. Bearden (jnb(AT)eller.arizona.edu), Sep 16 2003
Number of peaks at odd level in all Dyck paths of semilength n+1. Example: a(2)=6 because we have U*DU*DU*D, U*DUUDD, UUDDU*D, UUDUDD, UUU*DDD, where U=(1,1), D=(1,-1) and * indicates a peak at odd level. Number of ascents of length 1 in all Dyck paths of semilength n+1 (an ascent in a Dyck path is a maximal string of up steps). Example: a(2)=6 because we have uDuDuD, uDUUDD, UUDDuD, UUDuDD, UUUDDD, where an ascent of length 1 is indicated by a lower case letter. - Emeric Deutsch, Dec 05 2003 a(n-1) = number of subsets of 2n-1 distinct elements taken n at a time that contain a given element. E.g., n=4 -> a(3)=20 and if we consider the subsets of 7 taken 4 at a time with a 1 we get (1234, 1235, 1236, 1237, 1245, 1246, 1247, 1256, 1257, 1267, 1345, 1346, 1347, 1356, 1357, 1367, 1456, 1457, 1467, 1567) and there are 20 of them. - Jon Perry, Jan 20 2004 The dimension of a particular (necessarily existent) absolutely universal embedding of the unitary dual polar space DSU(2n,q^2) where q>2. - J. Taylor (jt_cpp(AT)yahoo.com), Apr 02 2004.
Number of standard tableaux of shape (n+1, 1^n). - Emeric Deutsch, May 13 2004 Erdős, Graham et al. conjectured that a(n) is never squarefree for sufficiently large n (cf. Graham, Knuth, Patashnik, Concrete Math., 2nd ed., Exercise 112). Sárközy showed that if s(n) is the square part of a(n), then s(n) is asymptotically (sqrt(2)-2)*(sqrt(n))*(Riemann Zeta Function(1/2)). Granville and Ramare proved that the only squarefree values are a(1)=2, a(2)=6 and a(4)=70. - Jonathan Vos Post, Dec 04 2004 [For more about this conjecture, see A261009. - N. J. A. Sloane, Oct 25 2015] The MathOverflow link contains the following comment (slightly edited): The Erdős squarefree conjecture (that a(n) is never squarefree for n>4) was proved in 1980 by Sárközy, A. (On divisors of binomial coefficients. I. J. Number Theory 20 (1985), no. 1, 70-80.) who showed that the conjecture holds for all sufficiently large values of n, and by A. Granville and O. Ramaré (Explicit bounds on exponential sums and the scarcity of squarefree binomial coefficients. Mathematika 43 (1996), no. 1, 73-107) who showed that it holds for all n>4. - Fedor Petrov, Nov 13 2010. [From N. J. A. Sloane, Oct 29 2015] p divides a((p-1)/2)-1= A030662(n) for prime p=5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97, ... = A002144(n) Pythagorean primes: primes of form 4n+1. - Alexander Adamchuk, Jul 04 2006 The number of direct routes from my home to Granny's when Granny lives n blocks south and n blocks east of my home in Grid City. To obtain a direct route, from the 2n blocks, choose n blocks on which one travels south. For example, a(2)=6 because there are 6 direct routes: SSEE, SESE, SEES, EESS, ESES and ESSE. - Dennis P. Walsh, Oct 27 2006 Inverse: With q = -log(log(16)/(pi a(n)^2)), ceiling((q + log(q))/log(16)) = n. - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Feb 26 2007
Number of partitions with Ferrers diagrams that fit in an n X n box (including the empty partition of 0). Example: a(2) = 6 because we have: empty, 1, 2, 11, 21 and 22. - Emeric Deutsch, Oct 02 2007 The number of walks of length 2n on an infinite linear lattice that begins and ends at the origin. - Stefan Hollos (stefan(AT)exstrom.com), Dec 10 2007
The number of lattice paths from (0,0) to (n,n) using steps (1,0) and (0,1). - Joerg Arndt, Jul 01 2011 Integral representation: C(2n,n)=1/Pi Integral [(2x)^(2n)/sqrt(1 - x^2),{x,-1, 1}], i.e., C(2n,n)/4^n is the moment of order 2n of the arcsin distribution on the interval (-1,1). - N-E. Fahssi, Jan 02 2008 Straub, Amdeberhan and Moll: "... it is conjectured that there are only finitely many indices n such that C_n is not divisible by any of 3, 5, 7 and 11." - Jonathan Vos Post, Nov 14 2008 Also, in sports, the number of ordered ways for a "Best of 2n-1 Series" to progress. For example, a(2) = 6 means there are six ordered ways for a "best of 3" series to progress. If we write A for a win by "team A" and B for a win by "team B" and if we list the played games chronologically from left to right then the six ways are AA, ABA, BAA, BB, BAB, and ABB. (Proof: To generate the a(n) ordered ways: Write down all a(n) ways to designate n of 2n games as won by team A. Remove the maximal suffix of identical letters from each of these.) - Lee A. Newberg, Jun 02 2009 Number of n X n binary arrays with rows, considered as binary numbers, in nondecreasing order, and columns, considered as binary numbers, in nonincreasing order. - R. H. Hardin, Jun 27 2009 Hankel transform is 2^n. - Paul Barry, Aug 05 2009 It appears that a(n) is also the number of quivers in the mutation class of twisted type BC_n for n>=2.
Number of words on {a,b} of length 2n such that no prefix of the word contains more b's than a's. - Jonathan Nilsson, Apr 18 2012 From Pascal's triangle take row(n) with terms in order a1,a2,..a(n) and row(n+1) with terms b1,b2,..b(n), then 2*(a1*b1 + a2*b2 + ... + a(n)*b(n)) to get the terms in this sequence. - J. M. Bergot, Oct 07 2012. For example using rows 4 and 5: 2*(1*(1) + 4*(5) + 6*(10) + 4*(10) + 1*(5)) = 252, the sixth term in this sequence. Take from Pascal's triangle row(n) with terms b1, b2, ..., b(n+1) and row(n+2) with terms c1, c2, ..., c(n+3) and find the sum b1*c2 + b2*c3 + ... + b(n+1)*c(n+2) to get A000984(n+1). Example using row(3) and row(5) gives sum 1*(5)+3*(10)+3*(10)+1*(5) = 70 = A000984(4). - J. M. Bergot, Oct 31 2012 a(n) == 2 mod n^3 iff n is a prime > 3. (See Mestrovic link, p. 4.) - Gary Detlefs, Feb 16 2013 Conjecture: For any positive integer n, the polynomial sum_{k=0}^n a(k)x^k is irreducible over the field of rational numbers. In general, for any integer m>1 and n>0, the polynomial f_{m,n}(x) = Sum_{k=0..n} (m*k)!/(k!)^m*x^k is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 23 2013 This comment generalizes the comment dated Oct 31 2012 and the second of the sequence's original comments. For j = 1 to n, a(n) = Sum_{k=0..j} C(j,k)* C(2n-j, n-k) = 2*Sum_{k=0..j-1} C(j-1,k)*C(2n-j, n-k). - Charlie Marion, Jun 07 2013 The differences between consecutive terms of the sequence of the quotients between consecutive terms of this sequence form a sequence containing the reciprocals of the triangular numbers. In other words, a(n+1)/a(n)-a(n)/a(n-1) = 2/(n*(n+1)). - Christian Schulz, Jun 08 2013 Number of distinct strings of length 2n using n letters A and n letters B. - Hans Havermann, May 07 2014 Expansion of G.f. A(x) = 1/(1+q*x*c(x)), where parameter q is positive or negative (except q=-1), and c(x) is the g.f. of A000108 for Catalan numbers. The case of q=-1 recovers the g.f. of A000108 as xA^2-A+1=0. The present sequence A000984 refers to q=-2. Recurrence: (1+q)*(n+2)*a(n+2) + ((q*q-4*q-4)*n + 2*(q*q-q-1))*a(n+1) - 2*q*q*(2*n+1)*a(n) = 0, a(0)=1, a(1)=-q. Asymptotics: a(n) ~ ((q+2)/(q+1))*(q^2/(-q-1))^n, q<=-3, a(n) ~ (-1)^n*((q+2)/(q+1))*(q^2/(q+1))^n, q>=5, and a(n) ~ -Kq*2^(2*n)/sqrt(Pi*n^3), where the multiplicative constant Kq is given by K1=1/9 (q=1), K2=1/8 (q=2), K3=3/25 (q=3), K4=1/9 (q=4). These formulas apply to existing sequences A126983 (q=1), A126984 (q=2), A126982 (q=3), A126986 (q=4), A126987 (q=5), A127017 (q=6), A127016 (q=7), A126985 (q=8), A127053 (q=9), and to A007854 (q=-3), A076035 (q=-4), A076036 (q=-5), A127628 (q=-6), A126694 (q=-7), A115970 (q=-8). (End) a(n)*(2^n)^(j-2) equals S(n), where S(n) is the n-th number in the self-convolved sequence which yields the powers of 2^j for all integers j, n>=0. For example, when n=5 and j=4, a(5)=252; 252*(2^5)^(4-2) = 252*1024 = 258048. The self-convolved sequence which yields powers of 16 is {1, 8, 96, 1280, 17920, 258048, ...}; i.e., S(5) = 258048. Note that the convolved sequences will be composed of numbers decreasing from 1 to 0, when j<2 (exception being j=1, where the first two numbers in the sequence are 1 and all others decreasing). - Bob Selcoe, Jul 16 2014 The variance of the n-th difference of a sequence of pairwise uncorrelated random variables each with variance 1. - Liam Patrick Roche, Jun 04 2015 Number of ordered trees with n edges where vertices at level 1 can be of 2 colors. Indeed, the standard decomposition of ordered trees leading to the equation C = 1 + zC^2 (C is the Catalan function), yields this time G = 1 + 2zCG, from where G = 1/sqrt(1-4z). - Emeric Deutsch, Jun 17 2015 Number of monomials of degree at most n in n variables. - Ran Pan, Sep 26 2015 Let V(n, r) denote the volume of an n-dimensional sphere with radius r, then V(n, 2^n) / Pi = V(n-1, 2^n) * a(n/2) for all even n. - Peter Luschny, Oct 12 2015 a(n) is the number of sets {i1,...,in} of length n such that n >= i1 >= i2 >= ... >= in >= 0. For instance, a(2) = 6 as there are only 6 such sets: (2,2) (2,1) (2,0) (1,1) (1,0) (0,0). - Anton Zakharov, Jul 04 2016 By analytic continuation to the entire complex plane there exist regularized values for divergent sums such as:
Sum_{k>=0} a(k)/(-2)^k = 1/sqrt(3).
Sum_{k>=0} a(k)/(-1)^k = 1/sqrt(5).
Sum_{k>=0} a(k)/(-1/2)^k = 1/3.
Sum_{k>=0} a(k)/(1/2)^k = -1/sqrt(7)i.
Sum_{k>=0} a(k)/(1)^k = -1/sqrt(3)i.
Sum_{k>=0} a(k)/2^k = -i. (End)
Number of sequences (e(1), ..., e(n+1)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) > e(j). [Martinez and Savage, 2.18] - Eric M. Schmidt, Jul 17 2017 The o.g.f. for the sequence equals the diagonal of any of the following the rational functions: 1/(1 - (x + y)), 1/(1 - (x + y*z)), 1/(1 - (x + x*y + y*z)) or 1/(1 - (x + y + y*z)). - Peter Bala, Jan 30 2018 Let s denote West's stack-sorting map. a(n) is the number of permutations pi of [n+1] such that s(pi) avoids the patterns 132, 231, and 321. a(n) is also the number of permutations pi of [n+1] such that s(pi) avoids the patterns 132, 312, and 321.
a(n) is the number of permutations of [n+1] that avoid the patterns 1342, 3142, 3412, and 3421. (End)
All binary self-dual codes of length 4n, for n>0, must contain at least a(n) codewords of weight 2n. More to the point, there will always be at least one, perhaps unique, binary self-dual code of length 4n that will contain exactly a(n) codewords that have a hamming weight equal to half the length of the code (2n). This code can be constructed by direct summing the unique binary self-dual code of length 2 (up to permutation equivalence) to itself an even number of times. A permutation equivalent code can be constructed by augmenting two identity matrices of length 2n together. - Nathan J. Russell, Nov 25 2018 Let [b/p] denote the Legendre symbol and 1/b denote the inverse of b mod p. Then, for m and n, where n is not divisible by p,
[(m+n)/p] == [n/p]*Sum_{k=0..(p-1)/2} (-m/(4*n))^k * a(k) (mod p).
Evaluating this identity for m = -1 and n = 1 demonstrates that, for all odd primes p, Sum_{k=0..(p-1)/2} (1/4)^k * a(k) is divisible by p. (End)
Number of vertices of the subgraph of the (2n-1)-dimensional hypercube induced by all bitstrings with n-1 or n many 1s. The middle levels conjecture asserts that this graph has a Hamilton cycle. - Torsten Muetze, Feb 11 2019 a(n) is the number of walks of length 2n from the origin with steps (1,1) and (1,-1) that stay on or above the x-axis. Equivalently, a(n) is the number of walks of length 2n from the origin with steps (1,0) and (0,1) that stay in the first octant. - Alexander Burstein, Dec 24 2019 Number of permutations of length n>0 avoiding the partially ordered pattern (POP) {3>1, 1>2} of length 4. That is, number of length n permutations having no subsequences of length 4 in which the first element is larger than the second element but smaller than the third elements. - Sergey Kitaev, Dec 08 2020 Also the number of integer compositions of 2n+1 with alternating sum 1, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. For example, the a(0) = 1 through a(2) = 6 compositions are:
(1) (2,1) (3,2)
(1,1,1) (1,2,2)
(2,2,1)
(1,1,2,1)
(2,1,1,1)
(1,1,1,1,1)
The following relate to these compositions:
- The unordered version is A000070. - The alternating sum 0 version is counted by A088218, ranked by A344619. - Including even indices gives A126869. - The complement is counted by A202736. Equivalently, a(n) counts binary numbers with 2n+1 digits and one more 1 than 0's. For example, the a(2) = 6 binary numbers are: 10011, 10101, 10110, 11001, 11010, 11100.
(End)
a(n) is the number of nx2 Young tableaux with a single horizontal wall between the first and second column. If there is a wall between two cells, the entries may be decreasing; see [Banderier, Wallner 2021].
Example for a(2)=6:
3 4 2 4 3 4 3|4 4|3 2|4
1|2, 1|3, 2|1, 1 2, 1 2, 1 3
a(n) is also the number of nx2 Young tableaux with n "walls" between the first and second column.
Example for a(2)=6:
3|4 2|4 4|3 3|4 4|3 4|2
1|2, 1|3, 1|2, 2|1, 2|1, 3|1 (End)
a(n)/4^n is the probability that a fair coin tossed 2n times will come up heads exactly n times and tails exactly n times, or that a random walk with steps of +-1 will return to the starting point after 2n steps (not necessarily for the first time). As n becomes large, this number asymptotically approaches 1/sqrt(n*Pi), using Stirling's approximation for n!.
a(n)/(4^n*(2n-1)) is the probability that a random walk with steps of +-1 will return to the starting point for the first time after 2n steps. The absolute value of the n-th term of A144704 is denominator of this fraction. Considering all possible random walks of exactly 2n steps with steps of +-1, a(n)/(2n-1) is the number of such walks that return to the starting point for the first time after 2n steps. See the absolute values of A002420 or A284016 for these numbers. For comparison, as mentioned by Stefan Hollos, Dec 10 2007, a(n) is the number of such walks that return to the starting point after 2n steps, but not necessarily for the first time. (End) Also the size of the shuffle product of two words of length n, such that the union of the two words consist of 2n distinct elements. - Robert C. Lyons, Mar 15 2023
REFERENCES
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 160.
Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 575, line -3, with a=b=n.
E. Deutsch and L. Shapiro, Seventeen Catalan identities, Bulletin of the Institute of Combinatorics and its Applications, 31, 31-38, 2001.
H. W. Gould, Combinatorial Identities, Morgantown, 1972, (3.66), page 30.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, Second Ed., see Exercise 112.
M. Griffiths, The Backbone of Pascal's Triangle, United Kingdom Mathematics Trust (2008), 3-124.
Leonard Lipshitz and A. van der Poorten. "Rational functions, diagonals, automata and arithmetic." In Number Theory, Richard A. Mollin, ed., Walter de Gruyter, Berlin (1990): 339-358.
J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
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).
LINKS
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy]. Dennis E. Davenport, Lara K. Pudwell, Louis W. Shapiro, and Leon C. Woodson, The Boundary of Ordered Trees, Journal of Integer Sequences, Vol. 18 (2015), Article 15.5.8. L. W. Shapiro, S. Getu, Wen-Jin Woan and L. C. Woodson, The Riordan Group, Discrete Appl. Maths. 34 (1991) 229-239. M. Wallner, Lattice Path Combinatorics, Diplomarbeit, Institut für Diskrete Mathematik und Geometrie der Technischen Universität Wien, 2013.
FORMULA
a(n)/(n+1) = A000108(n), the Catalan numbers. G.f.: A(x) = (1 - 4*x)^(-1/2) = 1F0(1/2;;4x).
D-finite with recurrence: n*a(n) + 2*(1-2*n)*a(n-1)=0.
a(n) = 2^n/n! * Product_{k=0..n-1} (2*k+1).
Using Stirling's formula in A000142 it is easy to get the asymptotic expression a(n) ~ 4^n / sqrt(Pi * n). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001 Integral representation as n-th moment of a positive function on the interval [0, 4]: a(n) = Integral_{x=0..4}(x^n*((x*(4-x))^(-1/2))/Pi), n=0, 1, ... This representation is unique. - Karol A. Penson, Sep 17 2001 a(n) = Max_{ (i+j)!/(i!j!) | 0<=i,j<=n }. - Benoit Cloitre, May 30 2002 E.g.f.: exp(2*x)*I_0(2x), where I_0 is Bessel function. - Michael Somos, Sep 08 2002 E.g.f.: I_0(2*x) = Sum a(n)*x^(2*n)/(2*n)!, where I_0 is Bessel function. - Michael Somos, Sep 09 2002 Determinant of n X n matrix M(i, j) = binomial(n+i, j). - Benoit Cloitre, Aug 28 2003 Given m = C(2*n, n), let f be the inverse function, so that f(m) = n. Letting q denote -log(log(16)/(m^2*Pi)), we have f(m) = ceiling( (q + log(q)) / log(16) ). - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Oct 30 2003
a(n+1) = Sum_{j=n..n*2+1} binomial(j, n). E.g., a(4) = C(7,3) + C(6,3) + C(5,3) + C(4,3) + C(3,3) = 35 + 20 + 10 + 4 + 1 = 70. - Jon Perry, Jan 20 2004 a(n) = (-1)^(n)*Sum_{j=0..(2*n)} (-1)^j*binomial(2*n, j)^2. - Helena Verrill (verrill(AT)math.lsu.edu), Jul 12 2004
a(n) = Sum_{k=0..n} binomial(2n+1, k)*sin((2n-2k+1)*Pi/2). - Paul Barry, Nov 02 2004 a(n-1) = (1/2)*(-1)^n*Sum_{0<=i, j<=n}(-1)^(i+j)*binomial(2n, i+j). - Benoit Cloitre, Jun 18 2005 G.f.: c(x)^2/(2*c(x)-c(x)^2) where c(x) is the g.f. of A000108. - Paul Barry, Feb 03 2006 G.f.: 1/(1-2x-2x^2/(1-2x-x^2/(1-2x-x^2/(1-2x-x^2/(1-... (continued fraction);
G.f.: 1/(1-2x/(1-x/(1-x/(1-x/(1-... (continued fraction). (End)
Let A(x) be the g.f. and B(x) = A(-x), then B(x) = sqrt(1-4*x*B(x)^2). - Vladimir Kruchinin, Jan 16 2011 a(n) = (-4)^n*sqrt(Pi)/(gamma((1/2-n))*gamma(1+n)). - Gerry Martens, May 03 2011 a(n) = upper left term in M^n, M = the infinite square production matrix:
2, 2, 0, 0, 0, 0, ...
1, 1, 1, 0, 0, 0, ...
1, 1, 1, 1, 0, 0, ...
1, 1, 1, 1, 1, 0, ...
1, 1, 1, 1, 1, 1, ....
a(n) = Hypergeometric([-n,-n],[1],1). - Peter Luschny, Nov 01 2011 G.f.: 1 + 2*x/(U(0)-2*x) where U(k) = 2*(2*k+1)*x + (k+1) - 2*(k+1)*(2*k+3)*x/U(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Jun 28 2012 a(n) = Sum_{k=0..n} binomial(n,k)^2*H(k)/(2*H(n)-H(2*n)), n>0, where H(n) is the n-th harmonic number. - Gary Detlefs, Mar 19 2013 G.f.: Q(0)*(1-4*x), where Q(k) = 1 + 4*(2*k+1)*x/( 1 - 1/(1 + 2*(k+1)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 11 2013 G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013 E.g.f.: E(0)/2, where E(k) = 1 + 1/(1 - 2*x/(2*x + (k+1)^2/(2*k+1)/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013 Special values of Jacobi polynomials, in Maple notation: a(n) = 4^n*JacobiP(n,0,-1/2-n,-1). - Karol A. Penson, Jul 27 2013 a(n) = 2^(4*n)/((2*n+1)*Sum_{k=0..n} (-1)^k*C(2*n+1,n-k)/(2*k+1)). - Mircea Merca, Nov 12 2013 a(n) = C(2*n-1,n-1)*C(4*n^2,2)/(3*n*C(2*n+1,3)), n>0. - Gary Detlefs, Jan 02 2014 0 = a(n)*(16*a(n+1) - 6*a(n+2)) + a(n+1)*(-2*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Sep 17 2014 G.f.: Sum_{n>=0} x^n/(1-x)^(2*n+1) * Sum_{k=0..n} C(n,k)^2 * x^k. - Paul D. Hanna, Nov 08 2014 a(n) = 4^n*hypergeom([-n,1/2],[1],1). - Peter Luschny, May 19 2015 a(n) = Sum_{k=0..floor(n/2)} C(n,k)*C(n-k,k)*2^(n-2*k). - Robert FERREOL, Aug 29 2015 a(n) ~ 4^n*(2-2/(8*n+2)^2+21/(8*n+2)^4-671/(8*n+2)^6+45081/(8*n+2)^8)/sqrt((4*n+1) *Pi). - Peter Luschny, Oct 14 2015 A(-x) = 1/x * series reversion( x*(2*x + sqrt(1 + 4*x^2)) ). Compare with the o.g.f. B(x) of A098616, which satisfies B(-x) = 1/x * series reversion( x*(2*x + sqrt(1 - 4*x^2)) ). See also A214377. - Peter Bala, Oct 19 2015 Sum_{n>=0} (-1)^n/a(n) = 4*(5 - sqrt(5)*log(phi))/25 = 0.6278364236143983844442267..., where phi is the golden ratio. - Ilya Gutkovskiy, Jul 04 2016 This sequence occurs as the closed-form expression for several binomial sums:
a(n) = Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n,k)*binomial(2*n + 1,k).
a(n) = 2*Sum_{k = 0..2*n-1} (-1)^(n+k)*binomial(2*n - 1,k)*binomial(2*n,k) for n >= 1.
a(n) = 2*Sum_{k = 0..n-1} binomial(n - 1,k)*binomial(n,k) for n >= 1.
a(n) = Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(x + k,n)*binomial(y + k,n) = Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(x - k,n)*binomial(y - k,n) for arbitrary x and y.
For m = 3,4,5,... both Sum_{k = 0..m*n} (-1)^k*binomial(m*n,k)*binomial(x + k,n)*binomial(y + k,n) and Sum_{k = 0..m*n} (-1)^k*binomial(m*n,k)*binomial(x - k,n)*binomial(y - k,n) appear to equal Kronecker's delta(n,0).
a(n) = (-1)^n*Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(x + k,n)*binomial(y - k,n) for arbitrary x and y.
For m = 3,4,5,... Sum_{k = 0..m*n} (-1)^k*binomial(m*n,k)*binomial(x + k,n)*binomial(y - k,n) appears to equal Kronecker's delta(n,0).
a(n) = Sum_{k = 0..2n} (-1)^k*binomial(2*n,k)*binomial(3*n - k,n)^2 = Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)* binomial(n + k,n)^2. (Gould, Vol. 7, 5.23).
a(n) = Sum_{k = 0..n} (-1)^(n+k)*binomial(2*n,n + k)*binomial(n + k,n)^2. (End)
Sum_{k>=0} a(k)/(p/q)^k = sqrt(p/(p-4q)) for q in N, p in Z/{-4q< (some p) <-2}.
...
Sum_{k>=0} a(k)/(-4)^k = 1/sqrt(2).
Sum_{k>=0} a(k)/(17/4)^k = sqrt(17).
Sum_{k>=0} a(k)/(18/4)^k = 3.
Sum_{k>=0} a(k)/5^k = sqrt(5).
Sum_{k>=0} a(k)/6^k = sqrt(3).
Sum_{k>=0} a(k)/8^k = sqrt(2).
...
Sum_{k>=0} a(k)/(p/q)^k = sqrt(p/(p-4q)) for p>4q.(End)
Boas-Buck recurrence: a(n) = (2/n)*Sum_{k=0..n-1} 4^(n-k-1)*a(k), n >= 1, a(0) = 1. Proof from a(n) = A046521(n, 0). See a comment there. - Wolfdieter Lang, Aug 10 2017 a(n) = Sum_{k = 0..n} (-1)^(n-k) * binomial(2*n+1, k) for n in N. - Rene Adad, Sep 30 2017 G.f. for {1/a(n)}: 4*(sqrt(4-x) + sqrt(x)*arcsin(sqrt(x)/2)) / (4-x)^(3/2).
E.g.f. for {1/a(n)}: 1 + exp(x/4)*sqrt(Pi*x)*erf(sqrt(x)/2)/2.
Sum_{n>=0} (-1)^n/a(n) = 4*(1/5 - arcsinh(1/2)/(5*sqrt(5))). (End)
a(n) = 2^(2*n)*Product_{k=1..2*n} k^((-1)^(k+1)) = A056040(2*n). a(n) = 4^n*binomial(n-1/2,-1/2) = 4^n*GegenbauerC(n,1/4,1). - Gerry Martens, Oct 19 2022 Occurs on the right-hand side of the binomial sum identities Sum_{k = -n..n} (-1)^k * (n + x - k) * binomial(2*n, n+k)^2 = (x + n)*a(n) and Sum_{k = -n..n} (-1)^k * (n + x - k)^2 * binomial(2*n, n+k)^3 = x*(x + 2*n)*a(n) (x arbitrary). Compare with the identity: Sum_{k = -n..n} (-1)^k * binomial(2*n, n+k)^2 = a(n). - Peter Bala, Jul 31 2023 4^n*a(n) = Sum_{k = 0..2*n} (-1)^k*a(k)*a(2*n-k).
16^n = Sum_{k = 0..2*n} a(k)*a(2*n-k). (End)
a(n) = Sum_{k=0..floor(n/2)} binomial(n,2k)*binomial(2*k,k)*2^(n-2*k). (H. W. Gould) - Gary Detlefs, May 28 2024 a(n) = Sum_{k=0..2*n} (-1)^k*binomial(2n,k)*binomial(2*n+2*k,n+k)*3^(2*n-k). (H. W. Gould) (End)
EXAMPLE
G.f.: 1 + 2*x + 6*x^2 + 20*x^3 + 70*x^4 + 252*x^5 + 924*x^6 + ...
For n=2, a(2) = 4!/(2!)^2 = 24/4 = 6, and this is the middle coefficient of the binomial expansion (a + b)^4 = a^4 + 4a^3b + 6a^2b^2 + 4ab^3 + b^4. - Michael B. Porter, Jul 06 2016
MAPLE
with(combstruct); [seq(count([S, {S=Prod(Set(Z, card=i), Set(Z, card=i))}, labeled], size=(2*i)), i=0..20)];
with(combstruct); [seq(count([S, {S=Sequence(Union(Arch, Arch)), Arch=Prod(Epsilon, Sequence(Arch), Z)}, unlabeled], size=i), i=0..25)];
with(combstruct):bin := {B=Union(Z, Prod(B, B))}: seq (count([B, bin, unlabeled], size=n)*n, n=1..25); # Zerinvary Lajos, Dec 05 2007 A000984List := proc(m) local A, P, n; A := [1, 2]; P := [1];
for n from 1 to m - 2 do P := ListTools:-PartialSums([op(P), 2*P[-1]]);
A := [op(A), 2*P[-1]] od; A end: A000984List(28); # Peter Luschny, Mar 24 2022
MATHEMATICA
CoefficientList[Series[1/Sqrt[1-4x], {x, 0, 25}], x] (* Harvey P. Dale, Mar 14 2011 *)
PROG
(Magma) a:= func< n | Binomial(2*n, n) >; [ a(n) : n in [0..10]];
(PARI) A000984(n)=binomial(2*n, n) \\ much more efficient than (2n)!/n!^2. \\ M. F. Hasler, Feb 26 2014 (PARI) fv(n, p)=my(s); while(n\=p, s+=n); s
(PARI) fv(n, p)=my(s); while(n\=p, s+=n); s
(Haskell)
(Python)
from __future__ import division
for n in range(10**3):
b = b*(4*n+2)//(n+1)
(GAP) List([1..1000], n -> Binomial(2*n, n)); # Muniru A Asiru, Jan 30 2018
CROSSREFS
Bisection of A001405 and of A226302. See also A025565, the same ordered partitions but without all in which are two successive zeros: 11110 (5), 11200 (18), 13000 (2), 40000 (0) and 22000 (1), total 26 and A025565(4)=26. See A261009 for a conjecture about this sequence. The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.
KEYWORD
nonn,easy,core,nice,walk,frac,changed
The Franel number a(n) = Sum_{k = 0..n} binomial(n,k)^3.
(Formerly M1971 N0781)
+10
137
1, 2, 10, 56, 346, 2252, 15184, 104960, 739162, 5280932, 38165260, 278415920, 2046924400, 15148345760, 112738423360, 843126957056, 6332299624282, 47737325577620, 361077477684436, 2739270870994736, 20836827035351596, 158883473753259752, 1214171997616258240
COMMENTS
Cusick gives a general method of deriving recurrences for the r-th order Franel numbers (this is the sequence of third-order Franel numbers), with floor((r+3)/2) terms.
This is the Taylor expansion of a special point on a curve described by Beauville. - Matthijs Coster, Apr 28 2004 An identity of V. Strehl states that a(n) = Sum_{k = 0..n} C(n,k)^2 * binomial(2*k,n). Zhi-Wei Sun conjectured that for every n = 2,3,... the polynomial f_n(x) = Sum_{k = 0..n} binomial(n,k)^2 * binomial(2*k,n) * x^(n-k) is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 21 2013 Conjecture: a(n) == 2 (mod n^3) iff n is prime. - Gary Detlefs, Mar 22 2013 a(p) == 2 (mod p^3) for any prime p since p | C(p,k) for all k = 1,...,p-1. - Zhi-Wei Sun, Aug 14 2013 a(n) is the maximal number of totally mixed Nash equilibria in games of 3 players, each with n+1 pure options. - Raimundas Vidunas, Jan 22 2014 This is one of the Apéry-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017 Diagonal of rational functions 1/(1 - x*y - y*z - x*z - 2*x*y*z), 1/(1 - x - y - z + 4*x*y*z), 1/(1 + y + z + x*y + y*z + x*z + 2*x*y*z), 1/(1 + x + y + z + 2*(x*y + y*z + x*z) + 4*x*y*z). - Gheorghe Coserea, Jul 04 2018 a(n) is the constant term in the expansion of ((1 + x) * (1 + y) + (1 + 1/x) * (1 + 1/y))^n. - Seiichi Manyama, Oct 27 2019 Diagonal of rational function 1 / ((1-x)*(1-y)*(1-z) - x*y*z). - Seiichi Manyama, Jul 11 2020 Named after the Swiss mathematician Jérôme Franel (1859-1939). - Amiram Eldar, Jun 15 2021 It appears that a(n) is equal to the coefficient of (x*y*z)^n in the expansion of (1 + x + y - z)^n * (1 + x - y + z)^n * (1 - x + y + z)^n. Cf. A036917. - Peter Bala, Sep 20 2021
REFERENCES
Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.
Jérôme Franel, On a question of Laisant, Intermédiaire des Mathématiciens, vol 1 1894 pp 45-47
H. W. Gould, Combinatorial Identities, Morgantown, 1972, (X.14), p. 56.
Murray Klamkin, ed., Problems in Applied Mathematics: Selections from SIAM Review, SIAM, 1990; see pp. 148-149.
John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 193.
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).
LINKS
P. Barrucand, A combinatorial identity, Problem 75-4, SIAM Rev., Vol. 17 (1975), p. 168. Solution by D. R. Breach, D. McCarthy, D. Monk and P. E. O'Neil, SIAM Rev., Vol. 18 (1976), p. 303. Juan Pla, Problem H-505, Advanced Problems and Solutions, The Fibonacci Quarterly, Vol. 33, No. 5 (1995), p. 473; Sum Formulae!, Solution to Problem H-505 by Paul S. Bruckman, ibid., Vol. 35, No. 1 (1997), pp. 93-95. Zhi-Wei Sun, Conjectures involving arithmetical sequences, arXiv:1208.2683v9 [math.CO] 2013; Number Theory: Arithmetic in Shangri-La (eds., S. Kanemitsu, H. Li and J. Liu), Proc. the 6th China-Japan Sem. (Shanghai, August 15-17, 2011), World Sci., Singapore, 2013, pp. 244-258.
FORMULA
A002893(n) = Sum_{m = 0..n} binomial(n, m)*a(m) [Barrucand].
Sum_{k = 0..n} C(n, k)^3 = (-1)^n*Integral_{x = 0..infinity} L_k(x)^3 exp(-x) dx. - from Askey's book, p. 43
D-finite with recurrence (n + 1)^2*a(n+1) = (7*n^2 + 7*n + 2)*a(n) + 8*n^2*a(n-1) [Franel]. - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jan 31 2001
a(n) ~ 2*3^(-1/2)*Pi^-1*n^-1*2^(3*n). - Joe Keane (jgk(AT)jgk.org), Jun 21 2002
O.g.f.: A(x) = Sum_{n >= 0} (3*n)!/n!^3 * x^(2*n)/(1 - 2*x)^(3*n+1). - Paul D. Hanna, Oct 30 2010 G.f.: hypergeom([1/3, 2/3], [1], 27 x^2 / (1 - 2x)^3) / (1 - 2x). - Michael Somos, Dec 17 2010 G.f.: Sum_{n >= 0} a(n)*x^n/n!^3 = [ Sum_{n >= 0} x^n/n!^3 ]^2. - Paul D. Hanna, Jan 19 2011 G.f.: A(x) = 1/(1-2*x)*(1+6*(x^2)/(G(0)-6*x^2)),
with G(k) = 3*(x^2)*(3*k+1)*(3*k+2) + ((1-2*x)^3)*((k+1)^2) - 3*(x^2)*((1-2*x)^3)*((k+1)^2)*(3*k+4)*(3*k+5)/G(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Dec 03 2011 In 2011 Zhi-Wei Sun found the formula Sum_{k = 0..n} C(2*k,n)*C(2*k,k)*C(2*(n-k),n-k) = (2^n)*a(n) and proved it via the Zeilberger algorithm. - Zhi-Wei Sun, Mar 20 2013 0 = a(n)*(a(n+1)*(-2048*a(n+2) - 3392*a(n+3) + 768*a(n+4)) + a(n+2)*(-1280*a(n+2) - 2912*a(n+3) + 744*a(n+4)) + a(n+3)*(+288*a(n+3) - 96*a(n+4))) + a(n+1)*(a(n+1)*(-704*a(n+2) - 1232*a(n+3) + 288*a(n+4)) + a(n+2)*(-560*a(n+2) - 1372*a(n+3) + 364*a(n+4)) + a(n+3)*(+154*a(n+3) - 53*a(n+4))) + a(n+2)*(a(n+2)*(+24*a(n+2) + 70*a(n+3) - 20*a(n+4)) + a(n+3)*(-11*a(n+3) + 4*a(n+4))) for all n in Z. - Michael Somos, Jul 16 2014 For r a nonnegative integer, Sum_{k = r..n} C(k,r)^3*C(n,k)^3 = C(n,r)^3*a(n-r), where we take a(n) = 0 for n < 0. - Peter Bala, Jul 27 2016 a(n) = (n!)^3 * [x^n] hypergeom([], [1, 1], x)^2. - Peter Luschny, May 31 2017 a(n) = Sum_{k=0..floor(n/2)} (n+k)!/(k!^3*(n-2*k)!) * 2^(n-2*k).
G.f. y=A(x) satisfies: 0 = x*(x + 1)*(8*x - 1)*y'' + (24*x^2 + 14*x - 1)*y' + 2*(4*x + 1)*y. (End)
a(n) = [x^n] (1 - x^2)^n*P(n,(1 + x)/(1 - x)), where P(n,x) denotes the n-th Legendre polynomial. See Gould, p. 56. - Peter Bala, Mar 24 2022 a(n) = (2^n/(4*Pi^2)) * Integral_{x,y=0..2*Pi} (1+cos(x)+cos(y)+cos(x+y))^n dx dy = (8^n/(Pi^2)) * Integral_{x,y=0..Pi} (cos(x)*cos(y)*cos(x+y))^n dx dy (Pla, 1995). - Amiram Eldar, Jul 16 2022 The g.f. T(x) obeys a period-annihilating ODE:
0=2*(1 + 4*x)*T(x) + (-1 + 14*x + 24*x^2)*T'(x) + x*(1 + x)*(-1 + 8*x)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = (4/243)*(1 - 8*x + 240*x^2 - 464*x^3 + 16*x^4);
g3 = -(8/19683)*(1 - 12*x - 480*x^2 + 3080*x^3 - 12072*x^4 + 4128*x^5 +
64*x^6);
which determine an elliptic surface with four singular fibers. (End)
For n >= 1, a(n) = 2 * Sum_{k = 0..n-1} binomial(n, k)^2 * binomial(n-1, k). Cf. A361716. For n >= 1, a(n) = 2 * hypergeom([-n, -n, -n + 1], [1, 1], -1). (End)
EXAMPLE
O.g.f.: A(x) = 1 + 2*x + 10*x^2 + 56*x^3 + 346*x^4 + 2252*x^5 + ...
O.g.f.: A(x) = 1/(1-2*x) + 3!*x^2/(1-2*x)^4 + (6!/2!^3)*x^4/(1-2*x)^7 + (9!/3!^3)*x^6/(1-2*x)^10 + (12!/4!^3)*x^8/(1-2*x)^13 + ... - Paul D. Hanna, Oct 30 2010 Let g.f. A(x) = Sum_{n >= 0} a(n)*x^n/n!^3, then
A(x) = 1 + 2*x + 10*x^2/2!^3 + 56*x^3/3!^3 + 346*x^4/4!^3 + ... where
A(x) = [1 + x + x^2/2!^3 + x^3/3!^3 + x^4/4!^3 + ...]^2. - Paul D. Hanna
MAPLE
add(binomial(n, k)^3, k=0..n) ;
end proc:
A000172_list := proc(len) series(hypergeom([], [1, 1], x)^2, x, len);
seq((n!)^3*coeff(%, x, n), n=0..len-1) end:
MATHEMATICA
Table[Sum[Binomial[n, k]^3, {k, 0, n}], {n, 0, 30}] (* Harvey P. Dale, Aug 24 2011 *) Table[ HypergeometricPFQ[{-n, -n, -n}, {1, 1}, -1], {n, 0, 20}] (* Jean-François Alcover, Jul 16 2012, after symbolic sum *) a[n_] := Sum[ Binomial[2k, n]*Binomial[2k, k]*Binomial[2(n-k), n-k], {k, 0, n}]/2^n; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Mar 20 2013, after Zhi-Wei Sun *) a[ n_] := SeriesCoefficient[ Hypergeometric2F1[ 1/3, 2/3, 1, 27 x^2 / (1 - 2 x)^3] / (1 - 2 x), {x, 0, n}]; (* Michael Somos, Jul 16 2014 *)
PROG
(PARI) {a(n)=polcoeff(sum(m=0, n, (3*m)!/m!^3*x^(2*m)/(1-2*x+x*O(x^n))^(3*m+1)), n)} \\ Paul D. Hanna, Oct 30 2010 (PARI) {a(n)=n!^3*polcoeff(sum(m=0, n, x^m/m!^3+x*O(x^n))^2, n)} \\ Paul D. Hanna, Jan 19 2011 (Haskell)
a000172 = sum . map a000578 . a007318_row
(Sage)
x, y, n = 1, 2, 1
while True:
yield x
n += 1
x, y = y, (8*(n-1)^2*x + (7*n^2-7*n + 2)*y) // n^2
(PARI) A000172(n)={sum(k=0, (n-1)\2, binomial(n, k)^3)*2+if(!bittest(n, 0), binomial(n, n\2)^3)} \\ M. F. Hasler, Sep 21 2015
CROSSREFS
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275- A291284 and A133370 respectively. Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.
Apery (Apéry) numbers: Sum_{k=0..n} (binomial(n,k)*binomial(n+k,k))^2.
(Formerly M4020)
+10
133
1, 5, 73, 1445, 33001, 819005, 21460825, 584307365, 16367912425, 468690849005, 13657436403073, 403676083788125, 12073365010564729, 364713572395983725, 11111571997143198073, 341034504521827105445, 10534522198396293262825, 327259338516161442321485
COMMENTS
Conjecture: For each n = 1,2,3,... the Apéry polynomial A_n(x) = Sum_{k = 0..n} binomial(n,k)^2*binomial(n+k,k)^2*x^k is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 21 2013 The expansions of exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 5*x + 49*x^2 + 685*x^3 + 11807*x^4 + 232771*x^5 + ... and exp( Sum_{n >= 1} a(n-1)*x^n/n ) = 1 + 3*x + 27*x^2 + 390*x^3 + 7038*x^4 + 144550*x^5 + ... both appear to have integer coefficients. See A267220. - Peter Bala, Jan 12 2016 Diagonal of the rational function R(x, y, z, w) = 1 / (1 - (w*x*y*z + w*x*y + w*z + x*y + x*z + y + z)); also diagonal of rational function H(x, y, z, w) = 1/(1 - w*(1+x)*(1+y)*(1+z)*(x*y*z + y*z + y + z + 1)). - Gheorghe Coserea, Jun 26 2018 Named after the French mathematician Roger Apéry (1916-1994). - Amiram Eldar, Jun 10 2021
REFERENCES
Julian Havil, The Irrationals, Princeton University Press, Princeton and Oxford, 2012, pp. 137-153.
Wolfram Koepf, Hypergeometric Identities. Ch. 2 in Hypergeometric Summation: An Algorithmic Approach to Summation and Special Function Identities. Braunschweig, Germany: Vieweg, pp. 55, 119 and 146, 1998.
Maxim Kontsevich and Don Zagier, Periods, pp. 771-808 of B. Engquist and W. Schmid, editors, Mathematics Unlimited - 2001 and Beyond, 2 vols., Springer-Verlag, 2001.
Leonard Lipshitz and Alfred van der Poorten, "Rational functions, diagonals, automata and arithmetic." In Number Theory, Richard A. Mollin, ed., Walter de Gruyter, Berlin (1990), pp. 339-358.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Roger Apéry, Irrationalité de zeta(2) et zeta(3), in Journées Arith. de Luminy. Colloque International du Centre National de la Recherche Scientifique (CNRS) held at the Centre Universitaire de Luminy, Luminy, Jun 20-24, 1978. Astérisque, Vol. 61 (1979), pp. 11-13. Frits Beukers, Consequences of Apéry's work on zeta(3), in "Zeta(3) irrationnel: les retombées", Rencontres Arithmétiques de Caen, June 2-3, 1995 [Mentions divisibility of a(n) by powers of 5 and powers of 11] Leonard Lipshitz and Alfred J. van der Poorten, Rational functions, diagonals, automata and arithmetic, in: Richard A. Mollin (ed.), Number theory, Proceedings of the First Conference of the Canadian Number Theory Association Held at the Banff Center, Banff, Alberta, April 17-27, 1988, de Gruyter, 2016, pp. 339-358; alternative link; Wayback Machine copy. Alfred van der Poorten, A proof that Euler missed ..., Math. Intelligencer, Vol. 1, No. 4 (December 1979), pp. 196-203, (b_n) after eq. (1.2), and Exercise 3.
FORMULA
D-finite with recurrence (n+1)^3*a(n+1) = (34*n^3 + 51*n^2 + 27*n + 5)*a(n) - n^3*a(n-1), n >= 1.
Representation as a special value of the hypergeometric function 4F3, in Maple notation: a(n)=hypergeom([n+1, n+1, -n, -n], [1, 1, 1], 1), n=0, 1, ... - Karol A. Penson Jul 24 2002 G.f.: (-1/2)*(3*x - 3 + (x^2-34*x+1)^(1/2))*(x+1)^(-2)*hypergeom([1/3,2/3],[1],(-1/2)*(x^2 - 7*x + 1)*(x+1)^(-3)*(x^2 - 34*x + 1)^(1/2)+(1/2)*(x^3 + 30*x^2 - 24*x + 1)*(x+1)^(-3))^2. - Mark van Hoeij, Oct 29 2011 Let g(x, y) = 4*cos(2*x) + 8*sin(y)*cos(x) + 5 and let P(n,z) denote the Legendre polynomial of degree n. Then G. A. Edgar posted a conjecture of Alexandru Lupas that a(n) equals the double integral 1/(4*Pi^2)*int {y = -Pi..Pi} int {x = -Pi..Pi} P(n,g(x,y)) dx dy. (Added Jan 07 2015: Answered affirmatively in Math Overflow question 178790) - Peter Bala, Mar 04 2012; edited by G. A. Edgar, Dec 10 2016 a(n) ~ (1+sqrt(2))^(4*n+2)/(2^(9/4)*Pi^(3/2)*n^(3/2)). - Vaclav Kotesovec, Nov 01 2012 a(n) = Sum_{k=0..n} C(n,k)^2 * C(n+k,k)^2. - Joerg Arndt, May 11 2013 0 = (-x^2+34*x^3-x^4)*y''' + (-3*x+153*x^2-6*x^3)*y'' + (-1+112*x-7*x^2)*y' + (5-x)*y, where y is g.f. - Gheorghe Coserea, Jul 14 2016 a(n) = Sum_{0 <= j, k <= n} (-1)^(n+j) * C(n,k)^2 * C(n+k,k)^2 * C(n,j) * C(n+k+j,k+j).
a(n) = Sum_{0 <= j, k <= n} C(n,k) * C(n+k,k) * C(k,j)^3 (see Koepf, p. 55).
a(n) = Sum_{0 <= j, k <= n} C(n,k)^2 * C(n,j)^2 * C(3*n-j-k,2*n) (see Koepf, p. 119).
Diagonal coefficients of the rational function 1/((1 - x - y)*(1 - z - t) - x*y*z*t) (Straub, 2014). (End)
a(n) = [x^n] 1/(1 - x)*( Legendre_P(n,(1 + x)/(1 - x)) )^m at m = 2. At m = 1 we get the Apéry numbers A005258. - Peter Bala, Dec 22 2020 a(n) = Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*binomial(n+k, k)* A108625(n, k). - Peter Bala, Jul 18 2024
EXAMPLE
G.f. = 1 + 5*x + 73*x^2 + 1445*x^3 + 33001*x^4 + 819005*x^5 + 21460825*x^6 + ...
a(2) = (binomial(2,0) * binomial(2+0,0))^2 + (binomial(2,1) * binomial(2+1,1))^2 + (binomial(2,2) * binomial(2+2,2))^2 = (1*1)^2 + (2*3)^2 + (1*6)^2 = 1 + 36 + 36 = 73. - Michael B. Porter, Jul 14 2016
MAPLE
a := proc(n) option remember; if n=0 then 1 elif n=1 then 5 else (n^(-3))* ( (34*(n-1)^3 + 51*(n-1)^2 + 27*(n-1) +5)*a((n-1)) - (n-1)^3*a((n-1)-1)); fi; end;
# Alternative:
a := n -> hypergeom([-n, -n, 1+n, 1+n], [1, 1, 1], 1):
MATHEMATICA
Table[HypergeometricPFQ[{-n, -n, n+1, n+1}, {1, 1, 1}, 1], {n, 0, 13}] (* Jean-François Alcover, Apr 01 2011 *) Table[Sum[(Binomial[n, k]Binomial[n+k, k])^2, {k, 0, n}], {n, 0, 30}] (* Harvey P. Dale, Oct 15 2011 *) a[ n_] := SeriesCoefficient[ SeriesCoefficient[ SeriesCoefficient[ SeriesCoefficient[ 1 / (1 - t (1 + x ) (1 + y ) (1 + z ) (x y z + (y + 1) (z + 1))), {t, 0, n}], {x, 0, n}], {y, 0, n}], {z, 0, n}]; (* Michael Somos, May 14 2016 *)
PROG
(Haskell)
a005259 n = a005259_list !! n
a005259_list = 1 : 5 : zipWith div (zipWith (-)
(tail $ zipWith (*) a006221_list a005259_list)
(zipWith (*) (tail a000578_list) a005259_list)) (drop 2 a000578_list)
(GAP) List([0..20], n->Sum([0..n], k->Binomial(n, k)^2*Binomial(n+k, k)^2)); # Muniru A Asiru, Sep 28 2018 (Magma) [&+[Binomial(n, k) ^2 *Binomial(n+k, k)^2: k in [0..n]]:n in [0..17]]; // Marius A. Burtea, Jan 20 2020 (Python)
m, g = 1, 0
for k in range(n+1):
g += m
m *= ((n+k+1)*(n-k))**2
m //=(k+1)**4
CROSSREFS
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275- A291284 and A133370 respectively.
Apéry numbers: a(n) = Sum_{k=0..n} binomial(n,k)^2 * binomial(n+k,k).
(Formerly M3057)
+10
114
1, 3, 19, 147, 1251, 11253, 104959, 1004307, 9793891, 96918753, 970336269, 9807518757, 99912156111, 1024622952993, 10567623342519, 109527728400147, 1140076177397091, 11911997404064793, 124879633548031009, 1313106114867738897, 13844511065506477501
COMMENTS
This is the Taylor expansion of a special point on a curve described by Beauville. - Matthijs Coster, Apr 28 2004 Conjecture: For each n=1,2,3,... the polynomial a_n(x) = Sum_{k=0..n} C(n,k)^2*C(n+k,k)*x^k is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 21 2013 Diagonal of rational functions 1/(1 - x - x*y - y*z - x*z - x*y*z), 1/(1 + y + z + x*y + y*z + x*z + x*y*z), 1/(1 - x - y - z + x*y + x*y*z), 1/(1 - x - y - z + y*z + x*z - x*y*z). - Gheorghe Coserea, Jul 07 2018
REFERENCES
Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.
S. Melczer, An Invitation to Analytic Combinatorics, 2021; p. 129.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
R. Apéry, Irrationalité de zeta(2) et zeta(3), in Journées Arith. de Luminy. Colloque International du Centre National de la Recherche Scientifique (CNRS) held at the Centre Universitaire de Luminy, Luminy, Jun 20-24, 1978. Astérisque, 61 (1979), 11-13.
FORMULA
D-finite with recurrence: (n+1)^2 * a(n+1) = (11*n^2+11*n+3) * a(n) + n^2 * a(n-1). - Matthijs Coster, Apr 28 2004 Let b(n) be the solution to the above recurrence with b(0) = 0, b(1) = 5. Then the b(n) are rational numbers with b(n)/a(n) -> zeta(2) very rapidly. The identity b(n)*a(n-1) - b(n-1)*a(n) = (-1)^(n-1)*5/n^2 leads to a series acceleration formula: zeta(2) = 5 * Sum_{n >= 1} 1/(n^2*a(n)*a(n-1)) = 5*(1/(1*3) + 1/(2^2*3*19) + 1/(3^2*19*147) + ...). Similar results hold for the constant e: see A143413. - Peter Bala, Aug 14 2008 G.f.: hypergeom([1/12, 5/12],[1], 1728*x^5*(1-11*x-x^2)/(1-12*x+14*x^2+12*x^3+x^4)^3) / (1-12*x+14*x^2+12*x^3+x^4)^(1/4). - Mark van Hoeij, Oct 25 2011 a(n) ~ ((11+5*sqrt(5))/2)^(n+1/2)/(2*Pi*5^(1/4)*n). - Vaclav Kotesovec, Oct 05 2012 1/Pi = 5*(sqrt(47)/7614)*Sum_{n>=0} (-1)^n a(n)*binomial(2n,n)*(682n+71)/15228^n. [Cooper, equation (4)] - Jason Kimberley, Nov 26 2012 a(-1 - n) = (-1)^n * a(n) if n>=0. a(-1 - n) = -(-1)^n * a(n) if n<0. - Michael Somos, Sep 18 2013 0 = a(n)*(a(n+1)*(+4*a(n+2) + 83*a(n+3) - 12*a(n+4)) + a(n+2)*(+32*a(n+2) + 902*a(n+3) - 147*a(n+4)) + a(n+3)*(-56*a(n+3) + 12*a(n+4))) + a(n+1)*(a(n+1)*(+17*a(n+2) + 374*a(n+3) - 56*a(n+4)) + a(n+2)*(+176*a(n+2) + 5324*a(n+3) - 902*a(n+4) + a(n+3)*(-374*a(n+3) + 83*a(n+4))) + a(n+2)*(a(n+2)*(-5*a(n+2) - 176*a(n+3) + 32*a(n+4)) + a(n+3)*(+17*a(n+3) - 4*a(n+4))) for all n in Z. - Michael Somos, Aug 06 2016 a(n) = binomial(2*n, n)*hypergeom([-n, -n, -n],[1, -2*n], 1). - Peter Luschny, Feb 10 2018 a(n) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*binomial(n+k,k)^2. - Peter Bala, Feb 10 2018 G.f. y=A(x) satisfies: 0 = x*(x^2 + 11*x - 1)*y'' + (3*x^2 + 22*x - 1)*y' + (x + 3)*y. - Gheorghe Coserea, Jul 01 2018 a(n) = Sum_{0 <= j, k <= n} (-1)^(j+k)*C(n,k)*C(n+k,k)^2*C(n,j)* C(n+k+j,k+j).
a(n) = Sum_{0 <= j, k <= n} (-1)^(n+j)*C(n,k)^2*C(n+k,k)*C(n,j)* C(n+k+j,k+j).
a(n) = Sum_{0 <= j, k <= n} (-1)^j*C(n,k)^2*C(n,j)*C(3*n-j-k,2*n). (End)
a(n) = [x^n] 1/(1 - x)*( Legendre_P(n,(1 + x)/(1 - x)) )^m at m = 1. At m = 2 we get the Apéry numbers A005259. - Peter Bala, Dec 22 2020 a(n) = (-1)^n*Sum_{j=0..n} (1 - 5*j*H(j) + 5*j*H(n - j))*binomial(n, j)^5, where H(n) denotes the n-th harmonic number, A001008/ A002805. (Paule/Schneider). - Peter Luschny, Jul 23 2021 The g.f. T(x) obeys a period-annihilating ODE:
0=(3 + x)*T(x) + (-1 + 22*x + 3*x^2)*T'(x) + x*(-1 + 11*x + x^2)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = 3*(1 - 12*x + 14*x^2 + 12*x^3 + x^4);
g3 = 1 - 18*x + 75*x^2 + 75*x^4 + 18*x^5 + x^6;
which determine an elliptic surface with four singular fibers. (End)
Conjecture: a(n)^2 = Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*binomial(n+k, k)* A143007(n, k). - Peter Bala, Jul 08 2024
EXAMPLE
G.f. = 1 + 3*x + 19*x^2 + 147*x^3 + 1251*x^4 + 11253*x^5 + 104959*x^6 + ...
MAPLE
with(combinat): seq(add((multinomial(n+k, n-k, k, k))*binomial(n, k), k=0..n), n=0..18); # Zerinvary Lajos, Oct 18 2006 a := n -> binomial(2*n, n)*hypergeom([-n, -n, -n], [1, -2*n], 1):
MATHEMATICA
a[n_] := HypergeometricPFQ[ {n+1, -n, -n}, {1, 1}, 1]; Table[ a[n], {n, 0, 18}] (* Jean-François Alcover, Jan 20 2012, after Vladeta Jovovic *) Table[Sum[Binomial[n, k]^2 Binomial[n+k, k], {k, 0, n}], {n, 0, 20}] (* Harvey P. Dale, Aug 25 2019 *)
PROG
(Haskell)
a005258 n = sum [a007318 n k ^ 2 * a007318 (n + k) k | k <- [0..n]]
(PARI) {a(n) = if( n<0, -(-1)^n * a(-1-n), sum(k=0, n, binomial(n, k)^2 * binomial(n+k, k)))} /* Michael Somos, Sep 18 2013 */ (GAP) a:=n->Sum([0..n], k->(-1)^(n-k)*Binomial(n, k)*Binomial(n+k, k)^2);;
(GAP) List([0..20], n->Sum([0..n], k->Binomial(n, k)^2*Binomial(n+k, k))); # Muniru A Asiru, Jul 29 2018 (Magma) [&+[Binomial(n, k)^2 * Binomial(n+k, k): k in [0..n]]: n in [0..25]]; // Vincenzo Librandi, Nov 28 2018 (Python)
m, g = 1, 0
for k in range(n+1):
g += m
m *= (n+k+1)*(n-k)**2
m //= (k+1)**3
CROSSREFS
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275- A291284 and A133370 respectively.
a(n) = Sum_{k=0..n} binomial(n,k)^2 * binomial(2*k,k).
(Formerly M2998 N1214)
+10
89
1, 3, 15, 93, 639, 4653, 35169, 272835, 2157759, 17319837, 140668065, 1153462995, 9533639025, 79326566595, 663835030335, 5582724468093, 47152425626559, 399769750195965, 3400775573443089, 29016970072920387, 248256043372999089
COMMENTS
This is the Taylor expansion of a special point on a curve described by Beauville. - Matthijs Coster, Apr 28 2004 a(n) is the 2n-th moment of the distance from the origin of a 3-step random walk in the plane. - Peter M. W. Gill (peter.gill(AT)nott.ac.uk), Feb 27 2004
a(n) is the number of Abelian squares of length 2n over a 3-letter alphabet. - Jeffrey Shallit, Aug 17 2010 Consider 2D simple random walk on honeycomb lattice. a(n) gives number of paths of length 2n ending at origin. - Sergey Perepechko, Feb 16 2011 Conjecture: For each n=1,2,3,... the polynomial g_n(x) = Sum_{k=0..n} binomial(n,k)^2*binomial(2k,k)*x^k is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 21 2013 This is one of the Apery-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017 a(n) is the sum of the squares of the coefficients of (x + y + z)^n. - Michael Somos, Aug 25 2018 a(n) is the constant term in the expansion of (1 + (1 + x) * (1 + y) + (1 + 1/x) * (1 + 1/y))^n. - Seiichi Manyama, Oct 28 2019
REFERENCES
Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.
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).
FORMULA
a(n) = Sum_{m=0..n} binomial(n, m) * A000172(m). [Barrucand] D-finite with recurrence: (n+1)^2 a(n+1) = (10*n^2+10*n+3) * a(n) - 9*n^2 * a(n-1). - Matthijs Coster, Apr 28 2004 Sum_{n>=0} a(n)*x^n/n!^2 = BesselI(0, 2*sqrt(x))^3. - Vladeta Jovovic, Mar 11 2003 a(n) = Sum_{p+q+r=n} (n!/(p!*q!*r!))^2 with p, q, r >= 0. - Michael Somos, Jul 25 2007 a(n) = hypergeom([1/2, -n, -n], [1, 1], 4). - Mark van Hoeij, Jun 02 2010 G.f.: 2*sqrt(2)/Pi/sqrt(1-6*z-3*z^2+sqrt((1-z)^3*(1-9*z))) * EllipticK(8*z^(3/2)/(1-6*z-3*z^2+sqrt((1-z)^3*(1-9*z)))). - Sergey Perepechko, Feb 16 2011 G.f.: Sum_{n>=0} (3*n)!/n!^3 * x^(2*n)*(1-x)^n / (1-3*x)^(3*n+1). - Paul D. Hanna, Feb 26 2012 G.f.: 1/(1-3*x)*(1-6*x^2*(1-x)/(Q(0)+6*x^2*(1-x))), where Q(k) = (54*x^3 - 54*x^2 + 9*x -1)*k^2 + (81*x^3 - 81*x^2 + 18*x -2)*k + 33*x^3 - 33*x^2 +9*x - 1 - 3*x^2*(1-x)*(1-3*x)^3*(k+1)^2*(3*k+4)*(3*k+5)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013 G.f.: G(0)/(2*(1-9*x)^(2/3)), where G(k) = 1 + 1/(1 - 3*(3*k+1)^2*x*(1-x)^2/(3*(3*k+1)^2*x*(1-x)^2 - (k+1)^2*(1-9*x)^2/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 31 2013 a(n) = [x^(2n)] 1/agm(sqrt((1-3*x)*(1+x)^3), sqrt((1+3*x)*(1-x)^3)). - Gheorghe Coserea, Aug 17 2016 0 = +a(n)*(+a(n+1)*(+729*a(n+2) -1539*a(n+3) +243*a(n+4)) +a(n+2)*(-567*a(n+2) +1665*a(n+3) -297*a(n+4)) +a(n+3)*(-117*a(n+3) +27*a(n+4))) +a(n+1)*(+a(n+1)*(-324*a(n+2) +720*a(n+3) -117*a(n+4)) +a(n+2)*(+315*a(n+2) -1000*a(n+3) +185*a(n+4)) +a(n+3)*(+80*a(n+3) -19*a(n+4))) +a(n+2)*(+a(n+2)*(-9*a(n+2) +35*a(n+3) -7*a(n+4)) +a(n+3)*(-4*a(n+3) +a(n+4))) for all n in Z. - Michael Somos, Oct 30 2017 G.f. y=A(x) satisfies: 0 = x*(x - 1)*(9*x - 1)*y'' + (27*x^2 - 20*x + 1)*y' + 3*(3*x - 1)*y. - Gheorghe Coserea, Jul 01 2018 Sum_{k>=0} binomial(2*k,k) * a(k) / 6^(2*k) = A086231 = (sqrt(3)-1) * (Gamma(1/24) * Gamma(11/24))^2 / (32*Pi^3). - Vaclav Kotesovec, Apr 23 2023 The g.f. T(x) obeys a period-annihilating ODE:
0=3*(-1 + 3*x)*T(x) + (1 - 20*x + 27*x^2)*T'(x) + x*(-1 + x)*(-1 + 9*x)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = (3/64)*(1 + 3*x)*(1 - 15*x + 75*x^2 + 3*x^3);
g3 = -(1/512)*(-1 + 6*x + 3*x^2)*(1 - 12*x + 30*x^2 - 540*x^3 + 9*x^4);
which determine an elliptic surface with four singular fibers. (End)
EXAMPLE
G.f.: A(x) = 1 + 3*x + 15*x^2 + 93*x^3 + 639*x^4 + 4653*x^5 + 35169*x^6 + ...
G.f.: A(x) = 1/(1-3*x) + 6*x^2*(1-x)/(1-3*x)^4 + 90*x^4*(1-x)^2/(1-3*x)^7 + 1680*x^6*(1-x)^3/(1-3*x)^10 + 34650*x^8*(1-x)^4/(1-3*x)^13 + ... - Paul D. Hanna, Feb 26 2012
MAPLE
series(1/GaussAGM(sqrt((1-3*x)*(1+x)^3), sqrt((1+3*x)*(1-x)^3)), x=0, 42) # Gheorghe Coserea, Aug 17 2016 A002893 := n -> hypergeom([1/2, -n, -n], [1, 1], 4):
MATHEMATICA
Table[Sum[Binomial[n, k]^2 Binomial[2k, k], {k, 0, n}], {n, 0, 20}] (* Harvey P. Dale, Aug 19 2011 *) a[ n_] := If[ n < 0, 0, HypergeometricPFQ[ {1/2, -n, -n}, {1, 1}, 4]]; (* Michael Somos, Oct 16 2013 *) a[n_] := SeriesCoefficient[BesselI[0, 2*Sqrt[x]]^3, {x, 0, n}]*n!^2; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Dec 30 2013 *) a[ n_] := If[ n < 0, 0, Block[ {x, y, z}, Expand[(x + y + z)^n] /. {t_Integer -> t^2, x -> 1, y -> 1, z -> 1}]]; (* Michael Somos, Aug 25 2018 *)
PROG
(PARI) {a(n) = if( n<0, 0, n!^2 * polcoeff( besseli(0, 2*x + O(x^(2*n+1)))^3, 2*n))};
(PARI) {a(n) = sum(k=0, n, binomial(n, k)^2 * binomial(2*k, k))}; /* Michael Somos, Jul 25 2007 */ (PARI) {a(n)=polcoeff(sum(m=0, n, (3*m)!/m!^3 * x^(2*m)*(1-x)^m / (1-3*x+x*O(x^n))^(3*m+1)), n)} \\ Paul D. Hanna, Feb 26 2012 (PARI) N = 42; x='x + O('x^N); v = Vec(1/agm(sqrt((1-3*x)*(1+x)^3), sqrt((1+3*x)*(1-x)^3))); vector((#v+1)\2, k, v[2*k-1]) \\ Gheorghe Coserea, Aug 17 2016 (Magma) [&+[Binomial(n, k)^2 * Binomial(2*k, k): k in [0..n]]: n in [0..25]]; // Vincenzo Librandi, Aug 26 2018 (SageMath)
def A002893(n): return simplify(hypergeometric([1/2, -n, -n], [1, 1], 4))
CROSSREFS
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275- A291284 and A133370 respectively.
Domb numbers: number of 2n-step polygons on diamond lattice.
(Formerly M3626 N1473)
+10
66
1, 4, 28, 256, 2716, 31504, 387136, 4951552, 65218204, 878536624, 12046924528, 167595457792, 2359613230144, 33557651538688, 481365424895488, 6956365106016256, 101181938814289564, 1480129751586116848, 21761706991570726096, 321401321741959062016
COMMENTS
a(n) is the (2n)th moment of the distance from the origin of a 4-step random walk in the plane. - Peter M.W. Gill (peter.gill(AT)nott.ac.uk), Mar 03 2004
Conjecture: Let D(n) be the (n+1) X (n+1) Hankel-type determinant with (i,j)-entry equal to a(i+j) for all i,j = 0..n. Then the number D(n)/12^n is always a positive odd integer. - Zhi-Wei Sun, Aug 14 2013 It appears that the expansions exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 4*x + 22*x^2 + 152*x^3 + 1241*x^4 + ... and exp( Sum_{n >= 1} 1/4*a(n)*x^n/n ) = 1 + x + 4*x^2 + 25*x^3 + 199*x^4 + ... have integer coefficients. See A267219. - Peter Bala, Jan 12 2016 This is one of the Apéry-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017 Named after the British-Israeli theoretical physicist Cyril Domb (1920-2012). - Amiram Eldar, Mar 20 2021
REFERENCES
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).
LINKS
Shaun Cooper, James G. Wan and Wadim Zudilin, Holonomic Alchemy and Series for 1/pi, in: G. Andrews and F. Garvan (eds.) Analytic Number Theory, Modular Forms and q-Hypergeometric Series, ALLADI60 2016, Springer Proceedings in Mathematics & Statistics, Vol 221. Springer, Cham, 2016; arXiv preprint, arXiv:1512.04608 [math.NT], 2015. Zhi-Wei Sun, Conjectures involving arithmetical sequences, in: S. Kanemitsu, H. Li and J. Liu (eds.), Number Theory: Arithmetic in Shangri-La, Proc. the 6th China-Japan Sem. Number Theory (Shanghai, August 15-17, 2011), World Sci., Singapore, 2013, pp. 244-258; alternative link.
FORMULA
a(n) = Sum_{k=0..n} binomial(n, k)^2 * binomial(2n-2k, n-k) * binomial(2k, k).
D-finite with recurrence: n^3*a(n) = 2*(2*n-1)*(5*n^2-5*n+2)*a(n-1) - 64*(n-1)^3*a(n-2). - Vladeta Jovovic, Jul 16 2004 Sum_{n>=0} a(n)*x^n/n!^2 = BesselI(0, 2*sqrt(x))^4. - Vladeta Jovovic, Aug 01 2006 G.f.: hypergeom([1/6, 1/3],[1],108*x^2/(1-4*x)^3)^2/(1-4*x). - Mark van Hoeij, Oct 29 2011 (1) 4^n*a(n) = Sum_{k = 0..n} (binomial(2k,k)*binomial(2(n-k),n-k))^3/ binomial(n,k)^2,
(2) a(n) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*binomial(2k,n)*binomial(2k,k)* binomial(2(n-k),n-k). (End)
G.f. y=A(x) satisfies: 0 = x^2*(4*x - 1)*(16*x - 1)*y''' + 3*x*(128*x^2 - 30*x + 1)*y'' + (448*x^2 - 68*x + 1)*y' + 4*(16*x - 1)*y. - Gheorghe Coserea, Jun 26 2018 a(n) = Sum_{p+q+r+s=n} (n!/(p!*q!*r!*s!))^2 with p,q,r,s >= 0. See Verrill, p. 5. - Peter Bala, Jan 06 2020 a(n) = 2*Sum_{k = 1..n} (k/n)*binomial(n, k)^2*binomial(2*n-2*k, n-k)* binomial(2*k, k) for n >= 1.
a(n-1) = (1/2)*Sum_{k = 1..n} (k/n)^3*binomial(n, k)^2*binomial(2*n-2*k, n-k)* binomial(2*k, k) for n >= 1. Cf. A081085. (End)
MAPLE
A002895 := n -> binomial(2*n, n)*hypergeom([1/2, -n, -n, -n], [1, 1, 1/2 - n], 1):
MATHEMATICA
Table[Sum[Binomial[n, k]^2 Binomial[2n-2k, n-k]Binomial[2k, k], {k, 0, n}], {n, 0, 30}] (* Harvey P. Dale, Aug 15 2011 *) a[n_] = Binomial[2*n, n]*HypergeometricPFQ[{1/2, -n, -n, -n}, {1, 1, 1/2-n}, 1]; (* or *) a[n_] := SeriesCoefficient[BesselI[0, 2*Sqrt[x]]^4, {x, 0, n}]*n!^2; Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Dec 30 2013, after Vladeta Jovovic *) max = 19; Total /@ MatrixPower[Table[Binomial[n, k]^2, {n, 0, max}, {k, 0, max}], 3] (* Jean-François Alcover, Mar 24 2015, after Peter Bala *)
PROG
(PARI) C=binomial;
a(n) = sum(k=0, n, C(n, k)^2 * C(2*n-2*k, n-k) * C(2*k, k) );
CROSSREFS
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692, A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.) For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275- A291284 and A133370 respectively.
Search completed in 0.077 seconds
|