login
a(n) = floor(n/5)*floor((n+1)/5).
2

%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_