OFFSET
1,2
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 Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). A partition with no two distinct parts divisible is said to be stable, and a partition with no two distinct parts relatively prime is said to be intersecting, so these are Heinz numbers of stable intersecting partitions.
EXAMPLE
The sequence of terms together with their prime indices begins:
1: {}
2: {1}
3: {2}
4: {1,1}
5: {3}
7: {4}
8: {1,1,1}
9: {2,2}
11: {5}
13: {6}
16: {1,1,1,1}
17: {7}
19: {8}
23: {9}
25: {3,3}
27: {2,2,2}
29: {10}
31: {11}
32: {1,1,1,1,1}
37: {12}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
stableQ[u_, Q_]:=!Apply[Or, Outer[#1=!=#2&&Q[#1, #2]&, u, u, 1], {0, 1}];
Select[Range[100], stableQ[Union[primeMS[#]], GCD[#1, #2]==1&]&&stableQ[Union[primeMS[#]], Divisible]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Nov 12 2019
STATUS
approved