OFFSET
1,1
COMMENTS
Sequence A158923, in which every term is the sum of the previous term and the average prime gap, is a "simulation" of the prime number sequence A000040. This sequence lists the terms in A158923 that match those in A000040 both in value and in position, or A158923(m) = A000040(m).
There are 68 matches found for m up to 1073741824 (prime(1073741824)=24563311309), with a(68) = 12496469849.
PROG
(Python)
from sympy import nextprime
from math import log
print(2)
a_last = b_last = m = 2
n = 1
while m >= 2:
a = a_last + int(log(a_last) + 0.5)
b = nextprime(b_last)
if a == b:
n += 1
print (m)
a_last = a
b_last = b
m += 1
(PARI) lista(nn) = {my(a = 2); for (n=1, nn, if (a == prime(n), print1(a, ", ")); a += round(log(a)); ); } \\ Michel Marcus, Oct 31 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ya-Ping Lu, Oct 24 2020
STATUS
approved