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.

Showing 1-4 of 4 results.

A023200 Primes p such that p + 4 is also prime.

Original entry on oeis.org

3, 7, 13, 19, 37, 43, 67, 79, 97, 103, 109, 127, 163, 193, 223, 229, 277, 307, 313, 349, 379, 397, 439, 457, 463, 487, 499, 613, 643, 673, 739, 757, 769, 823, 853, 859, 877, 883, 907, 937, 967, 1009, 1087, 1093, 1213, 1279, 1297, 1303, 1423, 1429, 1447, 1483
Offset: 1

Views

Author

Keywords

Comments

Smaller member p of cousin prime pairs (p, p+4).
A015913 contains the composite number 305635357, so it is different from both the present sequence and A029710. (305635357 is the only composite member of A015913 < 10^9.) - Jud McCranie, Jan 07 2001
Apart from the first term, all terms are of the form 6n + 1.
Complement of A067775 (primes p such that p + 4 is composite) with respect to A000040 (primes). With prime 2 also primes p such that q^2 + p is prime for some prime q (q = 3 if p = 2, q = 2 if p > 2). Subsequence of A232012. - Jaroslav Krizek, Nov 23 2013
Conjecture: The sequence is infinite and for every n, a(n+1) < a(n)^(1+1/n). Namely a(n)^(1/n) is a strictly decreasing function of n. - Jahangeer Kholdi and Farideh Firoozbakht, Nov 24 2014
From Alonso del Arte, Jan 12 2019: (Start)
If p splits in Z[sqrt(-2)], p + 4 is an inert prime in that domain. Likewise, if p splits in Z[sqrt(2)], p + 4 is an inert prime in that domain.
The only way for p or p + 4 to split in both domains is if it is congruent to 1 modulo 24, in which case the other prime is inert in both domains.
For example, 3 = (1 - sqrt(-2))*(1 + sqrt(-2)) but is inert in Z[sqrt(2)], while 7 = (3 - sqrt(2))*(3 + sqrt(2)) but is inert in Z[sqrt(-2)]. And also 11 = (3 - sqrt(-2))*(3 + sqrt(-2)) but 15 is composite in Z or any quadratic integer ring.
And 97 = (5 - 6*sqrt(-2))*(5 + 6*sqrt(-2)) = (1 - 7*sqrt(2))*(1 + 7*sqrt(2)), but 101 is inert in both Z[sqrt(-2)] and Z[sqrt(2)]. (End)

Crossrefs

Exactly the same as A029710 except for the exclusion of 3.

Programs

  • Haskell
    a023200 n = a023200_list !! (n-1)
    a023200_list = filter ((== 1) . a010051') $
                   map (subtract 4) $ drop 2 a000040_list
    -- Reinhard Zumkeller, Aug 01 2014
  • Magma
    [p: p in PrimesUpTo(1500) | NextPrime(p)-p eq 4]; // Bruno Berselli, Apr 09 2013
    
  • Maple
    A023200 := proc(n) option remember; if n = 1 then 3; else p := nextprime(procname(n-1)) ; while not isprime(p+4) do p := nextprime(p) ;  end do: p ; end if; end proc: # R. J. Mathar, Sep 03 2011
  • Mathematica
    Select[Range[10^2], PrimeQ[#] && PrimeQ[# + 4] &] (* Vladimir Joseph Stephan Orlovsky, Apr 29 2008 *)
    Select[Prime[Range[250]],PrimeQ[#+4]&] (* Harvey P. Dale, Oct 09 2023 *)
  • PARI
    print1(3);p=7;forprime(q=11,1e3,if(q-p==4,print1(", "p)); p=q) \\ Charles R Greathouse IV, Mar 20 2013
    

Formula

a(n) = A046132(n) - 4 = A087679(n) - 2.
a(n) >> n log^2 n via the Selberg sieve. - Charles R Greathouse IV, Nov 20 2016

Extensions

Definition modified by Vincenzo Librandi, Aug 02 2009
Definition revised by N. J. A. Sloane, Mar 05 2010

A067775 Primes p such that p + 4 is composite.

Original entry on oeis.org

2, 5, 11, 17, 23, 29, 31, 41, 47, 53, 59, 61, 71, 73, 83, 89, 101, 107, 113, 131, 137, 139, 149, 151, 157, 167, 173, 179, 181, 191, 197, 199, 211, 227, 233, 239, 241, 251, 257, 263, 269, 271, 281, 283, 293, 311, 317, 331, 337, 347, 353, 359, 367, 373, 383, 389
Offset: 1

Views

Author

Benoit Cloitre, Feb 06 2002

Keywords

Comments

Primes n such that n!*B(n+3) is an integer where B(k) are the Bernoulli numbers B(1) = -1/2, B(2) = 1/6, B(4) = -1/30, ..., B(2m+1) = 0 for m > 1.
If n is prime n!*B(n-1) is always an integer. Note that if Goldbach's conjecture (2n = p1 + p2 for all n >= 2) is false and K is the smallest value of n for which it fails, then for 2(K-2) = p3 + p4, the primes p3 and p4 must be taken from this list. See similar comment for A140555. - Keith Backman, Apr 06 2012
Complement of A023200 (primes p such that p + 4 is also prime) with respect to A000040 (primes). For p > 2: primes p such that there is no prime of the form r^2 + p where r is prime, subsequence of A232010. Example: the prime 7 is not in the sequence because 2^2 + 7 = 11 (prime). A232009(a(n)) = 0 for n > 1 . - Jaroslav Krizek, Nov 22 2013

Crossrefs

Programs

  • Mathematica
    A067775 = {}; Do[p = Prime@ n; If[ IntegerQ[ p! BernoulliB[p + 3]], AppendTo[A067775, p]], {n, 77}]; A067775 (* Robert G. Wilson v, Aug 19 2008 *)
    Select[Prime[Range[80]], Not[PrimeQ[# + 4]] &] (* Alonso del Arte, Apr 02 2014 *)
  • PARI
    lista(nn) = {forprime(p=1, nn, if (! isprime(p+4), print1(p, ", ")););} \\ Michel Marcus, Nov 22 2013

Formula

a(n) ~ n log n. - Charles R Greathouse IV, Nov 22 2013

Extensions

New name from Klaus Brockhaus at the suggestion of Michel Marcus, Nov 22 2013

A232009 a(n) = the smallest squarefree number (from A005117) of the form p*q with prime factors in a p^2 + n progression, or 0 if no such number exists.

Original entry on oeis.org

10, 33, 14, 39, 0, 155, 22, 51, 26, 57, 0, 185, 34, 69, 38, 205, 0, 215, 46, 87, 0, 93, 0, 511, 58, 0, 62, 111, 0, 553, 0, 123, 74, 129, 0, 305, 82, 141, 86, 623, 0, 335, 94, 159, 0, 355, 0, 365, 106, 177, 0, 183, 0, 395, 118, 0, 122, 201, 0, 763, 0, 213, 134
Offset: 1

Views

Author

Jaroslav Krizek, Nov 20 2013

Keywords

Comments

a(n) = the smallest squarefree number m of the form p*q with prime factors p and q = p^2 + n, or 0 if no such number exists; m = p^3 + p*n.
a(n) = 0 for numbers n from A232010.

Crossrefs

Formula

a(4) = 39 because 39 = 3 * 13 = 3 * (3^2 + 4).

A232010 Numbers k such that p^2 + k is composite for all prime p.

Original entry on oeis.org

5, 11, 17, 21, 23, 26, 29, 31, 35, 41, 45, 47, 51, 53, 56, 59, 61, 65, 68, 71, 73, 77, 81, 83, 86, 87, 89, 91, 95, 101, 107, 110, 111, 113, 115, 116, 117, 119, 121, 125, 129, 131, 134, 137, 139, 141, 143, 146, 149, 151, 152, 155, 157, 161, 165, 167, 171, 173, 176
Offset: 1

Views

Author

Jaroslav Krizek, Nov 20 2013

Keywords

Comments

Primes from this sequence are in A067775.
Complement of A232012.

Examples

			7 is not in the sequence because 2^2 + 7 = 11 (prime).
		

Crossrefs

Formula

A232009(a(n)) = 0.
Showing 1-4 of 4 results.