facs[n_]:=如果[n<=1,{{}},连接@@表[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}];
physim[n_]:=前缀[Join@@Table[Tuples[physim/@f],{f,Select[facs[n],Length[#]>1&]}],n];
表[Sum[Length[physemi[Times@@Prime/@m]],{m,IntegerPartitions[n]}],{n,8}]