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-5 of 5 results.

A156874 Number of Sophie Germain primes <= n.

Original entry on oeis.org

0, 1, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 18 2009

Keywords

Comments

a(n) = Sum_{k=1..n} A156660(k).
a(n) = A156875(2*n+1).
Hardy-Littlewood conjecture: a(n) ~ 2*C2*n/(log(n))^2, where C2=0.6601618158... is the twin prime constant (see A005597).
The truth of the above conjecture would imply that there exists an infinity of Sophie Germain primes (which is also conjectured).
a(n) ~ 2*C2*n/(log(n))^2 is also conjectured by Hardy-Littlewood for the number of twin primes <= n.

Examples

			a(120) = #{2, 3, 5, 11, 23, 29, 41, 53, 83, 89, 113} = 11.
		

Crossrefs

Cf. A005384 Sophie Germain primes p: 2p+1 is also prime.
Cf. A092816.

Programs

  • Mathematica
    Accumulate[Table[Boole[PrimeQ[n]&&PrimeQ[2n+1]],{n,1,200}]] (* Enrique Pérez Herrero, Apr 26 2012 *)
    Accumulate[Table[If[AllTrue[{n,2n+1},PrimeQ],1,0],{n,200}]]

Formula

a(10^n)= A092816(n). - Enrique Pérez Herrero, Apr 26 2012

Extensions

Edited and commented by Daniel Forgues, Jul 31 2009

A156875 Number of safe primes <= n.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 18 2009

Keywords

Comments

a(n) = SUM(A156659(k): 1<=k<=n);
a(n) = A156874(floor((n-1)\2)).

Examples

			a(120) = #{5, 7, 11, 23, 47, 59, 83, 107} = 8.
		

Crossrefs

A156877 Number of primes <= n that are safe primes and also Sophie Germain primes.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 18 2009

Keywords

Examples

			a(120) = #{5, 11, 23, 83} = 4.
		

Crossrefs

Formula

a(n) = A156874(n) + A156875(n) - A156876(n).
a(n) = Sum_{k=1..n} A156659(k)*A156660(k).

A156878 Number of primes <= n that are neither safe primes and nor Sophie Germain primes.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 18 2009

Keywords

Examples

			a(120) = #{13,17,19,31,37,43,61,67,71,73,79,97,10,103,109} = 15.
		

Crossrefs

Formula

a(n) = A000720(n) - A156876(n).

A156658 Primes p such that also 2*p+1 or (p-1)/2 is prime.

Original entry on oeis.org

2, 3, 5, 7, 11, 23, 29, 41, 47, 53, 59, 83, 89, 107, 113, 131, 167, 173, 179, 191, 227, 233, 239, 251, 263, 281, 293, 347, 359, 383, 419, 431, 443, 467, 479, 491, 503, 509, 563, 587, 593, 641, 653, 659, 683, 719, 743, 761, 809, 839, 863, 887, 911, 953, 983, 1013
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 13 2009

Keywords

Comments

Union of A005384 and A005385;
The intersection of A005384 and A005385 is given by A059455.
A156660(a(n)) + A156659(a(n)) > 0;
primes occurring in Cunningham chains of the first kind.
A156876 gives the number of these numbers <= n. [Reinhard Zumkeller, Feb 18 2009]

Crossrefs

Programs

  • Maple
    select(t -> isprime(t) and (isprime(2*t+1) or isprime((t-1)/2)), [2,seq(p,p=3..10000,2)]); # Robert Israel, May 03 2016
  • Mathematica
    Select[Prime@ Range@ 180, PrimeQ[2 # + 1] || PrimeQ[(# - 1)/2] &] (* Michael De Vlieger, Apr 06 2016 *)
  • PARI
    lista(nn) = {forprime(p=2, nn, if (isprime(2*p+1) || isprime((p-1)/2), print1(p, ", ")););} \\ Michel Marcus, Apr 06 2016
Showing 1-5 of 5 results.