[go: up one dir, main page]

login
Search: a258344 -id:a258344
     Sort: relevance | references | number | modified | created      Format: long | short | data
Expansion of Product_{m>0} (1+q^m)^(m(m+1)/2).
+10
34
1, 1, 3, 9, 19, 46, 100, 218, 460, 965, 1975, 3993, 7975, 15712, 30650, 59150, 113093, 214300, 402812, 751165, 1390714, 2557004, 4670770, 8479232, 15302657, 27462424, 49021252, 87057783, 153850769, 270614429, 473850031, 826125184, 1434286323, 2480145226
OFFSET
0,3
COMMENTS
Convolved with aerated A000294: [1, 0, 2, 0, 4, 0, 10, 0, 26, ...] = A000294. - Gary W. Adamson, Jun 13 2009
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -n*(n+1)/2, g(n) = -1. - Seiichi Manyama, Nov 14 2017
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Alois P. Heinz)
FORMULA
a(n) ~ 7^(1/8) * exp(2 * 7^(1/4) * Pi * n^(3/4) / (3^(5/4) * 5^(1/4)) + 3^(3/2) * 5^(1/2) * Zeta(3) * n^(1/2) / (2 * 7^(1/2) * Pi^2) - 3^(13/4) * 5^(5/4) * Zeta(3)^2 * n^(1/4) / (4 * 7^(5/4) * Pi^5) + 2025 * Zeta(3)^3 / (98*Pi^8)) / (2^(49/24) * 15^(1/8) * n^(5/8)), where Zeta(3) = A002117. - Vaclav Kotesovec, Mar 11 2015
a(0) = 1 and a(n) = (1/(2*n)) * Sum_{k=1..n} b(k)*a(n-k) where b(n) = Sum_{d|n} d^2*(d+1)*(-1)^(1+n/d). - Seiichi Manyama, Nov 14 2017
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k/(k*(1 - x^k)^3)). - Ilya Gutkovskiy, May 28 2018
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
add(binomial(i*(i+1)/2, j)*b(n-i*j, i-1), j=0..n/i)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..50); # Alois P. Heinz, Aug 03 2013
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, Sum[Binomial[i*(i+1)/2, j]*b[n-i*j, i-1], {j, 0, n/i}]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Oct 13 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
STATUS
approved
Expansion of Product_{k>=1} (1+x^k)^(k*(k+1)*(k+2)/6).
+10
21
1, 1, 4, 14, 36, 101, 260, 669, 1669, 4116, 9932, 23636, 55483, 128532, 294422, 667026, 1496232, 3324720, 7323570, 15998749, 34679966, 74622839, 159454379, 338472749, 713956569, 1496950669, 3120663129, 6469901522, 13343153563, 27379250529, 55907749171
OFFSET
0,3
LINKS
FORMULA
a(n) ~ (3*Zeta(5))^(1/10) / (2^(523/720) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401 * Pi^16 / (10497600000000 * Zeta(5)^3) + 49*Pi^8 * Zeta(3) / (16200000 * Zeta(5)^2) - Zeta(3)^2 / (150*Zeta(5)) + (343*Pi^12 / (2430000000 * 2^(3/5) * 15^(1/5) * Zeta(5)^(11/5)) - 7*Pi^4 * Zeta(3) / (4500 * 2^(3/5) * 15^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (1080000 * 2^(1/5) * 15^(2/5) * Zeta(5)^(7/5)) + Zeta(3) / (2^(6/5) * (15*Zeta(5))^(2/5))) * n^(2/5) + 7*Pi^4 / (180 * 2^(4/5) * (15*Zeta(5))^(3/5)) * n^(3/5) + 5*(15*Zeta(5))^(1/5) / 2^(12/5) * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k/(k*(1 - x^k)^4)). - Ilya Gutkovskiy, May 28 2018
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, add(
binomial(binomial(i+2, 3), j)*b(n-i*j, i-1), j=0..n/i)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..30); # Alois P. Heinz, May 28 2018
MATHEMATICA
nmax=40; CoefficientList[Series[Product[(1+x^k)^(k*(k+1)*(k+2)/6), {k, 1, nmax}], {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved
Expansion of Product(1+q^m)^(m(m-1)/2); m=1..inf.
+10
9
1, 0, 1, 3, 6, 13, 24, 49, 91, 181, 334, 632, 1163, 2138, 3880, 7006, 12531, 22279, 39369, 69078, 120597, 209282, 361405, 620829, 1061687, 1807014, 3062642, 5168784, 8688820, 14549659, 24274226, 40353748, 66854518, 110391391, 181695436, 298129605, 487706902
OFFSET
0,4
LINKS
FORMULA
a(n) ~ 7^(1/8) / (2^(47/24) * 15^(1/8) * n^(5/8)) * exp(-2025 * Zeta(3)^3 / (98*Pi^8) - 135*(15/7)^(1/4) * Zeta(3)^2 / (28*Pi^5) * n^(1/4) - 3*sqrt(15/7) * Zeta(3) / (2*Pi^2) * sqrt(n) + 2*(7/15)^(1/4) * Pi/3 * n^(3/4)), where Zeta(3) = A002117. - Vaclav Kotesovec, May 27 2015
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
add(binomial(binomial(i, 2), j)*b(n-i*j, i-1), j=0..n/i)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..50); # Alois P. Heinz, Aug 03 2013
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, Sum[Binomial[Binomial[i, 2], j]*b[n-i*j, i-1], {j, 0, n/i}]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Oct 13 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
Expansion of Product_{k>=1} (1+x^k)^(k*(k+1)).
+10
9
1, 2, 7, 24, 65, 184, 487, 1254, 3145, 7706, 18480, 43490, 100692, 229472, 515802, 1144416, 2508948, 5439642, 11671859, 24801738, 52221911, 109013538, 225718717, 463769652, 945915199, 1915895576, 3854803572, 7706786958, 15314564282, 30255672820, 59440488874
OFFSET
0,2
LINKS
FORMULA
a(n) ~ 7^(1/8) / (2^(47/24) * 15^(1/8) * n^(5/8)) * exp(2025*Zeta(3)^3 / (49*Pi^8) - 135*(15/14)^(1/4) * Zeta(3)^2 / (14*Pi^5) * n^(1/4) + 3*sqrt(15/14) * Zeta(3) / Pi^2 * sqrt(n) + 2*(14/15)^(1/4)*Pi/3 * n^(3/4)), where Zeta(3) = A002117.
MATHEMATICA
nmax=30; CoefficientList[Series[Product[(1+x^k)^(k*(k+1)), {k, 1, nmax}], {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved
Expansion of Product_{k>=1} (1+x^k)^(k*(k+1)*(k+2)).
+10
8
1, 6, 39, 224, 1131, 5412, 24411, 105078, 435048, 1740312, 6755877, 25533330, 94205738, 340064322, 1203313782, 4180514846, 14279610417, 48013553310, 159086287869, 519912616614, 1677331973910, 5345927500226, 16843574682291, 52494817082952, 161923200857711
OFFSET
0,2
LINKS
FORMULA
a(n) ~ 3^(1/5) * Zeta(5)^(1/10) / (2^(91/120) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401*Pi^16 / (1749600000000 * Zeta(5)^3) + 49*Pi^8 * Zeta(3) / (2700000 * Zeta(5)^2) - Zeta(3)^2 / (25*Zeta(5)) + (343*Pi^12/(405000000 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(11/5)) - 7*Pi^4 * Zeta(3) / (750 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (180000 * 2^(3/5) * 3^(4/5) * 5^(2/5) * Zeta(5)^(7/5)) + 3^(1/5) * Zeta(3) / (2^(3/5) * (5*Zeta(5))^(2/5))) * n^(2/5) + 7*Pi^4 / (180 * 2^(2/5) * 3^(1/5) * (5*Zeta(5))^(3/5)) * n^(3/5) + 5*3^(2/5) * (5*Zeta(5)/2)^(1/5)/4 * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.
MATHEMATICA
nmax=30; CoefficientList[Series[Product[(1+x^k)^(k*(k+1)*(k+2)), {k, 1, nmax}], {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved
Expansion of Product_{k>=1} (1+x^k)^(k*(k-1)*(k-2)).
+10
8
1, 0, 0, 6, 24, 60, 135, 354, 972, 2684, 6990, 17802, 44627, 111582, 277329, 684164, 1671984, 4050096, 9735209, 23238480, 55120950, 129940442, 304502583, 709464798, 1643920584, 3789158988, 8690016942, 19833550266, 45056952957, 101900481462, 229462378987
OFFSET
0,4
LINKS
FORMULA
a(n) ~ 3^(1/5) * Zeta(5)^(1/10) / (2^(91/120) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401 * Pi^16 / (1749600000000*Zeta(5)^3) + 49 * Pi^8 * Zeta(3) / (2700000 * Zeta(5)^2) - Zeta(3)^2 / (25*Zeta(5)) + (-343 * Pi^12 / (405000000 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(11/5)) + 7*Pi^4 * Zeta(3) / (750 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (180000 * 2^(3/5) * 3^(4/5) * 5^(2/5) * Zeta(5)^(7/5)) + 3^(1/5) * Zeta(3) / (2^(3/5) * (5*Zeta(5))^(2/5))) * n^(2/5) - 7*Pi^4 / (180 * 2^(2/5) * 3^(1/5) * (5*Zeta(5))^(3/5)) * n^(3/5) + 5*3^(2/5) * ((5*Zeta(5))/2)^(1/5)/4 * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.
MATHEMATICA
nmax=40; CoefficientList[Series[Product[(1+x^k)^(k*(k-1)*(k-2)), {k, 1, nmax}], {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved
Expansion of Product_{k>=1} (1+x^k)^(k*(k-1)*(k-2)/6).
+10
8
1, 0, 0, 1, 4, 10, 20, 39, 72, 144, 280, 567, 1112, 2187, 4204, 8073, 15309, 28986, 54548, 102286, 190881, 354717, 656194, 1208712, 2217624, 4052633, 7379630, 13390098, 24215587, 43649482, 78435884, 140513905, 250988186, 447037367, 794031641, 1406585604
OFFSET
0,5
LINKS
FORMULA
a(n) ~ (3*Zeta(5))^(1/10) / (2^(523/720) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401 * Pi^16 / (10497600000000 * Zeta(5)^3) + 49*Pi^8 * Zeta(3) / (16200000 * Zeta(5)^2) - Zeta(3)^2 / (150*Zeta(5)) + (-343*Pi^12 / (2430000000 * 2^(3/5) * 15^(1/5) * Zeta(5)^(11/5)) + 7*Pi^4 * Zeta(3) / (4500 * 2^(3/5) * 15^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (1080000 * 2^(1/5) * 15^(2/5) * Zeta(5)^(7/5)) + Zeta(3) / (2^(6/5) * (15*Zeta(5))^(2/5))) * n^(2/5) - 7*Pi^4 / (180 * 2^(4/5) * (15*Zeta(5))^(3/5)) * n^(3/5) + 5*(15*Zeta(5))^(1/5) / 2^(12/5) * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.
MATHEMATICA
nmax=50; CoefficientList[Series[Product[(1+x^k)^(k*(k-1)*(k-2)/6), {k, 1, nmax}], {x, 0, nmax}], x]
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved
Expansion of Product_{k>=1} 1/(1-x^k)^(k*(k-1)).
+10
8
1, 0, 2, 6, 15, 32, 79, 172, 397, 860, 1879, 3986, 8462, 17586, 36408, 74366, 150875, 303006, 604511, 1195872, 2350614, 4587484, 8898857, 17154278, 32883109, 62679852, 118858190, 224238730, 421021209, 786793776, 1463796383, 2711552690, 5002097398, 9190449808
OFFSET
0,3
LINKS
FORMULA
a(n) ~ 1 / (2^(3/2) * 15^(5/48) * Pi^(1/12) * n^(29/48)) * exp(-Zeta'(-1) - Zeta(3)/(4*Pi^2) - 75*Zeta(3)^3 / Pi^8 - 15^(5/4) * Zeta(3)^2 / (2*Pi^5) * n^(1/4) - sqrt(15) * Zeta(3) / Pi^2 * sqrt(n) + 4*Pi / (3*15^(1/4)) * n^(3/4)), where Zeta(3) = A002117, Zeta'(-1) = A084448 = 1/12 - log(A074962).
G.f.: exp(Sum_{k>=1} (sigma_3(k) - sigma_2(k))*x^k/k). - Ilya Gutkovskiy, Aug 22 2018
MATHEMATICA
nmax=40; CoefficientList[Series[Product[1/(1-x^k)^(k*(k-1)), {k, 1, nmax}], {x, 0, nmax}], x]
Clear[a]; a[n_]:= a[n] = 1/n*Sum[(DivisorSigma[3, k]-DivisorSigma[2, k])*a[n-k], {k, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 100}] (* Vaclav Kotesovec, Apr 11 2016, following a suggestion of George Beck *)
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, May 27 2015
STATUS
approved

Search completed in 0.011 seconds