A123274 a(0)=a(1)=a(2) = 1. a(n) = (a(n-1) +a(n-2)) /GCD(a(n-1)+a(n-2),a(n-3)), for n >= 3.

1, 1, 1, 2, 3, 5, 4, 3, 7, 5, 4, 9, 13, 11, 8, 19, 27, 23, 50, 73, 123, 98, 221, 319, 270, 589, 859, 724, 1583, 2307, 1945, 4252, 6197, 10449, 8323, 18772, 27095, 45867, 36481, 82348, 118829, 201177, 160003, 361180, 521183, 882363, 701773, 1584136, 2285909

Offset

0,4

Links

Table of n, a(n) for n=0..48.

Mathematica

f[l_List] := Append[l, (l[[ -1]] + l[[ -2]])/GCD[l[[ -1]] + l[[ -2]], l[[ -3]]]]; Nest[f, {1, 1, 1}, 46] (* Ray Chandler, Oct 17 2006 *)

Crossrefs

Sequence in context: A058981 A117339 A096016 * A214674 A185332 A023818

Adjacent sequences: A123271 A123272 A123273 * A123275 A123276 A123277

Keyword

nonn

Author

Leroy Quet, Oct 10 2006

Extensions

Extended by Ray Chandler, Oct 17 2006

Status

approved