|
%I #9 Jan 23 2017 15:05:42
%S 17,19,31,37,67,79,97,127,151,157,181,199,277,331,337,379,409,421,457,
%T 499,541,547,577,601,631,661,727,739,751,757,787,829,877,907,991,1009,
%U 1021,1087,1117,1171,1201,1249,1291,1381,1399,1459,1549,1597,1609,1669
%N Primes p such that 2p - 31 is prime.
%H Ivan Neretin, <a href="/A145485/b145485.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = 2*A145479(n) - 31.
%t aa = {}; k = 31; Do[If[PrimeQ[(k + Prime[n])/2], AppendTo[aa, k + Prime[n])/2]], {n, 1, 500}];aa
%o (PARI) list(lim)=my(v=List()); forprime(p=2,lim, if(isprime(2*p-31), listput(v,p))); Vec(v) \\ _Charles R Greathouse IV_, Jan 23 2017
%Y Cf. A063908-A063913, A092109, A145471-A145490.
%K easy,nonn
%O 1,1
%A _Artur Jasinski_, Oct 11 2008
|
