login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A295426
a(n) is the numerator of det(I+H) where H is the n X n Hilbert matrix.
2
2, 29, 2927, 659251, 46508430817, 616473989937916861, 3577562384224548869428843, 1314142513507030576449489451528961, 1914627150738259149750867704875720944260093, 59112836238579742851313392516538890376380560892536927, 21782568597204534349136837897139663659824535306651051308429796609
OFFSET
1,1
FORMULA
det(I+H) = Sum_{subsets S of {1,2,...,n}} Product_{1<=i<j<=|S|} (S_i-S_j)^2 / Product_{1<= i,j <= |S|} (S_i+S_j-1).
MAPLE
f := n -> numer(LinearAlgebra:-Determinant(LinearAlgebra:-IdentityMatrix(n)+LinearAlgebra:-HilbertMatrix(n))):
map(f, [$1..30]);
MATHEMATICA
a[n_] := Det[IdentityMatrix[n] + HilbertMatrix[n]] // Numerator;
Array[a, 11] (* Jean-François Alcover, Feb 26 2018 *)
PROG
(PARI) A295426(n) = numerator( matdet( matrix(n, n, i, j, 1/(i+j-1)+(i==j)) ) ); \\ Max Alekseyev, Feb 16 2018
CROSSREFS
Denominators in A295427.
Sequence in context: A059725 A112784 A252042 * A055559 A350858 A345041
KEYWORD
nonn,frac
AUTHOR
Robert Israel, Feb 12 2018
STATUS
approved