%I #13 Nov 17 2015 14:29:43
%S 0,4,32,104,250,492,845,1349,2021,2871,3949,5267,6830,8698,10878,
%T 13370,16244,19502,23139,27235,31787,36785,42319,48381,54956,62144,
%U 69932,78300,87358,97088,107465,118609,130497,143099,156545,170807,185850,201814,218666
%N Number of ordered triples (w,x,y) with all terms in {-n,...,-1,1,...,n} and w+3x+3y>0.
%C For a guide to related sequences, see A211422.
%H Colin Barker, <a href="/A211625/b211625.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_08">Index entries for linear recurrences with constant coefficients</a>, signature (2,-1,2,-4,2,-1,2,-1).
%F 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).
%F G.f.: x*(7*x^6+23*x^5+48*x^4+66*x^3+44*x^2+24*x+4) / ((x-1)^4*(x^2+x+1)^2). - _Colin Barker_, Nov 17 2015
%t Remove["Global`*"];
%t t = Compile[{{u, _Integer}},
%t Module[{s = 0}, (Do[If[w + 3 x + 3 y > 0,
%t s = s + 1], {w, #}, {x, #}, {y, #}] &[
%t Flatten[{Reverse[-#], #} &[Range[1, u]]]]; s)]];
%t Map[t[#] &, Range[0, 60]] (* A211625 *)
%t FindLinearRecurrence[%]
%t (* _Peter J. C. Moses_, Apr 13 2012 *)
%t LinearRecurrence[{2, -1, 2, -4, 2, -1, 2, -1},{0, 4, 32, 104, 250, 492, 845, 1349},36] (* _Ray Chandler_, Aug 02 2015 *)
%o (PARI) concat(0, Vec(x*(7*x^6+23*x^5+48*x^4+66*x^3+44*x^2+24*x+4)/((x-1)^4*(x^2+x+1)^2) + O(x^100))) \\ _Colin Barker_, Nov 17 2015
%Y Cf. A211422.
%K nonn,easy
%O 0,2
%A _Clark Kimberling_, Apr 17 2012