OFFSET
1,2
COMMENTS
1 is included somewhat arbitrarily. 1 has no nonzero exponents in its prime factorization, but it also has no prime factorization exponents that are not coprime to 1.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 9, 75, 746, 7433, 74270, 742714, 7427050, 74270567, 742705640, 7427055214, ... . Apparently, the asymptotic density of this sequence exists and equals 0.742705... . - Amiram Eldar, Feb 11 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
40 has the prime-factorization 2^3 * 5^1. The exponents are therefore 3 and 1. Since both 3 and 1 are coprime to 40, then 40 is included in the sequence.
MAPLE
filter:= proc(n) local E;
E:= map(t -> t[2], ifactors(n)[2]);
andmap(t -> igcd(t, n)=1, E)
end proc:
select(filter, [$1..200]); # Robert Israel, Oct 24 2019
MATHEMATICA
Select[Range[100], GCD[Times @@ Table[FactorInteger[ # ][[i, 2]], {i, 1, Length[FactorInteger[ # ]]}], # ] == 1 &] (* Stefan Steinerberger, Sep 15 2008 *)
PROG
(PARI) is(n) = {my(e = factor(n)[, 2]); for(i=1, #e, if(gcd(e[i], n) > 1, return(0))); 1; }; \\ Amiram Eldar, Feb 11 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Sep 11 2008
EXTENSIONS
More terms from Stefan Steinerberger, Sep 15 2008
STATUS
approved