login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A212903 Number of (w,x,y,z) with all terms in {0,...,n} and |w-x|<=|x-y|<=|y-z|. 2

%I #8 Jun 13 2015 00:54:14

%S 1,8,35,94,209,398,697,1130,1743,2568,3661,5060,6833,9024,11711,14946,

%T 18817,23382,28741,34958,42143,50368,59753,70376,82369,95816,110859,

%U 127590,146161,166670,189281,214098,241295,270984,303349,338508

%N Number of (w,x,y,z) with all terms in {0,...,n} and |w-x|<=|x-y|<=|y-z|.

%C For a guide to related sequences, see A211795.

%H <a href="/index/Rec">Index entries for linear recurrences with constant coefficients</a>, signature (2,1,-3,-1,1,3,-1,-2,1).

%F a(n) = 2*a(n-1)+a(n-2)-3*a(n-3)-a(n-4)+a(n-5)+3*a(n-6)-a(n-7)-2*a(n-8)+a(n-9).

%F G.f.: (1 + 6*x + 18*x^2 + 19*x^3 + 11*x^4 - 2*x^5 - 2*x^6 + x^7 )/(1 - 2*x - x^2 + 3*x^3 + x^4 - x^5 - 3*x^6 + x^7 + 2*x^8 - x^9).

%t t = Compile[{{n, _Integer}}, Module[{s = 0},(Do[If[Abs[w - x] <= Abs[x - y] <= Abs[y - z], s = s + 1],

%t {w, 0, n}, {x, 0, n}, {y, 0, n}, {z, 0, n}]; s)]];

%t m = Map[t[#] &, Range[0, 40]] (* A212903 *)

%t LinearRecurrence[{2, 1, -3, -1, 1, 3, -1, -2, 1}, {1, 8, 35, 94, 209, 398, 697, 1130, 1743}, 40]

%Y Cf. A211795.

%K nonn,easy

%O 0,2

%A _Clark Kimberling_, Jun 01 2012

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 18 17:36 EDT 2024. Contains 371781 sequences. (Running on oeis4.)