|
|
A203162
|
|
(n-1)-st elementary symmetric function of the first n terms of (1,2,3,1,2,3,1,2,3,...).
|
|
6
|
|
|
1, 3, 11, 17, 40, 132, 168, 372, 1188, 1404, 3024, 9504, 10800, 22896, 71280, 79056, 165888, 513216, 559872, 1166400, 3592512, 3872448, 8024832, 24634368, 26313984, 54307584, 166281984, 176359680, 362797056, 1108546560, 1169012736
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
Let esf abbreviate "elementary symmetric function". Then
0th esf of {1}: 1;
1st esf of {1,2}: 1+2=3;
2nd esf of {1,2,3} is 1*2+1*3+2*3=11.
|
|
MATHEMATICA
|
f[k_] := 1 + Mod[k + 2, 3]; t[n_] := Table[f[k], {k, n}]; a[n_] := SymmetricPolynomial[n - 1, t[n]]; Table[a[n], {n, 40}] (* A203162 *)
Rest[CoefficientList[Series[x*(1 + 3*x + 11*x^2 + 5*x^3 + 4*x^4)/(1 - 6*x^3)^2, {x, 0, 30}], x]] (* Vaclav Kotesovec, May 10 2019 *)
|
|
PROG
|
(PARI) my(x='x+O('x^40)); Vec(x*(1+3*x+11*x^2+5*x^3+4*x^4)/(1-6*x^3)^2) \\ G. C. Greubel, May 10 2019
(Magma) I:=[1, 3, 11, 17, 40, 132]; [n le 6 select I[n] else 12*Self(n-3) -36*Self(n-6): n in [1..40]]; // G. C. Greubel, May 10 2019
(Sage) a=(x*(1+3*x+11*x^2+5*x^3+4*x^4)/(1-6*x^3)^2).series(x, 40).coefficients(x, sparse=False); a[1:] # G. C. Greubel, May 10 2019
(GAP) a:=[1, 3, 11, 17, 40, 132];; for n in [7..40] do a[n]:=12*a[n-1]-36*a[n-2]-a[n-3]; od; a; # G. C. Greubel, May 10 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|