OFFSET
1,2
COMMENTS
First differs from A080197 at n = 28.
Erdős et al. (2008) proved that the asymptotic density of numbers k such that gcd(k, phi(k)) > (log(log(k)))^u for a real number u > 0 is equal to exp(-gamma) * Integral_{t=u..oo} rho(t) dt, where rho(t) is the Dickman-de Bruijn function and gamma is Euler's constant (A001620). For this sequence u = 1, and therefore its asymptotic density is 1 - exp(-gamma) = 0.43854... (A227242).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Paul Erdős, Florian Luca and Carl Pomerance, On the proportion of numbers coprime to a given integer, in: J.-M. De Koninck, A. Granville and F. Luca (eds.), Anatomy of Integers, AMS, 2008, pp. 47-64.
Wikipedia, Dickman function.
EXAMPLE
16 is a term since gcd(16, phi(16)) = gcd(16, 8) = 8 > log(log(16)) = 1.0197...
17 is not a term since gcd(17, phi(17)) = gcd(17, 16) = 1 < log(log(17)) = 1.0414...
MATHEMATICA
Select[Range[100], GCD[#, EulerPhi[#]] > Log[Log[#]] &]
PROG
(PARI) isok(k) = (k==1) || (gcd(k, eulerphi(k)) > log(log(k))); \\ Michel Marcus, Feb 19 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 18 2021
STATUS
approved