|
Search: id:A127941
|
|
|
| A127941 |
|
a(n) = numerator of b(n), where b(1)=1, b(n) = sum{1<=k<n,GCD(k,n)=1} 1/b(k). |
|
+0
2
|
|
| 1, 1, 2, 3, 19, 25, 12091, 942579, 802405765442, 9728923151534907, 40387361143971542769608346059516047956596, 16013133965687337121814734147004294320263129575
(list; graph; listen)
|
|
|
OFFSET
|
1,3
|
|
|
LINKS
|
Leroy Quet, Home Page (listed in lieu of email address)
|
|
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, then 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]] (*Chandler*)
|
|
CROSSREFS
|
Cf. A127942.
Adjacent sequences: A127938 A127939 A127940 this_sequence A127942 A127943 A127944
Sequence in context: A038957 A141750 A090476 this_sequence A140555 A058912 A040145
|
|
KEYWORD
|
frac,nonn
|
|
AUTHOR
|
Leroy Quet Feb 08 2007
|
|
EXTENSIONS
|
Extended by Ray Chandler (rayjchandler(AT)sbcglobal.net), Feb 09 2007
|
|
|
Search completed in 0.002 seconds
|
