cp's OEIS Frontend

A365883 Numbers k whose least prime divisor is equal to its exponent in the prime factorization of k.

%I A365883 #15 Jan 10 2024 05:05:52
%S A365883 4,12,20,27,28,36,44,52,60,68,76,84,92,100,108,116,124,132,135,140,
%T A365883 148,156,164,172,180,188,189,196,204,212,220,228,236,244,252,260,268,
%U A365883 276,284,292,297,300,308,316,324,332,340,348,351,356,364,372,380,388,396
%N A365883 Numbers k whose least prime divisor is equal to its exponent in the prime factorization of k.
%C A365883 Numbers k such that A020639(k) = A051904(k).
%C A365883 The asymptotic density of terms with least prime factor prime(n) (within all the positive integers) is d(n) = (1/prime(n)^prime(n) - 1/prime(n)^(prime(n)+1)) * Product_{k=1..(n-1)} (1-1/prime(k)). For example, for n = 1, 2, 3, 4 and 5, d(n) = 1/8, 1/81, 4/46875, 8/28824005 and 16/21968998637047.
%C A365883 The asymptotic density of this sequence is Sum_{n>=1} d(n) = 0.13743128989284883653... .
%H A365883 Amiram Eldar, <a href="/A365883/b365883.txt">Table of n, a(n) for n = 1..10000</a>
%e A365883 4 = 2^2 is a term since its least prime factor, 2, is equal to its exponent.
%p A365883 filter:= proc(n) local F;
%p A365883    F:= sort(ifactors(n)[2],(s,t) -> s[1]<t[1]);
%p A365883    F[1][1]=F[1][2]
%p A365883 end proc:
%p A365883 select(filter, [$2..1000]); # _Robert Israel_, Sep 22 2023
%t A365883 q[n_] := Equal @@ FactorInteger[n][[1]]; Select[Range[2, 400], q]
%o A365883 (PARI) is(n) = n > 1 && #Set(factor(n)[1,]) == 1;
%Y A365883 Cf. A020639, A051904.
%Y A365883 Subsequence of A100717 and A365889.
%Y A365883 Subsequences: A017113, A365884, A365885.
%K A365883 nonn,easy
%O A365883 1,1
%A A365883 _Amiram Eldar_, Sep 22 2023