|
|
A074089
|
|
Coefficient of q^3 in nu(n), where nu(0)=1, nu(1)=b and, for n >= 2, nu(n) = b*nu(n-1) + lambda*(1+q+q^2+...+q^(n-2))*nu(n-2) with (b,lambda)=(2,3).
|
|
20
|
|
|
0, 0, 0, 0, 0, 78, 501, 2574, 11757, 50034, 203229, 797316, 3046362, 11394774, 41885913, 151732722, 542840175, 1921208586, 6735519249, 23417342568, 80810560596, 277008392478, 943826398893, 3198199361910, 10783017814065
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
The coefficient of q^0 is A014983(n+1).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (78*x^5 -123*x^6 -498*x^7 +297*x^8 +1134*x^9 +567*x^10)/(1 -2*x -3*x^2)^4.
a(n) = 8*a(n-1) -12*a(n-2) -40*a(n-3) +74*a(n-4) +120*a(n-5) -108*a(n-6) -216*a(n-7) -81*a(n-8) for n>=11.
|
|
EXAMPLE
|
The first 6 nu polynomials are nu(0)=1, nu(1)=2, nu(2)=7, nu(3) = 20 + 6q, nu(4) = 61 + 33q + 21q^2, nu(5) = 182 + 144q + 120q^2 + 78q^3 + 18q^4, so the coefficients of q^3 are 0,0,0,0,0,78.
|
|
MATHEMATICA
|
b=2; lambda=3; expon=3; nu[0]=1; nu[1]=b; nu[n_] := nu[n]= Together[ b*nu[n-1]+lambda(1-q^(n-1))/(1-q)nu[n-2]]; a[n_] := Coefficient[nu[n], q, expon]
(* Second program: *)
CoefficientList[Series[(78*x^5-123*x^6-498*x^7+297*x^8+1134*x^9 + 567*x^10)/( 1-2*x-3*x^2)^4, {x, 0, 50}], x] (* G. C. Greubel, May 26 2018 *)
|
|
PROG
|
(PARI) x='x+O('x^30); concat([0, 0, 0, 0, 0], Vec((78*x^5 -123*x^6 -498*x^7 +297*x^8 +1134*x^9 +567*x^10)/(1 -2*x -3*x^2)^4)) \\ G. C. Greubel, May 26 2018
(Magma) m:=25; R<x>:=PowerSeriesRing(Integers(), m); [0, 0, 0, 0, 0] cat Coefficients(R!((78*x^5 -123*x^6 -498*x^7 +297*x^8 +1134*x^9 +567*x^10)/(1 -2*x -3*x^2)^4)); // G. C. Greubel, May 26 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 19 2002
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|