OFFSET
1,2
COMMENTS
The terms 70 and 90 show that sequence is strictly different from A042965. - Andrew S. Plewe, Jun 13 2007
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Range@94, Length@Union@Mod[Range[#]^2, #] <= EulerPhi[#] &] (* Ivan Neretin, Dec 14 2016 *)
f1[p_, e_] := Floor[p^(e+1)/(2p + 2)] + 1; f1[2, e_] := Floor[2^e/6] + 2; f[p_, e_] := f1[p, e]/((p-1) * p^(e-1)); q[1] = True; q[k_] := Times @@ f @@@ FactorInteger[k] <= 1; Select[Range[100], q] (* Amiram Eldar, Nov 11 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Max Alekseyev, Sep 18 2006
STATUS
approved