cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A057708 Numbers m such that 2^m reversed is prime.

This page as a plain text file.
%I A057708 #55 May 12 2025 10:13:10
%S A057708 1,4,5,7,10,17,24,37,45,55,70,77,107,137,150,271,364,1157,1656,2004,
%T A057708 2126,3033,3489,3645,4336,6597,7279,12690,13840,20108,21693,28888,
%U A057708 84155,102930
%N A057708 Numbers m such that 2^m reversed is prime.
%C A057708 a(35) > 105000. - _Giovanni Resta_, Feb 22 2013
%C A057708 From _Bernard Schott_, Jan 30 2022: (Start)
%C A057708 If m is an even term, then u = m/2 is a term of A350441, this because 2^m = 4^(m/2). In fact, terms of A350441 are half the even terms of this sequence here.
%C A057708 If m is a term multiple of 3, then k = m/3 is a term of A350442, this because 2^m = 8^(m/3). First examples: m = 24, 45, 150, 1656, ... and corresponding k = 8, 15, 50, 552, ... (End)
%C A057708 a(35) > 200000. - _Michael S. Branicky_, May 12 2025
%e A057708 4 is a term because 2^4 reversed is 61 and prime.
%p A057708 with(numtheory): myarray := []: for n from 1 to 4000 do it1 := convert(2^n, base, 10): it2 := sum(10^(nops(it1)-i)*it1[i], i=1..nops(it1)): if isprime(it2) then printf(`%d,`,n) fi: od:
%t A057708 Do[ If[ PrimeQ[ FromDigits[ Reverse[ IntegerDigits[2^n]] ]], Print[ n]], {n, 20000}] (* _Robert G. Wilson v_, Jan 29 2005 *)
%t A057708 Select[Range[4400],PrimeQ[IntegerReverse[2^#]]&] (* Requires Mathematica version 10 or later *) (* The program generates the first 25 terms of the sequence; to generate more, increase the Range constant, but the program will take longer to run. *) (* _Harvey P. Dale_, Aug 05 2020 *)
%o A057708 (Python)
%o A057708 from sympy import isprime
%o A057708 k, m, A057708_list = 1, 2,  []
%o A057708 while k <= 10**3:
%o A057708     if isprime(int(str(m)[::-1])):
%o A057708         A057708_list.append(k)
%o A057708     k += 1
%o A057708     m *= 2 # _Chai Wah Wu_, Mar 09 2021
%o A057708 (PARI) isok(m) = isprime(fromdigits(Vecrev(digits(2^m)))) \\ _Mohammed Yaseen_, Jul 20 2022
%Y A057708 Cf. A004094, A341713.
%Y A057708 Numbers m such that k^m reversed is prime: this sequence (k=2), A350441 (k=4), A058993 (k=5), A058994 (k=7), A350442 (k=8), A058995 (k=13).
%K A057708 base,nonn,more
%O A057708 1,2
%A A057708 _G. L. Honaker, Jr._, Oct 23 2000
%E A057708 More terms from Chris Nash (chris_nash(AT)prodigy.net), Oct 25 2000
%E A057708 Two more terms from _Robert G. Wilson v_, Jan 29 2001
%E A057708 3 more terms from _Farideh Firoozbakht_, Aug 05 2004
%E A057708 a(33)-a(34) from _Giovanni Resta_, Feb 22 2013