Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A126261
a(n) is the numerator of the sum of the reciprocals of the positive integers k, k<=n, where every positive integer <= k and coprime to k is also coprime to n.
1
1, 3, 11, 7, 137, 5, 363, 49, 19, 37, 83711, 7, 1145993, 167, 19, 1321, 42142223, 65, 275295799, 19, 43, 2887, 444316699, 133, 127, 1177, 18469, 85, 9227046511387, 23, 290774257297357, 3877999, 212, 22737, 971, 229, 2040798836801833, 233731, 1039
OFFSET
1,2
EXAMPLE
The positive integers k, k <= 6, where every positive integer <=k and coprime to k is also coprime to 6, are 1,2,6. So a(6) = 5 is the numerator of 1 +1/2 +1/6 = 5/3.
MATHEMATICA
f[n_] := Select[Range[n], GCD[ #, n] == 1 &]; g[n_] := Plus @@ (1/# &) /@ Select[Range[n], Times @@ GCD[f[ # ], n] == 1 &]; Table[Numerator[g[n]], {n, 40}] (* Ray Chandler, Dec 24 2006 *)
CROSSREFS
Sequence in context: A250034 A006495 A112286 * A050097 A330566 A361974
KEYWORD
frac,nonn
AUTHOR
Leroy Quet, Dec 22 2006
EXTENSIONS
Extended by Ray Chandler, Dec 24 2006
STATUS
approved