|
|
|
|
1, 2, 2, 2, 4, 4, 2, 4, 5, 4, 6, 4, 4, 8, 4, 4, 8, 6, 6, 8, 8, 4, 6, 8, 5, 12, 8, 4, 12, 8, 6, 8, 8, 8, 12, 10, 4, 12, 8, 8, 16, 8, 6, 12, 12, 8, 10, 8, 9, 14, 12, 8, 12, 16, 8, 16, 8, 4, 18, 8, 12, 16, 10, 8, 16, 16, 6, 16, 16, 8, 14, 12, 8, 20, 14, 12, 16, 8, 10, 16, 17, 8, 18, 16, 8, 20, 12
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (Sum_{k>0} x^(k^2 + k - 1) / (1 - x^(2*k - 1))^2) / (Sum_{k>0} x^(k*(k - 1))). - Michael Somos, Jul 05 2015
|
|
EXAMPLE
|
G.f. = x + 2*x^2 + 2*x^3 + 2*x^4 + 4*x^5 + 4*x^6 + 2*x^7 + 4*x^8 + 5*x^9 + ...
|
|
PROG
|
(PARI) {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( sum(k=1, sqrtint(4*n + 9)\2, x^(k^2 + k - 2) / (1 - x^(2*k - 1))^2, A) / sum(k=1, sqrtint(4*n + 1)\2 + 1, x^(k^2 - k), A), n))}; /* Michael Somos, Jul 05 2015 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|