|
|
A106803
|
|
Expansion of x*(1-x)/(1-2*x-x^2+x^3).
|
|
5
|
|
|
0, 1, 1, 3, 6, 14, 31, 70, 157, 353, 793, 1782, 4004, 8997, 20216, 45425, 102069, 229347, 515338, 1157954, 2601899, 5846414, 13136773, 29518061, 66326481, 149034250, 334876920, 752461609, 1690765888, 3799116465, 8536537209
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
COMMENTS
|
a(n) appears in the formula for the nonnegative powers of sigma, the ratio of the smaller diagonal in the heptagon to the side length s=2*sin(Pi/7), when expressed in the basis <1,rho,sigma>, with rho = 2*cos(Pi/7), the ratio of the smaller heptagon diagonal to the side length, as follows. sigma^n = a(n-1)*1 + B(n)*rho + a(n)*sigma, n>=0, with B(n)=A006054(n). Put a(-1):= 1. See the Steinbach reference, and a comment under A052547.
a(n-1) is the top left entry of the n-th power of the 3X3 matrix [0, 1, 0; 1, 1, 1; 0, 1, 1] or of the 3X3 matrix [0, 0, 1; 0, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 2*a(n-1) + a(n-2) - a(n-3), a(0)=0, a(1)=a(2)=1. - G. C. Greubel, Aug 14 2015
|
|
MATHEMATICA
|
m = {{0, 0, 1}, {1, 2, 0}, {1, 1, 0}}; v[0] = {0, 1, 1}; v[n_] := m.v[n - 1]; Table[v[n][[1]], {n, 0, 30}] (* Edited and corrected by L. Edson Jeffery, Oct 18 2017 *)
RecurrenceTable[{a[1]== 0, a[2]== 1, a[3]== 1, a[n]== 2*a[n-1] + a[n-2] - a[n-3]}, a, {n, 30}] (* G. C. Greubel, Aug 14 2015 *)
|
|
PROG
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,less
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|