OFFSET
1,5
COMMENTS
The empty subset is not allowed, otherwise we would get a(2)=1. - Alois P. Heinz, Sep 03 2009
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..100
FORMULA
a(n) is the coefficient of x^3y in product(x^(-2k)+x^k(y^k+y^(-k)), k=1..n) for n>2.
EXAMPLE
For n=5 we have splittings 4/23/15, 4/5/123, 13/5/24, so a(5)=3.
MAPLE
A113039:=proc(n) local i, j, p, t; t:= 0, 0; for j from 3 to n do p:=1; for i to j do p:=p*(x^(-2*i)+x^(i)*(y^i+y^(-i))); od; t:=t, coeff(coeff(p, x, 3), y, 1); od; t; end;
# second Maple program:
b:= proc() option remember; local i, j, t; `if` (args[1]=0, `if` (nargs=2, 1, b(args[t] $t=2..nargs)), add (`if` (args[j] -args[nargs] <0, 0, b(sort ([seq (args[i] -`if` (i=j, args[nargs], 0), i=1..nargs-1)])[], args[nargs]-1)), j=1..nargs-1)) end: a:= proc(n) local m; m:= n*(n+1)/2; `if` (n>2 and irem (m, 3)=0, b(m/3-1, m/3, m/3+1, n), 0) end: seq (a(n), n=1..42); # Alois P. Heinz, Sep 03 2009
MATHEMATICA
a[n_] := If[n <= 2, 0, Product[x^(-2k)+x^k(y^k+y^(-k)), {k, 1, n}] // SeriesCoefficient[#, {x, 0, 3}, {y, 0, 1}]&];
Table[an = a[n]; Print[n, " ", an]; an, {n, 1, 26}] (* Jean-François Alcover, Nov 17 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Floor van Lamoen, Oct 12 2005
EXTENSIONS
Extended beyond a(25) by Alois P. Heinz, Sep 03 2009
STATUS
approved