cp's OEIS Frontend

A059002 Primes whose reversal is a seventh power.

%I A059002 #10 Dec 22 2024 20:04:44
%S A059002 821,90367894271,188372457491,1938510215909,67717549154641,
%T A059002 95168980930291,772155318141637,1558489460499871,7505006230374799,
%U A059002 8237815094854781,23021614989689299,39712513595115047,919042243755035867,3573813286514042801,7597919404531928707
%N A059002 Primes whose reversal is a seventh power.
%H A059002 Robert Israel, <a href="/A059002/b059002.txt">Table of n, a(n) for n = 1..10000</a>
%p A059002 rev:= proc(n) local L,i;
%p A059002   L:= convert(n,base,10);
%p A059002   add(L[-i]*10^(i-1),i=1..nops(L))
%p A059002 end proc:
%p A059002 S:= {}: count:= 0:
%p A059002 for d from 1 while count < 60 do
%p A059002  for i from 10^(d-1) to 10^d do
%p A059002   if i mod 10 = 0 then next fi;
%p A059002   p:= rev(i^7);
%p A059002   if isprime(p) then S:= S union {p}; count:= count+1; fi;
%p A059002 od od:
%p A059002 sort(convert(S,list)); # _Robert Israel_, Dec 22 2024
%t A059002 Do[ If[ PrimeQ[ n ] && IntegerQ[ ToExpression[ StringReverse[ ToString[ n ] ] ]^(1/7) ], Print[ n ] ], {n, 1, 10^19} ]
%Y A059002 Cf. A007488.
%K A059002 nonn,base
%O A059002 1,1
%A A059002 _Robert G. Wilson v_, Jan 16 2001
%E A059002 More terms from _Sean A. Irvine_, Sep 09 2022