A384512 Record terms in A384698.

2, 3, 13, 17, 37, 41, 61, 613, 829, 1861, 2269, 7333, 35149, 1008229, 909889549, 1423665384101, 10341624100573, 440171836495742615578609, 471206109194322691633610979351605854911441181, 4466501842784976704198682186832272945270823914876207595593007001786562643495541

Offset

1,1

Links

Table of n, a(n) for n=1..20.

Formula

a(n) mod 4 = 1, n > 2.

Example

41 is a term because iterating the map on 12 results in a prime in 3 steps: 12 -> 2*12+1=25 -> 25-5=20 -> 2*20+1=41 and 41 is a record prime for starting integers <= 12.

Prog

(Python)
from sympy import isprime, primefactors; rec = 1
for n in range (2, 158163):
    while not isprime(n): n = n - min(primefactors(n)) if n%2 else 2*n + 1
    if n > rec: rec = n; print(n, end = ', ')

Crossrefs

Cf. A020639 (lpf), A383777, A384698.

Sequence in context: A144759 A215359 A115898 * A215350 A118134 A215386

Adjacent sequences: A384509 A384510 A384511 * A384513 A384514 A384515

Keyword

nonn

Author

Ya-Ping Lu, May 31 2025

Status

approved