|
%I #17 Apr 02 2023 08:17:35
%S 1,379,467,479,1325,7144,10311,26079
%N Numbers m such that 13^m reversed is prime.
%t Do[ If[ PrimeQ[ ToExpression[ StringReverse[ ToString[13^n] ] ] ], Print[n] ], {n, 1, 2300} ]
%t Select[Range[1400],PrimeQ[IntegerReverse[13^#]]&] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Jul 18 2020 *)
%o (Python)
%o from sympy import isprime
%o k, m, A058995_list = 1, 13, []
%o while k <= 10**3:
%o if isprime(int(str(m)[::-1])):
%o A058995_list.append(k)
%o k += 1
%o m *= 13 # _Chai Wah Wu_, Mar 09 2021
%o (PARI) isok(m) = isprime(fromdigits(Vecrev(digits(13^m)))) \\ _Mohammed Yaseen_, Jul 21 2022
%Y Numbers m such that k^m reversed is prime: A057708 (k=2), A350441 (k=4), A058993 (k=5), A058994 (k=7), A350442 (k=8), this sequence (k=13).
%K nonn,base,more
%O 1,2
%A _Robert G. Wilson v_, Jan 17 2001
%E a(6)-a(7) from _Chai Wah Wu_, Mar 09 2021
%E a(8) from _Michael S. Branicky_, Apr 01 2023
| 