%I #30 Jun 18 2017 02:16:30
%S 1,2,257,6562,65537,390626,1679617,5764802,16777217,43046722,
%T 100000001,214358882,429981697,815730722,1475789057,2562890626,
%U 4294967297,6975757442,11019960577,16983563042,25600000001,37822859362,54875873537
%N n^8 + 1.
%C a(n) = Phi_16(n) where Phi_k(x) is the k-th cyclotomic polynomial.
%H Harry J. Smith, <a href="/A060890/b060890.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Cy#CyclotomicPolynomialsValuesAtX">Index to values of cyclotomic polynomials of integer argument</a>
%H <a href="/index/Rec#order_09">Index entries for linear recurrences with constant coefficients</a>, signature (9,-36,84,-126,126,-84,36,-9,1).
%F a(0)=1, a(1)=2, a(2)=257, a(3)=6562, a(4)=65537, a(5)=390626, a(6)=1679617, a(7)=5764802, a(8)=16777217, a(n)=9*a(n-1)-36*a(n-2)+84*a(n-3)-126*a(n-4)+ 126*a(n-5)-84*a(n-6)+36*a(n-7)-9*a(n-8)+a(n-9). - _Harvey P. Dale_, Mar 12 2013
%F Sum_{n>=0} 1/a(n) = 1/2 + Pi*((sqrt(2 + sqrt(2)) * sin(sqrt(2 + sqrt(2))*Pi) + sqrt(2 - sqrt(2)) * sinh(sqrt(2 - sqrt(2))*Pi)) / (cosh(sqrt(2 - sqrt(2))*Pi) - cos(sqrt(2 + sqrt(2))*Pi)) + (sqrt(2 - sqrt(2)) * sin(sqrt(2 - sqrt(2))*Pi) + sqrt(2 + sqrt(2)) * sinh(sqrt(2 + sqrt(2))*Pi)) / (cosh(sqrt(2 + sqrt(2))*Pi) - cos(sqrt(2 - sqrt(2))*Pi))) / 8 = 1.5040621333147995112929... . - _Vaclav Kotesovec_, Feb 14 2015
%F Sum_{n>=0} (-1)^n/a(n) = 1/2 + Pi*((sqrt(2 - sqrt(2)) * sin(sqrt(2 - sqrt(2))*Pi/2) - sqrt(2 + sqrt(2)) * sinh(sqrt(2 + sqrt(2))*Pi/2)) / (cos(sqrt(2 - sqrt(2))*Pi/2) + cosh(sqrt(2 + sqrt(2))*Pi/2)) - (sqrt(2 - sqrt(2)) * sin(sqrt(2 - sqrt(2))*Pi/2) + sqrt(2 + sqrt(2)) * sinh(sqrt(2 + sqrt(2))*Pi/2)) / (cos(sqrt(2 - sqrt(2))*Pi/2) - cosh(sqrt(2 + sqrt(2))*Pi/2)) + (sqrt(2 + sqrt(2)) * sin(sqrt(2 + sqrt(2))*Pi/2) - sqrt(2 - sqrt(2)) * sinh(sqrt(2 - sqrt(2))*Pi/2)) / (cos(sqrt(2 + sqrt(2))*Pi/2) + cosh(sqrt(2 - sqrt(2))*Pi/2)) - (sqrt(2 + sqrt(2)) * sin(sqrt(2 + sqrt(2))*Pi/2) + sqrt(2 - sqrt(2)) * sinh(sqrt(2 - sqrt(2))*Pi/2)) / (cos(sqrt(2 + sqrt(2))*Pi/2) - cosh(sqrt(2 - sqrt(2))*Pi/2)))/16 = 0.5037518217314416642671664241... . - _Vaclav Kotesovec_, Feb 14 2015
%e G.f.: (1-7*x+275*x^2+4237*x^3+15689*x^4+15563*x^5+4321*x^6+239*x^7+2*x^8)/ (1-x)^9. - _Colin Barker_, Apr 21 2012
%p A060890 := proc(n)
%p numtheory[cyclotomic](16,n) ;
%p end proc:
%p seq(A060890(n),n=0..20) ; # _R. J. Mathar_, Feb 11 2014
%t Table[n^8+1,{n,0,40}] (* _Vladimir Joseph Stephan Orlovsky_, Apr 15 2011 *)
%t LinearRecurrence[{9,-36,84,-126,126,-84,36,-9,1},{1,2,257,6562,65537,390626,1679617,5764802,16777217},40] (* _Harvey P. Dale_, Mar 12 2013 *)
%o (PARI) for (n=0, 1000, write("b060890.txt", n, " ", n^8 + 1)) \\ _Harry J. Smith_, Jul 14 2009
%Y Cf. A002522, A001093, A002523, A002561, A002604.
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_, May 05 2001