OFFSET
1,2
COMMENTS
a(2) = 1+3. Every other a(n) of the sequence involves only adding integers that are <= n.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
12 is coprime to 1, 5, 7, 11,... Now 1 = 1, 1+5 = 6, 1+5+7 = 13, 1+5+7+11 = 24, ... 13 is the smallest such partial sum that is >= 12. So a(12) = 13.
MAPLE
printf("1, 4, ") ; for n from 3 to 60 do resul :=0 ; for m from 1 to 1000 do if gcd(n, m) <= 1 then resul := resul + m ; if resul >= n then printf("%a, ", resul) ; break ; fi ; fi ; od : od: # R. J. Mathar, Apr 02 2006
MATHEMATICA
f[n_] := Block[{k = 1, s}, While[s = Plus @@ Select[Range[k], GCD[ #, n] == 1 &]; s < n, k++ ]; s]; Table[f[n], {n, 75}] (* Ray Chandler, Dec 11 2006 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Mar 30 2006
EXTENSIONS
More terms from R. J. Mathar, Apr 02 2006
a(61) and later terms from Ray Chandler, Dec 11 2006
STATUS
approved