A215007 -id:A215007 - cp's OEIS frontend

A215455 a(n) = 6a(n-1) - 9a(n-2) + a(n-3), with a(0)=3, a(1)=6 and a(2)=18.

3, 6, 18, 57, 186, 621, 2109, 7251, 25146, 87726, 307293, 1079370, 3798309, 13382817, 47191491, 166501902, 587670810, 2074699233, 7325660010, 25869337773, 91359785781, 322660334739, 1139593274178, 4024976418198, 14216179376325, 50211881768346, 177350652641349

Offset: 0

Roman Witula, Aug 11 2012

The Berndt-type sequence number 1 for the argument 2*Pi/9 (see also A215007, A215008) is connected with the following trigonometric identities: f(n;x)=g(n;x)=const for n=1,2 (and are equal to 6 and 18 respectively), f(n;x)+g(n;x)=const for n=3,4,5 (and are equal to 120, 420 and 1512 respectively). Moreover each of the functions f(3;x), g(3;x) and f(6;x)+g(6;x) is not the constant function. Here f(n;x) := (2*cos(x))^(2n) + (2*cos(x-Pi/3))^(2n) + (2*cos(x+Pi/3))^(2n), and g(n;x) := (2*sin(x))^(2n) + (2*cos(x-Pi/6))^(2n) + (2*cos(x+Pi/6))^(2n), for every n=1,2,..., and x in R (see Witula-Slota paper for details).

			From the identity c(j)^2 = 2 + c(2*j) we deduce that a(1)=6 is equivalent with c(2) + c(4) + c(8) = 0, where c(j) := 2*cos(Pi*j/9).

Andrew Howroyd, Table of n, a(n) for n = 0..500
R. Witula and D. Slota, On modified Chebyshev polynomials, J. Math. Anal. Appl., 324 (2006), 321-343.
Index entries for linear recurrences with constant coefficients, signature (6,-9,1).

Cf. A094831, A215007, A215008.

LinearRecurrence[{6,-9,1}, {3,6,18}, 50]

Vec((3-12*x+9*x^2)/(1-6*x+9*x^2-x^3)+O(x^99)) \\ Charles R Greathouse IV, Sep 27 2012

a(n) = c(1)^(2*n) + c(2)^(2*n) + c(4)^(2*n), where c(j) = 2*cos(Pi*j/9).
G.f.: 3*(1 - x)*(1 - 3*x)/(1 - 6*x + 9*x^2 - x^3).
a(n) = 3*A094831(n). - Andrew Howroyd, Apr 28 2020

Terms a(22) and beyond from Andrew Howroyd, Apr 28 2020

A215512 a(n) = 5a(n-1) - 6a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

Original entry on oeis.org

1, 3, 8, 23, 70, 220, 703, 2265, 7327, 23748, 77043, 250054, 811760, 2635519, 8557089, 27784091, 90213440, 292919743, 951102166, 3088205812, 10027335807, 32558546329, 105716922615, 343260670908, 1114560365179, 3618954723062, 11750672095144, 38154192502527

Offset: 0

Roman Witula, Aug 14 2012

The Berndt-type sequence number 7 for the argument 2Pi/7 defined by the relation: sqrt(7)*a(n) = s(1)*c(4)^(2*n) + s(2)*c(1)^(2*n) + s(4)*c(2)^(2*n), where c(j):=2*cos(2*Pi*j/7) and s(j):=2*sin(2*Pi*j/7). If we additionally defined the following sequences:
sqrt(7)*b(n) = s(2)*c(4)^(2*n) + s(4)*c(1)^(2*n) + s(1)*c(2)^(2*n),
sqrt(7)*c(n) = s(4)*c(4)^(2*n) + s(1)*c(1)^(2*n) + s(2)*c(2)^(2*n), and
sqrt(7)*a1(n) = s(1)*c(4)^(2*n+1) + s(2)*c(1)^(2*n+1) + s(4)*c(2)^(2*n+1),
sqrt(7)*b1(n) = s(2)*c(4)^(2*n+1) + s(4)*c(1)^(2*n+1) + s(1)*c(2)^(2*n+1),
sqrt(7)*c1(n) = s(4)*c(4)^(2*n+1) + s(1)*c(1)^(2*n+1) + s(2)*c(2)^(2*n+1), then the following simple relationships between elements of these sequences hold true: a(n)=c1(n), c(n+1)=a1(n), -a(n)-b(n)=b1(n), which means that the sequences a1(n), b1(n), and c1(n) are completely and in very simple way determined by the sequences a(n), b(n) and c(n). However the last one's satisfy the following system of recurrence equations: a(n+1) = 2*a(n) + b(n), b(n+1) = a(n) + 2*b(n) - c(n), c(n+1) = c(n) - b(n). We have b(n)=A215694(n) and c(n)=A215695(n).
We note that a(n)=A000782(n) for every n=0,1,...,4 and A000782(5)-a(5)=2.
From general recurrence relation: a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3), i.e. a(n) = 5*(a(n-1)-a(n-2)) + (a(n-3)-a(n-2)) the following summation formula can be easily obtained: sum{k=3,..,n} a(k) = 5*a(n-1)-a(n-2)+a(0)-5*a(1). Hence in discussed sequence it follows that: sum{k=3,..,n} a(k) = 5*a(n-1) - a(n-2) - 14.

			We have a(6) = 10*a(4)+a(1), a(5) = 11*(a(3)-a(1)), a(10)-a(4)+a(3)+a(1)+a(0) = 77*10^3, and a(11)-a(4)+a(3)-a(2)+a(0) = 25*10^4 = (5^6)*(2^4).

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6
Index entries for linear recurrences with constant coefficients, signature (5, -6, 1).

Cf.A215694, A215695, A215007, A215008, A215143, A215493, A215494, A215510, A215575, A215455, A214683, A214699.

I:=[1,3,8]; [n le 3 select I[n] else 5*Self(n-1) - 6*Self(n-2) + Self(n-3): n in [1..30]]; // G. C. Greubel, Apr 23 2018

Mathematica
```
LinearRecurrence[{5,-6,1}, {1,3,8}, 50]
```

x='x+O('x^30); Vec((1-2*x-x^2)/(1-5*x+6*x^2-x^3)) \\ G. C. Greubel, Apr 23 2018

G.f.: (1-2*x-x^2)/(1-5*x+6*x^2-x^3).

A094429 Given the 3 X 3 matrix M = [0 1 0 / 0 0 1 / 7 -14 7], a(n) = (-) rightmost term of M^n * [1 1 1].

Original entry on oeis.org

0, 7, 42, 196, 833, 3381, 13377, 52136, 201341, 773122, 2958032, 11291903, 43042727, 163918671, 623875840, 2373568575, 9028148962, 34334213564, 130560389505, 496440779373, 1887579497489, 7176808297736, 27286630574917

Offset: 1

Gary W. Adamson, May 02 2004

M is derived from the Lucas polynomial: x^3 - 7*x^2 + 14*x - 7 with a root (and eigenvalue of the matrix): 3.801377358... = (2*sin(3*Pi/7))^2, the convergent of the sequence.
From Roman Witula, Sep 29 2012: (Start)
The Berndt-type sequence number 16 for the argument 2*Pi/7 (see Formula section and Crossrefs for other Berndt-type sequences for the argument 2*Pi/7 - for numbers from 1 to 18 without 16).
Note that all numbers of the form a(n)*7^(-1 - floor((n-1)/3)) are integers and even a(10) and a(11) are divisible by 7^5. (End)

			a(5) = 833. M^5 * [1 1 1] = [ -42 -196 -833].
We have 4*a(4) - a(5) = 4*a(5) - a(6) = 7*a(2) = 49, 88*a(10) = 23*a(11), and a(3) = 6*a(2), which implies the equalities c(4)*(s(1))^6 + c(2)*(s(4))^6 + c(1)*(s(2))^6 = 6*(c(4)*(s(1))^4 + c(2)*(s(4))^4 + c(1)*(s(2))^4) and
s(2)*(s(1))^8 + s(4)*(s(2))^8 + s(1)*(s(4))^8 = 6*( s(2)*(s(1))^6 + s(4)*(s(2))^6 + s(1)*(s(4))^6). - _Roman Witula_, Sep 29 2012

G. C. Greubel, Table of n, a(n) for n = 1..1700
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094430, A217274.

I:=[0,7,42]; [n le 3 select I[n] else 7*Self(n-1) -14*Self(n-2) +7*Self(n-3): n in [1..30]]; // G. C. Greubel, May 09 2018

Table[(MatrixPower[{{0, 1, 0}, {0, 0, 1}, {7, -14, 7}}, n].{-1, -1, -1})[[3]], {n, 23}] (* Robert G. Wilson v, May 08 2004 *)
LinearRecurrence[{7,-14,7}, {0,7,42}, 50] (* Roman Witula, Aug 13 2012 *)

x='x+O('x^30); concat([0], Vec(7*x^2*(1-x)/(1-7*x+14*x^2-7*x^3))) \\ G. C. Greubel, May 09 2018

a(n) = -(([0, 1, 0; 0, 0, 1; 7, -14, 7]^n)*[1,1,1]~)[3]; \\ Michel Marcus, May 10 2018

From Colin Barker, Jun 19 2012: (Start)
a(n) = 7*a(n-1) - 14*a(n-2) + 7*a(n-3).
G.f.: 7*x^2*(1-x)/(1 - 7*x + 14*x^2 - 7*x^3). (End)
a(n) = c(4)*(s(1))^(2*n) + c(2)*(s(4))^(2*n) + c(1)*(s(2))^(2*n) = (-1/sqrt(7))*(c(1)*(s(1))^(2*n+3) + c(2)*(s(2))^(2*n+3) + c(3)*(s(3))^(2*n+3)) = (-1/sqrt(7))*(s(2)*(s(1))^(2*n+2) + s(4)*(s(2))^(2*n+2) + s(1)*(s(4))^(2*n+2)), where c(j) := 2*cos(2*Pi*j/7) and s(j) := 2*sin(2*Pi*j/7) (for the sums of the respective odd powers see A217274, see also A215493 and comments to A215494). For the proof of these formulas see Witula-Slota's paper. - Roman Witula, Jul 24 2012

More terms from Robert G. Wilson v, May 08 2004

A094430 a(n) is the rightmost term of M^n * [1 0 0], where M is the 3 X 3 matrix [0 1 0 / 0 0 1 / 7 -14 7].

Original entry on oeis.org

7, 49, 245, 1078, 4459, 17836, 69972, 271313, 1044435, 4002467, 15294370, 58337097, 222255768, 846131608, 3219700183, 12247849145, 46582062709, 177142452214, 673583231587, 2561162729076, 9737971026812, 37024601601729

Offset: 1

Gary W. Adamson, May 02 2004

In A094429 the multiplier is [1 1 1] instead of [1 0 0]. The matrix M is derived from the 3rd-order Lucas polynomial x^3 - 7x^2 + 14x - 7, with a convergent of the series = 3.801937735... = (2 sin 3*Pi/7)^2; (an eigenvalue of the matrix and a root of the polynomial).
From Roman Witula, Sep 29 2012: (Start)
This sequence is the Berndt-type sequence number 17 for the argument 2*Pi/7 (see Formula section and Crossrefs for other Berndt-type sequences for the argument 2*Pi/7 - for numbers from 1 to 18 without 17).
Note that all numbers of the form a(n)*7^(-floor((n+4)/3)) are integers. (End)

			a(4) = 1078 since M^4 * [1 0 0] = [49 245 1078] = [a(2), a(3), a(4)].
We have a(2)=7*a(1), a(3)=5*a(2), 22*a(3)=5*a(4), and a(6)=4*a(5), which implies s(2)*s(1)^15 + s(4)*s(2)^15 + s(1)*s(4)^15 = 4*(s(2)*s(1)^13 + s(4)*s(2)^13 + s(1)*s(4)^13). - _Roman Witula_, Sep 29 2012

G. C. Greubel, Table of n, a(n) for n = 1..1000
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094429, A217274.

I:=[49,245,1078]; [7] cat [n le 3 select I[n] else 7*Self(n-1) -14*Self(n-2) + 7*Self(n-3): n in [1..30]]; // G. C. Greubel, May 09 2018

Table[(MatrixPower[{{0, 1, 0}, {0, 0, 1}, {7, -14, 7}}, n].{1, 0, 0})[[3]], {n, 22}] (* Robert G. Wilson v, May 08 2004 *)
Join[{7}, LinearRecurrence[{7,-14,7}, {49,245,1078}, 50]] (* Roman Witula, Aug 13 2012 *)(* corrected by G. C. Greubel, May 09 2018 *)

x='x+O('x^30); Vec(7*x/(1-7*x+14*x^2-7*x^3)) \\ G. C. Greubel, May 09 2018

a(n) = (([0, 1, 0; 0, 0, 1; 7, -14, 7]^n)*[1,0,0]~)[3]; \\ Michel Marcus, May 10 2018

From Colin Barker, Jun 19 2012: (Start)
a(n) = 7*a(n-1)-14*a(n-2)+7*a(n-3).
G.f.: 7*x/(1-7*x+14*x^2-7*x^3). (End)
-a(n) = s(2)*s(1)^(2*n+3) + s(4)*s(2)^(2*n+3) + s(1)*s(4)^(2*n+3), where s(j) := 2*sin(2*Pi*j/7); for the proof see A215494 and the Witula-Slota paper. This formula and the respective recurrence also give a(0)=a(-1)=0. - Roman Witula, Aug 13 2012

More terms from Robert G. Wilson v, May 08 2004

Name edited by Michel Marcus, May 10 2018

A215665 a(n) = 3*a(n-2) - a(n-3), with a(0)=0, a(1)=a(2)=-3.

Original entry on oeis.org

0, -3, -3, -9, -6, -24, -9, -66, -3, -189, 57, -564, 360, -1749, 1644, -5607, 6681, -18465, 25650, -62076, 95415, -211878, 348321, -731049, 1256841, -2541468, 4501572, -8881245, 16046184, -31145307, 57019797, -109482105, 202204698, -385466112, 716096199

Offset: 0

Roman Witula, Aug 20 2012

The Berndt-type sequence number 6 for the argument 2Pi/9 defined by the first relation from the section "Formula" below. Two sequences connected with a(n) (possessing the respective numbers 5 and 7) are discussed in A215664 and A215666 - for more details see comments to A215664 and Witula's reference. We have a(n) - a(n+1) = A215664(n).
From initial values and the recurrence formula we deduce that a(n)/3 are all integers.
We note that a(10) is the first element of a(n) which is positive integer and all (-1)^n*a(n+10) are positive integer, which can be obtained from the title recurrence relation.
The following decomposition holds (X - c(1)*c(2)^n)*(X - c(2)*c(4)^n)*(X - c(4)*c(1)^n) = X^3 - a(n)*X^2 - A215917(n-1)*X + (-1)^n.
If X(n) = 3*X(n-2) - X(n-3), n in Z, with X(n) = a(n) for every n=0,1,..., then X(-n) = abs(A215919(n)) = (-1)^n*A215919(n) for every n=0,1,...

			We have a(1)=a(2)=a(8)=-3, a(3)=a(6)=-9, a(4)+a(11)=-10*a(10), and 47*a(5)=2*a(11).

R. Witula, Ramanujan type formulas for arguments 2Pi/7 and 2Pi/9, Demonstratio Math., (in press, 2012).
D. Chmiela and R. Witula, Two parametric quasi-Fibonacci numbers of the nine order, (submitted, 2012).

Index entries for linear recurrences with constant coefficients, signature (0,3,-1).

Cf. A215455, A215634, A215635, A215636, A215664, A214699, A215007, A214683.

LinearRecurrence[{0,3,-1}, {0,-3,-3}, 50]

concat(0,Vec(-3*(1+x)/(1-3*x^2+x^3)+O(x^99))) \\ Charles R Greathouse IV, Oct 01 2012

a(n) = c(1)*c(2)^n + c(2)*c(4)^n + c(4)*c(1)^n, where c(j) := 2*cos(2*Pi*j/9).

G.f.: -3*x*(1+x)/(1-3*x^2+x^3).

A215666 a(n) = 3*a(n-2) - a(n-3), with a(0)=0, a(1)=-3, and a(2)=6.

Original entry on oeis.org

0, -3, 6, -9, 21, -33, 72, -120, 249, -432, 867, -1545, 3033, -5502, 10644, -19539, 37434, -69261, 131841, -245217, 464784, -867492, 1639569, -3067260, 5786199, -10841349, 20425857, -38310246, 72118920, -135356595, 254667006, -478188705, 899357613

Offset: 0

Roman Witula, Aug 20 2012

The Berndt-type sequence number 7 for the argument 2Pi/9 defined by the first relation from the section "Formula" below. Two sequences connected with a(n) (possessing the respective numbers 5 and 6) are discussed in A215664 and A215665 - for more details see comments to A215664 and Witula's reference. We have a(n) = A215664(n+2) - 2*A215664(n) and a(n+1) = A215664(n+1) - A215664(n).
From initial values and the title recurrence formula we deduce that a(n)/3 and a(3*n)/9 are all integers.
If we set X(n) = 3*X(n-2) - X(n-3), n in Z, with a(n) = X(n), for every n=0,1,..., then X(-n) = -abs(A215917(n)) = (-1)^n*A215917(n), for every n=0,1,...

			We have 8*a(3)+a(6)=5*a(6)+3*a(7)=0, a(5) + a(12) = 3000, and (a(30)-1000*a(10)-a(2))/10^5 is an integer. Further we obtain  c(4)*cos(4*Pi/7)^7 + c(1)*cos(8*Pi/7)^7 + c(2)*c(2*Pi/7)^7 = -15/16.

R. Witula, Ramanujan type formulas for arguments 2Pi/7 and 2Pi/9, Demonstratio Math., (in press, 2012).
D. Chmiela and R. Witula, Two parametric quasi-Fibonacci numbers of the nine order, (submitted, 2012).

Index entries for linear recurrences with constant coefficients, signature (0,3,-1).

Cf. A215455, A215634, A215635, A215636, A215664, A214699, A215007, A214683.

LinearRecurrence[{0,3,-1}, {0,-3,6}, 50]

concat(0,Vec(-3*(1-2*x)/(1-3*x^2+x^3)+O(x^99))) \\ Charles R Greathouse IV, Oct 01 2012

a(n) = c(4)*c(2)^n + c(1)*c(4)^n + c(2)*c(1)^n, where c(j) := 2*cos(2*Pi*j/9).

G.f.: -3*x*(1-2*x)/(1-3*x^2+x^3).

A217274 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(0)=0, a(1)=1, a(2)=7.

Original entry on oeis.org

0, 1, 7, 35, 154, 637, 2548, 9996, 38759, 149205, 571781, 2184910, 8333871, 31750824, 120875944, 459957169, 1749692735, 6654580387, 25306064602, 96226175941, 365880389868, 1391138718116, 5289228800247, 20109822277181, 76457523763621, 290689756066542

Offset: 0

Roman Witula, Sep 29 2012

This is the Berndt-type sequence number 18 for the argument 2*Pi/7 defined by the relation
a(n)*sqrt(7) = c(4)*s(1)^(2n+1) + c(2)*s(4)^(2n+1) + c(1)*s(2)^(2n+1) = (1/s(4))*s(1)^(2n+2) + (1/s(2))*s(4)^(2n+2) + (1/s(1))*s(2)^(2n+2), where c(j) := 2*cos(2*Pi*j/7) and s(j) := 2*sin(2*Pi*j/7) (for the sums of the respective even powers see A094429). For the proof of this formula see the Witula/Slota and Witula references.
The definitions of the other Berndt-type sequences for the argument 2*Pi/7 (with numbers from 1 to 17) are in the cross references.
We note that all numbers of the form a(n)*7^(-floor((n+1)/3)) = A217444(n) are integers.
It can be proved that Sum_{k=2..n}a(k) = 7*(a(n-1) - a(n-2)).

			Writing c(j) as cj and s(k) as sk,
we have 7*sqrt(7) = c4*s1^5 + c2*s4^5 + c1*s2^5
and c4*s1^13 + c2*s4^13 + c1*s2^13 = 4(c4*s1^11 + c2*s4^11 + c1*s2^11).
We note that a(9) = 87*a(3)*a(2)^2 and a(11) = 2*a(3)*a(5)*a(2)^2.

G. C. Greubel, Table of n, a(n) for n = 0..1000
R. Witula, Ramanujan type trigonometric formulas, Demonstratio Math., Vol. XLV, No. 4, 2012, pp. 789-796.
Roman Witula and Damian Slota, New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A033304, A094429, A094430, A094648, A108716, A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A215794, A215817, A215828, A215877, A217444.

I:=[0,1,7]; [n le 3 select I[n] else 7*Self(n-1)-14*Self(n-2)+7*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jul 26 2015

LinearRecurrence[{7,-14,7}, {0,1,7}, 30]
CoefficientList[Series[x/(1 - 7*x + 14*x^2 - 7*x^3), {x,0,50}], x] (* G. C. Greubel, Apr 16 2017 *)

a[0]:0$
a[1]:1$
a[2]:7$
a[n]:=7*a[n-1] - 14*a[n-2] + 7*a[n-3];
makelist(a[n], n, 0, 25); /* Martin Ettl, Oct 11 2012 */

concat(0, Vec(x/(1-7*x+14*x^2-7*x^3) + O(x^40))) \\ Michel Marcus, Jul 25 2015

G.f.: x/(1-7*x+14*x^2-7*x^3).

A215635 a(n) = - 12a(n-1) - 54a(n-2) - 112a(n-3) - 105a(n-4) -36a(n-5) - 2a(n-6), with a(0)=3, a(1)=-6, a(2)=18, a(3)=-60, a(4)=210, a(5)=-756.

Original entry on oeis.org

3, -6, 18, -60, 210, -756, 2772, -10296, 38610, -145860, 554268, -2116296, 8112462, -31201644, 120347532, -465328200, 1803025410, -6999149124, 27213719148, -105960069864, 413078158350, -1612098272460, 6297409350492, -24620247483624, 96324799842498, -377102656201956, 1477141800784668

Offset: 0

Roman Witula, Aug 18 2012

The Berndt-type sequence number 3 for the argument 2*Pi/9 defined by the relation: X(n) = a(n) + b(n)*sqrt(2), where X(n) := ((cos(Pi/24))^(2*n) + (cos(7*Pi/24))^(2*n) + (cos(3*Pi/8))^(2*n))*(-4)^n. We have b(n) = A215636(n).
We note that above formula is the Binet form of the following recurrence sequence: X(n+3) + 6*X(n+2) + 9*X(n+1) + (2 + sqrt(2))*X(n) = 0, which is a special type of the sequence X(n)=X(n;g) defined in the comments to A215634 for g:=Pi/24. The sequences a(n) and b(n) satisfy the following system of recurrence equations: a(n) = -b(n+3)-6*b(n+2)-9*b(n+1)-2*b(n), 2*b(n) =  -a(n+3)-6*a(n+2)-9*a(n+1)-2*a(n).
There exists an amazing relation: (-1)^n*a(n)=3*A000984(n) for every n=0,1,...,11 and 3*A000984(12)-a(12)=6.

Roman Witula and D. Slota, On modified Chebyshev polynomials, J. Math. Anal. Appl., 324 (2006), 321-343.
Index entries for linear recurrences with constant coefficients, signature (-12,-54,-112,-105,-36,-2).

Cf. A215455, A215634, A215636, A215007, A214699, A214683.

LinearRecurrence[{-12,-54,-112,-105,-36,-2}, {3,-6,18,-60,210,-756}, 50]

Vec((3+30*x+108*x^2+168*x^3+105*x^4+18*x^5) /(1+12*x+54*x^2+112*x^3+105*x^4+36*x^5+2*x^6)+O(x^99)) \\ Charles R Greathouse IV, Oct 01 2012

G.f.: (3+30*x+108*x^2+168*x^3+105*x^4+18*x^5) / (1+12*x+54*x^2+112*x^3+105*x^4+36*x^5+2*x^6).

A122068 Expansion of x(1-3x)(1-x)/(1-7x+14x^2-7x^3).

Original entry on oeis.org

1, 3, 10, 35, 126, 462, 1715, 6419, 24157, 91238, 345401, 1309574, 4970070, 18874261, 71705865, 272491891, 1035680954, 3936821259, 14965658694, 56893879910, 216295686467, 822315097387, 3126323230541, 11885921055638

Offset: 1

Gary W. Adamson, Oct 15 2006

G. C. Greubel, Table of n, a(n) for n = 1..1000
Peter Steinbach, Golden fields: a case for the heptagon, Math. Mag. Vol. 70, No. 1, Feb. 1997, 22-31.
R. Witula, P. Lorenc, M. Rozanski, and M. Szweda, Sums of the rational powers of roots of cubic polynomials, Zeszyty Naukowe Politechniki Slaskiej, Seria: Matematyka Stosowana z. 4, Nr. kol. 1920, 2014.
Index entries for linear recurrences with constant coefficients, signature (7,-14,7).

Cf. A087946, A081567.

Cf. A215007, A215008. - Roman Witula, May 16 2014

a:=[1,3,10];; for n in [4..30] do a[n]:=7*(a[n-1]-2*a[n-2]+a[n-3]); od; a; # G. C. Greubel, Oct 03 2019

I:=[1,3,10]; [n le 3 select I[n] else 7*(Self(n-1) -2*Self(n-2) + Self(n-3)): n in [1..30]]; // G. C. Greubel, Oct 03 2019

seq(coeff(series(x*(1-3*x)*(1-x)/(1-7*x+14*x^2-7*x^3), x, n+1), x, n), n =1..30); # G. C. Greubel, Oct 03 2019

M = {{2,1,0,0,0,0}, {1,2,1,0,0,0}, {0,1,2,1,0,0}, {0,0,1,2,1,0}, {0,0,0, 1,2,1}, {0,0,0,0,1,2}}; v[1] = {1,1,1,1,1,1}; v[n_]:= v[n] = M.v[n-1]; Table[v[n][[1]], {n,30}]
Rest@CoefficientList[Series[x*(1-3*x)*(1-x)/(1-7*x+14*x^2-7*x^3), {x, 0, 30}], x] (* G. C. Greubel, Oct 03 2019 *)
LinearRecurrence[{7,-14,7},{1,3,10},30] (* Harvey P. Dale, Mar 08 2020 *)

Vec(x*(1-3*x)*(1-x)/(1-7*x+14*x^2-7*x^3)+O(x^30)) \\ Charles R Greathouse IV, Sep 27 2012

def A122068_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(x*(1-3*x)*(1-x)/(1-7*x+14*x^2-7*x^3)).list()
a=A122068_list(30); a[1:] # G. C. Greubel, Oct 03 2019

From Roman Witula, May 16 2014: (Start)
a(n) = (1/2)*Sum_{k=0..2}(1 - 1/sqrt(7)*cot(2^k * alpha))* (2*sin(2^k * alpha))^(2n), where alpha := 2*Pi/7.
a(n) = (A215007(n) + A215008(n+1) - 2*A215008(n))/2. (End)
a(n) = binomial(2*n-1, n-1) + Sum_{k=1..n} (-1)^k*binomial(2*n, n+7*k). - Greg Dresden, Jan 28 2023

A217444 a(n) = A(n)7^(-floor(n+1)/3), where A(n) = 7A(n-1) - 14A(n-2) + 7A(n-3) with A(0)=0, A(1)=1, A(2)=7.

Original entry on oeis.org

0, 1, 1, 5, 22, 13, 52, 204, 113, 435, 1667, 910, 3471, 13224, 7192, 27367, 104105, 56563, 215098, 817909, 444276, 1689212, 6422529, 3488381, 13262821, 50424942, 27387681, 104126704, 395884336, 215018609, 817488295, 3108041875, 1688083894, 6417991803, 24400809980

Offset: 0

Roman Witula, Oct 03 2012

The Berndt-type sequence number 18a for the argument 2Pi/7, which is closely connected with the sequence A217274. Definitions other Berndt-type sequences for the argument 2Pi/7 like A215575, A215877, A033304 in sequences from Crossrefs are given.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,10,0,0,-17,0,0,1).

Cf. A215007, A215008, A215143, A215493, A215494, A215510, A215512, A215575, A215694, A215695, A108716, A215794, A215828, A215817, A215877, A094429, A094430, A217274, A094648, A033304.

i:=35; I:=[0, 1, 7]; A:=[m le 3 select I[m] else 7*Self(m-1)-14*Self(m-2)+7*Self(m-3): m in [1..i]]; [7^(-Floor(n/3))*A[n]: n in [1..i]]; // Bruno Berselli, Oct 03 2012

CoefficientList[Series[x*(1+x+5*x^2+12*x^3+3*x^4+2*x^5+x^6)/(1 - 10*x^3 + 17*x^6 - x^9), {x, 0, 40}], x] (* Vincenzo Librandi, Dec 15 2012 *)
LinearRecurrence[{0,0,10,0,0,-17,0,0,1}, {0, 1, 1, 5, 22, 13, 52, 204, 113}, 50] (* G. C. Greubel, Apr 23 2018 *)

x='x+O('x^30); concat([0], Vec(x*(1+x+5*x^2+12*x^3+3*x^4 +2*x^5 +x^6)/(1- 10*x^3+17*x^6-x^9))) \\ G. C. Greubel, Apr 23 2018

G.f.: x*(1+x+5*x^2+12*x^3+3*x^4+2*x^5+x^6)/(1-10*x^3+17*x^6-x^9). - Bruno Berselli, Oct 03 2012

a(n) = 10*a(n-3) - 17*a(n-6) + a(n-9). - G. C. Greubel, Apr 23 2018

cp's OEIS Frontend

A215455 a(n) = 6*a(n-1) - 9*a(n-2) + a(n-3), with a(0)=3, a(1)=6 and a(2)=18.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A215512 a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A094429 Given the 3 X 3 matrix M = [0 1 0 / 0 0 1 / 7 -14 7], a(n) = (-) rightmost term of M^n * [1 1 1].

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

Extensions

A094430 a(n) is the rightmost term of M^n * [1 0 0], where M is the 3 X 3 matrix [0 1 0 / 0 0 1 / 7 -14 7].

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

Extensions

A215665 a(n) = 3*a(n-2) - a(n-3), with a(0)=0, a(1)=a(2)=-3.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A215666 a(n) = 3*a(n-2) - a(n-3), with a(0)=0, a(1)=-3, and a(2)=6.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

A215455 a(n) = 6a(n-1) - 9a(n-2) + a(n-3), with a(0)=3, a(1)=6 and a(2)=18.

A215512 a(n) = 5a(n-1) - 6a(n-2) + a(n-3), with a(0)=1, a(1)=3, a(2)=8.

A217274 a(n) = 7a(n-1) - 14a(n-2) + 7*a(n-3) with a(0)=0, a(1)=1, a(2)=7.

A215635 a(n) = - 12a(n-1) - 54a(n-2) - 112a(n-3) - 105a(n-4) -36a(n-5) - 2a(n-6), with a(0)=3, a(1)=-6, a(2)=18, a(3)=-60, a(4)=210, a(5)=-756.

A122068 Expansion of x(1-3x)(1-x)/(1-7x+14x^2-7x^3).

A217444 a(n) = A(n)7^(-floor(n+1)/3), where A(n) = 7A(n-1) - 14A(n-2) + 7A(n-3) with A(0)=0, A(1)=1, A(2)=7.