|
%I #8 May 10 2013 12:45:08
%S 1,5,4,13,12,20,8,29,13,56,12,52,24,40,48,61,28,65,20,144,32,60,24,
%T 116,81,112,40,104,44,224,32,125,48,136,96,169,52,100,96,320,60,160,
%U 44,156,156,120,48,244,57,369,112,288,72,200,144,232,80,212,60,576,84,160,104
%N Let r+i*s be the sum of the distinct first-quadrant Gaussian integers dividing n; sequence gives r+s values.
%C A078910(n)-A078911(n) = A000203(n).
%C Also equal to the sum of all divisors of n (in Z[i]) having positive real part. - _M. F. Hasler_, Nov 22 2007
%H M. F. Hasler, <a href="/A078930/b078930.txt">Table of n, a(n) for n = 1..1000</a>.
%H Project Euler, <a href="http://projecteuler.net/index.php?section=problems&id=153">Problem 153: Investigating Gaussian Integers</a>.
%F a(n) = A078910(n)+A078911(n).
%F a(n) = A000203(n)+2*A078911(n). - M. F. Hasler, Nov 22 2007
%o (PARI) A078930(n,S=[])=sigma(n)+sumdiv(n*I,d,if(real(d)&imag(d)&!setsearch(S,d=vecsort(abs([real(d),imag(d)]))),S=setunion(S,[d]);(d[1]+d[2])<<(d[1]!=d[2]))) - _M. F. Hasler_, Nov 22 2007
%K nonn
%O 1,2
%A _Vladeta Jovovic_, Jan 11 2003
|
