OFFSET
1,4
COMMENTS
a(n) is also the number of ways to choose a perfect divisor of each factor in a factorization of n.
FORMULA
Dirichlet g.f.: 1/Product_{n > 1}(1 - A089723(n)/n^s).
EXAMPLE
The a(24) = 10 twice-factorizations are:
(2)*(2)*(2)*(3), (2)*(3)*(2*2), (3)*(2*2*2)
(2)*(2)*(6), (2*2)*(6),
(2)*(3)*(4),
(2)*(12),
(3)*(8),
(4)*(6),
(24).
MATHEMATICA
facs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Table[Sum[Product[Length[Divisors[GCD@@FactorInteger[d][[All, 2]]]], {d, f}], {f, facs[n]}], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Nov 29 2017
STATUS
approved