cp's OEIS Frontend

A317797 Sum of the norm of divisors of n over Gaussian integers, with associated divisors counted only once.

%I A317797 #28 Dec 06 2024 10:24:54
%S A317797 1,7,10,31,36,70,50,127,91,252,122,310,196,350,360,511,324,637,362,
%T A317797 1116,500,854,530,1270,961,1372,820,1550,900,2520,962,2047,1220,2268,
%U A317797 1800,2821,1444,2534,1960,4572,1764,3500,1850,3782,3276,3710,2210,5110,2451,6727
%N A317797 Sum of the norm of divisors of n over Gaussian integers, with associated divisors counted only once.
%C A317797 Equivalent of sigma (A000203) in the ring of Gaussian integers. Note that only norms are summed up.
%H A317797 Jianing Song, <a href="/A317797/b317797.txt">Table of n, a(n) for n = 1..10000</a>
%H A317797 Wikipedia, <a href="https://en.wikipedia.org/wiki/Gaussian_integer">Gaussian integer</a>
%F A317797 Multiplicative with a(2^e) = sigma(2^(2e)) = 2^(2e+1) - 1, a(p^e) = sigma(p^e)^2 = ((p^(e+1) - 1)/(p - 1))^2 if p == 1 (mod 4) and sigma_2(p^e) = A001157(p^e) = (p^(2e+2) - 1)/(p^2 - 1) if p == 3 (mod 4).
%e A317797 Let ||d|| denote the norm of d.
%e A317797 a(2) = ||1|| + ||1 + i|| + ||2|| = 1 + 2 + 4 = 7.
%e A317797 a(5) = ||1|| + ||2 + i|| + ||2 - i|| + ||5|| = 1 + 5 + 5 + 25 = 36. Note that 2 - i and 1 + 2i are associated so their norm is only counted once.
%t A317797 f[p_, e_] := If[p == 2, 2^(2*e + 1) - 1, Switch[Mod[p, 4], 1, ((p^(e + 1) - 1)/(p - 1))^2, 3, (p^(2 e + 2) - 1)/(p^2 - 1)]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Feb 12 2020 *)
%o A317797 (PARI)
%o A317797 a(n)=
%o A317797 {
%o A317797     my(r=1, f=factor(n));
%o A317797     for(j=1, #f[, 1], my(p=f[j, 1], e=f[j, 2]);
%o A317797         if(p==2, r*=(2^(2*e+1)-1));
%o A317797         if(Mod(p,4)==1, r*=((p^(e+1)-1)/(p-1))^2);
%o A317797         if(Mod(p,4)==3, r*=(p^(2*e+2)-1)/(p^2-1));
%o A317797     );
%o A317797     return(r);
%o A317797 }
%Y A317797 Cf. A001157.
%Y A317797 Equivalent of arithmetic functions in the ring of Gaussian integers (the corresponding functions in the ring of integers are in the parentheses): A062327 ("d", A000005), this sequence ("sigma", A000203), A079458 ("phi", A000010), A227334 ("psi", A002322), A086275 ("omega", A001221), A078458 ("Omega", A001222), A318608 ("mu", A008683).
%Y A317797 Equivalent in the ring of Eisenstein integers: A319449.
%K A317797 nonn,mult,look
%O A317797 1,2
%A A317797 _Jianing Song_, Aug 07 2018