OFFSET
1,2
MAPLE
A076920 := proc(nmax) local a, b, k; a := [1, 2] ; b := [1] ; while nops(b) < nmax do k := op(-1, a)+1 ; while gcd(k, op(-1, a)) <= 1 or gcd(k, op(-1, a)) = gcd(op(-1, a), op(-2, a)) do k := k+1 ; od ; b := [op(b), gcd(k, op(-1, a))] ; a := [op(a), k] ; od ; RETURN(b) ; end: A076920(80) ; # R. J. Mathar, Jul 01 2007
MATHEMATICA
f[1] = 1; f[2] = 2; (* f is A076919 *)
f[n_] := f[n] = Module[{k}, For[k = f[n-1] + 1, True, k++, If[GCD[f[n-1], f[n-2]] != GCD[k, f[n-1]] && GCD[k, f[n-1]] > 1, Return[k]]]];
GCD @@@ Partition[Table[f[n], {n, 1, 81}], 2, 1] (* Jean-François Alcover, Oct 25 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Oct 17 2002
EXTENSIONS
More terms from R. J. Mathar, Jul 01 2007
STATUS
approved