OFFSET
1,1
COMMENTS
Sequence extrapolated to n=1 using recurrence. - Andrew Howroyd, May 10 2017
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000 (terms 1..200 from Andrew Howroyd)
Eric Weisstein's World of Mathematics, Dominating Set
Eric Weisstein's World of Mathematics, Prism Graph
Index entries for linear recurrences with constant coefficients, signature (3,1,5,1,1,-1,-1).
FORMULA
From Andrew Howroyd, May 10 2017: (Start)
a(n) = 3*a(n-1) + a(n-2) + 5*a(n-3) + a(n-4) + a(n-5) - a(n-6) - a(n-7).
G.f.: x*(-7*x^6 - 6*x^5 + 5*x^4 + 4*x^3 + 15*x^2 + 2*x + 3)/((x^2 + 1)*(x^5 + x^4 - 2*x^3 - 2*x^2 - 3*x + 1)). (End)
MATHEMATICA
LinearRecurrence[{3, 1, 5, 1, 1, -1, -1}, {3, 11, 51, 183, 663, 2435,
8935}, 20] (* Eric W. Weisstein, May 17 2017 *)
Rest[CoefficientList[Series[x (-7 x^6 - 6 x^5 + 5 x^4 + 4 x^3 + 15 x^2 + 2 x + 3)/((x^2 + 1) (x^5 + x^4 - 2 x^3 - 2 x^2 - 3 x + 1)), {x, 0, 20}], x]] (* G. C. Greubel, May 17 2017 *)
Table[2 Cos[n Pi/2] + RootSum[1 + #1 - 2 #1^2 - 2 #1^3 - 3 #1^4 + #1^5 &, #^n &], {n, 20}] (* Eric W. Weisstein, May 26 2017 *)
PROG
(PARI)
Vec((-7*x^6-6*x^5+5*x^4+4*x^3+15*x^2+2*x+3)/((x^2+1)*(x^5+x^4-2*x^3-2*x^2-3*x+1))+O(x^15)) \\ Andrew Howroyd, May 10 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Eric W. Weisstein, Apr 01 2017
EXTENSIONS
a(1)-a(2) and a(16)-a(25) from Andrew Howroyd, May 10 2017
STATUS
approved