OFFSET
1,1
COMMENTS
Conjecture: a(2*n+1) == 2 (mod 4), a(2*n) == 0 (mod 4).
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..730 (first 180 terms from Paul D. Hanna)
FORMULA
a(n) ~ c * exp(Pi*n) / n, where c = 0.107862... . - Vaclav Kotesovec, Nov 29 2015
EXAMPLE
L.g.f.: L(x) = 2*x + 24*x^2/2 + 386*x^3/3 + 6832*x^4/4 + 128442*x^5/5 + 2505720*x^6/6 + 50153770*x^7/7 + 1022997344*x^8/8 +...
such that exponentiation yields the g.f. of A219430:
exp(L(x)) = 1 + 2*x + 14*x^2 + 154*x^3 + 2062*x^4 + 31066*x^5 + 504886*x^6 + 8652402*x^7 + 154208270*x^8 +...+ A219430(n)*x^n +...
MATHEMATICA
nmax = 20; A015128 = Rest[CoefficientList[Series[Exp[Sum[(DivisorSigma[1, 2*n] - DivisorSigma[1, n])*(x^n/n), {n, 1, nmax^2}]], {x, 0, nmax^2}], x]]; A219430 = Table[A015128[[n^2]], {n, 1, nmax}]; Rest[CoefficientList[Series[Log[1 + Sum[A219430[[k]]*x^k, {k, 1, nmax}]], {x, 0, nmax}], x] * Range[0, nmax]] (* Vaclav Kotesovec, Nov 28 2015 *)
nmax = 20; Rest[CoefficientList[Series[Log[1 + Sum[Sum[PartitionsP[k^2 - j]*PartitionsQ[j], {j, 0, k^2}] * x^k, {k, 1, nmax}]], {x, 0, nmax}], x] * Range[0, nmax]] (* Vaclav Kotesovec, Nov 28 2015 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 19 2012
STATUS
approved