|
Search: id:A118679
|
|
|
| A118679 |
|
Absolute value of numerator of determinant of n X n matrix with M(i,j) = i/(i+1) if i=j otherwise 1. |
|
+0
6
|
|
| 1, 2, 1, 13, 19, 13, 17, 43, 53, 1, 19, 89, 103, 59, 67, 151, 13, 47, 1, 229, 251, 137, 149, 1, 349, 47, 101, 433, 463, 1, 263, 43, 593, 157, 83, 701, 739, 389, 409, 859, 53, 59, 1, 1033, 83, 563, 587, 1223, 67, 331, 1, 1429, 1483, 769, 797, 127, 1709, 1, 457, 1889
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
COMMENT
|
Some a(n) are equal to 1 (n=1,3,10,19,24,30,43..). It appears that all other a(n) are primes that belong to A038889 (17 is a square mod p).
Numbers n such that a(n) = 1 are listed in A127852 = {1,3,10,19,24,30,43,51,58,62,73,75,82,94,101,106,115,116,...}. All a(n)>1 are prime belonging to A038889 (17 is a square mod p).
|
|
FORMULA
|
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
|
Cf. A038889.
Cf. A118680, A127852, A127853.
Adjacent sequences: A118676 A118677 A118678 this_sequence A118680 A118681 A118682
Sequence in context: A074808 A113097 A032001 this_sequence A087451 A063558 A013020
|
|
KEYWORD
|
frac,nonn
|
|
AUTHOR
|
Alexander Adamchuk (alex(AT)kolmogorov.com), May 19 2006, Feb 03 2007
|
|
|
Search completed in 0.002 seconds
|
