cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A360264 Sum of mass(k/n) for all k, 1 <= k <= n, that are relatively prime to n.

Original entry on oeis.org

1, 2, 6, 8, 18, 12, 34, 26, 42, 32, 74, 36, 98, 56, 80, 78, 150, 64, 178, 92, 140, 116, 238, 100, 238, 148, 222, 160, 338, 112, 374, 214, 280, 220, 348, 180, 486, 260, 356, 248, 562, 192, 602, 316, 388, 344, 682, 264, 662, 328, 528, 404, 810, 308, 688, 424
Offset: 1

Views

Author

Jeffrey Shallit, Jan 31 2023

Keywords

Comments

The mass of a rational k/n is the sum of the partial quotients in the continued fraction for k/n. Alternatively, it is the number of steps in the "subtractive algorithm" to compute gcd(k,n).

Examples

			For n = 4 the two numbers relatively prime to n are 1 and 3; 1/4 = [0,4] and 3/4 = [0,1,3].  So the sum of all these is a(3) = 8.
		

Crossrefs

Programs

  • Maple
    a:= n-> add(`if`(igcd(n, k)=1, add(i, i=convert(k/n, confrac)), 0), k=1..n):
    seq(a(n), n=1..60);  # Alois P. Heinz, Jan 31 2023
  • Mathematica
    a[n_] := Total@ Flatten@ (ContinuedFraction[#/n] & /@ Select[Range[n], CoprimeQ[#, n] &]); Array[a, 100] (* Amiram Eldar, Dec 13 2024 *)

Formula

Panov (1982) and Liehl (1983) independently proved that a(n) is asymptotically (6/Pi)^2*n*(log n)^2.