OFFSET
1,2
COMMENTS
FORMULA
det(M) = (-1)^(n+1)*(n^2+3*n-2)/(2*(n+1)!), implying that a(n)=p, where p=A006530(n^2+3*n-2) is the largest prime divisor of (n^2+3*n-2), if p>n+1 or p=sqrt((n^2+3*n-2)/2); otherwise a(n)=1.
a(n) = Numerator[(-1)^(n+1) Det[ DiagonalMatrix[ Table[ i/(i+1) - 1, {i, 1, n} ] ] + 1 ]].
a(n) = Numerator[ (n^2+3n-2)/(2(n+1)!) ] = Numerator[ ((2n+3)^2-17)/(4(n+1)!) ].
MATHEMATICA
Numerator[Table[(-1)^(n+1) Det[ DiagonalMatrix[ Table[ i/(i+1) - 1, {i, 1, n} ] ] + 1 ], {n, 1, 70} ]]
Table[ Numerator[ (n^2+3n-2)/(2(n+1)!) ], {n, 1, 100} ]
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
Alexander Adamchuk, May 19 2006, Feb 03 2007
EXTENSIONS
Edited by Max Alekseyev, Jun 02 2009
STATUS
approved