a(n)=sum(k=1..n-1, (sum(j=0..k, binomial(k,j)*sum(i=j..n-k+j-1, binomial(j,i-j)*binomial(k-j,n-3*(k-j)-i-1))))*binomial(n+k-1,n-1)/n, n>1. a(1)=1. [_Vladimir Kruchinin kru(AT)ie.tusur.ru, _, May 10 2011]
a(n):=if n<2 then n else sum((sum(binomial(k, j)*sum(binomial(j, i-j)*binomial(k-j, n-3*(k-j)-i-1), i, j, n-k+j-1), j, 0, k))*binomial(n+k-1, n-1), k, 1, n-1)/n; [_Vladimir Kruchinin kru(AT)ie.tusur.ru, _, May 10 2011]
Vec(serreverse(x-sum(j=2, 6, x^j))) /* show terms */ /* _Joerg Arndt, _, May 12 2011 */
approved
editing