OFFSET
1,1
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. The multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}. This sequence lists all MM-numbers of not necessarily strict connected antichains of multisets spanning an initial interval of positive integers.
EXAMPLE
The sequence of multisystems whose MM-numbers belong to the sequence begins:
2: {{}}
3: {{1}}
7: {{1,1}}
9: {{1},{1}}
13: {{1,2}}
19: {{1,1,1}}
27: {{1},{1},{1}}
37: {{1,1,2}}
49: {{1,1},{1,1}}
53: {{1,1,1,1}}
61: {{1,2,2}}
81: {{1},{1},{1},{1}}
89: {{1,1,1,2}}
91: {{1,1},{1,2}}
113: {{1,2,3}}
131: {{1,1,1,1,1}}
151: {{1,1,2,2}}
169: {{1,2},{1,2}}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
normQ[sys_]:=Or[Length[sys]==0, Union@@sys==Range[Max@@Max@@sys]];
zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]], 2], And[Less@@#, GCD@@s[[#]]]>1&]}, If[c=={}, s, zsm[Sort[Append[Delete[s, List/@c[[1]]], LCM@@s[[c[[1]]]]]]]]];
stableQ[u_, Q_]:=!Apply[Or, Outer[#1=!=#2&&Q[#1, #2]&, u, u, 1], {0, 1}];
Select[Range[200], And[normQ[primeMS/@primeMS[#]], stableQ[primeMS[#], Divisible], Length[zsm[primeMS[#]]]==1]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Dec 16 2018
STATUS
approved