cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A102506 Numbers n such that for some positive number k, z=n+ik is a complex multiperfect number; that is, z divides sigma(z), where sigma is the sum of divisors function extended to the complex numbers.

Original entry on oeis.org

1, 5, 6, 10, 12, 28, 60, 72, 100, 108, 120, 140, 150, 204, 263, 300, 526, 600, 672, 720, 912, 1200, 1470, 1520, 1704, 3600, 4560, 4680, 4828, 5584, 5880, 6240, 6312
Offset: 1

Views

Author

T. D. Noe, Jan 12 2005

Keywords

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 Product_{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?

Examples

			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.
		

Crossrefs

Cf. A102507. Note that A101367 and A101366 use Mathematica's Divisors function, the sum of the first-quadrant divisors, which does not enjoy the nice multiplicative properties of Spira's sigma function.

Programs

  • 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]

A102531 Real part of absolute Gaussian perfect numbers, in order of increasing magnitude.

Original entry on oeis.org

3, 15, 6, 19, 111, 91, 159, 72, 472, 904, 2584, 1616, 999, 4328, 702, 4424, 7048, 7328, 2474, 9352, 7144, 7240, 5117, 739, 6327, 15128, 13168, 1263, 14280, 3224, 21704, 15160, 21992, 14044, 23132, 9135, 23656, 24614, 7272, 15464, 9040, 28424, 30956, 14728, 32399
Offset: 1

Views

Author

T. D. Noe, Jan 13 2005

Keywords

Comments

An absolute Gaussian perfect number z satisfies abs(sigma(z)-z) = abs(z), where sigma(z) is sum of the divisors of z, as defined by Spira for Gaussian integers.

Examples

			For z=3+7i, we have sigma(z)-z = 7+3i, which has the same magnitude as z.
		

Crossrefs

See A102532 for the imaginary part.
Cf. A102506 and A102507 (Gaussian multiperfect numbers). See also A101366, A101367.

Programs

  • Mathematica
    lst={}; nn=1000; Do[z=a+b*I; If[Abs[z]<=nn && Abs[(DivisorSigma[1, z]-z)] == Abs[z], AppendTo[lst, {Abs[z]^2, z}]], {a, nn}, {b, nn}]; Re[Transpose[Sort[lst]][[2]]]

Extensions

a(22)-a(45) from Amiram Eldar, Feb 10 2020

A102532 Imaginary part of absolute Gaussian perfect numbers, in order of increasing magnitude.

Original entry on oeis.org

7, 42, 57, 82, 78, 189, 341, 549, 664, 1048, 1016, 3776, 4072, 1672, 5049, 3816, 1128, 368, 7097, 2504, 6816, 7912, 10743, 12605, 12606, 248, 10392, 16853, 15208, 20824, 344, 15688, 392, 18076, 1340, 22302, 5912, 1573, 24392, 21328, 24896, 7048, 1220, 29288, 9678
Offset: 1

Views

Author

T. D. Noe, Jan 13 2005

Keywords

Comments

See A102531 for the real part.

Crossrefs

Programs

  • Mathematica
    lst={}; nn=1000; Do[z=a+b*I; If[Abs[z]<=nn && Abs[(DivisorSigma[1, z]-z)] == Abs[z], AppendTo[lst, {Abs[z]^2, z}]], {a, nn}, {b, nn}]; Im[Transpose[Sort[lst]][[2]]]

Extensions

a(22)-a(45) from Amiram Eldar, Feb 10 2020

A101367 Perfect Abs. Real part of complex z such that Abs[(Total[Divisors[z]]-z)]=Abs[z].

Original entry on oeis.org

5, 3, 19, 15, 29, 6, 74, 19, 111, 147, 185, 91, 197, 269, 122, 159, 72, 827, 1487, 2903, 968, 999, 702, 5803, 326, 2474, 7871
Offset: 0

Views

Author

Ed Pegg Jr, Jan 13 2005

Keywords

Comments

Having Perfect Abs is not as good as being Perfect. A complex number can also have Abundant Abs or Deficient Abs.

Examples

			The divisors for 269+92i are: 1, 2+I, 3+4i, 6+5i, 7+2i, 7+16i, 12+11i, 13+34i, 17+126i, 32+47i, 39+2i, 269+92i. The (sum - k) is 139+248i. Abs[139+248i] == Abs[269+92i]
		

Crossrefs

Programs

  • Mathematica
    Re[Sort[Select[Flatten[Table[a + b I, {a, 1, 500}, {b, 1, 500}]], Abs[Total[Divisors[ # ]] - # ] == Abs[ # ] &], Abs[ #1] < Abs[ #2] &]]

Extensions

Ten more terms from Hans Havermann, Jan 15 2005
Showing 1-4 of 4 results.