|
|
A206303
|
|
Expansion of e.g.f.: Product_{n>=1} (1 - x^(2*n-1))^(-1/(2*n-1)).
|
|
18
|
|
|
1, 1, 2, 8, 32, 184, 1264, 9568, 79232, 816128, 8769536, 101867776, 1322831872, 18122579968, 268425347072, 4436611211264, 73309336469504, 1303024044310528, 25235367455752192, 497968598916333568, 10431118327503650816, 234674470003955204096, 5359992446798535852032
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
a(n)/n! is the Euler transform of [1, 0, 1/3, 0, 1/5, 0, 1/7, 0, ...].
E.g.f.: A(x) = B(x) / sqrt(B(x^2)), where B(x) = e.g.f. of A028342.
E.g.f. A(x) satisfies: Product_{n>=0} A(x^(2^n))^(1/2^n) = e.g.f. of A028342.
E.g.f.: exp(Sum_{k>=1} (d(2*k) - d(k))*x^k/k), where d(k) = number of divisors of k (A000005). - Ilya Gutkovskiy, Sep 17 2018
|
|
EXAMPLE
|
G.f.: A(x) = 1 + x + 2*x^2/2! + 8*x^3/3! + 32*x^4/4! + 184*x^5/5! + ...
The e.g.f. equals the product:
A(x) = (1-x)^(-1) * (1-x^3)^(-1/3) * (1-x^5)^(-1/5) * (1-x^7)^(-1/7) * (1-x^9)^(-1/9) * (1-x^11)^(-1/11) * ...
|
|
MAPLE
|
with(numtheory):
b:= proc(n) option remember; `if`(n=0, 1, add(add(
`if`(d::odd, 1, 0), d=divisors(j))*b(n-j), j=1..n)/n)
end:
a:= n-> b(n)*n!:
|
|
MATHEMATICA
|
b[n_]:= b[n]= If[n==0, 1, Sum[Sum[If[OddQ[d], 1, 0], {d, Divisors[j]}]* b[n-j], {j, n}]/n];
a[n_]:= b[n]*n!;
|
|
PROG
|
(PARI) {a(n)=n!*polcoeff(prod(m=1, n, (1-x^(2*m-1)+x*O(x^n))^(-1/(2*m-1))), n)}
for(n=0, 31, print1(a(n), ", "))
(Magma)
m:=40;
f:= func< x | (&*[1/(1 - x^(2*n-1))^(1/(2*n-1)) : n in [1..m+10]]) >;
R<x>:=PowerSeriesRing(Rationals(), m);
(SageMath)
m=40
def f(x): return 1/product( (1 - x^(2*n-1))^(1/(2*n-1)) for n in range(1, m+11) )
P.<x> = PowerSeriesRing(QQ, prec)
return P( f(x) ).egf_to_ogf().list()
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|