%I #5 Jul 31 2015 21:38:51
%S 1,1,84,1397520,5314794912000,4855173934730716800000,
%T 1090093558153665322315192780800000,
%U 60907190511553979457004412118419080463155200000
%N Permanent of the inverse Hilbert matrix.
%t Permanent[m_List] := With[{v = Array[x, Length[m]]}, Coefficient[Times @@ (m . v), Times @@ v]]; f[n_] := Block[{i = Inverse[Table[1/(i + j - 1), {i, n}, {j, n}]]}, Permanent[i]]; Table[ f[n], {n, 7}] (* _Robert G. Wilson v_, Oct 24 2005 *)
%o J (www.jsoftware.com)
%o H =: % @: >: @: (+/~) @: i. @ x:
%o perm=: +/ .*
%o perm@%.@H n
%Y Cf. A005249 = determinant of inverse Hilbert matrix; and A092326 = (permanent/determinant) of inverse Hilbert matrix.
%K nonn
%O 0,3
%A _Roger Hui_, Oct 22 2005
|