%I #32 Sep 08 2022 08:44:36
%S 1,65,793,4825,19721,62281,164305,379793,793585,1531441,2771561,
%T 4757545,7812793,12356345,18920161,28167841,40914785,58149793,
%U 81058105,111045881,149766121,199146025,261415793,339138865,435243601,553056401,696336265,869310793,1076713625,1323823321
%N 6-dimensional centered cube numbers.
%C These are never prime, as a(n) = (2*n^2 + 2*n + 1) * (n^4 + 2*n^3 + 5*n^2 + 4*n + 1). - _Jonathan Vos Post_, Aug 17 2011
%H Vincenzo Librandi, <a href="/A008516/b008516.txt">Table of n, a(n) for n = 0..10000</a>
%H <a href="/index/Rec#order_07">Index entries for linear recurrences with constant coefficients</a>, signature (7,-21,35,-35,21,-7,1).
%F From _Colin Barker_, Jul 09 2012: (Start)
%F G.f.: (1 + 58*x + 359*x^2 + 604*x^3 + 359*x^4 + 58*x^5 + x^6)/(1-x)^7.
%F a(n) = 1 + 6*n + 15*n^2 + 20*n^3 + 15*n^4 + 6*n^5 + 2*n^6. (End)
%F E.g.f.: (1 +64*x +332*x^2 +440*x^3 +205*x^4 +36*x^5 +2*x^6)*exp(x). - _G. C. Greubel_, Nov 09 2019
%p seq(n^6+(n+1)^6, n=0..35);
%t Table[n^6 + (n+1)^6, {n,0,35}] (* _Alonso del Arte_, Aug 17 2011 *)
%t LinearRecurrence[{7,-21,35,-35,21,-7,1},{1,65,793,4825,19721,62281,164305},30] (* _Harvey P. Dale_, Jun 19 2021 *)
%o (Magma) [(n+1)^6+n^6: n in [0..35]]; // _Vincenzo Librandi_, Aug 27 2011
%o (PARI) vector(36, n, n^6+(n-1)^6) \\ _G. C. Greubel_, Nov 09 2019
%o (Sage) [n^6+(n+1)^6 for n in (0..35)] # _G. C. Greubel_, Nov 09 2019
%o (GAP) List([0..35], n-> n^6+(n+1)^6); # _G. C. Greubel_, Nov 09 2019
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_
|