OFFSET
1,1
LINKS
Sebastian F. Orellana, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = gpf(a(2) + prime(3)) = gpf(5+5) = 5.
MATHEMATICA
gpf[n_] := FactorInteger[n][[-1, 1]]; a[1] = 2; a[n_] := a[n] = gpf[a[n - 1] + Prime[n]]; Array[a, 100] (* Amiram Eldar, Jun 15 2023 *)
PROG
(Python)
from sympy import factorint, prime
list=[2]
num=1
k=50
for i in range(0, k):
list.append(max(factorint(list[i]+prime(i+1))))
print(list)
CROSSREFS
KEYWORD
nonn
AUTHOR
Sebastian F. Orellana, Jun 12 2023
STATUS
approved