|
%I #12 Jun 25 2018 22:55:26
%S 0,1,1,5,1,21,1,21,10,50,1,85,1,85,50,85,1,210,1,210,91,210,1,341,26,
%T 250,91,341,1,610,1,341,170,455,122,850,1,546,260,850,1,1300,1,850,
%U 500,850,1,1365,50,1300,362,1050,1,1911,260,1365,500,1300,1,2562,1,1365
%N Sum_{d |C(n)} d^2, where C(n) is the Cototient function n - phi(n) (A051953).
%H Harry J. Smith, <a href="/A063476/b063476.txt">Table of n, a(n) for n = 1..1000</a>
%o (PARI) C(n)=n-eulerphi(n); j=[]; for(n=1,200,j=concat(j,sumdiv(C(n),d,d^2))); j
%o (PARI) { for (n=1, 1000, if (n>1, a=sumdiv(n-eulerphi(n), d, d^2), a=0); write("b063476.txt", n, " ", a) ) } \\ _Harry J. Smith_, Aug 22 2009
%Y Cf. A051953.
%K nonn
%O 1,4
%A _Jason Earls_, Jul 27 2001
|
