login
Numbers divisible by prime(k)^k for each prime index k.
18

%I #20 Nov 23 2020 02:05:48

%S 1,2,4,8,9,16,18,27,32,36,54,64,72,81,108,125,128,144,162,216,243,250,

%T 256,288,324,432,486,500,512,576,625,648,729,864,972,1000,1024,1125,

%U 1152,1250,1296,1458,1728,1944,2000,2048,2187,2250,2304,2401,2500,2592

%N Numbers divisible by prime(k)^k for each prime index k.

%C These are a kind of self-describing numbers (cf. A001462, A304679).

%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. The prime signature of a number is the multiset of multiplicities (or exponents) in its prime factorization.

%C Also Heinz numbers of integer partitions where the multiplicity of k is at least k (A117144). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

%H Amiram Eldar, <a href="/A324525/b324525.txt">Table of n, a(n) for n = 1..10000</a>

%F Closed under multiplication.

%F Sum_{n>=1} 1/a(n) = Product_{k>=1} 1 + 1/(prime(k)^(k-1) * (prime(k)-1)) = 2.35782843100111139159... - _Amiram Eldar_, Nov 23 2020

%e The sequence of terms together with their prime indices begins as follows. For example, 36 = prime(1) * prime(1) * prime(2) * prime(2) is a term because the prime multiplicities are {2,2}, which are greater than or equal to the prime indices {1,2}.

%e 1: {}

%e 2: {1}

%e 4: {1,1}

%e 8: {1,1,1}

%e 9: {2,2}

%e 16: {1,1,1,1}

%e 18: {1,2,2}

%e 27: {2,2,2}

%e 32: {1,1,1,1,1}

%e 36: {1,1,2,2}

%e 54: {1,2,2,2}

%e 64: {1,1,1,1,1,1}

%e 72: {1,1,1,2,2}

%e 81: {2,2,2,2}

%e 108: {1,1,2,2,2}

%e 125: {3,3,3}

%e 128: {1,1,1,1,1,1,1}

%p q:= n-> andmap(i-> i[2]>=numtheory[pi](i[1]), ifactors(n)[2]):

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

%t Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k>=PrimePi[p]]&]

%t seq[max_] := Module[{ps = {2}, p, s = {1}, s1, s2, emax}, While[ps[[-1]]^Length[ps] < max, AppendTo[ps, NextPrime[ps[[-1]]]]]; Do[p = ps[[k]]; emax = Floor[Log[p, max]]; s1 = Join[{1}, p^Range[k, emax]]; s2 = Select[Union[Flatten[Outer[Times, s, s1]]], # <= max &]; s = Union[s, s2], {k, 1, Length[ps]}]; s]; seq[3000] (* _Amiram Eldar_, Nov 23 2020 *)

%Y Cf. A001156, A033461, A056239, A062457, A112798, A117144, A118914 (prime signature), A124010, A181819, A276078, A304679.

%Y Cf. A109298, A324524, A324571, A324572, A324587, A324588.

%Y Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

%K nonn

%O 1,2

%A _Gus Wiseman_, Mar 08 2019