A085000 and related seqs

all at abouthugo.de all at abouthugo.de
Tue Sep 30 23:22:01 CEST 2003


SeqFans,

%S A085000 1,10,412,40800,6839492,1865999570
%N A085000 Maximal determinant of an n X n matrix using the integers 1
to n^2.

is a sequence I've found extremely interesting after stumbling upon it
more or less by chance when browsing through "bref" sequences.
At that time it only had the first 3 terms.

10 days ago I started a thread in sci.math
"Maximal determinant of matrix with elements 1..n^2"
You can locate the thread in Google Groups by searching for
"pfoertner determinant" or read it at
http://mathforum.org/discuss/sci.math/t/539004

This is my last post to sci.math:
(had planned to submit it as A088238, but only one non-zero term
seems a little bit poor)
Hugo Pfoertner wrote:
> 
> Kent Paul Dolan wrote:
> >
> > "Hugo Pfoertner" <nothing at abouthugo.de> wrote;
> > > For n=7 I'm currently ... improved to
> > > a(7)=761973807872 for matrix
> >
> > > 37  6 47 35 18  9 23
> > > 49 30 12  3 26 19 36
> > > 25 40 33 21  1 43 11
> > >  2 20 39  8 34 31 41
> > > 17 22  5 46 13 27 44
> > > 28 10 15 29 45 42  7
> > > 16 48 24 32 38  4 14
> >
> > And yet since you suspect the largest terms should be on the main
> > diagonal, and only one of them is, it seems you have a very long way to
> > go indeed.
> >
> 
> [...]
> 
> > xanthian.
> 
> Using row permutations the matrix can also be written
> (did not check the sign of the determinant)
> 
> 49 30 12  3 26 19 36
> 16 48 24 32 38  4 14
> 37  6 47 35 18  9 23
> 17 22  5 46 13 27 44
> 28 10 15 29 45 42  7
> 25 40 33 21  1 43 11
>  2 20 39  8 34 31 41
> 
> which is not too far from the "max diagonal" form. My earlier
> conjecture that the optimal solutions might have the n largest
> elements on the diagonal does probably not hold for n>5
> because also my n=6 solution, which I couldn't improve has
> different diagonal elements in the lower right corner.
> 
> Hugo

In the meantime the OEIS entry
http://www.research.att.com/projects/OEIS?Anum=A085000 has
been updated twice. There are two minor glitches in the updated page:
The estimate given by Tim Paulden for a(6) was
"a(6) is at least 1862125166." and not 1865999570 (which is the value
I've found a little bit later with my program). The example given for
the
3X3 matrix says:
"The following 3 X 3 matrix is one of 72 whose determinant is 412:"
There are only 36 matrices with det(A_3)=412, the other 36 matrices have
det(A_3)=-412.
Currently I'm searching n=7 with best result found so far
det(A_7)=762074188050 with A_7=
49  5 14 19 23 24 42
25 48 33 15 10  7 37
21 18 47  4 26 43 16
 8  9 39 46 29 13 32
38 31 22 27 45 11  1
28 30 17 44  2 41 12
 6 34  3 20 40 36 35

To test formulas for the asymptotic behavior I've als run the
same program for n=8,9,10:
a(8)>=4.40799*10^14
a(9)>=3.46119*10^17
a(10)>=3.568545*10^20

The heuristic formula a(n)~=0.44 * n^(2.06*n) doesn't fit the data
for n>=6. Can we find a more appropriate formula?

An observation I made for the "optimal" matrices: All row and
column sums are very close or equal to the average value
n*(n^2+1)/2. For A_7 given above all row or column sums are
either 174,175 or 176.

There are closely related problems, that are even harder to solve.
See OEIS:
http://www.research.att.com/projects/OEIS?Anum=A088237
http://www.research.att.com/projects/OEIS?Anum=A088216
http://www.research.att.com/projects/OEIS?Anum=A088217

Maybe the toughest of all:
How many different of these matrices have determinant 0?
I had planned to submit it as A088238, but only one non-zero term
seems a little bit poor :-(

I only know that 0 2X2 matrices and 36*76=2736 3X3 have
det(A_n)=0. (See http://www.research.att.com/projects/OEIS?Anum=A088215
)

Hugo Pfoertner





More information about the SeqFan mailing list