A287405 a(n) = (7*(cot(1*Pi/7))^2)^n + (7*(cot(2*Pi/7))^2)^n + (7*(cot(4*Pi/7))^2)^n.

3, 35, 931, 27587, 830403, 25054435, 756187747, 22824258947, 688917131651, 20793986742179, 627637106311971, 18944339609269571, 571808137046942019, 17259221092289630307, 520945214725090792931, 15723995613526902256387, 474606601742375424297731

Offset: 0

Author

Kai Wang, May 24 2017

Keywords

• nonn
• easy

Comments

a(n) = x1^n + x2^n + x3^n, where x1, x2, x3 are the roots of x^3 - 35*x^2 + 147*x - 49, x1 = 7*(cot(1*Pi/7))^2, x2 = 7*(cot(2*Pi/7))^2, x3 = 7*(cot(4*Pi/7))^2.

Links

• Colin Barker, Table of n, a(n) for n = 0..650
• 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.
• Index entries for linear recurrences with constant coefficients, signature (35,-147,49).

Crossrefs

Cf. A108716, A033304, A094648, A274220, A215076, A274075 A274032, A275195, A215575, A274975.

Programs

Mathematica

LinearRecurrence[{35,-147,49},{3,35,931},30] (* Harvey P. Dale, Mar 15 2018 *)

PARI

Vec((3 - 7*x)*(1 - 21*x) / (1 - 35*x + 147*x^2 - 49*x^3) + O(x^30)) \\ Colin Barker, May 26 2017

PARI

polsym(x^3 - 35*x^2 + 147*x - 49, 20) \\ Joerg Arndt, May 26 2017

Formula

a(n) = 35*a(n-1) - 147*a(n-2) + 49*a(n-3), a(0) = 3, a(1) = 35, a(2) = 931.

Bisection of A215575: a(n) = A215575(2*n).

G.f.: (3 - 7*x)*(1 - 21*x) / (1 - 35*x + 147*x^2 - 49*x^3). - Colin Barker, May 26 2017