OFFSET
1,3
EXAMPLE
The sequence's terms, among terms a(1) through a(7), which are coprime to 8 are a(1)=1, a(2)=1, a(5)=3 and a(7) = 23. So a(8) is the numerator of 1 +1 +1/3 +1/23 = 164/69, which is 164.
MATHEMATICA
f[l_List] := Sum[1/l[[k]], {k, Length[l]}]; g[l_List] := Block[{n = Length[l] + 1}, Append[l, Numerator@f[Select[l, GCD[ #, n] == 1 &]]]]; Nest[g, {1}, 13] (* Ray Chandler, Jan 04 2007 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Jan 02 2007
EXTENSIONS
Extended by Ray Chandler, Jan 04 2007
STATUS
approved