|
|
A213393
|
|
Number of (w,x,y) with all terms in {0,...,n} and 2*max(w,x,y) > 3*min(w,x,y).
|
|
3
|
|
|
0, 6, 24, 54, 108, 192, 300, 450, 648, 882, 1176, 1536, 1944, 2430, 3000, 3630, 4356, 5184, 6084, 7098, 8232, 9450, 10800, 12288, 13872, 15606, 17496, 19494, 21660, 24000, 26460, 29106, 31944, 34914, 38088, 41472, 45000, 48750, 52728
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Every term is even.
For a guide to related sequences, see A212959.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 2*a(n-1)-a(n-2)+2*a(n-3)-4*a(n-4)+2*a(n-5)-a(n-6)+2*a(n-7)-a(n-8).
G.f.: 6*x*(x^2+1)*(x+1)^2 / ((x^2+x+1)^2*(x-1)^4).
a(n) = (8*n^2+16*n+8-8*n*((2*n+2) mod 3)-8*((2*n+2) mod 3)+2*((2*n+2) mod 3)^2)/3. (End)
E.g.f.: 2*exp(-x/2)*(6*exp(3*x/2)*(1 + x*(13 + 2*x*(6 + x))) - 3*(2 + x)*cos(sqrt(3)*x/2) - sqrt(3)*(2 - 3*x)*sin(sqrt(3)*x/2))/27. - Stefano Spezia, Feb 25 2023
|
|
MATHEMATICA
|
t = Compile[{{n, _Integer}}, Module[{s = 0},
(Do[If[2*Max[w, x, y] > 3*Min[w, x, y], s = s + 1],
{w, 0, n}, {x, 0, n}, {y, 0, n}]; s)]];
m = Map[t[#] &, Range[0, 45]] (* A213393 *)
m/2 (* integers *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|