OFFSET
1,2
COMMENTS
This sequence uses a number-theoretic extension of the sigma function that is due to Spira. A nonzero Gaussian integer has a unique factorization as u q1^e1 q2^e2..qn^en, where u is a unit (1,-1,i,-i), the qk are Gaussian primes in the first quadrant and the ek are positive integers.
Then Spira defines the sum of divisors to be prod_{k=1..n) (qk^(ek+1)-1)/(qk-1). This appears to be the natural number-theoretic extension. Spira's definition preserves the multiplicative property: if GCD(x,y)=1, then sigma(x*y)=sigma(x)*sigma(y). (Mathematica's DivisorSigma function uses this formula.)
It appears that the value of k, A102507, is unique for each n. The sum of divisors function, as defined by Spira, is implemented in Mathematica for complex z as the DivisorSigma[1,z]. For the z=n+ik given here, sigma(z)/z is usually a small Gaussian integer. The first instance of a positive integral value of sigma(z)/z is z=600+3800i, in which case the ratio is 3. The complex multiperfect numbers can be arranged into classes according to the value of sigma(z)/z. Does each class have a finite number of members?
LINKS
R. Spira, The Complex Sum Of Divisors, American Mathematical Monthly, 1961 Vol. 68, pp. 120-124.
Eric Weisstein's World of Mathematics, Multiperfect Number
EXAMPLE
For n=1, we have z=1+3i. The divisors of z are 1, 1+i, 1+3i and 2+i. Hence sigma(z)=5+5i and sigma(z)/z = 2-i.
MATHEMATICA
lst={}; Do[z=n+k*I; s=DivisorSigma[1, z]; If[Mod[s, z]==0, AppendTo[lst, z]; Print[{z, s, s/z}]], {n, 1200}, {k, 10000}]; Re[lst]
CROSSREFS
KEYWORD
nice,nonn,more
AUTHOR
T. D. Noe, Jan 12 2005
STATUS
approved