OFFSET
0,2
COMMENTS
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Wikipedia, Quasi-polynomial.
Index entries for linear recurrences with constant coefficients, signature (3,-6,10,-12,12,-10,6,-3,1).
FORMULA
a(n) = +3*a(n-1) -6*a(n-2) +10*a(n-3) -12*a(n-4) +12*a(n-5) -10*a(n-6) +6*a(n-7) -3*a(n-8) +a(n-9).
a(n) = 3*a(n-4) -3*a(n-8) +a(n-12).
G.f.: x*(-3+2*x-3*x^2-3*x^3+x^7) / ( (x-1)^3*(x^2+1)^3 ).
a(n) = n*(n+5)*((-1)^((2*n-(-1)^n-3)/4)+3)/8 = n*(n+5)*(3-i^(n*(n+1)))/8, where i=sqrt(-1); also a(n) = a(n-4)*A028557(n)/A028557(n-4) for n>4. - Bruno Berselli, Dec 30 2010
From Peter Bala, Aug 07 2022: (Start)
a(n) = numerator of n*(n+5)/4.
a(n) is quasi-polynomial in n: a(4*n) = n*(4*n+5) = A343560(n+1); a(4*n+1) = (2*n+3)*(4*n+1); a(4*n+2) = (2*n+1)*(4*n+7); a(4*n+3) = (n+2)*(4*n+3) = A180863(n+2). (End)
Sum_{n>=1} 1/a(n) = 112/75 - Pi/10. - Amiram Eldar, Aug 16 2022
EXAMPLE
The reduced fractions are 0, 3/5, 7/9, 6/7, 9/10, 25/27, 33/35, 21/22, 26/27, 63/65, 75/77, 44/45, ..
MAPLE
A178242 := proc(n) n*(5+n)/(n+1)/(n+4) ; numer(%) ; end proc:
seq(A178242(n), n=0..80) ; # R. J. Mathar, Dec 20 2010
MATHEMATICA
f[n_] := n/GCD[n, 4]; Array[f[#] f[# + 5] &, 50, 0]
Table[Numerator[n*(5+n)/((n+1)*(n+4))], {n, 0, 50}] (* G. C. Greubel, Sep 21 2018 *)
PROG
(Magma) [Floor(n*(n+5)*((-1)^((2*n-(-1)^n-3)/4)+3)/8) : n in [0..50]]; // Vincenzo Librandi, Oct 08 2011
(PARI) vector(50, n, n--; numerator(n*(5+n)/((n+1)*(n+4)))) \\ G. C. Greubel, Sep 21 2018
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
Paul Curtz, Dec 20 2010
STATUS
approved