login
Primes prime(n) such that 10*prime(n)+prime(n+2) is prime.
1

%I #6 Apr 17 2017 14:37:43

%S 3,5,7,11,13,17,41,67,89,101,109,181,193,199,223,227,283,313,337,347,

%T 367,419,421,449,461,479,491,587,641,647,709,727,773,811,823,827,857,

%U 859,863,881,887,911,953,1021,1049,1063,1091,1181,1217,1361,1399,1447

%N Primes prime(n) such that 10*prime(n)+prime(n+2) is prime.

%H Harvey P. Dale, <a href="/A130122/b130122.txt">Table of n, a(n) for n = 1..1000</a>

%e a(7)=41 because 10*41+47=457 and 457 is prime,

%e a(8)=67 because 10*67+73=743 and 743 is prime,

%e a(9)=89 because 10*89+101=991 and 991 is prime, ...

%p ts_p6_01:=proc(n) local a,b,i,ans; ans := [ ]: for i from 1 to n do a := 10*ithprime(i)+ithprime(i+2): if (isprime(a)=true) then ans := [ op(ans), ithprime(i) ]: fi od; RETURN(ans) end: ts_p6_01(1000);

%t Select[Partition[Prime[Range[300]],3,1],PrimeQ[10#[[1]]+#[[3]]]&][[All,1]] (* _Harvey P. Dale_, Apr 17 2017 *)

%Y Cf. A129900.

%K nonn

%O 1,1

%A _Jani Melik_, Aug 01 2007