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.

A079340 Absolute value of determinant of n X n matrix whose entries are the integers from 1 to n^2 spiraling outward, ending in a corner.

Original entry on oeis.org

1, 5, 72, 1380, 31920, 861840, 26611200, 925404480, 35805369600, 1526139014400, 71066912716800, 3590219977344000, 195589552648089600, 11430978821982720000, 713448513897799680000, 47363888351558338560000
Offset: 1

Views

Author

Kit Vongmahadlek (kit119(AT)yahoo.com), Jan 03 2003

Keywords

Comments

If n == 0 or 1 (mod 4), the sign of the determinant will be independent of the orientation of the spiral. For n == 2 or 3 (mod 4), the sign will be reversed when the orientation is rotated by 1/4 or flipped on the horizontal or vertical axis. - Franklin T. Adams-Watters, Dec 31 2013
This distribution of the integers is sometimes known as Ulam's spiral, although that is sometimes reserved for when the primes are marked out in some way. - Franklin T. Adams-Watters, Dec 31 2013

Examples

			n=2, det=-5: {1 2 / 4 3 }.
n=3, det=72: {7 8 9 / 6 1 2 / 5 4 3 }.
n=4, det=-1380: { 7 8 9 10 / 6 1 2 11 /5 4 3 12 / 16 15 14 13 }.
n=5, det=31920: { 21 22 23 24 25 / 20 7 8 9 10 / 19 6 1 2 11 /18 5 4 3 12 / 17 16 15 14 13 }
		

Crossrefs

Programs

  • Mathematica
    M[0, 0] = 1;
    M[i_, j_] := If[i <= j,
      If[i + j >= 0, If[i != j, M[i + 1, j] + 1, M[i, j - 1] + 1],
       M[i, j + 1] + 1],
      If[i + j > 1, M[i, j - 1] + 1, M[i - 1, j] + 1]
      ]
    M[n_] := If[EvenQ[n],
      Table[M[i, j], {j, n/2, -n/2 + 1, -1}, {i, -n/2 + 1, n/2}],
      Table[M[i, j], {j, (n - 1)/2, -(n - 1)/2, -1}, {i, -(n - 1)/2, (n - 1)/2}]]
    a[n_]:=Det[M[n]] (* Christian Krattenthaler, Apr 19 2017 *)
  • Maxima
    A079340(n):=if n=1 then 1 else (2*n^2-3*n+3)*(2*n-2)!/(2*(n-1)!)$
    makelist(A079340(n),n,1,30); /* Martin Ettl, Nov 05 2012 */

Formula

a(n) = (2*n^2-3*n+3) (2n-2)!/(2 (n-1)!) = A096376(n-1)*A000407(n-2), n>1. - Conjectured by Dean Hickerson, Jan 30 2003. Proved in the article by Bhatnagar and Krattenthaler.
D-finite with recurrence (2*n^2-7*n+8)*a(n) -2*(2*n-3)*(2*n^2-3*n+3)*a(n-1)=0. - R. J. Mathar, May 03 2019

Extensions

Extended by Robert G. Wilson v, Jan 25 2003