OFFSET
2,2
COMMENTS
Bisection of A053298.
LINKS
Colin Barker, Table of n, a(n) for n = 2..1000
Index entries for linear recurrences with constant coefficients, signature (8,-26,45,-45,26,-8,1).
FORMULA
G.f.: x^2*(1+26*x+35*x^2-12*x^3-2*x^4) / ((1-x)^5*(1-3*x+x^2)). - Colin Barker, Feb 19 2016
From G. C. Greubel, Jul 22 2025: (Start)
a(n) = Lucas(2*n+7) - (12*n^4 + 20*n^3 + 81*n^2 + 169*n + 174)/6.
E.g.f.: exp(3*x/2)*(29*cosh(p*x) + 26*p*sinh(p*x)) - (1/6)*(174 + 282*x + 225*x^2 + 92*x^3 + 12*x^4)*exp(x), where 2*p = sqrt(5). (End)
MATHEMATICA
A027006[n_]:= LucasL[2*n+7] -(12*n^4 +20*n^3 +81*n^2 +169*n +174)/6;
Table[A027006[n], {n, 2, 42}] (* G. C. Greubel, Jul 22 2025 *)
PROG
(PARI) Vec(x^2*(1+26*x+35*x^2-12*x^3-2*x^4) / ((1-x)^5*(1-3*x+x^2)) + O(x^40)) \\ Colin Barker, Feb 19 2016
(Magma)
A027006:= func< n | Lucas(2*n+7) -(12*n^4 +20*n^3 +81*n^2 +169*n +174)/6 >;
[A027006(n): n in [2..40]]; // G. C. Greubel, Jul 22 2025
(SageMath)
def A027006(n): return lucas_number2(2*n+7, 1, -1) -(12*n^4 +20*n^3 +81*n^2 +169*n +174)//6
print([A027006(n) for n in range(2, 41)]) # G. C. Greubel, Jul 22 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
