login
Start from the sequence of primes, keep the 1st, then delete 2 primes, keep the next, delete 3 primes, keep the next, delete 5 primes, etc ...
0

%I #8 Dec 18 2021 14:52:22

%S 2,7,19,43,79,139,223,317,443,601,809,1021,1291,1601,1949,2311,2729,

%T 3251,3727,4283,4937,5563,6263,6983,7817,8713,9623,10597,11657,12641,

%U 13723,14957,16217,17581,19031,20479,21997,23567,25247,26927,28711,30671,32531

%N Start from the sequence of primes, keep the 1st, then delete 2 primes, keep the next, delete 3 primes, keep the next, delete 5 primes, etc ...

%D J.-P. Delahaye, Des suites fractales d’entiers, Pour la Science, No. 531 January 2022. Sequence c) p. 82.

%o (Python)

%o from sympy import nextprime

%o from itertools import islice

%o def primes(p=2):

%o while True: yield p; p = nextprime(p)

%o def agen():

%o seq1, seq2 = primes(), primes()

%o while True:

%o p, q = next(seq1), next(seq2)

%o yield p

%o for i in range(q):

%o p = next(seq1)

%o print(list(islice(agen(), 43))) # _Michael S. Branicky_, Dec 18 2021

%Y Cf. A000040.

%K nonn

%O 1,1

%A _Michel Marcus_, Dec 18 2021