OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..17948
EXAMPLE
a(4) = 13 is in the sequence because 19^12 == 1 (mod 13^2), and 1 < 19 < 13^2-1.
MAPLE
Aeven:= remove(t -> igcd(t-1, numtheory:-phi(t^2))=1, {seq(i, i=2..1000, 2)}):
Aodd:= {seq(i, i=3..1000, 2)} minus {seq(3^i, i=0..floor(log[3](1000)))}:
sort(convert(Aeven union Aodd, list));
MATHEMATICA
Aeven = DeleteCases[Range[2, 1000, 2], t_ /; GCD[t-1, EulerPhi[t^2]] == 1];
Aodd = Complement[Range[3, 1000, 2], Table[3^i, {i, 0, Floor[Log[3, 1000]]} ]];
Union[Aeven, Aodd] (* Jean-François Alcover, Apr 24 2019, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel and Thomas Ordowski, Dec 28 2016
STATUS
approved