OFFSET
1,2
COMMENTS
Denominators are A118392. Fractions are: 1/1, 5/4, 27/20, 7/5, 10/7, 81/56, 35/24, 22/15, 81/55, 65/44, 77/52, 135/91, 52/35, 119/80, 405/272, 76/51, 85/57, 567/380, 209/140, 115/77, 378/253, 275/184, 299/200, 486/325, 175/117, 377/252, 1215/812, 217/145, 232/155, 1485/992.
2n+3 divides a(2n). 2n-1 divides a(2n-1). p divides a(p) for prime p>2. The only primes in a(n) are a(2) = 5 and a(4) = 7. - Alexander Adamchuk, May 08 2007
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
FORMULA
EXAMPLE
a(1) = 1 = numerator of 1/1.
a(2) = 5 = numerator of 5/4 = 1/1 + 1/4.
a(3) = 27 = numerator of 27/20 = 1/1 + 1/4 + 1/10.
a(4) = 7 = numerator of 7/5 = 1/1 + 1/4 + 1/10 + 1/20.
a(5) = 10 = numerator of 10/7 = 1/1 + 1/4 + 1/10 + 1/20 + 1/35.
a(20) = 115 = numerator of 115/77 = 1/1 + 1/4 + 1/10 + 1/20 + 1/35 + 1/56 + 1/84 + 1/120 + 1/165 + 1/220 + 1/286 + 1/364 + 1/455 + 1/560 + 1/680 + 1/816 + 1/969 + 1/1140 + 1/1330 + 1/1540.
MAPLE
A118391:= n-> numer(3*n*(n+3)/(2*(n+1)*(n+2))); seq(A118391(n), n=1..60) # G. C. Greubel, Feb 18 2021
MATHEMATICA
Table[ Numerator[3n(n+3)/(2(n+1)(n+2))], {n, 1, 100} ] (* Alexander Adamchuk, May 08 2007 *)
Accumulate[1/Binomial[Range[60]+2, 3]]//Numerator (* Harvey P. Dale, Aug 31 2023 *)
PROG
(PARI) s=0; for(i=3, 50, s+=1/binomial(i, 3); print(numerator(s))) /* Phil Carmody, Mar 27 2012 */
(Sage) [numerator(3*n*(n+3)/(2*(n+1)*(n+2))) for n in (1..60)] # G. C. Greubel, Feb 18 2021
(Magma) [Numerator(3*n*(n+3)/(2*(n+1)*(n+2))): n in [1..60]]; // G. C. Greubel, Feb 18 2021
CROSSREFS
KEYWORD
easy,frac,nonn
AUTHOR
Jonathan Vos Post, Apr 27 2006
EXTENSIONS
More terms from Alexander Adamchuk, May 08 2007
STATUS
approved