OFFSET
1,18
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..101025
FORMULA
G.f. A(x) satisfies: A(x) = x + Sum_{k>=2} A(x^(k*(k+1)/2)).
a(1) = 1; a(n) = Sum_{d|n, d<n} A010054(n/d)*a(d).
MATHEMATICA
terms = 100; A[_] = 0; Do[A[x_] = x + Sum[A[x^(k (k + 1)/2)], {k, 2, terms}] + O[x]^(terms + 1) // Normal, terms + 1]; Rest[CoefficientList[A[x], x]]
a[n_] := If[n == 1, n, Sum[If[(d < n && IntegerQ[Sqrt[8 n/d + 1]]), a[d], 0], {d, Divisors[n]}]]; Table[a[n], {n, 1, 100}]
PROG
(PARI)
A010054(n) = issquare(1 + 8*n);
A308264(n) = if(1==n, n, sumdiv(n, d, if(d<n, A010054(n/d)*A308264(d), 0))); \\ Antti Karttunen, Oct 05 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 17 2019
EXTENSIONS
Data section extended to 105 terms by Antti Karttunen, Oct 05 2023
STATUS
approved