login
Characteristic polynomials of a square matrix based on A051731 where A051731(1,N)=1 and A051731(N,N)=0 and where N=size of matrix, analogous to the Redheffer matrix.
2

%I #42 Feb 16 2025 08:33:36

%S 1,1,-1,-1,-1,1,-1,0,2,-1,0,0,2,-3,1,-1,2,1,-5,4,-1,1,-3,5,-8,9,-5,1,

%T -1,4,-4,-5,15,-14,6,-1,0,-1,6,-17,29,-31,20,-7,1,0,0,2,-13,36,-55,50,

%U -27,8,-1,1,-7,23,-50,84,-112,112,-78,35,-9,1

%N Characteristic polynomials of a square matrix based on A051731 where A051731(1,N)=1 and A051731(N,N)=0 and where N=size of matrix, analogous to the Redheffer matrix.

%C From _Mats Granvik_, Sep 30 2017: (Start)

%C Conjecture: The largest absolute value of the eigenvalues of these characteristic polynomials appear to have the same prime signature in the factorization of the matrix sizes N.

%C In other words: Let b(N) equal the sequence of the largest absolute values of the eigenvalues of the characteristic polynomials of the matrices of size N. b(N) is then a sequence of truncated eigenvalues starting:

%C b(N=1..infinity)

%C = 1.00000, 1.61803, 1.61803, 2.00000, 1.61803, 2.20557, 1.61803, 2.32472, 2.00000, 2.20557, 1.61803, 2.67170, 1.61803, 2.20557, 2.20557, 2.61803, 1.61803, 2.67170, 1.61803, 2.67170, 2.20557, 2.20557, 1.61803, 3.08032, 2.00000, 2.20557, 2.32472, 2.67170, 1.61803, 2.93796, 1.61803, 2.89055, 2.20557, 2.20557, 2.20557, 3.21878, 1.61803, 2.20557, 2.20557, 3.08032, 1.61803, 2.93796, 1.61803, 2.67170, 2.67170, 2.20557, 1.61803, 3.45341, 2.00000, 2.67170, 2.20557, 2.67170, 1.61803, 3.08032, 2.20557, 3.08032, 2.20557, 2.20557, 1.61803, 3.53392, 1.61803, 2.20557, 2.67170, ...

%C It then appears that for n = 1,2,3,4,5,...,infinity we have the table:

%C Prime signature: b(Axxxxxx(n)) = Largest abs(eigenvalue):

%C p^0 : b(1) = 1.0000000000000000000000000000...

%C p : b(A000040(n)) = 1.6180339887498949025257388711...

%C p^2 : b(A001248(n)) = 2.0000000000000000000000000000...

%C p*q : b(A006881(n)) = 2.2055694304005917238953315973...

%C p^3 : b(A030078(n)) = 2.3247179572447480566665944934...

%C p^2*q : b(A054753(n)) = 2.6716998816571604358216518448...

%C p^4 : b(A030514(n)) = 2.6180339887498917939012699207...

%C p^3*q : b(A065036(n)) = 3.0803227214906021558249449299...

%C p*q*r : b(A007304(n)) = 2.9379558827528557962693867011...

%C p^5 : b(A050997(n)) = 2.8905508875432590620846440288...

%C p^2*q^2 : b(A085986(n)) = 3.2187765853016649941764626419...

%C p^4*q : b(A178739(n)) = 3.4534111136673804054453285061...

%C p^2*q*r : b(A085987(n)) = 3.5339198574905377192578725953...

%C p^6 : b(A030516(n)) = 3.1478990357047909043330946587...

%C p^3*q^2 : b(A143610(n)) = 3.7022736187975437971431347250...

%C p^5*q : b(A178740(n)) = 3.8016448153137023524550386355...

%C p^3*q*r : b(A189975(n)) = 4.0600260453688532535920785448...

%C p^7 : b(A092759(n)) = 3.3935083220984414431597997463...

%C p^4*q^2 : b(A189988(n)) = 4.1453038440113498808159420150...

%C p^2*q^2*r: b(A179643(n)) = 4.2413382309993874486053755390...

%C p^6*q : b(A189987(n)) = 4.1311805192254587026923218218...

%C p*q*r*s : b(A046386(n)) = 3.8825338629275134572083061357...

%C ...

%C b(Axxxxxx(1)) in the sequences above, is given by A025487.

%C (End)

%C First column in the coefficients of the characteristic polynomials is the Möbius function A008683.

%C Row sums of coefficients start: 0, -1, 0, 0, 0, 0, 0, 0, 0, ...

%C Third diagonal is a signed version of A000096.

%C Most of the eigenvalues are equal to 1. The number of eigenvalues equal to 1 are given by A075795 for n>1.

%C The first three of the eigenvalues above can be calculated as nested radicals. The fourth eigenvalue 2.205569430400590... minus 1 = 1.205569430400590... is also a nested radical.

%H Mats Granvik, <a href="/A275345/a275345_1.txt">Mathematica program to verify that eigenvalues determine prime signature</a>

%H OEIS Wiki, <a href="http://oeis.org/wiki/Prime_signatures">Prime signatures</a>

%H Eric Weisstein, <a href="https://mathworld.wolfram.com/PrimeSignature.html">Prime signature</a>

%H <a href="https://oeis.org/wiki/Index_to_OEIS:_Section_Pri#prime_signature">Index to sequences related to prime signature</a>

%e {

%e { 1},

%e { 1, -1},

%e {-1, -1, 1},

%e {-1, 0, 2, -1},

%e { 0, 0, 2, -3, 1},

%e {-1, 2, 1, -5, 4, -1},

%e { 1, -3, 5, -8, 9, -5, 1},

%e {-1, 4, -4, -5, 15, -14, 6, -1},

%e { 0, -1, 6, -17, 29, -31, 20, -7, 1},

%e { 0, 0, 2, -13, 36, -55, 50, -27, 8, -1},

%e { 1, -7, 23, -50, 84, -112, 112, -78, 35, -9, 1}

%e }

%t Clear[x, AA, nn, s]; Monitor[AA = Flatten[Table[A = Table[Table[If[Mod[n, k] == 0, 1, 0], {k, 1, nn}], {n, 1, nn}]; MatrixForm[A]; a = A[[1, nn]]; A[[1, nn]] = A[[nn, nn]]; A[[nn, nn]] = a; CoefficientList[CharacteristicPolynomial[A, x], x], {nn, 1, 10}]], nn]

%Y Cf. A051731, A008683, A000040, A001248, A006881, A030078, A030514, A054753, A000096, A001622, A272874, A075795.

%Y Cf. A025487. - _Mats Granvik_, Sep 30 2017

%K sign,tabl

%O 0,9

%A _Mats Granvik_, Jul 24 2016