A069042 Numbers k such that A000364(k) == 1 (mod k^2).
1, 2, 17, 37, 1153, 1303
Offset: 1
Programs
-
Mathematica
Select[Range[1, 10000], Divisible[Abs[EulerE[2*#]] - 1, #^2] &] (* Amiram Eldar, Jun 03 2017 *)
-
PARI
is(k) = Mod(abs(eulerfrac(2*k)), k^2) == 1; \\ Jinyuan Wang, Mar 13 2025
Extensions
More terms from Hans Havermann, Apr 07 2003
a(1) prepended by Jinyuan Wang, Mar 13 2025