Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #16 Oct 24 2020 16:51:55
%S 2,3,5,7,11,13,19,31,109,139,149,179,199,229,239,683,739,809,823,907,
%T 977,991,1019,1033,1061,1103,1117,1187,1201,1229,1327,1439,1453,1481,
%U 1523,1579,1607,1621,1663,1733,1747,1789,4931,4967,5003,5021,5039,5147,5237
%N The prime numbers in A158923.
%C This sequence is consists of the prime numbers in A158923, in which A158923(m) is the sum of the previous term A158923(m-1) and the average prime gap log(A158923(m-1)) rounded to the nearest integer with A158923(1) = 2.
%o (Python)
%o from sympy import isprime
%o from math import log
%o print(2)
%o a_last = p_last = m = 2
%o n = 1
%o while m >= 2:
%o a = a_last + int(log(a_last) + 0.5)
%o if isprime(a) == 1:
%o print(a)
%o n += 1
%o a_last = a
%o m += 1
%Y Cf. A000040, A158923.
%K nonn
%O 1,1
%A _Ya-Ping Lu_, Oct 23 2020