|
%I #19 Jun 30 2018 02:27:02
%S 5,13,31,157,239,683,1697,3067,4273,5563,6323,6781,9043,13841,16427,
%T 22247,22901,23767,27109,32687,33721,45833,58217,61933,83341,90583,
%U 94811,115979,158161,172751,193619,202877,221567,284233,297967,308773
%N Primes of the form 2*k*prime(k) + 1.
%H Harry J. Smith, <a href="/A062403/b062403.txt">Table of n, a(n) for n = 1..1000</a>
%p select(isprime,[seq(2*i*ithprime(i)+1,i=1..200)]); # _Muniru A Asiru_, Jun 29 2018
%o (PARI) je=[]; for(n=0,200, if(isprime(2*n*prime(n)+1), je=concat(je,2*n*prime(n)+1),)); je
%o (PARI) { n=c=0; forprime (p=2, 5*10^5, q=2*(c++)*p + 1; if (isprime(q), write("b062403.txt", n++, " ", q); if (n==1000, break)) ) } \\ _Harry J. Smith_, Aug 07 2009
%K easy,nonn
%O 1,1
%A _Jason Earls_, Jul 08 2001
%E Definition corrected by _Harry J. Smith_, Aug 07 2009
%E Name simplified by _Jon E. Schoenfield_, Jun 29 2018
|
