OFFSET
1,3
EXAMPLE
{b(n)}: 1, 1, 2, 3/2, 19/6, 25/19, 12091/2850, ... Since 1 and 5 are the positive integers which are coprime to 6 and are < 6, b(6) = 1/b(1) + 1/b(5) = 1 + 6/19 = 25/19.
MATHEMATICA
f[l_List] := Block[{n = Length[l] + 1, d}, d = Select[Range[n - 1], GCD[ #, n] == 1 &]; Append[l, Sum[1/l[[d[[i]]]], {i, Length[d]}]]]; Numerator[Nest[f, {1}, 12]] (* Ray Chandler, Feb 09 2007 *)
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
Leroy Quet, Feb 08 2007
EXTENSIONS
Extended by Ray Chandler, Feb 09 2007
STATUS
approved