OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..800
FORMULA
a(n) = A255704(2*n+1,n+1).
a(n) ~ c * d^n / sqrt(n), where d = A051491 = 2.955765285651994974714817524... and c = 0.70755335886284109851526791506579... . - Vaclav Kotesovec, Feb 28 2016
MAPLE
with(numtheory):
g:= proc(n, k) option remember; `if`(n=0, 1, add(add(d*(g(d-1, k)-
`if`(d=k, 1, 0)), d=divisors(j))*g(n-j, k), j=1..n)/n)
end:
a:= a-> g(2*n, n+1) -`if`(n=0, 0, g(2*n, n)):
seq(a(n), n=0..40);
MATHEMATICA
g[n_, k_] := g[n, k] = If[n == 0, 1, Sum[DivisorSum[j, #*(g[# - 1, k] - If[# == k, 1, 0]) &]*g[n - j, k], {j, 1, n}]/n];
a[n_] := g[2n, n+1] - If[n == 0, 0, g[2n, n]];
Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 24 2017, translated from Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Mar 02 2015
STATUS
approved