|
%I #16 Sep 09 2014 20:48:22
%S 83,167,227,503,677,1217,1487,4973,5843,10337,12437,13163,15683,15923,
%T 21803,22067,23057,24107,24113,36587,40883,47717,53633,57287,65183,
%U 73013,95633,104123,107777,109583,113147,122267,128987,129533,132173,148403,148727
%N Primes p such that p - m^2, m = 2, 4, 6, 8, are all (positive) primes.
%H Chai Wah Wu, <a href="/A246874/b246874.txt">Table of n, a(n) for n = 1..1000</a>
%e 83, 79, 67, 47 and 19 are primes so 83 is a term in this sequence.
%o (PARI) select(p->isprime(p-4) && isprime(p-16) && isprime(p-36) && isprime(p-64), primes(200000)) \\ _Colin Barker_, Sep 08 2014
%o (Python)
%o from sympy import prime, isprime
%o A246874_list = [p for p in (prime(n) for n in range(1,3*10**6)) if all([isprime(p-m*m) for m in range(2,10,2)])]
%o # _Chai Wah Wu_, Sep 08 2014
%Y Subsequence of A246873.
%K nonn
%O 1,1
%A _Zak Seidov_, Sep 06 2014
%E More terms from _Colin Barker_, Sep 08 2014
|
