# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a214699 Showing 1-1 of 1 %I A214699 #53 Jan 10 2024 04:58:48 %S A214699 0,3,0,9,-3,27,-18,84,-81,270,-327,891,-1251,3000,-4644,10251,-16932, %T A214699 35397,-61047,123123,-218538,430416,-778737,1509786,-2766627,5308095, %U A214699 -9809667,18690912,-34737096,65882403,-122902200,232384305,-434589003,820055115,-1536151314 %N A214699 a(n) = 3*a(n-2) - a(n-3) with a(0)=0, a(1)=3, a(2)=0. %C A214699 All a(n) are divisible by 3. %C A214699 The Ramanujan-type sequence number 1 for the argument 2*Pi/9 defined by the following identity: %C A214699 3^(1/3)*a(n) = (c(1)/c(2))^(1/3)*c(1)^n + (c(2)/c(4))^(1/3)*c(2)^n + (c(4)/c(1))^(1/3)*c(4)^n = -( (c(1)/c(2))^(1/3)*c(2)^(n+1) + (c(2)/c(4))^(1/3)*c(4)^(n+1) + (c(4)/c(1))^(1/3)*c(1)^(n+1) ), where c(j) := 2*cos(2*Pi*j/9). %C A214699 The definitions of other Ramanujan-type sequences, for the argument of 2*Pi/9 in one's, are given in the Crossrefs section. %D A214699 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 %H A214699 G. C. Greubel, Table of n, a(n) for n = 0..1000 %H A214699 Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7. %H A214699 Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5. %H A214699 Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796. %H A214699 Index entries for linear recurrences with constant coefficients, signature (0,3,-1). %F A214699 G.f.: 3*x/(1 - 3*x^2 + x^3). %F A214699 From _Roman Witula_, Oct 06 2012: (Start) %F A214699 a(n+1) = 3*(-1)^n*A052931(n), which from recurrence relations for a(n) and A052931 can easily be proved inductively. %F A214699 a(n) = -A214779(n+1) - A214779(n). (End) %e A214699 We have a(2) = a(1) + a(4) = a(4) + a(7) + a(8) = -a(3) + a(5) + a(6) = 0, which implies %e A214699 (c(1)/c(2))^(1/3)*c(1)^2 + (c(2)/c(4))^(1/3)*c(2)^2 + (c(4)/c(1))^(1/3)*c(4)^2 = (c(1)/c(2))^(1/3)*(c(1) + c(1)^4) + (c(2)/c(4))^(1/3)*(c(2) + c(2)^4) + (c(4)/c(1))^(1/3)*(c(4) + c(4)^4) = (c(1)/c(2))^(1/3)*(c(1)^4 + c(1)^7 + c(1)^8) + (c(2)/c(4))^(1/3)*(c(2)^4 + c(2)^7 + c(2)^8) + (c(4)/c(1))^(1/3)*(c(4)^4 + c(4)^7 + c(4)^8) = 0. %e A214699 Moreover we have 3000*3^(1/3) = (c(1)/c(2))^(1/3)*c(1)^13 + (c(2)/c(4))^(1/3)*c(2)^13 + (c(4)/c(1))^(1/3)*c(4)^13. - _Roman Witula_, Oct 06 2012 %t A214699 LinearRecurrence[{0,3,-1}, {0,3,0}, 30] %t A214699 CoefficientList[Series[3*x/(1 - 3*x^2 + x^3),{x,0,34}],x] (* _James C. McMahon_, Jan 09 2024 *) %o A214699 (Magma) [n le 3 select 3*(1+(-1)^n)/2 else 3*Self(n-2) - Self(n-3): n in [1..40]]; // _G. C. Greubel_, Jan 08 2024 %o A214699 (SageMath) %o A214699 def a(n): # a=A214699 %o A214699 if (n<3): return 3*(n%2) %o A214699 else: return 3*a(n-2) - a(n-3) %o A214699 [a(n) for n in range(41)] # _G. C. Greubel_, Jan 08 2024 %Y A214699 Cf. A006053, A052931, A214683, A214778, A214779. %Y A214699 Cf. A214951, A214954, A217053, A217052, A217069. %K A214699 sign,easy %O A214699 0,2 %A A214699 _Roman Witula_, Jul 26 2012 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE