|
|
|
|
0, 1, 2, 5, 9, 15, 27, 42, 65, 99, 148, 214, 308, 435, 605, 839, 1145, 1548, 2080, 2769, 3659, 4812, 6278, 8145, 10518, 13506, 17257, 21961, 27821, 35095, 44117, 55243, 68928, 85735, 106285, 131357, 161893, 198944, 243817, 298060, 363446
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
a(n) ~ 3^(1/4) * log(2) * exp(Pi*sqrt(2*n/3)) / (2^(7/4) * Pi * n^(1/4)). - Vaclav Kotesovec, Oct 09 2018
|
|
MAPLE
|
b:= proc(n) option remember; `if`(n=0, 1, add(add(d*[0, 1][1+
irem(d, 2)], d=numtheory[divisors](j))*b(n-j), j=1..n)/n)
end:
g:= proc(n, i) option remember; `if`(i*(i+1)/2<n, 0, `if`(n=0, [1, 0],
add((l->[l[1], l[2]+l[1]*j])(g(n-i*j, i-1)), j=0..min(n/i, 1))))
end:
a:= n-> add(b(n-j)*g(j$2)[2], j=0..n):
|
|
MATHEMATICA
|
Table[Sum[Count[#, _?OddQ] - Count[#, _?EvenQ] &@ Divisors@ k SeriesCoefficient[QPochhammer[q, q^2]^-2, {q, 0, #}] &[n - k], {k, n}], {n, 0, 40}] (* Michael De Vlieger, Jun 18 2016, after Michael Somos at A022567 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|