|
|
A335283
|
|
a(n) = 1 + Sum_{d|n, n/d odd, d < n} a(d).
|
|
2
|
|
|
1, 1, 2, 1, 2, 2, 2, 1, 4, 2, 2, 2, 2, 2, 6, 1, 2, 4, 2, 2, 6, 2, 2, 2, 4, 2, 8, 2, 2, 6, 2, 1, 6, 2, 6, 4, 2, 2, 6, 2, 2, 6, 2, 2, 16, 2, 2, 2, 4, 4, 6, 2, 2, 8, 6, 2, 6, 2, 2, 6, 2, 2, 16, 1, 6, 6, 2, 2, 6, 6, 2, 4, 2, 2, 16, 2, 6, 6, 2, 2, 16, 2, 2, 6, 6, 2, 6, 2, 2, 16, 6, 2, 6, 2, 6, 2, 2, 4, 16, 4, 2, 6, 2, 2, 26
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
LINKS
|
|
|
FORMULA
|
G.f. A(x) satisfies: A(x) = x / (1 - x) + Sum_{k>=2} A(x^(2*k-1)).
G.f.: x / (1 - x) + Sum_{n>=1} a(n) * x^(3*n) / (1 - x^(2*n)).
a(1) = 1; a(2*n) = a(n), a(2*n+1) = 2 * A074206(2*n+1).
|
|
MATHEMATICA
|
a[n_] := 1 + Sum[If[d < n && OddQ[n/d], a[d], 0], {d, Divisors[n]}]; Table[a[n], {n, 85}]
terms = 85; A[_] = 0; Do[A[x_] = x/(1 - x) + Sum[A[x^(2 k - 1)], {k, 2, terms}] + O[x]^(terms + 1) // Normal, terms + 1]; CoefficientList[A[x], x] // Rest
|
|
PROG
|
(PARI)
up_to = 20000;
A335283lista(up_to) = {my(v = vector(up_to)); for(n=1, up_to, v[n] = 1 + sumdiv(n, d, if(d<n&&((n/d)%2), v[d], 0))); (v); };
v335283 = A335283lista(up_to);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|