A323306 Heinz numbers of integer partitions that can be arranged into a matrix with equal row-sums and equal column-sums.
1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 36, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 100, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193
Offset: 1
Keywords
Examples
6480 belongs to the sequence because it is the Heinz number of (3,2,2,2,2,1,1,1,1), which can be arranged in the following ways: [1 1 3] [1 2 2] [1 2 2] [1 3 1] [2 1 2] [2 1 2] [2 2 1] [2 2 1] [3 1 1] [2 2 1] [1 2 2] [3 1 1] [2 1 2] [1 3 1] [2 1 2] [1 1 3] [2 2 1] [1 2 2] [2 2 1] [3 1 1] [1 2 2] [2 1 2] [2 1 2] [1 3 1] [2 2 1] [1 1 3] [1 2 2]
Crossrefs
Programs
-
Mathematica
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]]; ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&]; Select[Range[100],!Select[ptnmats[#],And[SameQ@@Total/@#,SameQ@@Total/@Transpose[#]]&]=={}&]
Comments