cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A103807 Primes p such that 2*p-27, 2*p+27, 2*p-33 and 2*p+33 are primes or -1 times primes.

Original entry on oeis.org

2, 5, 7, 23, 37, 103, 313, 457, 733, 863, 2053, 2063, 2917, 4723, 7187, 7817, 8017, 9007, 9473, 9973, 10687, 11527, 11923, 13477, 13883, 15787, 26833, 31477, 34897, 36097, 36353, 36493, 39937, 44417, 46447, 47623, 52103, 53377, 55813, 60737
Offset: 1

Views

Author

Zak Seidov, Feb 16 2005

Keywords

Comments

Intersection of A103805 and A103806.

Crossrefs

Programs

  • Magma
    [ p: p in PrimesUpTo(61000) | IsPrime(2*p-27) and IsPrime(2*p+27) and IsPrime(2*p-33) and IsPrime(2*p+33) ];
    
  • Mathematica
    Intersection[Select[Range[100000], PrimeQ[ # ]&&PrimeQ[2#+33]&&PrimeQ[2#-33]&&PrimeQ[ # ]&&PrimeQ[2#+27]&&PrimeQ[2#-27]&]]
    okQ[n_]:=Module[{x=2n},And@@PrimeQ[{x-27,x+27,x-33,x+33}]]; Select[Prime[Range[7000]],okQ]  (* Harvey P. Dale, Jan 23 2011 *)
  • PARI
    {forprime(p=2, 61000, if(isprime(abs(2*p-27))&&isprime(2*p+27)&&isprime(abs(2*p-33))&&isprime(2*p+33), print1(p, ", ")))}

Extensions

Definition clarified, comment adjusted, MAGMA and PARI programs added by Klaus Brockhaus, Mar 21 2010