This site is supported by donations to The OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A215794 a(n) = -7^n*A(2*n+1), where A(n) = A(n-1) + A(n-2) + A(n-3)/7, with A(0)=3, A(1)=1, A(2)=3. 16
 -1, -31, -609, -11711, -224833, -4315871, -82846113, -1590286719, -30526618241, -585978870687, -11248256653025, -215917815567167, -4144686996149441, -79560041170858591, -1527208244431770145, -29315784501060168447, -562736106255347592449 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,2 COMMENTS The Berndt-type sequence number 12 for the argument 2Pi/7   defined by the relation sqrt(7)*a(n) = t(1)^(2*n+1) + t(2)^(2*n+1) + t(4)^(2*n+1) = (-sqrt(7) + 4*s(1))^(2*n+1) + (-sqrt(7) + 4*s(2))^(2*n+1) + (-sqrt(7) + 4*s(4))^(2*n+1), where t(j) := tan(2*Pi*j/7) and s(j) := sin(2*Pi*j/7) (the respective sum with even powers in A108716 are given, see also A215828). We note that sqrt(7)*a(n) = B(2*n+1), where B(n) is defined in the comments to A215575. From Witula-Slota's (Section 6) and Witula's (Remark 11) papers it follows that B(n) is equal to the product (-sqrt(7))^n by the value of big omega function with index n for the argument 2*i/sqrt(7). The last value is equal to A(n). The respective recurrence relation for A(n) from the following decomposition follow (see Witula-Slota's paper for details): (X-1-2*i*d*s(1))*(X-1-2*i*d*s(2))*(X-1- 2*i*d*s(4)) = X^3 - (3+i*sqrt(7))*X^2 + (3+i*2*sqrt(7)*d)*X - (1+i*sqrt(7)*d + i*sqrt(7)*d^3), since the big omega function with index n for the argument d is equal to the sum: (1 + 2*i*d*s(1))^n + (1 + 2*i*d*s(2))^n + (1 + 2*i*d*s(4))^n and it is equal to 3 for n=0, 3 + i*sqrt(7)*d for n=1, and at last 3 + 2*i*sqrt(7)*d - 7*d^2 for n=2. The sequence a(n+1)/a(n) is decreasing and convergent to (t(2))^2 = 19,195669... Moreover we have floor(a(n+1)/a(n)) = 19 for every n=1,2,... LINKS Vincenzo Librandi, Table of n, a(n) for n = 0..200 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, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5 Index entries for linear recurrences with constant coefficients, signature (21,-35,7). FORMULA G.f.: -(1+10*x-7*x^2)/(1-21*x+35*x^2-7*x^3). [Bruno Berselli, Aug 30 2012] a(n) = -A275195(2*n-1)/(7^n). - Kai Wang, Aug 02 2016 EXAMPLE We have -31*sqrt(7) = t(1)^3 + t(2)^3 + t(4)^3. MATHEMATICA LinearRecurrence[{21, -35, 7}, {-1, -31, -609}, 17] (* Bruno Berselli, Aug 30 2012 *) PROG (MAGMA) m:=17; R:=PowerSeriesRing(Integers(), m); Coefficients(R!(-(1+10*x-7*x^2)/(1-21*x+35*x^2-7*x^3)));  // Bruno Berselli, Aug 30 2012 (MAGMA) I:=[-1, -31, -609]; [n le 3 select I[n] else 21*Self(n-1)-35*Self(n-2)+7*Self(n-3): n in [1..20]]; // Vincenzo Librandi, Mar 19 2013 CROSSREFS Cf. A215575, A108716, A215828, A275195. Sequence in context: A005462 A028201 A028184 * A028182 A028155 A028175 Adjacent sequences:  A215791 A215792 A215793 * A215795 A215796 A215797 KEYWORD sign,easy AUTHOR Roman Witula, Aug 23 2012 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified November 13 17:34 EST 2019. Contains 329106 sequences. (Running on oeis4.)