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

%I #39 Aug 03 2016 00:17:32

%S 3,-9,83,-753,6851,-62329,567059,-5159009,46935811,-427014249,

%T 3884905043,-35344223825,321555905219,-2925462465753,26615373873171,

%U -242142271419073,2202970354179075,-20042260085157577,182341168849178195,-1658909809373582257

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

%C A Berndt-type sequence for tan(2*Pi/7).

%C a(n) is always an integer.

%C a(n) is x1^n + x2^n + x3^n, where x1, x2, x3 are the roots of the polynomial

%C x^3 + 9*x^2 - x - 1.

%C x1 = tan(Pi/7)/tan(2*Pi/7),

%C x2 = tan(2*Pi/7)/tan(4*Pi/7),

%C x3 = tan(4*Pi/7)/tan(Pi/7).

%C This is a two sided sequence. The other half is A274075. - _Kai Wang_, Aug 02 2016

%H Seiichi Manyama, <a href="/A274032/b274032.txt">Table of n, a(n) for n = 0..1000</a>

%H B. C. Berndt, L.-C. Zhang, <a href="http://dx.doi.org/10.1007/BF01444636">Ramanujan's identities for eta-functions</a>, Math. Ann. 292 (1992), 561-573.

%H Roman Witula, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL12/Witula/witula17.html">Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2Pi/7</a>, J. Integer Seq., 12 (2009), Article 09.8.5.

%H Roman Witula and Damian Slota, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL10/Slota/witula13.html">New Ramanujan-Type Formulas and Quasi-Fibonacci Numbers of Order 7</a>, Journal of Integer Sequences, Vol. 10 (2007), Article 07.5.6

%H Roman Witula and Damian Slota, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL10/Slota2/slota99.html">Quasi-Fibonacci Numbers of Order 11</a>, Journal of Integer Sequences, Vol. 10 (2007), Article 07.8.5

%H Roman Witula, Damian Slota and Adam Warzynski, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL9/Slota/slota57.html">Quasi-Fibonacci Numbers of the Seventh Order</a>, J. Integer Seq., 9 (2006), Article 06.4.3.

%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (-9,1,1).

%F 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.

%F From _Colin Barker_, Jun 07 2016: (Start)

%F a(n) = -9*a(n-1)+a(n-2)+a(n-3) for n>2.

%F G.f.: (3+18*x-x^2) / (1+9*x-x^2-x^3).

%F (End)

%o (PARI) Vec((3+18*x-x^2)/(1+9*x-x^2-x^3) + O(x^30)) \\ _Colin Barker_, Jun 07 2016

%o (PARI) polsym(x^3 + 9*x^2 - x - 1, 30) \\ _Charles R Greathouse IV_, Jul 20 2016

%Y Cf. A033304, A094648, A215076, A274075.

%K sign,easy

%O 0,1

%A _Kai Wang_, Jun 07 2016

%E Edited by _N. J. A. Sloane_, Jun 07 2016