%I #19 Dec 15 2023 06:19:16
%S 1,3,1,7,1,3,1,15,10,3,1,7,1,3,1,31,1,30,1,7,1,3,1,15,26,3,10,7,1,3,1,
%T 63,1,3,1,70,1,3,1,15,1,3,1,7,10,3,1,31,50,78,1,7,1,30,1,15,1,3,1,7,1,
%U 3,10,127,1,3,1,7,1,3,1,150,1,3,26,7,1,3,1,31,91,3,1,7,1,3,1,15,1,30,1
%N Sum of squares and twice squares dividing n.
%H Antti Karttunen, <a href="/A065745/b065745.txt">Table of n, a(n) for n = 1..65537</a>
%F Multiplicative with a(2^e) = 2^(e+1)-1, a(p^e) = (p^(e+2)-1)/(p-1)/(p+1) for odd p and even e and a(p^e) = (p^(e+1)-1)/(p-1)/(p+1) for odd p and odd e.
%F From _Amiram Eldar_, Dec 15 2023: (Start)
%F Dirichlet g.f.: (1 + 1/2^(s-1)) * zeta(2*s-2) * zeta(s).
%F Sum_{k=1..n} a(k) ~ c * n^(3/2), where c = ((2+sqrt(2))/6) * zeta(3/2) = 1.4865345575818562471... . (End)
%t f[2, e_] := 2^(e+1) - 1; f[p_, e_] := If[OddQ[e], (p^(e+1)-1)/(p^2-1), (p^(e+2)-1)/(p^2-1)]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Sep 13 2020 *)
%o (PARI) a(n) = sumdiv(n, d, issquare(d)*d + (1 - d%2)*issquare(d/2)*d) \\ _Michel Marcus_, Jun 17 2013
%Y Cf. A035316, A065704, A078434.
%K mult,nonn,easy
%O 1,2
%A _Vladeta Jovovic_, Dec 04 2001
%E More terms from _Matthew Conroy_, Jan 19 2002