OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
Index entries for linear recurrences with constant coefficients, signature (15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,-120).
FORMULA
G.f.: (t^16 + 2*t^15 + 2*t^14 + 2*t^13 + 2*t^12 + 2*t^11 + 2*t^10 + 2*t^9 + 2*t^8 + 2*t^7 + 2*t^6 + 2*t^5 + 2*t^4 + 2*t^3 + 2*t^2 + 2*t + 1)/( 120*t^16 - 15*t^15 - 15*t^14 - 15*t^13 - 15*t^12 - 15*t^11 - 15*t^10 - 15*t^9 - 15*t^8 - 15*t^7 - 15*t^6 - 15*t^5 - 15*t^4 - 15*t^3 - 15*t^2 - 15*t + 1).
From G. C. Greubel, Sep 10 2023: (Start)
G.f.: (1+t)*(1-t^16)/(1 - 16*t + 135*t^16 - 120*t^17).
a(n) = 15*Sum_{j=1..15} a(n-j) - 120*a(n-16). (End)
MATHEMATICA
coxG[{16, 120, -15}] (* The coxG program is at A169452 *) (* Harvey P. Dale, Nov 27 2015 *)
CoefficientList[Series[(1+t)*(1-t^16)/(1-16*t+135*t^16-120*t^17), {t, 0, 50}], t] (* G. C. Greubel, Jul 01 2016; Sep 10 2023 *)
PROG
(Magma) R<x>:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+x)*(1-x^16)/(1-16*x+135*x^16-120*x^17) )); // G. C. Greubel, Sep 10 2023
(SageMath)
def A167926_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1+x)*(1-x^16)/(1-16*x+135*x^16-120*x^17) ).list()
A167926_list(40) # G. C. Greubel, Sep 10 2023
(PARI) Vec((1+x^2)*(1+x^4)*(1+x^8)*(1+x)^2/(1-15*x-15*x^2-15*x^3-15*x^4-15*x^5-15*x^6-15*x^7-15*x^8-15*x^9-15*x^10-15*x^11-15*x^12-15*x^13-15*x^14-15*x^15+120*x^16)+O(x^99)) \\ Charles R Greathouse IV, May 16 2026
CROSSREFS
KEYWORD
nonn,easy,changed
AUTHOR
John Cannon and N. J. A. Sloane, Dec 03 2009
STATUS
approved
