login
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.
0
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
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
KEYWORD
nonn
AUTHOR
Leroy Quet, Oct 10 2006
EXTENSIONS
Extended by Ray Chandler, Oct 17 2006
STATUS
approved