A145480 - OEIS

%I #10 Jan 23 2017 13:06:29

%S 37,97,109,157,181,241,277,349,409,421,577,661,709,757,829,877,937,

%T 1009,1117,1201,1249,1381,1429,1609,1621,1669,1777,1801,2029,2089,

%U 2137,2221,2269,2389,2437,2521,2557,2617,2857,3049,3061,3121,3181,3217,3301,3361

%N Primes p such that (p+37)/2 is prime.

%C All these primes are congruent to 1 mod 12

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

%F a(n) >> n log^2 n. - _Charles R Greathouse IV_, Jan 23 2017

%t aa = {}; k = 37; Do[If[PrimeQ[(k + Prime[n])/2], AppendTo[aa, Prime[n]]], {n, 1, 500}];aa

%t Select[Prime[Range[500]],PrimeQ[(37+#)/2]&] (* _Harvey P. Dale_, Jun 23 2016 *)

%o (PARI) list(n)=my(t=1, p, i=1); while(i<n, p=prime(i); i=i+1; if(p>2&&isprime((37+p)/2), print(p,", "))) \\_Anders Hellström_, Jan 23 2017

%o (PARI) list(lim)=my(v=List()); forprime(p=3,lim, if(isprime((p+37)/2), listput(v,p))); Vec(v) \\ _Charles R Greathouse IV_, Jan 23 2017

%Y A092109, A145471-A145480.

%K nonn,easy

%O 1,1

%A _Artur Jasinski_, Oct 11 2008