OFFSET
1,1
COMMENTS
Sequence extrapolated to n=1 using recurrence. - Andrew Howroyd, May 10 2017
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
Eric Weisstein's World of Mathematics, Dominating Set
Eric Weisstein's World of Mathematics, Moebius Ladder
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) for n>7.
G.f.: x*(1-x)*(1+x)*(3*x^4+2*x^3+6*x^2+6*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, 15, 51, 179, 663, 2439,
8935}, 20] (* _Eric W. Weisstein, May 17 2017 *)
Rest[CoefficientList[Series[x*(1 - x)*(1 + x)*(3*x^4 + 2*x^3 + 6*x^2 + 6*x + 3)/((x^2 + 1)*(x^5 + x^4 - 2*x^3 - 2*x^2 - 3*x + 1)), {x, 0, 50}], x]] (* G. C. Greubel, May 17 2017 *)
Table[RootSum[1 + # - 2 #^2 - 2 #^3 - 3 #^4 + #^5 &, #^n &] - 2 Cos[(n Pi)/2], {n, 20}] (* _Eric W. Weisstein, Jun 14 2017 *)
PROG
(PARI)
Vec((1-x)*(1+x)*(3*x^4+2*x^3+6*x^2+6*x+3)/((x^2+1)*(x^5+x^4-2*x^3-2*x^2-3*x+1))+O(x^50)) \\ Andrew Howroyd, May 10 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Mar 31 2017
EXTENSIONS
a(1)-(2) and a(16)-a(25) from Andrew Howroyd, May 10 2017
STATUS
approved