login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A088746 a(n) = determinant of upper left n X n matrix A_n of matrix A088745, where A_n has elements 1..n^2 with abs(det(A_n)) = max by adding elements on the lower and right border. 1

%I

%S 1,-10,205,-6300,276363,-15615642,861774366

%N a(n) = determinant of upper left n X n matrix A_n of matrix A088745, where A_n has elements 1..n^2 with abs(det(A_n)) = max by adding elements on the lower and right border.

%C Terms were computed by Hugo Pfoertner (see A088745 for program link).

%e The upper left 7x7 matrix of A088745 is:

%e [ 1, 3, 8, 11, 24, 26, 48]

%e [ 4, 2, 7, 16, 18, 36, 39]

%e [ 6, 9, 5, 13, 22, 30, 40]

%e [15, 12, 14, 10, 21, 34, 42]

%e [17, 25, 23, 20, 19, 31, 47]

%e [35, 27, 29, 33, 32, 28, 45]

%e [37, 43, 46, 41, 49, 44, 38]

%t (* first do *) Needs["LinearAlgebra`MatrixManipulation`"]; Needs["DiscreteMath`Combinatorica`"] (* then *) f[m_List] := Block[{k = 1, lmt, mxm = 0, mtx = m, mtxm, nbr = Sqrt@Length@Flatten@m, pm}, lmt = (2nbr + 1)!; pm = Reverse@ Range[nbr^2 + 1, (nbr + 1)^2]; While[k <= lmt, mtxr = AppendRows[mtx, Partition[Take[pm, nbr], 1]]; mtxc = AppendColumns[mtxr, {Take[pm, -(nbr + 1)]}]; ab = Det@mtxc; If[Abs@ab >= mxm, mxm = Abs@ab; mtxm = mtxc; Print[{ab, mtxc // MatrixForm}]]; k++; pm = NextPermutation@pm]; mtxm]; NestList[f, {{1}}, 6] (* _Robert G. Wilson v_ *)

%Y Cf. A088745.

%K sign

%O 1,2

%A _Paul D. Hanna_, Oct 17 2003

%E Entries checked by _Robert G. Wilson v_, Sep 22 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 28 22:08 EDT 2022. Contains 354122 sequences. (Running on oeis4.)