A381991 Numbers whose prime indices have a unique multiset partition into constant multisets with distinct sums.
1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79
Offset: 1
Keywords
Examples
The prime indices of 270 are {1,2,2,2,3}, and there are two multiset partitions into constant multisets with distinct sums: {{1},{2},{3},{2,2}} and {{1},{3},{2,2,2}}, so 270 is not in the sequence. The prime indices of 300 are {1,1,2,3,3}, of which there are no multiset partitions into constant multisets with distinct sums, so 300 is not in the sequence. The prime indices of 360 are {1,1,1,2,2,3}, of which there is only one multiset partition into constant multisets with distinct sums: {{1},{1,1},{3},{2,2}}, so 360 is in the sequence. The terms together with their prime indices begin: 1: {} 2: {1} 3: {2} 4: {1,1} 5: {3} 6: {1,2} 7: {4} 9: {2,2} 10: {1,3} 11: {5} 13: {6} 14: {1,4} 15: {2,3} 17: {7} 18: {1,2,2} 19: {8} 20: {1,1,3} 21: {2,4} 22: {1,5} 23: {9} 24: {1,1,1,2} 25: {3,3}
Crossrefs
Programs
-
Mathematica
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]]; pfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[pfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]]; Select[Range[100],Length[Select[pfacs[#],UnsameQ@@hwt/@#&]]==1&]
Comments