OFFSET
0,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,1,-1).
FORMULA
G.f.: 5*(197 + 10*x + 10*x^2 - 182*x^3)/((1-x)^2*(1+x+x^2)). - R. J. Mathar, Nov 14 2007
From G. C. Greubel, Sep 08 2025: (Start)
a(n) = (5/3)*(35*n + 586 + 5*A057078(n)).
E.g.f.: (5/3)*( 5*exp(-x/2)*( cos((sqrt(3)*x)/2) + (1/sqrt(3))*sin((sqrt(3)*x)/2)) + (586 + 35*x)*exp(x) ). (End)
MAPLE
A131640 := proc(n) option remember ; if n =0 then 985 ; elif n mod 3 = 0 then A131640(n-1)+75 ; else A131640(n-1)+50 ; fi ; end: seq(A131640(n), n=0..80) ; # R. J. Mathar, Oct 24 2007
MATHEMATICA
LinearRecurrence[{1, 0, 1, -1}, {985, 1035, 1085, 1160}, 50] (* Ray Chandler, Aug 25 2015 *)
Table[5*(35*n +591 -5*Mod[n, 3])/3, {n, 0, 50}] (* G. C. Greubel, Sep 08 2025 *)
PROG
(PARI) Vec(5*(197+10*x+10*x^2-182*x^3)/((1-x)^2*(1+x+x^2)) + O(x^40)) \\ Andrew Howroyd, Feb 20 2018
(Magma)
A131640:= func< n | (5/3)*(35*n + 591 - 5*(n mod 3) ) >;
[A131640(n): n in [0..50]]; // G. C. Greubel, Sep 08 2025
(SageMath)
def A131640(n): return 5*(35*n + 591 - 5*(n%3))//3
print([A131640(n) for n in range(51)]) # G. C. Greubel, Sep 08 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Eric M. Adler (eadler(AT)simi.k12.ca.us), Sep 05 2007
EXTENSIONS
Definition supplied by N. J. A. Sloane, Sep 14 2007
More terms from R. J. Mathar, Oct 24 2007
STATUS
approved
