OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: exp( Sum_{n>=1} x^n * Sum_{d|n} (1+x)^d / d ).
G.f.: Product {n >= 1} 1/(1 - (1 + x)*x^n). - Peter Bala, Jan 20 2015
a(n) ~ phi^(n+1) / (sqrt(5)* A276987), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Jul 16 2019
EXAMPLE
G.f.: A(x) = 1 + x + 3*x^2 + 6*x^3 + 12*x^4 + 23*x^5 + 43*x^6 + 79*x^7 +...
where
log(A(x)) = x*(1+x)/(1-x) + x^2*(1+x)^2/(2*(1-x^2)) + x^3*(1+x)^3/(3*(1-x^3)) + x^4*(1+x)^4/(4*(1-x^4)) + x^5*(1+x)^5/(5*(1-x^5)) +...
Explicitly,
log(A(x)) = x + 5*x^2/2 + 10*x^3/3 + 17*x^4/4 + 26*x^5/5 + 38*x^6/6 + 57*x^7/7 + 81*x^8/8 + 118*x^9/9 + 180*x^10/10 +...
MATHEMATICA
nmax = 50; CoefficientList[Series[Product[1/(1 - x^k*(1 + x)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jul 31 2021 *)
PROG
(PARI) {a(n)=polcoeff(exp(sum(m=1, n+1, x^m/m*(1+x)^m/(1-x^m +x*O(x^n)) )), n)}
for(n=0, 50, print1(a(n), ", "))
(PARI) {a(n)=polcoeff(exp(sum(m=1, n+1, x^m*sumdiv(m, d, (1+x +x*O(x^n))^d/d) )), n)}
for(n=0, 50, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 19 2013
STATUS
approved