%I #24 Feb 26 2024 14:27:06
%S 1,2,65537,43046722,4294967297,152587890626,2821109907457,
%T 33232930569602,281474976710657,1853020188851842,10000000000000001,
%U 45949729863572162,184884258895036417,665416609183179842,2177953337809371137,6568408355712890626
%N n^16 + 1.
%C a(n) = Phi_32(n) where Phi_k is the k-th cyclotomic polynomial.
%H Harry J. Smith, <a href="/A060895/b060895.txt">Table of n, a(n) for n = 0..1000</a>
%H Hisanori Mishima, <a href="http://www.asahi-net.or.jp/~KC2H-MSM/mathland/matha1/matha103.htm">Factorizations of many number sequences</a>
%H <a href="/index/Cy#CyclotomicPolynomialsValuesAtX">Index to values of cyclotomic polynomials of integer argument</a>
%H <a href="/index/Rec#order_17">Index entries for linear recurrences with constant coefficients</a>, signature (17,-136,680,-2380,6188,-12376,19448,-24310,24310,-19448,12376,-6188,2380,-680,136,-17,1).
%F G.f.: -(2*x^16 +65503*x^15 +41932865*x^14 +3572084695*x^13 +85383240369*x^12 +782115513931*x^11 +3207483186165*x^10 +6382798914035*x^9 +6382798938345*x^8 +3207483166717*x^7 +782115526307*x^6 +85383234181*x^5 +3572087075*x^4 +41932185*x^3 +65639*x^2 -15*x +1)/(x -1)^17. [_Colin Barker_, Oct 29 2012]
%p A060895 := proc(n)
%p numtheory[cyclotomic](32,n) ;
%p end proc:
%p seq(A060895(n),n=0..20) ; # _R. J. Mathar_, Feb 11 2014
%t Range[0, 25]^16 + 1 (* or *)
%t Cyclotomic[32, Range[0, 25]] (* _Paolo Xausa_, Feb 26 2024 *)
%o (PARI) { for (n=0, 1000, write("b060895.txt", n, " ", n^16 + 1); ) } \\ _Harry J. Smith_, Jul 14 2009
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_, May 05 2001