OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
MAPLE
b:= proc(n, k) option remember; `if`(n<2, k*n, `if`(n::odd, 0,
(t-> t*(1-t)/2)(b(n/2, k)))+add(b(i, k)*b(n-i, k), i=1..n/2))
end:
a:= n-> add(add((-1)^i*binomial(k, i)*b(n, k-i), i=0..k), k=0..n):
seq(a(n), n=1..23); # Alois P. Heinz, Sep 07 2019
MATHEMATICA
b[n_, k_] := b[n, k] = If[n < 2, k n, If[OddQ[n], 0, Function[t, t(1 - t)/2][b[n/2, k]]] + Sum[b[i, k] b[n - i, k], {i, 1, n/2}]];
a[n_] := Sum[Sum[(-1)^i Binomial[k, i] b[n, k - i], {i, 0, k}], {k, 0, n}];
Array[a, 23] (* Jean-François Alcover, Apr 10 2020, after Alois P. Heinz *)
PROG
(PARI) \\ here R(n, k) is k-th column of A319539 as a vector.
R(n, k)={my(v=vector(n)); v[1]=k; for(n=2, n, v[n]=sum(j=1, (n-1)\2, v[j]*v[n-j]) + if(n%2, 0, binomial(v[n/2]+1, 2))); v}
seq(n)={sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) )}
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Sep 23 2018
STATUS
approved