OFFSET
1,2
COMMENTS
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
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.
FORMULA
a(n) = n iff n belongs to A319630.
EXAMPLE
The prime indices of 12 are {1,1,2}, with sparse submultisets {{},{1},{2},{1,1}}, with Heinz numbers {1,2,3,4}, so a(12) = 4.
The prime indices of 36 are {1,1,2,2}, with sparse submultisets {{},{1},{2},{1,1},{2,2}}, with Heinz numbers {1,2,3,4,9}, so a(36) = 9.
The prime indices of 462 are {1,2,4,5}, with sparse submultisets {{},{1},{2},{4},{5},{1,4},{2,4},{1,5},{2,5}}, with Heinz numbers {1,2,3,7,11,14,21,22,33}, so a(462) = 33.
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Max@@Select[Divisors[n], FreeQ[Differences[prix[#]], 1]&], {n, 100}]
CROSSREFS
A000005 counts divisors (or submultisets of prime indices).
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 05 2025
STATUS
approved
