OFFSET
1,1
COMMENTS
Subsequence of A068407. - Giovanni Resta, Jul 25 2013
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..100
EXAMPLE
193 is a prime and 193^5=267785184193 ends in 193, hence 193 is in the sequence.
MAPLE
with(numtheory):with(StringTools):KD := proc() local a, b, d, e, f; a:= ithprime(n); b:= a^5; d:=length(a); e:=floor(b/(10^d))*10^d; f:=b-e; if a=f then RETURN (a) fi:end:seq(KD(), n=1..500000);
MATHEMATICA
d[n_] := Block[{x}, Select[x /. List@ ToRules@ Reduce[x^5 == x, {x}, Modulus -> 10^n], # > 10^(n-1) && PrimeQ@# &]]; Union @@ d /@ Range@ 9 (* Giovanni Resta, Jul 25 2013 *)
PROG
(PARI) is(n)=isprime(n) && Mod(n, 10^#digits(n))^5==n \\ Charles R Greathouse IV, Jul 26 2013
CROSSREFS
KEYWORD
nonn,base
AUTHOR
K. D. Bajpai, Jul 25 2013
EXTENSIONS
a(35)-a(39) from Giovanni Resta, Jul 25 2013
STATUS
approved