OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..3217
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(d*
max(1, 2^(d-1)-1), d=numtheory[divisors](j)), j=1..n)/n)
end:
a:= n-> b(n)-max(1, 2^(n-1)-1):
seq(a(n), n=1..35); # Alois P. Heinz, Jan 05 2022
MATHEMATICA
b[n_] := b[n] = If[n == 0, 1, Sum[b[n - j]*Sum[d*
Max[1, 2^(d-1) - 1], {d, Divisors[j]}], {j, 1, n}]/n];
a[n_] := b[n] - Max[1, 2^(n-1)-1];
Table[a[n], {n, 1, 35}] (* Jean-François Alcover, Mar 11 2022, Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Salah Uddin Mohammad, Nov 19 2021
STATUS
approved