OFFSET
0,2
COMMENTS
a(n) is Euler transform of A010734. - Alois P. Heinz, Oct 17 2008
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Alois P. Heinz)
P. Nataf, M. Lajkó, A. Wietek, K. Penc, F. Mila, A. M. Läuchli, Chiral spin liquids in triangular lattice SU (N) fermionic Mott insulators with artificial gauge fields, arXiv preprint arXiv:1601.00958 [cond-mat.quant-gas], 2016.
N. J. A. Sloane, Transforms
FORMULA
a(n) ~ 3^(5/2) * exp(Pi * sqrt(6*n)) / (256 * n^3). - Vaclav Kotesovec, Feb 28 2015
a(0) = 1, a(n) = (9/n)*Sum_{k=1..n} A000203(k)*a(n-k) for n > 0. - Seiichi Manyama, Mar 27 2017
G.f.: exp(9*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 06 2018
MAPLE
with(numtheory): a:= proc(n) option remember; `if`(n=0, 1, add(add(d*9, d=divisors(j)) *a(n-j), j=1..n)/n) end: seq(a(n), n=0..40); # Alois P. Heinz, Oct 17 2008
MATHEMATICA
nmax=50; CoefficientList[Series[Product[1/(1-x^k)^9, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Feb 28 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved