login
A228647
a(n) = A001609(n^2) for n>=1, where g.f. of A001609 is x*(1+3*x^2)/(1-x-x^3).
2
1, 5, 31, 453, 14131, 946781, 136250983, 42115660581, 27961563559891, 39874307297033165, 122134599693975367423, 803522677430288749340325, 11354589189995520431547851761, 344634362031276605039944979868611, 22467750416780812361715214948922598721, 3146114090698891414621617889648190060326821
OFFSET
1,2
COMMENTS
A001609 forms the logarithmic derivative of Narayana's cows sequence A000930.
FORMULA
Equals the logarithmic derivative of A228648.
EXAMPLE
L.g.f.: L(x) = x + 5*x^2/2 + 31*x^3/3 + 453*x^4/4 + 14131*x^5/5 +...
where
exp(L(x)) = 1 + x + 3*x^2 + 13*x^3 + 128*x^4 + 2974*x^5 + 161048*x^6 + 19632276*x^7 +...+ A228648(n)*x^n +...
PROG
(PARI) {A001609(n)=n*polcoeff(-log(1-x-x^3 +x*O(x^n)), n)}
{a(n)=A001609(n^2)}
for(n=1, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 28 2013
STATUS
approved