login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A046657 a(n) = A002088(n)/2. 7
1, 2, 3, 5, 6, 9, 11, 14, 16, 21, 23, 29, 32, 36, 40, 48, 51, 60, 64, 70, 75, 86, 90, 100, 106, 115, 121, 135, 139, 154, 162, 172, 180, 192, 198, 216, 225, 237, 245, 265, 271, 292, 302, 314, 325, 348, 356, 377, 387, 403, 415, 441, 450, 470, 482 (list; graph; refs; listen; history; text; internal format)
OFFSET
2,2
COMMENTS
a(n) = |{(x,y) : 1 <= x <= y <= n, x + y <= n, 1 = gcd(x,y)}| = |{(x,y) : 1 <= x <= y <= n, x + y > n, 1 = gcd(x,y)}|. - Steve Butler, Mar 31 2006
Brousseau proved that if the starting numbers of a generalized Fibonacci sequence are <= n (for n > 1) then the number of such sequences with relatively prime successive terms is a(n). - Amiram Eldar, Mar 31 2017
LINKS
Alfred Brousseau, A Note on the Number of Fibonacci Sequences, The Fibonacci Quarterly, Vol. 10, No. 6 (1972), pp. 657-658.
FORMULA
a(n) = 1/2 + Sum_{i<j<=n, gcd(i, j) = 1} i/j. - Joseph Wheat, Feb 22 2018
MAPLE
a:=n->sum(numtheory[phi](k), k=1..n): seq(a(n)/2, n=2..60); # Muniru A Asiru, Mar 05 2018
MATHEMATICA
Rest@ Accumulate[EulerPhi@ Range@ 56]/2 (* Michael De Vlieger, Apr 02 2017 *)
PROG
(PARI) a(n) = sum(k=1, n, eulerphi(k))/2; \\ Michel Marcus, Apr 01 2017
(GAP) List([2..60], n->Sum([1..n], k->Phi(k)/2)); # Muniru A Asiru, Mar 05 2018
(Python)
from functools import lru_cache
@lru_cache(maxsize=None)
def A046657(n): # based on second formula in A018805
if n == 0:
return 0
c, j = 0, 2
k1 = n//j
while k1 > 1:
j2 = n//k1 + 1
c += (j2-j)*(4*A046657(k1)-1)
j, k1 = j2, n//j2
return (n*(n-1)-c+j)//4 # Chai Wah Wu, Mar 25 2021
CROSSREFS
Cf. A002088.
Partial sums of A023022.
Sequence in context: A104214 A349523 A355616 * A102825 A070991 A225527
KEYWORD
nonn
AUTHOR
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 21:54 EDT 2024. Contains 371254 sequences. (Running on oeis4.)