|
%I #19 Apr 07 2021 19:54:55
%S 1,6,120,5040,554400,86486400,23524300800,8045310873600,
%T 4070927302041600,3305592969257779200,3074201461409734656000,
%U 4094836346597766561792000,6715531608420337161338880000,12128250084807128913378017280000
%N Sum of the integers in the reduced residue system of A002110(n).
%C Sum of the integers up to A002110(n) and coprime to A002110(n).
%C The sequence gives the sum of row n of A286941(n).
%H Michael De Vlieger, <a href="/A335334/b335334.txt">Table of n, a(n) for n = 1..196</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Totative.html">Totative</a>.
%F a(n) = A023896(A002110(n)).
%F a(n) = A002110(n)*A005867(n)/2 = A070826(n)*A005867(n).
%F a(n) = (A002110(n)*A038110(n+1)/2)*A058250(n).
%e For n = 3: A002110(3) = 30, the reduced residue system of 30 is {1, 7, 11, 13, 17, 19, 23, 29}. The sum is a(3) = 120.
%t n = 15;
%t A002110 = Drop[FoldList[Times, 1, Prime[Range[n]]], 1];
%t A005867 = Drop[EulerPhi@FoldList[Times, 1, Prime@Range@n], 1];
%t A002110*A005867/2
%t (* Second program: *)
%t Map[# EulerPhi[#]/2 &, FoldList[Times, Prime@ Range@ 14]] (* _Michael De Vlieger_, Apr 07 2021 *)
%o (PARI) a(n) = my(P=factorback(primes(n))); P*eulerphi(P)/2; \\ _Michel Marcus_, Jun 02 2020
%Y Cf. A023896, A002110, A005867, A070826, A038110, A058250, A286941.
%K nonn
%O 1,2
%A _Jamie Morken_, Jun 02 2020
| 