OFFSET
1,1
COMMENTS
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.
The asymptotic density of this sequence is Sum_{n>=1} d(n) = 0.13743128989284883653... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
4 = 2^2 is a term since its least prime factor, 2, is equal to its exponent.
MAPLE
filter:= proc(n) local F;
F:= sort(ifactors(n)[2], (s, t) -> s[1]<t[1]);
F[1][1]=F[1][2]
end proc:
select(filter, [$2..1000]); # Robert Israel, Sep 22 2023
MATHEMATICA
q[n_] := Equal @@ FactorInteger[n][[1]]; Select[Range[2, 400], q]
PROG
(PARI) is(n) = n > 1 && #Set(factor(n)[1, ]) == 1;
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Sep 22 2023
STATUS
approved