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.

A167706 The single or isolated numbers. The union of single (or isolated or non-twin) primes and single (or isolated or average of twin prime pairs) nonprimes.

Original entry on oeis.org

2, 4, 6, 12, 18, 23, 30, 37, 42, 47, 53, 60, 67, 72, 79, 83, 89, 97, 102, 108, 113, 127, 131, 138, 150, 157, 163, 167, 173, 180, 192, 198, 211, 223, 228, 233, 240, 251, 257, 263, 270, 277, 282, 293, 307, 312, 317, 331, 337, 348, 353, 359, 367, 373, 379, 383, 389
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Nov 10 2009, Nov 14 2009

Keywords

Comments

Equals A007510 U A014574.

Crossrefs

Programs

  • Mathematica
    With[{nn = 78}, {2}~Join~Union[Transpose[Select[Partition[Prime@ Range@ nn, 3, 1], And[#[[2]] - #[[1]] != 2, #[[3]] - #[[2]] != 2] &]][[2]], Map[Mean, Select[Partition[Prime@ Range@ nn, 2, 1], Differences@ # == {2} &]]]] (* Michael De Vlieger, Feb 22 2017, after Harvey P. Dale at A007510 and A014574 *)
  • PARI
    is(n)=if(n%6, (isprime(n) && !isprime(n-2) && !isprime(n+2)) || n==4, isprime(n-1) && isprime(n+1)) \\ Charles R Greathouse IV, Apr 29 2015
    
  • PARI
    lista(pmax) = {my(p1 = 2, p2 = 3); print1(2, ", "); forprime(p3 = 5, pmax, if(p2 == p1 + 2, print1(p1 + 1, ", ")); if(p2 != p1 + 2 && p2 != p3 - 2, print1(p2, ", ")); p1 = p2; p2 = p3);} \\ Amiram Eldar, May 17 2024

Formula

a(n) ~ n log n. - Charles R Greathouse IV, Apr 29 2015

Extensions

Corrected (97 inserted) by R. J. Mathar, Nov 16 2009