%I #10 Feb 02 2021 04:35:45
%S 1,3,5,7,11,13,15,17,19,21,23,25,29,31,33,35,37,39,41,43,47,49,51,53,
%T 55,57,59,61,65,67,69,71,73,75,77,79,83,85,87,89,91,93,95,97,101,103,
%U 105,107,109,111,113,115,119,121,123,127,129,131,133,137,139,141
%N Numbers in whose prime factorization the exponent of prime(k) is less than k for all prime indices k.
%C 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.
%C The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions where each part k appears fewer than k times. Such partitions are counted by A087153.
%C The asymptotic density of this sequence is Product_{k>=1} (1 - 1/prime(k)^k) = 0.44070243286030291209... - _Amiram Eldar_, Feb 02 2021
%H Amiram Eldar, <a href="/A325128/b325128.txt">Table of n, a(n) for n = 1..10000</a>
%e The sequence of terms together with their prime indices begins:
%e 1: {}
%e 3: {2}
%e 5: {3}
%e 7: {4}
%e 11: {5}
%e 13: {6}
%e 15: {2,3}
%e 17: {7}
%e 19: {8}
%e 21: {2,4}
%e 23: {9}
%e 25: {3,3}
%e 29: {10}
%e 31: {11}
%e 33: {2,5}
%e 35: {3,4}
%e 37: {12}
%e 39: {2,6}
%e 41: {13}
%e 43: {14}
%e 47: {15}
%e 49: {4,4}
%t Select[Range[100],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k<PrimePi[p]]&]
%Y Cf. A056239, A062457, A087153, A109298, A112798, A115584, A118914, A276078.
%Y Cf. A324524, A324525, A324571, A325127, A325130.
%K nonn
%O 1,2
%A _Gus Wiseman_, Apr 01 2019