A330836 Numbers of the form 2^(2*p-1)*3*M_p^2, where p > 2 is a Mersenne exponent, A000043, and M_p is the corresponding Mersenne prime, A000668.

4704, 1476096, 396386304, 6753750274277376, 442715102395357986816, 113336363243719574421504, 31901471869127420013759771876790370304, 42404329554681223873219247037048711787234652848116929825491652260298489856

Offset: 1

Author

Walter Kehowski, Jan 12 2020

Keywords

• nonn

Comments

Also numbers with power-spectral basis {M_p^2*(M_p+2)^2, M_p^2*(M_p+1)^2, (M_p^2-1)^2}. The first element of the spectral basis of a(n) is A330819(n+1), the second element is A330837(n+1), and the third element is A330820(n+1). Generally, a power-spectral basis is a spectral basis that consists of primes and powers.

The spectral sum of a(n), that is, the sum of the elements of its spectral basis, is 2*a(n)+1. In this case, we say that a(n) has index 2.

a(n) is also isospectral with A330838(n), that is, a(n) and A330838(n) have the same spectral basis, but A330838(n) has index 1. Thus, A330838(n) and a(n) form an isospectral pair.

Subsequence of Zumkeller numbers (A083207), since a(n) = 2^r * 3 * s, where s is relatively prime to 6. - Ivan N. Ianakiev, Feb 03 2020

Examples

			If p = 3, then a(1) = 2^(2*3-1)*3*7^2 = 4704, and the spectral basis of 4704 is {63^2, 56^2, 48^2}, consisting of powers. The spectral sum of a(1), that is, the sum of the elements of its spectral basis, is 2*4704+1 = 9409. In this case, we say that a(1) has index 2. The number A330838(1) = 9704 has the same spectral basis as a(1), but with index 1. We say that A330838(1) and a(1) are isospectral and form an isospectral pair.
		

Links

• G. Sobczyk, The Missing Spectral Basis in Algebra and Number Theory, The American Mathematical Monthly 108(4), April 2001.
• Wikipedia, Idempotent (ring theory)
• Wikipedia, Peirce decomposition

Crossrefs

Cf. A000043, A000668, A133049, A139306, A330819, A330820, A330837, A330838.

Programs

Maple

a := proc(n::posint)
  local p, m;
  p:=NumberTheory[IthMersenne](n+1);
  m:=2^p-1;
  return 2^(2*p-1)*3*m^2;
end:

Mathematica

f[p_] := 2^(2*p - 1)*3*(2^p - 1)^2; f /@ MersennePrimeExponent /@ Range[2, 9] (* Amiram Eldar, Jan 12 2020 *)

Formula

a(n) = A139306(n+1) * 3 * A133049(n+1).