OFFSET
1,2
COMMENTS
EXAMPLE
For n = 600 the unique multiset partition is {{1},{1,3},{1,2,3}}. The unique factorization is 2*10*30.
MATHEMATICA
hwt[n_]:=Total[Cases[FactorInteger[n], {p_, k_}:>PrimePi[p]*k]];
sfacs[n_]:=If[n<=1, {{}}, Join@@Table[(Prepend[#, d]&)/@Select[sfacs[n/d], Min@@#>=d&], {d, Select[Rest[Divisors[n]], SquareFreeQ]}]];
Select[Range[100], Length[Select[sfacs[#], UnsameQ@@hwt/@#&]]==1&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 12 2025
STATUS
approved
