OFFSET
1,1
MAPLE
G:= proc(nv)
local n, v, k, np;
n:= nv[1]; v:= nv[2];
k:= min(seq((-n-1) mod p, p = numtheory:-factorset(v-n-1)));
np:= n+k+1;
[np, v+k+igcd(np, v+k)^2]
end proc:
R:= 3: t:= [3, 12]:
for i from 2 to 24 do
t:= G(t);
R:= R, t[1];
od:
R;
MATHEMATICA
2 + Position[Differences @ Nest[Append[#1, #1[[-1]] + GCD[#2, #1[[-1]]]^2] & @@ {#, Length[#] + 2} &, {3}, 10000], _?(# > 1 &)] // Flatten (* Amiram Eldar, Apr 14 2021 after Michael De Vlieger at A260815 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Apr 13 2021
STATUS
approved