OFFSET
1,1
COMMENTS
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). A prime index of n is a number m such that prime(m) divides n. A multiset m whose distinct elements are m_1, m_2, ..., m_k with multiplicities y_1, y_2, ..., y_k is reducible if either m is of size 1 or gcd(m_1,...,m_k) = 1 and the multiset {y_1,...,y_k} is also reducible.
EXAMPLE
60 has relatively prime prime indices {1,1,2,3} with multiplicities {1,1,2} corresponding to A181819(90) = 12. 12 has relatively prime prime indices {1,1,2} with multiplicities {1,2} corresponding to A181819(12) = 6. 6 has relatively prime prime indices {1,2} with multiplicities {1,1} corresponding to A181819(6) = 4. 4 has relatively prime prime indices {1,1} with multiplicities {2} corresponding to A181819(4) = 3. 3 is prime, so we conclude that 60 belongs to the sequence.
MATHEMATICA
rdzQ[n_]:=And[n>1, Or[PrimeQ[n], And[rdzQ[Times@@Prime/@FactorInteger[n][[All, 2]]], GCD@@PrimePi/@FactorInteger[n][[All, 1]]==1]]];
Select[Range[50], rdzQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jun 22 2018
STATUS
approved
