|
|
A051634
|
|
Strong primes: prime(k) > (prime(k-1) + prime(k+1))/2.
|
|
40
|
|
|
11, 17, 29, 37, 41, 59, 67, 71, 79, 97, 101, 107, 127, 137, 149, 163, 179, 191, 197, 223, 227, 239, 251, 269, 277, 281, 307, 311, 331, 347, 367, 379, 397, 419, 431, 439, 457, 461, 479, 487, 499, 521, 541, 557, 569, 587, 599, 613, 617, 631, 641, 659, 673, 701
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The inequality above is false. The least counterexample is a(19799) = 496283 < A051635(19799) = 496291. - Amiram Eldar, Nov 26 2023
Conjecture: Limit_{N->oo} Sum_{n=1..N} (NextPrime(a(n))-a(n)) / a(N) = 1/4. [A heuristic proof is available at www.primepuzzles.net - Conjecture 91] - Alain Rocchelli, Nov 14 2022
|
|
REFERENCES
|
A. Murthy, Smarandache Notions Journal, Vol. 11 N. 1-2-3 Spring 2000.
|
|
LINKS
|
|
|
FORMULA
|
Conjecture: Limit_{n->oo} n / PrimePi(a(n)) = 1/2. - Alain Rocchelli, Mar 17 2024
|
|
EXAMPLE
|
11 belongs to the sequence because 11 > (7 + 13)/2.
|
|
MAPLE
|
q:= n-> isprime(n) and 2*n>prevprime(n)+nextprime(n):
|
|
MATHEMATICA
|
p=Prime[Range[200]]; p[[Flatten[1+Position[Sign[Differences[p, 2]], -1]]]]
|
|
PROG
|
(Haskell)
a051634 n = a051634_list !! (n-1)
a051634_list = f a000040_list where
f (p:qs@(q:r:ps)) = if 2 * q > (p + r) then q : f qs else f qs
(Python)
from sympy import nextprime
def aupto(limit):
alst, p, q, r = [], 2, 3, 5
while q <= limit:
if 2*q > p + r: alst.append(q)
p, q, r = q, r, nextprime(r)
return alst
|
|
CROSSREFS
|
|
|
KEYWORD
|
nice,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|