A215076 -id:A215076 - cp's OEIS frontend

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

0, 0, 1, 1, 3, 4, 9, 14, 28, 47, 89, 155, 286, 507, 924, 1652, 2993, 5373, 9707, 17460, 31501, 56714, 102256, 184183, 331981, 598091, 1077870, 1942071, 3499720, 6305992, 11363361, 20475625, 36896355, 66484244, 119801329, 215873462, 388991876, 700937471

Offset: 0

N. J. A. Sloane

a(n+1) = S(n) for n>=1, where S(n) is the number of 01-words of length n, having first letter 1, in which all runlengths of 1's are odd. Example: S(4) counts 1000, 1001, 1010, 1110. See A077865. - Clark Kimberling, Jun 26 2004
For n>=1, number of compositions of n into floor(j/2) kinds of j's (see g.f.). - Joerg Arndt, Jul 06 2011
Counts walks of length n between the first and second nodes of P_3, to which a loop has been added at the end. Let A be the adjacency matrix of the graph P_3 with a loop added at the end. A is a 'reverse Jordan matrix' [0,0,1; 0,1,1; 1,1,0]. a(n) is obtained by taking the (1,2) element of A^n. - Paul Barry, Jul 16 2004
Interleaves A094790 and A094789. - Paul Barry, Oct 30 2004
a(n) appears in the formula for the nonnegative powers of rho:= 2*cos(Pi/7), 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 sigma:=rho^2-1, the ratio of the larger heptagon diagonal to the side length, as follows. rho^n = C(n)*1 + C(n+1)*rho + a(n)*sigma, n>=0, with C(n) = A052547(n-2). See the Steinbach reference, and a comment under A052547. - Wolfdieter Lang, Nov 25 2010
If with the above notations the power basis <1,rho,rho^2> of Q(rho) is used, nonnegative powers of rho are given by rho^n  = -a(n-1)*1 + A052547(n-1)*rho + a(n)*rho^2. For negative powers see A006054. - Wolfdieter Lang, May 06 2011
-a(n-1) also appears in the formula for the nonpositive powers of sigma (see the above comment for the definition, and the Steinbach basis <1,rho,sigma>) as follows: sigma^(-n) = A(n)*1 -a(n+1)*rho -A(n-1)*sigma, with A(n) = A052547(n), A(-1):=0. - Wolfdieter Lang, Nov 25 2010

			G.f. = x^2 + x^3 + 3*x^4 + 4*x^5 + 9*x^6 + 14*x^7 + 28*x^8 + 47*x^9 + ...
Regarding the description "number of compositions of n into floor(j/2) kinds of j's," the a(6)=9 compositions of 6 are (2a, 2a, 2a), (3a, 3a), (2a, 4a), (2a, 4b), (4a, 2a), (4b, 2a), (6a), (6b), (6c). - _Bridget Tenner_, Feb 25 2022

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
R. Witula, E. Hetmaniok and D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the 15th International Conference on Fibonacci Numbers and Their Applications (2012).

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Robin Chapman and Nicholas C. Singer, Eigenvalues of a bidiagonal matrix, Amer. Math. Monthly, 111 (2004), p. 441.
Nachum Dershowitz, Between Broadway and the Hudson: A Bijection of Corridor Paths, arXiv:2006.06516 [math.CO], 2020.
Tomislav Došlić, Mate Puljiz, Stjepan Šebek, and Josip Žubrinić, On a variant of Flory model, arXiv:2210.12411 [math.CO], 2022.
Jia Huang, A coin flip game and generalizations of Fibonacci numbers, arXiv:2501.07463 [math.CO], 2025. See p. 9.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 433
Mohammed L. Nadji, Moussa Ahmia, Daniel F. Checa, and José L. Ramírez, Arndt Compositions with Restricted Parts, Palindromes, and Colored Variants, J. Int. Seq. (2025) Vol. 28, Issue 3, Article 25.3.6. See p. 11.
László Németh and Dragan Stevanović, Graph solution of system of recurrence equations, Research Gate, 2023. See Table 2 at p. 6.
László Németh and László Szalay, Sequences Involving Square Zig-Zag Shapes, J. Int. Seq., Vol. 24 (2021), Article 21.5.2.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
R. Sachdeva and A. K. Agarwal, Combinatorics of certain restricted n-color composition functions, Discrete Mathematics, 340, (2017), 361-372.
Genki Shibukawa, New identities for some symmetric polynomials and their applications, arXiv:1907.00334 [math.CA], 2019.
P. Steinbach, Golden fields: a case for the heptagon, Math. Mag. 70 (1997), no. 1, 22-31.
Kai Wang, Fibonacci Numbers And Trigonometric Functions Outline, (2019).
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
Index entries for linear recurrences with constant coefficients, signature (1,2,-1).

Cf. A006054, A052547, A077865, A094789, A094790, A096975, A096976.

Cf. A120757, A187065, A187066, A187067, A214683, A215076, A215100.

a006053 n = a006053_list !! n
a006053_list = 0 : 0 : 1 : zipWith (+) (drop 2 a006053_list)
   (zipWith (-) (map (2 *) $ tail a006053_list) a006053_list)
-- Reinhard Zumkeller, Oct 14 2011

[ n eq 1 select 0 else n eq 2 select 0 else n eq 3 select 1 else Self(n-1) +2*Self(n-2) -Self(n-3): n in [1..40] ]; // Vincenzo Librandi, Aug 19 2011

a[0]:=0: a[1]:=0: a[2]:=1: for n from 3 to 40 do a[n]:=a[n-1]+2*a[n-2]-a[n-3] od:seq(a[n], n=0..40); # Emeric Deutsch
A006053:=z**2/(1-z-2*z**2+z**3); # conjectured by Simon Plouffe in his 1992 dissertation

LinearRecurrence[{1,2,-1}, {0,0,1}, 50]  (* Vladimir Joseph Stephan Orlovsky, May 25 2011 *)

{a(n) = if( n<0, n = -1-n; polcoeff( -1 / (1 - 2*x - x^2 + x^3) + x * O(x^n), n), polcoeff( x^2 / (1 - x - 2*x^2 + x^3) + x * O(x^n), n))}; /* Michael Somos, Nov 30 2014 */

@CachedFunction
def a(n): # a = A006053
    if (n<3): return (n//2)
    else: return a(n-1) + 2*a(n-2) - a(n-3)
[a(n) for n in range(41)] # G. C. Greubel, Feb 12 2023

G.f.: x^2/(1 - x - 2*x^2 + x^3). - Emeric Deutsch, Dec 14 2004
a(n) = c^(n-2) - a(n-1)*(c-1) + (1/c)*a(n-2) for n > 3 where c = 2*cos(Pi/7). Example: a(7) = 14 = c^5 - 9*(c-1) + 4/c = 18.997607... - 7.21743962... + 2.219832528... - Gary W. Adamson, Jan 24 2010
G.f.: -1 + 1/(1 - Sum_{j>=1} floor(j/2)*x^j). - Joerg Arndt, Jul 06 2011
a(n+2) = A094790(n/2+1)*(1+(-1)^n)/2 + A094789((n+1)/2)*(1-(-1)^n)/2. - Paul Barry, Oct 30 2004
First differences of A028495. - Floor van Lamoen, Nov 02 2005
a(n) = A187065(2*n+1); a(n+1) = A187066(2*n+1) = A187067(2*n). - L. Edson Jeffery, Mar 16 2011
a(n) = 2^n*(c(1)^(n-1)*(c(1)+c(2)) + c(3)^(n-1)*(c(3)+c(6)) + c(5)^(n-1)*(c(5)+c(4)) )/7, with c(j):=cos(Pi*j/7). - Herbert Kociemba, Dec 18 2011
a(n+1)*(-1)^n*49^(1/3) = (c(1)/c(4))^(1/3)*(2*c(1))^n + (c(2)/c(1))^(1/3)*(2*c(2))^n + (c(4)/c(2))^(1/3)*(2c(4))^n = (c(2)/c(1))^(1/3)*(2*c(1))^(n+1) + (c(4)/c(2))^(1/3)*(c(2))^(n+1) + (c(1)/c(4))^(1/3)*(2*c(4))^(n+1), where c(j) := cos(2Pi*j/7); for the proof, see Witula et al.'s papers. - Roman Witula, Jul 21 2012
The previous formula connects the sequence a(n) with A214683, A215076, A215100, A120757. We may call a(n) the Ramanujan-type sequence number 2 for the argument 2*Pi/7. - Roman Witula, Aug 02 2012
a(n) = -A006054(1-n) for all n in Z. - Michael Somos, Nov 30 2014
G.f.: x^2 / (1 - x / (1 - 2*x / (1 + 5*x / (2 - x / (5 - 2*x))))). - Michael Somos, Jan 20 2017
a(n) ~ r*c^n, where r=0.241717... is one of the roots of 49*x^3-7*x+1, and c=2*cos(Pi/7) (as in Gary W. Adamson's formula). - Daniel Checa, Nov 04 2022
a(2n-1) = 2*a(n+1)*a(n) - a(n)^2 - a(n-1)^2. - Richard Peterson, May 25 2023

More terms from Emeric Deutsch, Dec 14 2004

Typo in definition fixed by Reinhard Zumkeller, Oct 14 2011

A120757 Expansion of x^2(2+x)/(1-3x-4*x^2-x^3).

Original entry on oeis.org

0, 2, 7, 29, 117, 474, 1919, 7770, 31460, 127379, 515747, 2088217, 8455018, 34233669, 138609296, 561217582, 2272323599, 9200450421, 37251863241, 150829715006, 610697048403, 2472661868474, 10011603514040, 40536155064419

Offset: 1

Gary W. Adamson & Roger L. Bagula, Jul 01 2006

The (1,1)-entry of the matrix M^n, where M is the 3 X 3 matrix [0,1,1; 1,1,2; 1,2,2].
a(n)/a(n-1) tends to 4.0489173...an eigenvalue of M and a root to the characteristic polynomial x^3 - 3x^2 - 4x - 1.
C(n):=a(n), with a(0):=1 (hence the o.g.f. for C(n) is (1-3*x-2*x^2)/(1-3*x-4*x^2-x^3)), appears in the following formula for the nonnegative powers of rho*sigma, where rho:=2*cos(Pi/7) and sigma:=sin(3*Pi/7)/sin(Pi/7) = rho^2-1 are the ratios of the smaller and larger diagonal length to the side length in a regular 7-gon (heptagon). See the Steinbach reference where the basis <1,rho,sigma> is used in an extension of the rational field. (rho*sigma)^n = C(n) + B(n)*rho + A(n)*sigma,n>=0, with B(n)= |A122600(n-1)|, B(0)=0, and A(n)= A181879(n). For the nonpositive powers see A085810(n)*(-1)^n, A181880(n-2)*(-1)^n and A116423(n+1)*(-1)^(n+1), respectively. See also a comment under A052547.
We have a(n)=cs(3n+1), where the sequence cs(n) and its two conjugate sequences as(n) and bs(n) are defined in the comments to the sequence A214683 (see also A215076, A215100, A006053). We call the sequence a(n) the Ramanujan-type sequence number 5 for the argument 2Pi/7. Since as(3n+1)=bs(3n+1)=0, we obtain the following relation: 49^(1/3)*a(n) = (c(1)/c(4))^(n + 1/3) + (c(4)/c(2))^(n + 1/3) + (c(2)/c(1))^(n + 1/3), where c(j) := Cos(2Pi/7) (for more details and proofs see Witula et al.'s papers). - Roman Witula, Aug 02 2012

			a(7)=1919 because M^7= [1919,3458,4312;3458,6231,7770;4312,7770,9689].

R. Witula, E. Hetmaniok and D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 1..1000
P. Steinbach, Golden fields: a case for the heptagon, Math. Mag. 70 (1997), no. 1, 22-31, MR 1439165
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,4,1).

Cf. A006053, A214683, A215076, A215100.

a:=[0,2,7]; [ n le 3 select a[n] else 3*Self(n-1) + 4*Self(n-2) + Self(n-3): n in [1..25]]; // Marius A. Burtea, Oct 03 2019

with(linalg): M[1]:=matrix(3,3,[0,1,1,1,1,2,1,2,2]): for n from 2 to 25 do M[n]:=multiply(M[1],M[n-1]) od: seq(M[n][1,1],n=1..25);

LinearRecurrence[{3,4,1},{0,2,7},40] (* Roman Witula, Aug 02 2012 *)

a(n)=([0,1,0; 0,0,1; 1,4,3]^(n-1)*[0;2;7])[1,1] \\ Charles R Greathouse IV, Jun 22 2016

@CachedFunction
def a(n): # a = A120757
    if (n<3): return (0,2,7)[n]
    else: return 3*a(n-1) + 4*a(n-2) + a(n-3)
[a(n) for n in range(40)] # G. C. Greubel, Nov 25 2022

a(n) = 3*a(n-1) + 4*a(n-2) + a(n-3) (follows from the minimal polynomial of the matrix M). See also the o.g.f. given in the name.

Edited by N. J. A. Sloane, Dec 03 2006

New name, old name as comment; o.g.f.; reference.

A122600 Expansion of 1/(1 + 3x - 4x^2 + x^3).

Original entry on oeis.org

1, -3, 13, -52, 211, -854, 3458, -14001, 56689, -229529, 929344, -3762837, 15235416, -61686940, 249765321, -1011279139, 4094585641, -16578638800, 67125538103, -271785755150, 1100438056662, -4455582728689, 18040286167865, -73043627475013, 295747609825188, -1197457625543481

Offset: 0

Roger L. Bagula and Gary W. Adamson, Sep 20 2006

Suggested by the Steinbach heptagon polynomial p^3 - 2*p^2*(1 - p) - p(1 - p)^2 + (1 - p)^3 = (1 - 4 p + 3 p^2 + p^3).

B(n):=|a(n-1)| = a(n-1)*(-1)^(n-1) with B(0):=0 (hence the o.g.f. for B(n) is x/(1 + 3*x - 4*x^2 + x^3)) appears in the following formula for the nonnegative powers of rho*sigma, where rho:=2*cos(Pi/7) and sigma:=sin(3*Pi/7)/sin(Pi/7)= rho^2-1 are the ratios of the smaller and larger diagonal length to the side length in a regular 7-gon (heptagon). See the Steinbach reference where the basis <1,rho,sigma> is used in an extension of the rational field. (rho*sigma)^n = C(n) + B(n)*rho + A(n)*sigma,n>=0, with C(n)= A120757(n) with C(0):=1, and A(n)= A181879(n). For the nonpositive powers see A085810*(-1)^n, A181880(n) and A116423(n)*(-1)^n, respectively. See also a comment under A052547.

Toufik Mansour and Mark Shattuck, Avoidance of vincular patterns by flattened derangements, arXiv:2504.14713 [math.CO], 2025. See p. 6.
Peter Steinbach, Golden fields: a case for the heptagon, Math. Mag. 70 (1997), no. 1, 22-31.
Index entries for linear recurrences with constant coefficients, signature (-3,4,-1).

Cf. A065941.

p[x_] := 1 - 4 x + 3x^2 + x^3; q[x_] := ExpandAll[x^3*p[1/x]]; Table[ SeriesCoefficient[ Series[x/q[x], {x, 0, 30}], n], {n, 0, 30}]
CoefficientList[Series[1/(1 + 3*x - 4*x^2 + x^3), {x, 0, 50}], x] (* or *) LinearRecurrence[{-3, 4, -1}, {1, -3, 13}, 40] (* Vladimir Joseph Stephan Orlovsky, Jan 31 2012 *)

a(n)= -3*a(n-1) + 4*a(n-2) - a(n-3), n>=2, a(-1):=0, a(1)=0, a(1)=-3 (from the o.g.f. given in the name).
a(n) = (-1)^n*Sum_{k=0..n} binomial(n+k+2,3*k+2)*7^k. - Emanuele Munarini, Aug 27 2017
From Kai Wang, Jul 05 2020: (Start)
a(n) =  Sum_{i+2j+3k=n} (-1)^(i+k)*3^i*4^j*((i+j+k)!)/(i!*j!*k!).
a(n) = (-1)^n*(6*A215076(n+4) - 21*A215076(n+3) - 13*A215076(n+2))/7. (End)

Edited by N. J. A. Sloane, Feb 01 2007

A215100 a(n) = 3a(n-1) + 4a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

Original entry on oeis.org

2, 5, 22, 88, 357, 1445, 5851, 23690, 95919, 388368, 1572470, 6366801, 25778651, 104375627, 422608286, 1711106017, 6928126822, 28051412820, 113577851765, 459867333397, 1861964820071, 7538941645566, 30524551550379, 123591386053472, 500411306007498, 2026124013786761

Offset: 0

Roman Witula, Aug 03 2012

Ramanujan-type sequence number 4 for the argument 2*Pi/7. We have a(n)=bs(3n+2), where the sequence bs(n) and its two conjugate sequences as(n) and cs(n) are defined in the comments to A214683 (see also A215076, A120757, A006053). Since we also have as(3n+2)=cs(3n+2)=0 from the formula for S(n) (see Comments at A214683) we obtain the relation 7^(1/3)*a(n)= (c(1)/c(4))^(n + 2/3) + (c(4)/c(2))^(n + 2/3) + (c(2)/c(1))^(n + 2/3).

			From 4*a(2) = a(3) = 88 we get 88*7^(1/3) = 4*((c(1)/c(4))^(8/3) + (c(4)/c(2))^(8/3) + (c(2)/c(1))^(8/3))=(c(1)/c(4))^(11/3) + (c(4)/c(2))^(11/3) + (c(2)/c(1))^(11/3).

R. Witula, E. Hetmaniok and D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,4,1).

Cf. A006053, A214683, A215076, A120757.

I:=[2,5,22]; [n le 3 select I[n] else 3*Self(n-1) +4*Self(n-2) +Self(n-3): n in [1..41]]; // G. C. Greubel, Nov 25 2022

LinearRecurrence[{3, 4, 1}, {2, 5, 22}, 40]

Vec((2-x-x^2)/(1-3*x-4*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

@CachedFunction
def a(n): # a = A215100
    if (n<3): return (2,5,22)[n]
    else: return 3*a(n-1) + 4*a(n-2) + a(n-3)
[a(n) for n in range(41)] # G. C. Greubel, Nov 25 2022

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

More terms from Michel Marcus, Apr 20 2016

A274032 Sum of n-th powers of the roots of x^3 + 9*x^2 - x - 1.

Original entry on oeis.org

3, -9, 83, -753, 6851, -62329, 567059, -5159009, 46935811, -427014249, 3884905043, -35344223825, 321555905219, -2925462465753, 26615373873171, -242142271419073, 2202970354179075, -20042260085157577, 182341168849178195, -1658909809373582257

Offset: 0

Kai Wang, Jun 07 2016

A Berndt-type sequence for tan(2*Pi/7).
a(n) is always an integer.
a(n) is x1^n + x2^n + x3^n, where x1, x2, x3 are the roots of the polynomial
x^3 + 9*x^2 - x - 1.
x1 = tan(Pi/7)/tan(2*Pi/7),
x2 = tan(2*Pi/7)/tan(4*Pi/7),
x3 = tan(4*Pi/7)/tan(Pi/7).
This is a two sided sequence. The other half is A274075. - Kai Wang, Aug 02 2016

Seiichi Manyama, Table of n, a(n) for n = 0..1000
B. C. Berndt, L.-C. Zhang, Ramanujan's identities for eta-functions, Math. Ann. 292 (1992), 561-573.
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
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
Roman Witula and Damian Slota, Quasi-Fibonacci Numbers of Order 11, Journal of Integer Sequences, Vol. 10 (2007), Article 07.8.5
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (-9,1,1).

Cf. A033304, A094648, A215076, A274075.

Vec((3+18*x-x^2)/(1+9*x-x^2-x^3) + O(x^30)) \\ Colin Barker, Jun 07 2016

polsym(x^3 + 9*x^2 - x - 1, 30) \\ Charles R Greathouse IV, Jul 20 2016

a(n) = (tan(Pi/7)/tan(2*Pi/7))^n + (-tan(2*Pi/7)/tan(3*Pi/7))^n + (-tan(3*Pi/7)/tan(Pi/7))^n.
From Colin Barker, Jun 07 2016: (Start)
a(n) = -9*a(n-1)+a(n-2)+a(n-3) for n>2.
G.f.: (3+18*x-x^2) / (1+9*x-x^2-x^3).
(End)

Edited by N. J. A. Sloane, Jun 07 2016

A215560 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=a(1)=3, a(2)=101.

Original entry on oeis.org

3, 3, 101, 444, 5981, 38468, 390974, 2948431, 26868565, 216624495, 1888775906, 15657923053, 134074085330, 1124375492334, 9556192325235, 80523923708399, 682280993578341, 5760499663646612, 48746948619251921, 411906111379078256, 3483838470286469746, 29447943482916260935

Offset: 0

Roman Witula, Aug 16 2012

The Ramanujan-type sequence number 6 for the argument 2Pi/7 (see also A214683, A215112, A006053, A006054, A215076, A215100, A120757 for the numbers: 1, 1a, 2, 2a, 3, 4 and 5 respectively).
The sequence a(n) is one of the three special sequences (the remaining two are A215569 and A215572) connected with the following recurrence relation: T(n):=49^(1/3)*T(n-2)+T(n-3), with T(0)=3, T(1)=0, and T(2)=2*49^(1/3) - see the comments to A214683.
It can be proved that
T(n) = (c(1)^4/c(2))^(n/3) + (c(2)^4/c(4))^(n/3) + (c(4)^4/c(1))^(n/3), where c(j):=2*cos(2*Pi*j/7), and the following decomposition hold true:
T(n) = at(n) + bt(n)*7^(1/3) + ct(n)*49^(1/3), where sequences at(n), bt(n), and ct(n) satisfy the following system of recurrence equations: at(n)=7*bt(n-2)+at(n-3),
bt(n)=7*ct(n-2)+bt(n-3), ct(n)=at(n-2)+ct(n-3), with at(0)=3, at(1)=at(2)=bt(0)=bt(1)=bt(2)=ct(0)=ct(1)=0, ct(2)=2 - for details see the first Witula reference.
It follows that a(n)=at(3*n), bt(3*n)=ct(3*n)=0.
Every difference of the form a(n)-a(n-2)-a(n-3) is divisible by 3. Because the difference a(n+1)-a(n) is congruent to the difference a(n-4)-a(n-2) modulo 3 we easily deduce that a(6)-a(5) and a(7)-a(6)-2 are both divisible by 3.

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012

Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,46,1).

Cf. A214683, A215112, A006053, A006054, A215076, A215100, A120757, A214699.

LinearRecurrence[{3, 46, 1}, {3, 3, 101}, 50]

Vec((3-6*x-46*x^2)/(1-3*x-46*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

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

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

More terms from Michel Marcus, Apr 20 2016

A215572 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=106.

Original entry on oeis.org

2, 5, 106, 550, 6531, 44999, 435973, 3384404, 30252969, 246877464, 2135653370, 17793576423, 151867661753, 1276243154087, 10832435479322, 91356359187721, 773637352766062, 6534137016412674, 55281085635664595, 467187197014742851, 3951025667301212597, 33398969150217473532

Offset: 0

Roman Witula, Aug 16 2012

The Ramanujan-type sequence number 8 for the argument 2Pi/7 (see also A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560, A215569 for the numbers: 1, 1a, 2, 2a, 3-7 respectively). The sequence a(n) is one of the three special sequences (the remaining two are A215560 and A215569) connected with the following recurrence relation:
(c(1)^4/c(2))^(n/3) + (c(2)^4/c(4))^(n/3) + (c(4)^4/c(1))^(n/3) = at(n) + bt(n)*7^(1/3) + ct(n)*49^(1/3), where c(j):=2*cos(2*Pi*j/7), and the sequences at(n), bt(n), and ct(n) are defined in comments to A215560 (see also A215569). It follows that a(n)=ct(3*n+2), at(3*n+2)=bt(3*n+2)=0, which implies the first formula below.
We note that if a(n), a(n+1) and a(n+2) are all odd for some n in N then a(n+3) is even, a(n+4) is odd, a(n+5) and a(n+6) are both even, and the numbers a(n+7), a(n+8), a(n+9) are all odd again. In consequence, this situation hold for every n of the form 7*k+4, k=0,1,..., in the other words cyclical through all sequence a(n), n=4,5,... (from n=1 whenever we start from odd-even-even sequence).

			From 4*a(1)+5*a(2)=a(3) we obtain 4*((c(1)^4/c(2))^(5/3) + (c(2)^4/c(4))^(5/3) + (c(4)^4/c(1))^(5/3)) + 5*((c(1)^4/c(2))^(8/3) + (c(2)^4/c(4))^(8/3) + (c(4)^4/c(1))^(8/3)) = (4 + 5*c(1)^4/c(2))*((c(1)^4/c(2))^(5/3) + (4 + 5*c(2)^4/c(4))*((c(2)^4/c(4))^(5/3) + (4 + 5*c(4)^4/c(1))*((c(4)^4/c(1))^(5/3) = (c(1)^4/c(2))^(11/3) + (c(2)^4/c(4))^(11/3) + (c(4)^4/c(1))^(11/3) = 550*49^(1/3).

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,46,1).

Cf. A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560, A215569.

LinearRecurrence[{3,46,1}, {2,5,106}, 50]
CoefficientList[Series[(2 - x - x^2)/(1 - 3*x - 46*x^2 - x^3), {x,0,50}], x] (* G. C. Greubel, Apr 16 2017 *)

Vec((2-x-x^2)/(1-3*x-46*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

49^(1/3)*a(n) = (c(1)^4/c(2))^(n+2/3) + (c(2)^4/c(4))^(n+2/3) + (c(4)^4/c(1))^(n+2/3) = (c(1)*(c(1)/c(2))^(1/3))^(3*n+2) + (c(2)*(c(2)/c(4))^(1/3))^(3*n+2) + (c(4)*(c(4)/c(1))^(1/3))^(3*n+2).

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

More terms from Michel Marcus, Apr 20 2016

A274075 Sum of n-th powers of the roots of x^3 + x^2 - 9*x - 1.

Original entry on oeis.org

3, -1, 19, -25, 195, -401, 2131, -5545, 24323, -72097, 285459, -910009, 3407043, -11311665, 41065043, -139462985, 497736707, -1711838529, 6052005907, -20960815961, 73717030595, -256312368337, 898804827731, -3131899112169, 10964830193411, -38253117375201

Offset: 0

Kai Wang, Jun 09 2016

a(n) is always an integer.
This is the other half of A274032.
a(n) is x1^n + x2^n + x3^n, where x1, x2, x3 are the roots of the polynomial
x^3 + x^2 - 9*x - 1.
x1 = tan(Pi/7)/tan(4*Pi/7),
x2 = tan(4*Pi/7)/tan(2*Pi/7),
x3 = tan(2*Pi/7)/tan(Pi/7).

Colin Barker, Table of n, a(n) for n = 0..1000
B. C. Berndt, L.-C. Zhang, Ramanujan's identities for eta-functions, Math. Ann. 292 (1992), 561-573.
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
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
Roman Witula and Damian Slota, Quasi-Fibonacci Numbers of Order 11, Journal of Integer Sequences, Vol. 10 (2007), Article 07.8.5
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (-1,9,1).

Cf. A033304, A094648, A215076, A274032.

FullSimplify[Table[(Tan[Pi/7]/Tan[4*Pi/7])^n + (Tan[4*Pi/7]/Tan[2*Pi/7])^n + (Tan[2*Pi/7]/Tan[Pi/7])^n, {n, 0, 12}]] (* Wesley Ivan Hurt, Jun 11 2016 *)

Vec((3+2*x-9*x^2)/(1+x-9*x^2-x^3) + O(x^30)) \\ Colin Barker, Jun 11 2016

polsym(x^3 + x^2 - 9*x - 1, 30) \\ Charles R Greathouse IV, Jul 20 2016

a(n) = (tan(Pi/7)/tan(4*Pi/7))^n + (tan(4*Pi/7)/tan(2*Pi/7))^n + (tan(2*Pi/7)/tan(Pi/7))^n.
a(n) = -a(n-1) + 9*a(n-2) + a(n-3) for n>2.
G.f.: (3+2*x-9*x^2) / (1+x-9*x^2-x^3). - Colin Barker, Jun 11 2016

A274220 a(n) = (-cos(Pi/7)/cos(2Pi/7))^n + (-cos(2Pi/7)/cos(3Pi/7))^n + (cos(3Pi/7)/cos(Pi/7))^n.

Original entry on oeis.org

3, -4, 10, -25, 66, -179, 493, -1369, 3818, -10672, 29865, -83626, 234237, -656205, 1838483, -5151080, 14432666, -40438941, 113306686, -317477255, 889550021, -2492461633, 6983719214, -19567941936, 54828148469, -153625048854, 430447808073, -1206087937261, 3379383275971, -9468821484028

Offset: 0

Kai Wang, Jun 14 2016

a(n) is an integer.
This is other half of A215076.
a(n) is the sum of n-th powers of the roots of x^3 + 4*x^2 + 3*x - 1. - Greg Dresden, Mar 11 2020

			a(0) = 3, a(1) = -4, a(2) = 10, a(3) = -25.

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

Colin Barker, Table of n, a(n) for n = 0..1000
B. C. Berndt, L.-C. Zhang, Ramanujan's identities for eta-functions, Math. Ann. 292 (1992), 561-573.
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7, J. Integer Seq., 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-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.
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (-4,-3,1).

Cf. A215076, A033304, A094648, A274032.

CoefficientList[Series[(3 + 8 x + 3 x^2)/(1 + 4 x + 3 x^2 - x^3), {x, 0, 29}], x] (* Michael De Vlieger, Jun 14 2016 *)

Vec((3+8*x+3*x^2)/(1+4*x+3*x^2-x^3) + O(x^30)) \\ Colin Barker, Jun 14 2016

polsym(x^3 + 4*x^2 + 3*x - 1,33) \\ Joerg Arndt, Mar 12 2020

a(n) = -4*a(n-1)-3*a(n-2)+a(n-3).

G.f.: (3+8*x+3*x^2) / (1+4*x+3*x^2-x^3). - Colin Barker, Jun 14 2016

Many terms corrected by Colin Barker, Jun 14 2016

A215569 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=0, a(1)=14, a(2)=49.

Original entry on oeis.org

0, 14, 49, 791, 4641, 50358, 365351, 3417162, 27107990, 238878773, 1967021021, 16916594611, 141471629572, 1204545261843, 10138247340452, 85965295695706, 725459810009753, 6140921279372187, 51879880394260905, 438847479843913070, 3709157858947113027

Offset: 0

Roman Witula, Aug 16 2012

The Ramanujan-type sequence number 7 for the argument 2Pi/7 (see also A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560 for the numbers: 1, 1a, 2, 2a, 3-6 respectively). The sequence a(n) is one of the three special sequences (the remaining two are A215560 and A215572) connected with the following recurrence relation:
(c(1)^4/c(2))^(n/3) + (c(2)^4/c(4))^(n/3) + (c(4)^4/c(1))^(n/3) = at(n) + bt(n)*7^(1/3) + ct(n)*49^(1/3), where c(j):=2*cos(2*Pi*j/7), and the sequences at(n), bt(n), and ct(n) satisfy the following system of recurrence equations: at(n)=7*bt(n-2)+at(n-3),
bt(n)=7*ct(n-2)+bt(n-3), ct(n)=at(n-2)+ct(n-3), with at(0)=3, at(1)=at(2)=bt(0)=bt(1)=bt(2)=ct(0)=ct(1)=0, ct(2)=2 - for details see the Witula's first paper (see also A215560). It follows that a(n)=bt(3*n+1), at(3*n+1)=ct(3*n+1)=0, which implies the first formula below.
We note that all numbers a(n) are divided by 7.

			We have  (c(1)^4/c(2))^(4/3) + (c(2)^4/c(4))^(4/3) + (c(4)^4/c(1))^(4/3) = (2/7)*(c(1)^4/c(2))^(7/3) + (c(2)^4/c(4))^(7/3) + (c(4)^4/c(1))^(7/3)).

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,46,1).

Cf. A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560.

LinearRecurrence[{3,46,1},{0,14,49},30] (* Harvey P. Dale, Jan 12 2015 *)

7^(1/3)*a(n) = (c(1)^4/c(2))^(n+1/3) + (c(2)^4/c(4))^(n+1/3) + (c(4)^4/c(1))^(n+1/3) = (c(1)*(c(1)/c(2))^(1/3))^(3*n+1) + (c(2)*(c(2)/c(4))^(1/3))^(3*n+1) + (c(4)*(c(4)/c(1))^(1/3))^(3*n+1).

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

More terms from Harvey P. Dale, Jan 12 2015

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A120757 Expansion of x^2*(2+x)/(1-3*x-4*x^2-x^3).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A122600 Expansion of 1/(1 + 3*x - 4*x^2 + x^3).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A215100 a(n) = 3*a(n-1) + 4*a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

Extensions

A274032 Sum of n-th powers of the roots of x^3 + 9*x^2 - x - 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Formula

Extensions

A215560 a(n) = 3*a(n-1) + 46*a(n-2) + a(n-3) with a(0)=a(1)=3, a(2)=101.

Views

Author

Keywords

Comments

A120757 Expansion of x^2(2+x)/(1-3x-4*x^2-x^3).

A122600 Expansion of 1/(1 + 3x - 4x^2 + x^3).

A215100 a(n) = 3a(n-1) + 4a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

A215560 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=a(1)=3, a(2)=101.

A215572 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=106.

A274220 a(n) = (-cos(Pi/7)/cos(2Pi/7))^n + (-cos(2Pi/7)/cos(3Pi/7))^n + (cos(3Pi/7)/cos(Pi/7))^n.

A215569 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=0, a(1)=14, a(2)=49.