A333463 a(n) = Sum_{k=1..n} floor(n/k) * gcd(n,k).

1, 4, 7, 13, 14, 27, 22, 38, 39, 53, 39, 85, 49, 81, 93, 106, 68, 143, 78, 165, 144, 143, 98, 243, 147, 176, 189, 252, 131, 338, 143, 281, 251, 243, 279, 440, 178, 278, 308, 470, 200, 514, 212, 438, 488, 350, 234, 660, 339, 522, 427, 538, 271, 670, 487, 714, 489, 462, 307, 1028

Offset

1,2

Links

Chai Wah Wu, Table of n, a(n) for n = 1..10000

Formula

a(n) = Sum_{d|n} phi(n/d) * A006218(d).

a(n) = Sum_{k=1..n} Sum_{d|k} gcd(n,d).

Mathematica

Table[Sum[Floor[n/k] GCD[n, k], {k, 1, n}], {n, 1, 60}]
Table[Sum[EulerPhi[n/d] Sum[DivisorSigma[0, k], {k, 1, d}], {d, Divisors[n]}], {n, 1, 60}]

Prog

(PARI) a(n) = sum(k=1, n, (n\k)*gcd(n, k)); \\ Michel Marcus, Mar 23 2020
(Python)
from math import isqrt
from sympy import divisors, totient
def A333463(n): return sum((2*sum(d//k for k in range(1, isqrt(d)+1))-isqrt(d)**2)*totient(n//d) for d in divisors(n, generator=True)) # Chai Wah Wu, Oct 07 2021

Crossrefs

Cf. A000005, A000010, A006218, A018804, A333465.

Sequence in context: A266568 A244436 A310794 * A310795 A310796 A175091

Adjacent sequences: A333460 A333461 A333462 * A333464 A333465 A333466

Keyword

nonn

Author

Ilya Gutkovskiy, Mar 22 2020

Status

approved