cp's OEIS Frontend

A137324 a(n) = Sum_{prime p < n} gcd(n,p).

%I A137324 #23 Sep 08 2022 08:45:32
%S A137324 1,3,2,6,3,5,6,9,4,8,5,13,12,7,6,10,7,13,16,19,8,12,13,22,11,16,9,17,
%T A137324 10,12,23,28,21,14,11,31,26,17,12,22,13,25,20,37,14,18,21,20,33,28,15,
%U A137324 19,30,23,36,45,16,24,17,49,26,19,34,31,18,36,43,30,19,23,20,58,27,40,37
%N A137324 a(n) = Sum_{prime p < n} gcd(n,p).
%F A137324 a(p) = A000720(p) - 1 for prime p. - _R. J. Mathar_, Apr 09 2008
%F A137324 a(n) = A048865(n) + A105221(n). - _Wesley Ivan Hurt_, Nov 21 2021
%e A137324 a(10) = 9 because gcd(10,2) = 2, gcd(10,3) = 1, gcd(10,5) = 5, gcd(10,7) = 1; 2 + 1 + 5 + 1 = 9.
%e A137324 The underlying irregular table of gcd(n,2), gcd(n,3), gcd(n,5), gcd(n,7), etc., for which a(n) provides row sums, is obtained by deleting columns from A050873(n,k) and looks as follows for n=3,4,5,...:
%e A137324   1
%e A137324   2 1
%e A137324   1 1
%e A137324   2 3 1
%e A137324   1 1 1
%e A137324   2 1 1 1
%e A137324   1 3 1 1
%e A137324   2 1 5 1
%e A137324   1 1 1 1
%e A137324   2 3 1 1 1
%e A137324   1 1 1 1 1
%e A137324   2 1 1 7 1 1
%e A137324   1 3 5 1 1 1
%e A137324   2 1 1 1 1 1
%e A137324   1 1 1 1 1 1
%e A137324   2 3 1 1 1 1 1
%e A137324   1 1 1 1 1 1 1
%e A137324   2 1 5 1 1 1 1 1
%p A137324 A137324 := proc(n) local a,i; a :=0 ; for i from 1 to numtheory[pi](n-1) do a := a+gcd(n,ithprime(i)) ; od: a; end: seq(A137324(n),n=3..80) ; # _R. J. Mathar_, Apr 09 2008
%t A137324 Table[Plus @@ GCD[n, Select[Range[n - 1], PrimeQ[ # ] &]], {n, 3, 70}] (* _Stefan Steinerberger_, Apr 09 2008 *)
%o A137324 (PARI) a(n) = sum(k=1, n-1, gcd(n,k)*isprime(k)); \\ _Michel Marcus_, Nov 07 2014
%o A137324 (Magma) [&+[Gcd(n,p):p in PrimesInInterval(1,n-1)]:n in [3..77]]; // _Marius A. Burtea_, Aug 07 2019
%o A137324 (Python)
%o A137324 from math import gcd
%o A137324 from sympy import primerange
%o A137324 def a(n): return sum(gcd(n, p) for p in primerange(1, n))
%o A137324 print([a(n) for n in range(3, 78)]) # _Michael S. Branicky_, Nov 21 2021
%Y A137324 Cf. A000720, A006579, A048865, A105221.
%K A137324 easy,nonn
%O A137324 3,2
%A A137324 _Max Sills_, Apr 06 2008
%E A137324 Corrected and extended by _R. J. Mathar_ and _Stefan Steinerberger_, Apr 09 2008