cp's OEIS Frontend

A243409 Primes p such that 100p-1, 100p-3, 100p-7, and 100p-9 are all prime.

%I A243409 #17 Feb 06 2019 18:10:05
%S A243409 2,797,1193,6803,15773,28793,35507,41579,53189,53279,57347,60161,
%T A243409 70457,77549,81839,140549,143387,150779,151241,164447,170627,201011,
%U A243409 255083,285287,293831,300317,316073,336671,343661,449921,470087,486947,488603,518801,556289,569243,602087
%N A243409 Primes p such that 100p-1, 100p-3, 100p-7, and 100p-9 are all prime.
%H A243409 K. D. Bajpai, <a href="/A243409/b243409.txt">Table of n, a(n) for n = 1..5500</a>
%e A243409 2 is prime, 100*2-1 = 199 is prime, 100*2-3 = 197 is prime, 100*2-7 = 193 is prime, and 100*2-9 = 191 is prime. Thus 2 is a member of this sequence.
%t A243409 Select[Prime[Range[50000]],PrimeQ[100# -1]&&PrimeQ[100# -3]&&PrimeQ[100# -7] &&PrimeQ[100# -9] &] (* _K. D. Bajpai_, Jun 13 2014 *)
%t A243409 Select[Prime[Range[50000]],AllTrue[100#-{1,3,7,9},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Feb 06 2019 *)
%o A243409 (Python)
%o A243409 import sympy
%o A243409 from sympy import isprime
%o A243409 from sympy import prime
%o A243409 {print(prime(n),end=', ') for n in range(1,10**5) if isprime(100*prime(n)-1) and isprime(100*prime(n)-3) and isprime(100*prime(n)-7) and isprime(100*prime(n)-9)}
%o A243409 (PARI) for(n=1,10^5,if(ispseudoprime(100*prime(n)-1)&& ispseudoprime(100*prime(n)-3)&& ispseudoprime(100*prime(n)-7)&& ispseudoprime(100*prime(n)-9),print1(prime(n),", ")))
%Y A243409 Cf. A236042, A064976.
%K A243409 nonn
%O A243409 1,1
%A A243409 _Derek Orr_, Jun 04 2014