%I #20 Aug 23 2023 08:37:19
%S 13,31,59,191,887,1019,1931,2903,5471,8087,9587,19031,23099,33119,
%T 57587,80651,129587,168083,188351,327179,359987,414731,678971,846383,
%U 898691,910103,984047,1040387,1044479,1132091,1331711,1411331,1444787,1517819,1669259,1909907
%N Primes of the form prime(k)^2 + 2*prime(k-1) where prime(k) is the k-th prime number.
%H Robert Israel, <a href="/A155820/b155820.txt">Table of n, a(n) for n = 1..10000</a>
%e prime(4)=7, prime(3)=5; 7^2+2*5=59, a prime. Hence 59 is a term.
%p count:= 0: q:= 2: R:= NULL:
%p while count < 100 do
%p p:= q; q:= nextprime(q);
%p v:= q^2 + 2*p;
%p if isprime(v) then count:= count+1; R:= R,v; fi;
%p od:
%p R; # _Robert Israel_, Aug 22 2023
%t list = {}; Do[m = Prime[k]^2 + 2*Prime[k - 1]; If[PrimeQ[m], AppendTo[list, m]], {k, 2, 300}]; list (* _Vaclav Kotesovec_, Feb 14 2019 *)
%Y Cf. A000040.
%K nonn
%O 1,1
%A Avik Roy (avik_3.1416(AT)yahoo.co.in), Jan 28 2009
%E More terms from _Vaclav Kotesovec_, Feb 14 2019
%E Changed offset to 1 by _Vaclav Kotesovec_, Feb 14 2019
|