A053176 -id:A053176 - cp's OEIS frontend

A059500 Primes p such that both q=(p-1)/2 and 2p + 1 = 4q + 3 are composite numbers. Intersection of A059456 and A053176.

13, 17, 19, 31, 37, 43, 61, 67, 71, 73, 79, 97, 101, 103, 109, 127, 137, 139, 149, 151, 157, 163, 181, 193, 197, 199, 211, 223, 229, 241, 257, 269, 271, 277, 283, 307, 311, 313, 317, 331, 337, 349, 353, 367, 373, 379, 389, 397, 401, 409, 421, 433, 439, 449

Offset: 1

Labos Elemer, Feb 05 2001

nonn

Primes which are neither safe nor of Sophie Germain type.
Primes not in Cunningham chains of the first kind. - Alonso del Arte, Jun 30 2005
A010051(a(n))*(1-A156660(a(n)))*(1-A156659(a(n))) = 1; A156878 gives numbers of these numbers <= n. - Reinhard Zumkeller, Feb 18 2009

			Prime p=17 is here because both 35 and 8 are composite numbers. Such primes fall "out of" any Cunningham chain of first kind (or generate Cunningham chains of 0-length).

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from _Reinhard Zumkeller)
C. K. Caldwell, Cunningham Chains

Cf. A005384, A005385, A053176, A059452-A059456, A007700, A005602, A023272, A023302, A023330, A156658.

Complement[Prime[Range[100]], Select[Prime[Range[100]], PrimeQ[2# + 1] &], Select[Prime[Range[100]], PrimeQ[(# - 1)/2] &]] (Delarte)
Select[Prime[Range[100]],!PrimeQ[q=2#+1]&&!PrimeQ[(#-1)/2]&] (* Zak Seidov, Mar 09 2013 *)

is(n)=isprime(n)&&!isprime(n\2)&&!isprime(2*n+1) \\ Charles R Greathouse IV, Jan 16 2013

a(n) ~ n log n. - Charles R Greathouse IV, Jan 16 2013

A156543 Multiplicative closure of primes that are not Sophie Germain primes (A053176).

Original entry on oeis.org

1, 7, 13, 17, 19, 31, 37, 43, 47, 49, 59, 61, 67, 71, 73, 79, 91, 97, 101, 103, 107, 109, 119, 127, 133, 137, 139, 149, 151, 157, 163, 167, 169, 181, 193, 197, 199, 211, 217, 221, 223, 227, 229, 241, 247, 257, 259, 263, 269, 271, 277, 283, 289, 301, 307, 311

Offset: 1

Reinhard Zumkeller, Feb 09 2009

nonn

a(n) mod 6 = 1 or = 5; subsequence of A007310;
A156542(a(n)) = 0;
A045979 is a subsequence.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000

Cf. A007310, A045979, A053176, A156541, A156542.

A079699 Duplicate of A053176.

Original entry on oeis.org

7, 13, 17, 19, 31, 37, 43, 47, 59, 61, 67, 71, 73, 79, 97, 101, 103, 107, 109, 127, 137

Offset: 1

dead

A010051 Characteristic function of primes: 1 if n is prime, else 0.

Original entry on oeis.org

0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0

Offset: 1

N. J. A. Sloane

The following sequences all have the same parity (with an extra zero term at the start of a(n)): a(n), A061007, A035026, A069754, A071574. - Jeremy Gardiner, Aug 09 2002
Hardy and Wright prove that the real number 0.011010100010... is irrational. See Nasehpour link. - Michel Marcus, Jun 21 2018
The spectral components (excluding the zero frequency) of the Fourier transform of the partial sequences {a(j)} with j=1..n and n an even number, exhibit a remarkable symmetry with respect to the central frequency component at position 1 + n/4. See the Fourier spectrum of the first 2^20 terms in Links, Comments in A289777, and Conjectures in A001223 of Sep 01 2019. It also appears that the symmetry grows with n. - Andres Cicuttin, Aug 23 2020

J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 3.
V. Brun, Über das Goldbachsche Gesetz und die Anzahl der Primzahlpaare, Arch. Mat. Natur. B, 34, no. 8, 1915.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Oxford University Press, London, 1975.
Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 65.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 132.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from N. J. A. Sloane)
Andres Cicuttin, Fourier spectrum of the first 2^20 terms of the characteristic function of primes.
William Craig, Jan-Willem van Ittersum and Ken Ono, Integer partitions detect the primes, PNAS, Vol. 121, No. 39 (2024), e2409417121.
Yoichi Motohashi, An overview of the Sieve Method and its History, arXiv:math/0505521 [math.NT], 2005-2006.
Peyman Nasehpour, A Computational Criterion for the Irrationality of Some Real Numbers, Journal of Algorithms and Computation, Vol. 52, No. 1 (2020), pp. 97-104, preprint, arXiv:1806.07560 [math.AC], 2018.
José L. Ramírez and Gustavo N. Rubiano, Properties and Generalizations of the Fibonacci Word Fractal, The Mathematica Journal, Vol. 16 (2014).
Eric Weisstein's World of Mathematics, Prime Number.
Eric Weisstein's World of Mathematics, Prime Constant.
Eric Weisstein's World of Mathematics, Prime zeta function primezeta(s).
Index entries for characteristic functions

Cf. A051006 (constant 0.4146825... (base 10) = 0.01101010001010001010... (base 2)), A001221 (inverse Moebius transform), A143519, A156660, A156659, A156657, A059500, A053176, A059456, A072762.
First differences of A000720, so A000720 gives partial sums.
Column k=1 of A117278.
Characteristic function of A000040.
Cf. A008683.

import Data.List (unfoldr)
a010051 :: Integer -> Int
a010051 n = a010051_list !! (fromInteger n-1)
a010051_list = unfoldr ch (1, a000040_list) where
   ch (i, ps'@(p:ps)) = Just (fromEnum (i == p),
                              (i + 1, if i == p then ps else ps'))
-- Reinhard Zumkeller, Apr 17 2012, Sep 15 2011

s:=[]; for n in [1..100] do if IsPrime(n) then s:=Append(s,1); else s:=Append(s,0); end if; end for; s;

[IsPrime(n) select 1 else 0: n in [1..100]];  // Bruno Berselli, Mar 02 2011

A010051:= n -> if isprime(n) then 1 else 0 fi;

Table[ If[ PrimeQ[n], 1, 0], {n, 105}] (* Robert G. Wilson v, Jan 15 2005 *)
Table[Boole[PrimeQ[n]], {n, 105}] (* Alonso del Arte, Aug 09 2011 *)
Table[PrimePi[n] - PrimePi[n-1], {n,50}] (* G. C. Greubel, Jan 05 2017 *)

a(n)=isprime(n) \\ Charles R Greathouse IV, Apr 16 2011

from sympy import isprime
def A010051(n): return int(isprime(n)) # Chai Wah Wu, Jan 20 2022

a(n) = floor(cos(Pi*((n-1)! + 1)/n)^2) for n >= 2. - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 07 2002
Let M(n) be the n X n matrix m(i, j) = 0 if n divides ij + 1, m(i, j) = 1 otherwise; then for n > 0 a(n) = -det(M(n)). - Benoit Cloitre, Jan 17 2003
n >= 2, a(n) = floor(phi(n)/(n - 1)) = floor(A000010(n)/(n - 1)). - Benoit Cloitre, Apr 11 2003
a(n) = Sum_{d|gcd(n, A034386(n))} mu(d). [Brun]
a(m*n) = a(m)*0^(n - 1) + a(n)*0^(m - 1). - Reinhard Zumkeller, Nov 25 2004
a(n) = 1 if n has no divisors other than 1 and n, and 0 otherwise. - Jon Perry, Jul 02 2005
Dirichlet generating function: Sum_{n >= 1} a(n)/n^s = primezeta(s), where primezeta is the prime zeta function. - Franklin T. Adams-Watters, Sep 11 2005
a(n) = (n-1)!^2 mod n. - Franz Vrabec, Jun 24 2006
a(n) = A047886(n, 1). - Reinhard Zumkeller, Apr 15 2008
Equals A051731 (the inverse Möbius transform) * A143519. - Gary W. Adamson, Aug 22 2008
a(n) = A051731((n + 1)! + 1, n) from Wilson's theorem: n is prime if and only if (n + 1)! is congruent to -1 mod n. - N-E. Fahssi, Jan 20 2009, Jan 29 2009
a(n) = A166260/A001477. - Mats Granvik, Oct 10 2009
a(n) = 0^A070824, where 0^0=1. - Mats Granvik, Gary W. Adamson, Feb 21 2010
It appears that a(n) = (H(n)*H(n + 1)) mod n, where H(n) = n!*Sum_{k=1..n} 1/k = A000254(n). - Gary Detlefs, Sep 12 2010
Dirichlet generating function: log( Sum_{n >= 1} 1/(A112624(n)*n^s) ). - Mats Granvik, Apr 13 2011
a(n) = A100995(n) - sqrt(A100995(n)*A193056(n)). - Mats Granvik, Jul 15 2011
a(n) * (2 - n mod 4) = A151763(n). - Reinhard Zumkeller, Oct 06 2011
(n - 1)*a(n) = ( (2*n + 1)!! * Sum_{k=1..n}(1/(2*k + 1))) mod n, n > 2. - Gary Detlefs, Oct 07 2011
For n > 1, a(n) = floor(1/A001222(n)). - Enrique Pérez Herrero, Feb 23 2012
a(n) = mu(n) * Sum_{d|n} mu(d)*omega(d), where mu is A008683 and omega A001222 or A001221 indistinctly. - Enrique Pérez Herrero, Jun 06 2012
a(n) = A003418(n+1)/A003418(n) - A217863(n+1)/A217863(n) = A014963(n) - A072211(n). - Eric Desbiaux, Nov 25 2012
For n > 1, a(n) = floor(A014963(n)/n). - Eric Desbiaux, Jan 08 2013
a(n) = ((abs(n-2))! mod n) mod 2. - Timothy Hopper, May 25 2015
a(n) = abs(F(n)) - abs(F(n)-1/2) - abs(F(n)-1) + abs(f(n)-3/2), where F(n) = Sum_{m=2..n+1} (abs(1 - (n mod m)) - abs(1/2 - (n mod m)) + 1/2), n > 0. F(n) = 1 if n is prime, > 1 otherwise, except F(1) = 0. a(n) = 1 if F(n) = 1, 0 otherwise. - Timothy Hopper, Jun 16 2015
For n > 4, a(n) = (n-2)! mod n. - Thomas Ordowski, Jul 24 2016
From Ilya Gutkovskiy, Jul 24 2016: (Start)
G.f.: A(x) = Sum_{n>=1} x^A000040(n) = B(x)*(1 - x), where B(x) is the g.f. for A000720.
a(n) = floor(2/A000005(n)), for n>1. (End)
a(n) = pi(n) - pi(n-1) = A000720(n) - A000720(n-1), for n>=1. - G. C. Greubel, Jan 05 2017
Decimal expansion of Sum_{k>=1} (1/10)^prime(k) = 9 * Sum_{k>=1} pi(k)/10^(k+1), where pi(k) = A000720(k). - Amiram Eldar, Aug 11 2020
a(n) = 1 - ceiling((2/n) * Sum_{k=2..floor(sqrt(n))} floor(n/k)-floor((n-1)/k)), n>1. - Gary Detlefs, Sep 08 2023
a(n) = Sum_{d|n} mu(d)*omega(n/d), where mu = A008683 and omega = A001221. - Ridouane Oudra, Apr 12 2025
a(n) = 0 if (n^2 - 3*n + 2) * A000203(n) - 8 * A002127(n) > 0 else 1 (n>2, see Craig link). - Bill McEachen, Jul 04 2025

A045979 Bernoulli number B_{2n} has denominator 6.

Original entry on oeis.org

1, 7, 13, 17, 19, 31, 37, 43, 47, 49, 59, 61, 67, 71, 73, 79, 91, 97, 101, 103, 107, 109, 127, 133, 137, 139, 149, 151, 157, 163, 167, 169, 181, 193, 197, 199, 211, 217, 223, 227, 229, 241, 247, 257, 259, 263, 269, 271, 277, 283, 289

Offset: 1

N. J. A. Sloane

All primes in A053176 are in the sequence. If n is in the sequence, its factorization contains only primes in A053176. - Benoit Cloitre, Oct 19 2002
B(2n) has denominator 6 iff (n^2-1)*B(2n) is an integer. - Benoit Cloitre, Feb 15 2004
Subsequence of A156543. - Reinhard Zumkeller, Feb 10 2009

B. C. Berndt, Ramanujan's Notebooks Part IV, Springer-Verlag, see p. 76.

Cf. A051222, A053176, A156543.

[n: n in [0..400] | Denominator(Bernoulli(2*n)) eq 6]; // Vincenzo Librandi, Feb 06 2016

ok[n_] := IntegerQ[(n^2 - 1)*BernoulliB[2n]]; Select[Range[300], ok] (* Jean-François Alcover, Jun 27 2012, after Benoit Cloitre *)
result = {}; Do[count = 0;
Do[If[Not[PrimeQ[2*Divisors[n][[i]] + 1]], count++],
{i, 2, DivisorSigma[0, n]}]; If[count == DivisorSigma[0, n] - 1, AppendTo[result, n]], {n, 1, 10000}]; result  (* Richard R. Forberg, Aug 06 2016 *)
Position[BernoulliB[2 Range[300]],?(Denominator[#]==6&)]//Flatten (* _Harvey P. Dale, Jan 28 2017 *)

isok(n) = denominator(bernfrac(2*n)) == 6; \\ Michel Marcus, Feb 06 2016

a(n) seems to be asymptotic to c*n, 5 < c < 6. - Benoit Cloitre, Oct 19 2002

A059455 Safe primes which are also Sophie Germain primes.

Original entry on oeis.org

5, 11, 23, 83, 179, 359, 719, 1019, 1439, 2039, 2063, 2459, 2819, 2903, 2963, 3023, 3623, 3779, 3803, 3863, 4919, 5399, 5639, 6899, 6983, 7079, 7643, 7823, 10163, 10799, 10883, 11699, 12203, 12263, 12899, 14159, 14303, 14699, 15803, 17939

Offset: 1

Labos Elemer, Feb 02 2001

nonn

Primes p such that both (p-1)/2 and 2*p+1 are prime.
Except for 5, all are congruent to 11 modulo 12.
Primes "inside" Cunningham chains of first kind.
Infinite under Dickson's conjecture. - Charles R Greathouse IV, Jul 18 2012
See A162019 for the subset of a(n) that are "reproduced" by the application of the transformations (a(n)-1)/2 and 2*a(n)+1 to the set a(n). - Richard R. Forberg, Mar 05 2015

			83 is a term because it is prime and 2*83+1 = 167 and (83-1)/2 = 41 are both primes.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
Chris K. Caldwell, Cunningham Chains.

Intersection of A005384 and A005385.

Cf. A053176, A059452, A059453, A059455, A059456, A007700, A005602, A023272, A023302, A023330, A156659, A156660, A156877, A162019.

[p: p in PrimesUpTo(20000) |IsPrime((p-1) div 2) and IsPrime(2*p+1)]; // Vincenzo Librandi, Oct 31 2014

lst={}; Do[p=Prime[n]; If[PrimeQ[(p-1)/2]&&PrimeQ[2*p+1], AppendTo[lst, p]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 02 2008 *)
Select[Prime[Range[1000]], AllTrue[{(# - 1)/2, 2 # + 1}, PrimeQ] &] (* requires Mathematica 10+; Feras Awad, Dec 19 2018 *)

forprime(p=2,1e5,if(isprime(p\2)&&isprime(2*p+1),print1(p", "))) \\ Charles R Greathouse IV, Jul 15 2011

from itertools import count, islice
from sympy import isprime, prime
def A059455_gen(): # generator of terms
    return filter(lambda p:isprime(p>>1) and isprime(p<<1|1),(prime(i) for i in count(1)))
A059455_list = list(islice(A059455_gen(),10)) # Chai Wah Wu, Jul 12 2022

A156660(a(n))*A156659(a(n)) = 1; A156877 gives numbers of these numbers <= n. - Reinhard Zumkeller, Feb 18 2009

A119480 Numbers n such that the Bernoulli number B_{4n} has denominator 30.

Original entry on oeis.org

1, 2, 17, 19, 31, 38, 47, 59, 61, 62, 71, 94, 101, 103, 107, 109, 118, 122, 137, 149, 151, 157, 167, 181, 197, 206, 211, 218, 223, 227, 229, 241, 257, 263, 269, 271, 283, 289, 302, 311, 313, 314, 317, 331, 334, 337, 347, 349, 353, 361, 362, 367, 379

Offset: 1

Alexander Adamchuk, Jul 26 2006

nonn

Most a(n) are primes from A043297(n) except for a(1) = 1 and composite a(n) for n=6,10,12,17,18,26,28,38,39,42,45,50,51, ... a(6) = 38 = 2*19, a(10) = 62 = 2*31, a(12) = 94 = 2*47, a(17) = 118 = 2*59, a(18) = 122 = 2*61, a(26) = 206 = 2*103, a(28) = 218 = 2*109, a(38) = 289 = 17*17, a(39) = 302 = 2*151, a(42) = 314 = 2*157, a(45) = 334 = 2*167, a(50) = 361 = 19*19, a(51) = 362 = 2*181, ... It appears that most composite a(n) are the doubles of some primes from A043297(n) belonging to A081092[n] and A045404[n] - Primes congruent to {3, 4, 5, 6} mod 7. The rest of composite a(n) are the squares of the primes from A043297(n).

Some a(n) are the products of different primes from A043297(n), for example a(77) = 527 = 17*31. a(n) belong to A045402 Primes congruent to {1, 3, 4, 5, 6} mod 7. a(n) is a subset of A053176 Primes p such that 2p+1 is composite, A045979 Bernoulli number B_{2n} has denominator 6, A090863 Numbers n such that F(n+1)*F(n-1)*B(2n) is an integer, where F(k)=k-th Fibonacci number and B(2k)=2k-th Bernoulli number. - Alexander Adamchuk, Jul 27 2006

E. Pérez Herrero, Table of n, a(n) for n=1..50000

Cf. A043297, A051225, A051222, A051230, A045404.

Cf. A045402, A053176, A045979, A090863.

Select[Range@ 400, Denominator@ BernoulliB[4 #] == 30 &] (* Michael De Vlieger, Aug 09 2017 *)

a(n) = A051225[n]/2.

A059452 Safe primes (A005385) that are not Sophie Germain primes.

Original entry on oeis.org

7, 47, 59, 107, 167, 227, 263, 347, 383, 467, 479, 503, 563, 587, 839, 863, 887, 983, 1187, 1283, 1307, 1319, 1367, 1487, 1523, 1619, 1823, 1907, 2027, 2099, 2207, 2447, 2579, 2879, 2999, 3119, 3167, 3203, 3467, 3947, 4007, 4079, 4127, 4139, 4259, 4283

Offset: 1

Labos Elemer, Feb 02 2001

nonn

Except for 7, these primes are congruent to 11 modulo 12.

Terminal primes in complete Cunningham chains of first kind, i.e., the chains cannot be continued from these primes.

			347 is a term because 173 is a prime but 695 is not.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Chris K. Caldwell, Cunningham Chain.
Warut Roonguthai, Yves Gallot's Proth.exe and Cunningham Chains. [Wayback Machine link]

Intersection of A005385 and A053176.

Cf. A005384, A059452, A059453, A059454, A059455, A059456, A007700, A005602, A023272, A023302, A023330, A156659, A156660.

lst={};Do[p=Prime[n];If[PrimeQ[(p-1)/2],If[ !PrimeQ[2*p+1],AppendTo[lst,p]]],{n,6!}];lst (* Vladimir Joseph Stephan Orlovsky, Jun 24 2009 *)

is(p) = p > 2 && isprime(p) && isprime((p-1)/2) && !isprime(2*p+1); \\ Amiram Eldar, Jul 15 2024

from itertools import count, islice
from sympy import isprime, prime
def A059452_gen(): # generator of terms
    return filter(lambda p:isprime(p>>1) and not isprime(p<<1|1),(prime(i) for i in count(1)))
A059452_list = list(islice(A059452_gen(),10)) # Chai Wah Wu, Jul 12 2022

A156659(a(n))*(1-A156660(a(n))) = 1. - Reinhard Zumkeller, Feb 18 2009

Broken link updated by R. J. Mathar, Apr 12 2010

A156660 Characteristic function of Sophie Germain primes.

Original entry on oeis.org

0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 0

Reinhard Zumkeller, Feb 13 2009

nonn

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000
Wikipedia, Sophie Germain prime
Index entries for characteristic functions

Cf. A156659.

Cf. A005384, A156874, A092816.

a156660 n = fromEnum $ a010051 n == 1 && a010051 (2 * n + 1) == 1
-- Reinhard Zumkeller, May 01 2012

a(n)=isprime(n)&&isprime(2*n+1) \\ Felix Fröhlich, Aug 11 2014

a(n) = if n and also 2*n+1 is prime then 1 else 0.
a(A005384(n)) = 1; a(A138887(n)) = 0; a(A053176(n)) = 0.
A156874(n) = Sum_{k=1..n} a(k). - Reinhard Zumkeller, Feb 18 2009
a(n) = A010051(n)*A010051(2*n+1).
For n>1 a(n) = floor((floor(phi(n)/(n-1)) + floor(phi(2*n+1)/(2*n)))/2). - Enrique Pérez Herrero, Apr 28 2012
For n>1 a(n) = floor(phi(2*n^2+n)/(2*n^2-2*n)). - Enrique Pérez Herrero, May 02 2012

Definition corrected by Daniel Forgues, Aug 04 2009

A059456 Unsafe primes: primes not in A005385.

Original entry on oeis.org

2, 3, 13, 17, 19, 29, 31, 37, 41, 43, 53, 61, 67, 71, 73, 79, 89, 97, 101, 103, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 173, 181, 191, 193, 197, 199, 211, 223, 229, 233, 239, 241, 251, 257, 269, 271, 277, 281, 283, 293, 307, 311, 313, 317, 331, 337

Offset: 1

Labos Elemer, Feb 02 2001

nonn

A010051(a(n))*(1-A156659(a(n))) = 1; subsequence of A156657. - Reinhard Zumkeller, Feb 18 2009
Also, primes p such that p-1 is a non-semiprime. - Juri-Stepan Gerasimov, Apr 28 2010
Conjecture: From the sequence of prime numbers, let 2 and remove the first data iteration of 2*p+1; leave 3 and remove the prime data by the iteration 2*p+1 and we get the sequence. Example for p=2, remove(5,11,23,47); p=3, remove(7); p=13, p=17, p=19, p=23, remove(47); and so on. - Vincenzo Librandi, Aug 07 2010

			31 is here because (31-1)/2=15 is not prime. 2 and 3 are here because 1/2 and 1 are not prime numbers.

Michael De Vlieger, Table of n, a(n) for n = 1..10000

Cf. A005384, A005385, A053176, A059452-A059456, A007700, A005602, A023272, A023302, A023330.

Initial terms for groups in A075712.

Complement[Prime@ Range@ PrimePi@ Max@ #, #] &@ Select[Prime@ Range@ 90, PrimeQ[(# - 1)/2] &] (* Michael De Vlieger, May 01 2016 *)
Select[Prime[Range[100]],PrimeOmega[#-1]!=2&] (* Harvey P. Dale, May 13 2018 *)

is(n)=isprime(n) && !isprime(n\2) \\ Charles R Greathouse IV, May 02 2016

a(n) ~ n log n. - Charles R Greathouse IV, Dec 29 2024

cp's OEIS Frontend

A059500 Primes p such that both q=(p-1)/2 and 2p + 1 = 4q + 3 are composite numbers. Intersection of A059456 and A053176.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A156543 Multiplicative closure of primes that are not Sophie Germain primes (A053176).

Views

Author

Keywords

Comments

Links

Crossrefs

A079699 Duplicate of A053176.

Views

Author

Keywords

A010051 Characteristic function of primes: 1 if n is prime, else 0.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

Python

Formula

A045979 Bernoulli number B_{2n} has denominator 6.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A059455 Safe primes which are also Sophie Germain primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Formula

A119480 Numbers n such that the Bernoulli number B_{4n} has denominator 30.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A059452 Safe primes (A005385) that are not Sophie Germain primes.

Views