|
%I #35 Sep 08 2022 08:44:36
%S 0,0,1,2,3,5,7,10,13,16,20,24,29,34,39,45,51,58,65,72,80,88,97,106,
%T 115,125,135,146,157,168,180,192,205,218,231,245,259,274,289,304,320,
%U 336,353,370,387,405,423,442,461,480,500,520,541,562,583,605,627,650,673
%N a(n) = floor((n^2 + 1)/5).
%C Without initial zeros, Molien series for 3-dimensional group [2+,n] = 2*(n/2).
%H Vincenzo Librandi, <a href="/A008738/b008738.txt">Table of n, a(n) for n = 0..5000</a>
%H <a href="/index/Mo#Molien">Index entries for Molien series</a>
%H <a href="/index/Rec#order_07">Index entries for linear recurrences with constant coefficients</a>, signature (2,-1,0,0,1,-2,1).
%F G.f.: x^2*(1+x^3)/((1-x)^2*(1-x^5)) = x^2*(1+x)*(1-x+x^2)/( (1-x)^3 *(1+x+x^2+x^3+x^4) ).
%F a(n+2)= A249020(n) + A249020(n-1). - _R. J. Mathar_, Aug 11 2021
%t Floor[(Range[0,60]^2 + 1)/5] (* _G. C. Greubel_, Aug 03 2019 *)
%o (PARI) a(n)=(n^2+1)\5;
%o (Magma) [(n^2+1) div 5: n in [0..60]]; // _Bruno Berselli_, Oct 28 2011
%o (Sage) [floor((n^2+1)/5) for n in (0..60)] # _G. C. Greubel_, Aug 03 2019
%o (GAP) List([0..60], n-> Int((n^2 + 1)/5)); # _G. C. Greubel_, Aug 03 2019
%Y Cf. A011858. Partial sums of A288156.
%K nonn,easy
%O 0,4
%A _N. J. A. Sloane_
%E More terms from Philip Mummert (s1280900(AT)cedarville.edu), May 10 2000
|
