[go: up one dir, main page]

login
n*R(n)-1 is prime, where R(n) is reverse of n.
1

%I #8 Jul 05 2020 17:27:29

%S 2,12,21,30,36,60,63,90,114,132,150,162,174,192,198,204,231,237,240,

%T 246,255,261,264,291,306,330,360,378,390,402,411,420,438,447,456,462,

%U 471,477,495,510,552,588,594,603,609,627,630,642,654,669,690,726,732

%N n*R(n)-1 is prime, where R(n) is reverse of n.

%C From _Robert Israel_, Jul 25 2019: (Start)

%C If n is in the sequence and is not divisible by 10, then A004086(n) is in the sequence.

%C All terms except 2 are divisible by 3. (End)

%H Robert Israel, <a href="/A073041/b073041.txt">Table of n, a(n) for n = 1..10000</a>

%e 12 is a term because 12*21-1=251 is prime.

%p revdigs:= proc(n) local L,i;

%p L:= convert(n,base,10);

%p add(L[-i]*10^(i-1),i=1..nops(L))

%p end proc:

%p select(t -> isprime(t*revdigs(t)-1), [$1..1000]); # _Robert Israel_, Jul 25 2019

%t Select[Range[800],PrimeQ[# IntegerReverse[#]-1]&] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Jul 05 2020 *)

%Y Cf. A004086, A061205.

%K easy,base,nonn

%O 1,1

%A _Shyam Sunder Gupta_, Aug 23 2002