login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(n) = 3*a(n-1) + 6*a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=26.
7

%I #35 Oct 03 2019 04:04:30

%S 2,5,26,110,491,2159,9533,42044,185489,818264,3609770,15924383,

%T 70250033,309906167,1367143082,6031116281,26606113502,117372181274,

%U 517784341115,2284192224491,10076654901437,44452902392372,196102828810229,865102555686356,3816377542312814

%N a(n) = 3*a(n-1) + 6*a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=26.

%C Ramanujan-type sequence number 4 for the argument 2*Pi/9 is defined by the following relation: 9^(1/3)*a(n)=(c(1)/c(2))^(n - 1/3) + (c(2)/c(4))^(n - 1/3) + (c(4)/c(1))^(n - 1/3), where c(j) := Cos(2Pi*j/9) - for the proof see Witula et al.'s papers. We have a(n)=bx(3n-1), where the sequence bx(n) and its two conjugate sequences ax(n) and cx(n) are defined in the comments to the sequence A214779. We note that ax(3n-1)=cx(3n-1)=0. Moreover we have ax(3n)=A214778(n), bx(3n)=cx(3n)=0 and cx(3n+1)=A214954(n), ax(3n+1)=bx(3n+1)=0.

%D 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. (in review)

%H Roman Witula, <a href="https://doi.org/10.1515/dema-2013-0418">Ramanujan Type Trigonometric Formulae</a>, Demonstratio Math. 45 (2012) 779-796.

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

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

%F a(n+1) - a(n) = A214778(n+1). - _Roman Witula_, Oct 06 2012

%e We have 2*9^(1/3) = (c(2)/c(1))^(1/3) + (c(4)/c(2))^(1/3) + (c(1)/c(4))^(1/3), 5*9^(1/3) = (c(1)/c(2))^(2/3) + (c(2)/c(4))^(2/3) + (c(4)/c(1))^(2/3), and 110*9^(1/3)=(c(1)/c(2))^(8/3) + (c(2)/c(4))^(8/3) + (c(4)/c(1))^(8/3). Moreover we obtain a(6)-a(2)-a(1)-a(0)=9500, a(12)-a(2)-a(1)-a(0)=70250000 and a(12)-a(6)=3^3*43*a(1)*a(3)^2. - _Roman Witula_, Oct 06 2012

%t LinearRecurrence[{3, 6, 1}, {2, 5, 26}, 40] (* _T. D. Noe_, Jul 30 2012 *)

%o (PARI) Vec((2-x-x^2)/(1-3*x-6*x^2-x^3)+O(x^99)) \\ _Charles R Greathouse IV_, Oct 01 2012

%Y Cf. A214778, A214779, A214954, A214699, A217053, A217052, A217069.

%K nonn,easy

%O 0,1

%A _Roman Witula_, Jul 30 2012