OFFSET
1,32
COMMENTS
A set partition is inseparable iff the underlying set has no permutation whose adjacent elements all belong to different blocks. Note that inseparability only depends on the sizes of the blocks.
FORMULA
a(2^n) = A111133(n).
EXAMPLE
The prime indices of 9216 are {1,1,1,1,1,1,1,1,1,1,2,2}, with a(9216) = 2 choices: {{2},{1,4,5}} and {{2},{1,3,6}}. The other 4 disjoint families {{2},{10}}, {{2},{4,6}}, {{2},{3,7}}, {{2},{1,9}} are separable.
The prime indices of 15552 are {1,1,1,1,1,1,2,2,2,2,2}, with a(15552) = 1 choice: {{5},{1,2,3}}. The other 5 disjoint families {{5},{6}}, {{5},{2,4}}, {{6},{2,3}}, {{6},{1,4}}, {{1,5},{2,3}} are separable.
MATHEMATICA
disjointFamilies[y_]:=Union[Sort/@Select[Tuples[IntegerPartitions/@Length/@Split[y]], UnsameQ@@Join@@#&]];
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
seps[ptn_, fir_]:=If[Total[ptn]==1, {{fir}}, Join@@Table[Prepend[#, fir]&/@seps[MapAt[#-1&, ptn, fir], nex], {nex, Select[DeleteCases[Range[Length[ptn]], fir], ptn[[#]]>0&]}]];
seps[ptn_]:=If[Total[ptn]==0, {{}}, Join@@(seps[ptn, #]&/@Range[Length[ptn]])];
Table[Length[Select[disjointFamilies[prix[n]], seps[Length/@#]=={}&]], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 31 2025
STATUS
approved
