cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A338396 The prime numbers in A158923.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 19, 31, 109, 139, 149, 179, 199, 229, 239, 683, 739, 809, 823, 907, 977, 991, 1019, 1033, 1061, 1103, 1117, 1187, 1201, 1229, 1327, 1439, 1453, 1481, 1523, 1579, 1607, 1621, 1663, 1733, 1747, 1789, 4931, 4967, 5003, 5021, 5039, 5147, 5237
Offset: 1

Views

Author

Ya-Ping Lu, Oct 23 2020

Keywords

Comments

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.

Crossrefs

Programs

  • Python
    from sympy import isprime
    from math import log
    print(2)
    a_last = p_last = m = 2
    n = 1
    while m >= 2:
        a = a_last + int(log(a_last) + 0.5)
        if isprime(a) == 1:
            print(a)
            n += 1
        a_last = a
        m += 1