OFFSET
0,3
FORMULA
a(n) / A046161(n) ~ exp(sqrt(n/6)*Pi) / (4*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Apr 12 2021
EXAMPLE
1, 1/2, 3/8, 13/16, 67/128, 239/256, 1031/1024, 2501/2048, 36579/32768, 109915/65536, 468653/262144, 1043851/524288, ...
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(
`if`(d::odd, d/2, 0), d=numtheory[divisors](j)), j=1..n)/n)
end:
a:= n-> numer(b(n)):
seq(a(n), n=0..25); # Alois P. Heinz, Apr 12 2021
MATHEMATICA
nmax = 25; CoefficientList[Series[Product[(1 + x^k)^(1/2), {k, 1, nmax}], {x, 0, nmax}], x] // Numerator
a[n_] := a[n] = If[n == 0, 1, (1/(2 n)) Sum[Sum[Mod[d, 2] d, {d, Divisors[k]}] a[n - k], {k, 1, n}]]; Table[a[n], {n, 0, 25}] // Numerator
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Ilya Gutkovskiy, Apr 07 2021
STATUS
approved