OFFSET
1,2
COMMENTS
MATHEMATICA
hwt[n_]:=Total[Cases[FactorInteger[n], {p_, k_}:>PrimePi[p]*k]];
facs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Table[Sum[Length[Select[IntegerPartitions[n], Select[facs[Times@@Prime/@#], Length[#]==k&&SameQ@@hwt/@#&]!={}&]], {k, Divisors[n]}], {n, 1, 10}]
PROG
(PARI) T(n, d) = my(v=partitions(n/d), w=List([])); forvec(s=vector(d, i, [1, #v]), listput(w, vecsort(concat(vector(d, i, v[s[i]])))), 1); #Set(w);
a(n) = sumdiv(n, d, T(n, d)); \\ Jinyuan Wang, Feb 13 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Apr 20 2024
EXTENSIONS
More terms from Jinyuan Wang, Feb 13 2025
STATUS
approved