A177005 Numbers k such that k^k = k (mod prime(k)).
1, 4, 169, 391, 1546, 16761, 18278, 20201, 21775, 31120, 126882, 178465, 9502273, 10553442, 24677776, 56923413, 422766345, 1759518201, 4152696703, 6800832991, 14421293461, 106195400697, 667339219893, 915091102299
Offset: 1
Programs
-
Mathematica
Do[If[PowerMod[n, n, Prime[n]] == n, Print[n]], {n, 350000000}] -
PARI
is(k) = Mod(k, prime(k))^k == k; \\ Jinyuan Wang, May 09 2020
Extensions
a(17) from Jinyuan Wang, May 09 2020
a(18)-a(24) from Giovanni Resta, May 10 2020
Comments