login
Numbers where the sum of distinct prime indices (A066328) is equal to the number of prime factors counted with multiplicity (A001222).
13

%I #11 Apr 13 2019 09:02:31

%S 1,2,9,12,18,40,100,112,125,240,250,352,360,392,405,540,600,672,675,

%T 810,832,900,1008,1125,1350,1372,1500,1512,1701,1875,1936,2112,2176,

%U 2240,2250,2268,2352,2401,3168,3402,3528,3750,3969,4752,4802,4864,4992,5292

%N Numbers where the sum of distinct prime indices (A066328) is equal to the number of prime factors counted with multiplicity (A001222).

%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. For example, 540 = prime(1)^2 * prime(2)^3 * prime(3)^1 has sum of distinct prime indices 1 + 2 + 3 = 6, while the number of prime factors counted with multiplicity is 2 + 3 + 1 = 6, so 540 belongs to the sequence.

%C Also Heinz numbers of the integer partitions counted by A114638. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

%F A066328(a(n)) = A001222(a(n)).

%e The sequence of terms together with their prime indices begins:

%e 1: {}

%e 2: {1}

%e 9: {2,2}

%e 12: {1,1,2}

%e 18: {1,2,2}

%e 40: {1,1,1,3}

%e 100: {1,1,3,3}

%e 112: {1,1,1,1,4}

%e 125: {3,3,3}

%e 240: {1,1,1,1,2,3}

%e 250: {1,3,3,3}

%e 352: {1,1,1,1,1,5}

%e 360: {1,1,1,2,2,3}

%e 392: {1,1,1,4,4}

%e 405: {2,2,2,2,3}

%e 540: {1,1,2,2,2,3}

%e 600: {1,1,1,2,3,3}

%e 672: {1,1,1,1,1,2,4}

%p with(numtheory):

%p q:= n-> is(add(pi(p), p=factorset(n))=bigomega(n)):

%p select(q, [$1..5600])[]; # _Alois P. Heinz_, Mar 07 2019

%t Select[Range[1000],Total[PrimePi/@First/@FactorInteger[#]]==PrimeOmega[#]&]

%Y Cf. A001221, A001222, A056239, A066328, A112798, A114638, A117144, A276078.

%Y Cf. A109298, A324524, A324525, A324570, A324571, A324572.

%K nonn

%O 1,2

%A _Gus Wiseman_, Mar 07 2019