|
Search: id:A018804
|
|
|
| A018804 |
|
Sum of gcd(k,n) for 1 <= k <= n. |
|
+0
36
|
|
| 1, 3, 5, 8, 9, 15, 13, 20, 21, 27, 21, 40, 25, 39, 45, 48, 33, 63, 37, 72, 65, 63, 45, 100, 65, 75, 81, 104, 57, 135, 61, 112, 105, 99, 117, 168, 73, 111, 125, 180, 81, 195, 85, 168, 189, 135, 93, 240, 133, 195, 165, 200, 105, 243, 189, 260, 185, 171, 117, 360
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
COMMENT
|
a(n) is the number of times the number 1 appears in the character table of the cyclic group C_n. - Ahmed Fares (ahmedfares(AT)my-deja.com), Jun 02 2001
a(n) is the number of ways to express all fractions f/g whereby each product (f/g)*n is a natural number between 1 and n (using fractions of the form f/g with 1 <= f,g <= n). For example, for n=4 there are 8 such fractions: 1/1, 1/2, 2/2, 3/3, 1/4, 2/4, 3/4 and 4/4. - Ron Lalonde (ronronronlalonde(AT)hotmail.com), Oct 03 2002
a(n) is the number of non-congruent solutions to xy = 0 mod n. - Yuval Dekel (dekelyuval(AT)hotmail.com), Oct 06 2003
Equals row sums of triangle A127375 [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 27 2008]
|
|
REFERENCES
|
J. O. Shallit, Problem E 2821, American Mathematical Monthly 87 (1980), 220. Solution in American Mathematical Monthly, 88 (1981), 444-445.
|
|
LINKS
|
T. D. Noe, Table of n, a(n) for n=1..2000
|
|
FORMULA
|
a(n)=Sum_{d|n} d*phi(n/d), where phi(n) is Euler totient function (cf. A000010) - Vladeta Jovovic (vladeta(AT)eunet.rs), Apr 04 2001
Multiplicative; for prime p, a(p^e) = p^(e-1)*((p-1)e+p).
Dirichlet g.f.: zeta(s-1)^2/zeta(s).
a(n)=Sum_{d|n} d*tau(d)*mu(n/d) - Benoit Cloitre (benoit7848c(AT)orange.fr), Oct 23 2003
Equals A054523 * [1,2,3,...]. Equals row sums of triangle A010766. - Gary W. Adamson (qntmpkt(AT)yahoo.com), May 20 2007
Equals Mobius transform of A029935 = A054525 * (1, 2, 4, 5, 8, 8, 12, 12,...). [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 02 2008]
Equals row sums of triangle A127478 [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 03 2008]
|
|
MAPLE
|
a:=n->sum(igcd(n, j), j=1..n): seq(a(n), n=1..60); - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Nov 05 2006
|
|
PROGRAM
|
(PARI) a(n)=direuler(p=2, n, (1-X)/(1-p*X)^2)[n]
|
|
CROSSREFS
|
Cf. A080997, A080998 for rankings of the positive integers in terms of centrality, defined to be the average fraction of an integer that it shares with the other integers as a gcd, or A018804(n)/n^2, also A080999, a permutation of this sequence (A080999(n) = A018804(A080997(n))).
Cf. A010766, A054523.
A127468 [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 03 2008]
A127375 [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 27 2008]
Sequence in context: A050094 A137319 A138808 this_sequence A032682 A022769 A067241
Adjacent sequences: A018801 A018802 A018803 this_sequence A018805 A018806 A018807
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
David W. Wilson (davidwwilson(AT)comcast.net)
|
|
|
Search completed in 0.003 seconds
|
