cp's OEIS Frontend

A322888 Chebyshev T-polynomials T_n(16).

%I A322888 #22 Sep 08 2022 08:46:23
%S A322888 1,16,511,16336,522241,16695376,533729791,17062657936,545471324161,
%T A322888 17438019715216,557471159562751,17821639086292816,569734979601807361,
%U A322888 18213697708171542736,582268591681887560191,18614381236112230383376,595077930963909484707841
%N A322888 Chebyshev T-polynomials T_n(16).
%H A322888 Colin Barker, <a href="/A322888/b322888.txt">Table of n, a(n) for n = 0..600</a>
%H A322888 Wikipedia, <a href="https://en.wikipedia.org/wiki/Chebyshev_polynomials">Chebyshev polynomials</a>.
%H A322888 <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%H A322888 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (32, -1).
%F A322888 a(0) = 1, a(1) = 16 and a(n) = 32*a(n-1) - a(n-2) for n > 1.
%F A322888 From _Colin Barker_, Dec 30 2018: (Start)
%F A322888 G.f.: (1 - 16*x) / (1 - 32*x + x^2).
%F A322888 a(n) = ((16+sqrt(255))^(-n) * (1+(16+sqrt(255))^(2*n))) / 2.
%F A322888 (End)
%p A322888 seq(coeff(series((1-16*x)/(1-32*x+x^2),x,n+1), x, n), n = 0 .. 20); # _Muniru A Asiru_, Dec 31 2018
%t A322888 Array[ChebyshevT[#, 16] &, 17, 0] (* or *)
%t A322888 With[{k = 16}, CoefficientList[Series[(1 - k x)/(1 - 2 k x + x^2), {x, 0, 16}], x]] (* _Michael De Vlieger_, Jan 01 2019 *)
%o A322888 (PARI) {a(n) = polchebyshev(n, 1, 16)}
%o A322888 (PARI) Vec((1 - 16*x) / (1 - 32*x + x^2) + O(x^20)) \\ _Colin Barker_, Dec 30 2018
%o A322888 (GAP) a:=[1,16];; for n in [3..20] do a[n]:=32*a[n-1]-a[n-2]; od; Print(a); # _Muniru A Asiru_, Dec 31 2018
%o A322888 (Magma)  I:=[1, 16]; [n le 2 select I[n] else 32*Self(n-1)-Self(n-2): n in [1..20]]; // _Vincenzo Librandi_, Jan 02 2019
%Y A322888 Column 16 of A322836.
%K A322888 nonn,easy
%O A322888 0,2
%A A322888 _Seiichi Manyama_, Dec 29 2018