|
%I #21 Sep 08 2022 08:44:35
%S 0,0,0,0,0,1,1,1,1,2,4,4,4,4,6,9,9,9,9,12,16,16,16,16,20,25,25,25,25,
%T 30,36,36,36,36,42,49,49,49,49,56,64,64,64,64,72,81,81,81,81,90,100,
%U 100,100,100,110,121,121,121,121
%N a(n) = floor(n/5)*floor((n+1)/5).
%H G. C. Greubel, <a href="/A008497/b008497.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_11">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,0,0,2,-2,0,0,0,-1,1).
%F From _R. J. Mathar_, Apr 16 2010: (Start)
%F a(n) = A002266(n)*A002266(n+1).
%F a(n)= a(n-1) + 2*a(n-5) - 2*a(n-6) - a(n-10) + a(n-11).
%F G.f.: x^5*(1+x^4)/ ((x^4+x^3+x^2+x+1)^2 * (1-x)^3). (End)
%p seq( mul(floor((n+j)/5), j=0..1), n=0..55); # _G. C. Greubel_, Nov 08 2019
%t Times@@@Partition[Floor[Range[0,60]/5],2,1] (* or *) LinearRecurrence[ {1,0,0,0,2,-2,0,0,0,-1,1},{0,0,0,0,0,1,1,1,1,2,4},60] (* _Harvey P. Dale_, Feb 01 2015 *)
%t Product[Floor[(Range[55] +j-1)/5], {j,0,1}] (* _G. C. Greubel_, Nov 08 2019 *)
%o (PARI) a(n) = (n\5)*((n+1)\5); \\ _Michel Marcus_, Jan 06 2017
%o (PARI) vector(56, n, prod(j=0,1, (n+j-1)\5) ) \\ _G. C. Greubel_, Nov 08 2019
%o (Magma) [&*[Floor((n+j)/5): j in [0..1]]: n in [0..55]]; // _G. C. Greubel_, Nov 08 2019
%o (Sage) [product(floor((n+j)/5) for j in (0..1)) for n in (0..55)] # _G. C. Greubel_, Nov 08 2019
%o (GAP) List([0..55], n-> Int(n/5)*Int((n+1)/5) ); # _G. C. Greubel_, Nov 08 2019
%Y Cf. A002266.
%K nonn
%O 0,10
%A _N. J. A. Sloane_
|
