OFFSET
1,4
COMMENTS
a(n) is also the number of ways to choose a perfect divisor d|n and then a sequence of log_d(n) factorizations of d.
EXAMPLE
The a(16) = 10 twice-factorizations are (2*2*2*2), (2*2*4), (2*8), (4*4), (16), (2*2)*(2*2), (2*2)*(4), (4)*(2*2), (4)*(4), (2)*(2)*(2)*(2).
MATHEMATICA
postfacs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[postfacs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
a[n_]:=Sum[Length[postfacs[n^(1/g)]]^g, {g, Divisors[GCD@@FactorInteger[n][[All, 2]]]}];
Array[a, 50]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Nov 30 2017
STATUS
approved