A215917 -id:A215917 - cp's OEIS frontend

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

3, -6, 18, -63, 234, -891, 3429, -13257, 51354, -199098, 772173, -2995218, 11619045, -45073827, 174857211, -678335958, 2631522330, -10208681991, 39603398850, -153636822171, 596016389349, -2312177133105, 8969825761002

Offset: 0

Roman Witula, Aug 18 2012

The Berndt-type sequence number 2 for the argument 2Pi/9 . Similarly like the respective sequence number 1 -- see A215455 -- the sequence a(n) is connected with the following general recurrence relation: X(n+3) + 6*X(n+2) + 9*X(n+1) + ((2*cos(3*g))^2)*X(n) = 0, X(0)=3, X(1)=-6, X(2)=18. The Binet formula for this one has the form: X(n) = (-4)^n*((cos(g))^(2*n) + cos(g+Pi/3))^(2*n) + cos(g-Pi/3))^(2*n)) - for details see Witula-Slota's reference and comments to A215455.

The characteristic polynomial of a(n) has the form x^3 + 6*x^2 + 9*x + 3 = (x + (2*cos(Pi/18))^2)*(x+(2*cos(5*Pi/18))^2)*(x+(2*cos(7*Pi/18))^2). We note that (2*cos(Pi/18))^2 = 2 - c(4), (2*cos(5*Pi/18))^2 = 2 - c(2), and (2*cos(7*Pi/18))^2 = 2 - c(1), where c(j) = 2*cos(2*Pi*j/9) - see trigonometric relations for A215455. Furthermore all numbers a(n)*3^(-ceiling((n+1)/3)) are integers.

R. Witula, On some applications of formulas for sums of the unimodular complex numbers, Wyd. Pracowni Komputerowej Jacka Skalmierskiego, Gliwice 2011 (in Polish).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
R. Witula, 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,-3).

Cf. A215455, A215635, A215636, A215664, A215885, A215665, A215666, A215829, A215831, A215917, A215919, A215945, A216034, A215948, A216757.

I:=[3,-6,18]; [n le 3 select I[n] else -6*Self(n-1)-9*Self(n-2)-3*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Aug 30 2017

LinearRecurrence[{-6,-9,-3}, {3,-6,18}, 50]
CoefficientList[Series[(3 + 12 x + 9 x^2)/(1 + 6 x + 9 x^2 + 3 x^3), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 30 2017 *)

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

a(n) = (-4)^n*((cos(Pi/18))^(2*n) + (cos(5*Pi/18))^(2*n) + (cos(7*Pi/18))^(2*n)).
G.f.: (3 + 12*x + 9*x^2)/(1 + 6*x + 9*x^2 + 3*x^3).
a(n)*(-1)^n = s(1)^(2*n) + s(2)^(2*n) + s(4)^(2*n), where s(j) := 2*sin(2*Pi*j/9) -- for the proof see Witula's book. The respective sums with odd powers of sines in A216757 are given. - Roman Witula, Sep 15 2012

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).

A215636 a(n) = - 12a(n-1) - 54a(n-2) - 112a(n-3) - 105a(n-4) - 36a(n-5) - 2a(n-6) with a(0)=a(1)=a(2)=0, a(3)=-3, a(4)=24, a(5)=-135.

Original entry on oeis.org

0, 0, 0, -3, 24, -135, 660, -3003, 13104, -55689, 232500, -958617, 3916440, -15890355, 64127700, -257698347, 1032023136, -4121456625, 16421256420, -65301500577, 259259758056, -1027901275131, 4070632899300, -16104283594083, 63657906293520, -251447560563465, 992593021410900

Offset: 0

Roman Witula, Aug 18 2012

The Berndt-type sequence number 4 for the argument 2*Pi/9 defined by the relation: X(n) = b(n) + a(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) = A215635(n) (see also section "Example" below). For more details - see comments to A215635, A215634 and Witula-Slota's reference.

			We have X(1)=-6, X(2)=18 and X(3)=-60-3*sqrt(2), which implies the equality: (cos(Pi/24))^6 + (cos(7*Pi/24))^6 + (cos(3*Pi/8))^6 = (60+3*sqrt(2))/64.

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, A215635, A215664, A215885, A215665, A215666, A215829, A215831, A215917, A215919, A215945, A216034, A215948, A216757.

LinearRecurrence[{-12,-54,-112,-105,-36,-2}, {0,0,0,-3,24,-135}, 50]

G.f.: (-3*x^3-12*x^4-9*x^5)/(1+12*x+54*x^2+112*x^3+105*x^4+36*x^5+2*x^6).

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

Original entry on oeis.org

0, -3, 12, -36, 105, -303, 873, -2514, 7239, -20844, 60018, -172815, 497601, -1432785, 4125540, -11879019, 34204272, -98487276, 283582809, -816544155, 2351145189, -6769852758, 19493014119, -56127897168, 161613838746, -465348502119, 1339917609189, -3858138988821

Offset: 0

Roman Witula, Aug 27 2012

The Berndt-type sequence number 10 for the argument 2Pi/9 defined by the first trigonometric relation from the section "Formula" below. The sequence a(n) is connected with sequences A215917 and A215885 - see the respective formula.

We have A035045(n)=abs(a(n+1)/3) for every n=0,1,...,5 and A035045(7) + a(7)/3 = 1, A035045(8) - a(8)/3 = 10, A035045(9) + a(9)/3 = 63, and A035045(10) - a(10)/3 = 320 - all these four results-numbers are in A069269.

			We have a(2)=-4*a(1), a(3)=-3*a(2), a(6)/a(3) = -24.25, and a(9) = 579*a(3).

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

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

Cf. A215917, A215885, A215664.

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

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).
a(n) = A215917(n+1) + A215917(n) - 2*(-1)^n*A215885(n).
G.f.: -3*x*(1-x)/(1+3*x-x^3).

A217336 a(n) = 3^(-1+floor(n/2))A(n), where A(n) = 3A(n-1) + A(n-2) - A(n-3)/3 with A(0)=A(1)=3, A(2)=11.

Original entry on oeis.org

1, 1, 11, 35, 345, 1129, 11091, 36315, 356721, 1168017, 11473371, 37567443, 369023049, 1208298105, 11869049763, 38863020555, 381749439969, 1249968331809, 12278374244523, 40203278289027, 394914722339385, 1293075627640713

Offset: 0

Roman Witula, Oct 01 2012

The Berndt-type sequence number 14 for the argument 2Pi/9 defined by the relation: A(n)*(-sqrt(3))^n = t(1)^n + (-t(2))^n + t(4)^n = (-sqrt(3) + 4*s(1))^n + (-sqrt(3) - 4*s(2))^n + (-sqrt(3) + 4*s(4))^n, where s(j) := sin(2*Pi*j/9) and t(j) := tan(2*Pi*j/9).
The definitions of the other Berndt-type sequences for the argument 2Pi/9 like A215945, A215948, A216034 in Crossrefs are given.
We note that all a(2*n), n=2,3,..., are divisible by 3, and it is only when n=5 that a(2*n) is divisible by 9.

			Note that A(0)=A(1)=3, a(0)=a(1)=1, A(2)=a(2)=11, A(3)=a(3)=35, A(4)=115, a(4)=345 and A(5) = 1129/3, which implies the equality  3387*sqrt(3) = -t(1)^5 + t(2)^5 - t(4)^5.

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

Paolo Xausa, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,33,0,-27,0,3).

Cf. A215455, A215634-A215636, A215664, A215885, A215665, A215666, A215829, A215831, A215917, A215919, A215945, A216034, A215948, A216757.

/* By definition: */ i:=22; I:=[3,3,11]; A:=[m le 3 select I[m] else 3*Self(m-1)+Self(m-2)-Self(m-3)/3: m in [1..i]]; [3^(-1+Floor((n-1)/2))*A[n]: n in [1..i]]; // Bruno Berselli, Oct 02 2012

LinearRecurrence[{0, 33, 0, -27, 0, 3}, {1, 1, 11, 35, 345, 1129},25] (* Paolo Xausa, Feb 23 2024 *)

G.f.: (1+x-22*x^2+2*x^3+9*x^4+x^5)/(1-33*x^2+27*x^4-3*x^6). - Bruno Berselli, Oct 01 2012

A218489 The sequence of coefficients of cubic polynomials p(x+n), where p(x) = x^3 - 3*x + 1.

Original entry on oeis.org

1, 0, -3, 1, 1, 3, 0, -1, 1, 6, 9, 3, 1, 9, 24, 19, 1, 12, 45, 53, 1, 15, 72, 111, 1, 18, 105, 199, 1, 21, 144, 323, 1, 24, 189, 489, 1, 27, 240, 703, 1, 30, 297, 971, 1, 33, 360, 1299, 1, 36, 429, 1693, 1, 39, 504, 2159, 1, 42, 585, 2703, 1, 45, 672, 3331

Offset: 0

Roman Witula, Oct 30 2012

sign

We note that p(x) =  (x - s(1))*(x + c(1))*(x - c(2)),
p(x+1) = x^3 + 3*x^2 -1 = (x + s(1)*c(1))*(x - s(1)*c(2))*(x + c(1)*c(2)), p(x+2) = x^3 + 6*x^2 + 9*x + 3 = (x + c(1/2)^2)*(x + s(2)^2)*(x + s(4)^2), and p(x + n) = (x + n - 2 + c(1/2)^2)*(x + n - 2 + s(2)^2)*(x + n - 2 + s(4)^2), n = 2,3,..., where c(j) := 2*cos(Pi*j/9) and s(j) := 2*sin(Pi*j/18). These one's are characteristic polynomials many sequences A... - see crossrefs.
A218332 is the sequence of coefficients of polynomials p(x-n).

Cf. A214699, A214779, A215634, A215664, A215666, A215917, A215919.

We have a(4*k) = 1, a(4*k + 1) = 3*k, a(4*k + 2) = 3*k^2 - 3, and a(4*k + 3) = k^3 - 3*k + 1. Moreover we obtain
b(k+1) = b(k) + 3, c(k+1) = 2*b(k) + c(k) + 3, d(k+1) = b(k) + c(k) + d(k) + 1, where p(x + k) = x^3 + b(k)*x^2 + c(k)*x + d(k).
Empirical g.f.: -(3*x^15-3*x^13+x^12-13*x^11+9*x^10+6*x^9-3*x^8+5*x^7-12*x^6-3*x^5+3*x^4-x^3+3*x^2-1) / ((x-1)^4*(x+1)^4*(x^2+1)^4). - Colin Barker, May 17 2013

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

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

Mathematica

PARI

Formula

A215636 a(n) = - 12*a(n-1) - 54*a(n-2) - 112*a(n-3) - 105*a(n-4) - 36*a(n-5) - 2*a(n-6) with a(0)=a(1)=a(2)=0, a(3)=-3, a(4)=24, a(5)=-135.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

A217336 a(n) = 3^(-1+floor(n/2))*A(n), where A(n) = 3*A(n-1) + A(n-2) - A(n-3)/3 with A(0)=A(1)=3, A(2)=11.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

Formula

A218489 The sequence of coefficients of cubic polynomials p(x+n), where p(x) = x^3 - 3*x + 1.

Views

Author

Keywords

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

A215636 a(n) = - 12a(n-1) - 54a(n-2) - 112a(n-3) - 105a(n-4) - 36a(n-5) - 2a(n-6) with a(0)=a(1)=a(2)=0, a(3)=-3, a(4)=24, a(5)=-135.

A217336 a(n) = 3^(-1+floor(n/2))A(n), where A(n) = 3A(n-1) + A(n-2) - A(n-3)/3 with A(0)=A(1)=3, A(2)=11.