A059456 Unsafe primes: primes not in A005385.
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
Keywords
Examples
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.
Links
- Michael De Vlieger, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
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 *)
-
PARI
is(n)=isprime(n) && !isprime(n\2) \\ Charles R Greathouse IV, May 02 2016
Formula
a(n) ~ n log n. - Charles R Greathouse IV, Dec 29 2024
Comments