[go: up one dir, main page]

login
Numbers m such that 2^m reversed is prime.
14

%I #51 Jul 20 2022 15:19:47

%S 1,4,5,7,10,17,24,37,45,55,70,77,107,137,150,271,364,1157,1656,2004,

%T 2126,3033,3489,3645,4336,6597,7279,12690,13840,20108,21693,28888,

%U 84155,102930

%N Numbers m such that 2^m reversed is prime.

%C a(35) > 105000. - _Giovanni Resta_, Feb 22 2013

%C From _Bernard Schott_, Jan 30 2022: (Start)

%C 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 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)

%e 4 is a term because 2^4 reversed is 61 and prime.

%p 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 Do[ If[ PrimeQ[ FromDigits[ Reverse[ IntegerDigits[2^n]] ]], Print[ n]], {n, 20000}] (* _Robert G. Wilson v_, Jan 29 2005 *)

%t 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 (Python)

%o from sympy import isprime

%o k, m, A057708_list = 1, 2, []

%o while k <= 10**3:

%o if isprime(int(str(m)[::-1])):

%o A057708_list.append(k)

%o k += 1

%o m *= 2 # _Chai Wah Wu_, Mar 09 2021

%o (PARI) isok(m) = isprime(fromdigits(Vecrev(digits(2^m)))) \\ _Mohammed Yaseen_, Jul 20 2022

%Y Cf. A004094, A341713.

%Y 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 base,nonn,more

%O 1,2

%A _G. L. Honaker, Jr._, Oct 23 2000

%E More terms from Chris Nash (chris_nash(AT)prodigy.net), Oct 25 2000

%E Two more terms from _Robert G. Wilson v_, Jan 29 2001

%E 3 more terms from _Farideh Firoozbakht_, Aug 05 2004

%E a(33)-a(34) from _Giovanni Resta_, Feb 22 2013