login
A090995
Number of meaningful differential operations of the n-th order on the space R^10.
13
10, 18, 32, 58, 104, 188, 338, 610, 1098, 1980, 3566, 6428, 11580, 20870, 37602, 67762, 122096, 220018, 396448, 714388, 1287266, 2319594, 4179738, 7531660, 13571542, 24455124, 44066548, 79405254, 143083226, 257827186, 464588384
OFFSET
1,1
COMMENTS
Also (starting 6,10,...) the number of zig-zag paths from top to bottom of a rectangle of width 6. - Joseph Myers, Dec 23 2008
Number of walks of length n on the path graph P_6. - Andrew Howroyd, Apr 17 2017
FORMULA
Equals 2 * A090990.
a(k+6) = 5*a(k+4) - 6*a(k+2) + a(k).
From Colin Barker, May 03 2012: (Start)
a(n) = a(n-1) + 2*a(n-2) - a(n-3).
G.f.: 2*x*(5+4*x-3*x^2)/(1-x-2*x^2+x^3). (End)
MAPLE
NUM := proc(k :: integer) local i, j, n, Fun, Identity, v, A; n := 10; # <- DIMENSION Fun := (i, j)->piecewise(((j=i+1) or (i+j=n+1)), 1, 0); Identity := (i, j)->piecewise(i=j, 1, 0); v := matrix(1, n, 1); A := piecewise(k>1, (matrix(n, n, Fun))^(k-1), k=1, matrix(n, n, Identity)); return(evalm(v&*A&*transpose(v))[1, 1]); end:
MATHEMATICA
a[n_ /; n <= 6] := {10, 18, 32, 58, 104, 188}[[n]]; a[n_] := a[n] = 5*a[n-2] - 6*a[n-4] + a[n-6]; Array[a, 31] (* Jean-François Alcover, Oct 07 2017 *)
2*LinearRecurrence[{1, 2, -1}, {5, 9, 16}, 40] (* G. C. Greubel, Feb 02 2019 *)
PROG
(PARI) my(x='x+O('x^40)); Vec(2*x*(5+4*x-3*x^2)/(1-x-2*x^2+x^3)) \\ G. C. Greubel, Feb 02 2019
(Magma) m:=40; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( 2*x*(5+4*x-3*x^2)/(1-x-2*x^2+x^3) )); // G. C. Greubel, Feb 02 2019
(Sage) a=(2*x*(5+4*x-3*x^2)/(1-x-2*x^2+x^3)).series(x, 40).coefficients(x, sparse=False); a[1:] # G. C. Greubel, Feb 02 2019
(GAP) a:=[10, 18, 32];; for n in [4..30] do a[n]:=a[n-1]+2*a[n-2]-a[n-3]; od; a; # G. C. Greubel, Feb 02 2019
CROSSREFS
Column 6 of A220062.
Sequence in context: A190038 A233695 A014006 * A363769 A153360 A189323
KEYWORD
nonn,easy
AUTHOR
Branko Malesevic, Feb 29 2004
EXTENSIONS
More terms from Joseph Myers, Dec 23 2008
STATUS
approved