%I #15 Sep 08 2022 08:45:51
%S 1,3,12,32,-625,-24624,-705894,-19922944,-588305187,-18500000000,
%T -622498190424,-22414085849088,-862029149531797,-35320307409809408,
%U -1537494104003906250,-70904672533321089024,-3454944623172347662151,-177423154932124201844736
%N Determinant of the symmetric n X n matrix M_n where M_n(j,k) = n for j = k, M_n(j,n) = n-j, M_n(n,k) = n-k, M_n(j,k) = 0 otherwise.
%D J.-M. Monier, Algèbre et géometrie, exercices corrigés. Dunod, 1997, p. 78.
%H Vincenzo Librandi, <a href="/A174963/b174963.txt">Table of n, a(n) for n = 1..100</a>
%F a(n) = n^n - ((n-1)*n*(2*n-1)/6)*n^(n-2).
%e a(5) = determinant(M_5) = -625 where M_5 is the matrix
%e [5 0 0 0 4]
%e [0 5 0 0 3]
%e [0 0 5 0 2]
%e [0 0 0 5 1]
%e [4 3 2 1 5]
%p with(numtheory):for n from 1 to 25 do:x:=n^n -((n-1)*n*(2*n-1)/6)*n^(n-2):print(x):od:
%o (Magma) [ n^n -((n-1)*n*(2*n-1)/6)*n^(n-2): n in [1..18] ]; // _Klaus Brockhaus_, Apr 11 2010
%o (Magma) [ Determinant( SymmetricMatrix( &cat[ [ i lt j select 0 else n: i in [1..j] ]: j in [1..n-1] ] cat [ 1+((n-1-k) mod n): k in [1..n] ] ) ): n in [1..18] ]; // _Klaus Brockhaus_, Apr 11 2010
%Y Cf. A174962.
%K sign
%O 1,2
%A _Michel Lagneau_, Apr 02 2010
%E Edited by _Klaus Brockhaus_, Apr 11 2010
|