login
a(n) = n^3*(n^10 + 1)/2.
1

%I #16 Sep 08 2022 08:45:48

%S 0,1,4100,797175,33554464,610351625,6530347116,48444505375,

%T 274877907200,1270932914529,5000000000500,17261356072631,

%U 53496602690400,151437553297225,396857386628444,973097534181375,2251799813687296,4952289016455425,10411482432838500

%N a(n) = n^3*(n^10 + 1)/2.

%H Vincenzo Librandi and T. D. Noe, <a href="/A168191/b168191.txt">Table of n, a(n) for n = 0..1000</a> (First 313 terms from Vincenzo Librandi)

%F G.f.: x (1 + 4086*x + 739866*x^2 + 22766750*x^3 + 211640655*x^4 + 752810988*x^5 + 1137585708*x^6 + 752810988*x^7 + 211640655*x^8 + 22766750*x^9 + 739866*x^10 + 4086*x^11 + x^12)/(1-x)^14. - _Vincenzo Librandi_, Jul 16 2016

%t Table[n^3*(n^10 + 1)/2, {n,0,50}] (* _G. C. Greubel_, Jul 15 2016 *)

%t CoefficientList[Series[x (1 + 4086 x + 739866 x^2 + 22766750 x^3 + 211640655 x^4 + 752810988 x^5 + 1137585708 x^6 + 752810988 x^7 + 211640655 x^8 + 22766750 x^9 + 739866 x^10 + 4086 x^11 + x^12)/(1-x)^14, {x, 0, 30}], x] (* _Vincenzo Librandi_, Jul 16 2016 *)

%o (Magma) [n^3*(n^10+1)/2: n in [0..50]]; // _Vincenzo Librandi_, Apr 25 2011

%K nonn

%O 0,3

%A _N. J. A. Sloane_, Dec 11 2009