OFFSET
1,1
COMMENTS
By definition, each term of this sequence is a positive integer.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
MAPLE
bb:= proc(n) option remember; false end: b:= proc(n) option remember; local k, m; if n<3 then bb(n):= true; n else m:= denom(b(n-1) /b(n-2)); for k from m by m while bb(k) do od; bb(k):= true; k fi end: a:= n-> b(n+2) *b(n+1) /b(n): seq(a(n), n=1..100); # Alois P. Heinz, May 18 2009
MATHEMATICA
bb[_] = False;
b[n_] := b[n] = Module[{k, m}, If[n<3, bb[n] = True; n, m = Denominator[ b[n-1]/b[n-2]]; For[k = m, bb[k], k += m]; bb[k] = True; k]];
a[n_] := b[n+2] b[n+1]/b[n];
Array[a, 100] (* Jean-François Alcover, Nov 12 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Leroy Quet, May 06 2009
EXTENSIONS
More terms from Alois P. Heinz, May 18 2009
STATUS
approved