zykov, благодарю за ссылки, буду разбираться!
Впрочем

- это перебор.
В какую память это влезет???
У меня на PARI/GP для

и

считается за две минуты. К сожалению мой алгоритм вычисляет не отдельные значения, а первые

значений, а также работает только для последовательностей

у которых

. Вот и он:
Код:
a(n)=n--; binomial(2*n, n)/(n+1)
f(n)=my(v, v1, v2); v=vector(n+2,i,if(i==1,vector(n+2,j,a(j)),vector(n-i+3,j,(j==1)))); v1=vector(n+1,i,0); v1[1]=a(2)-1; for(i=2,n+1,for(j=2,i,v[j][i-j+2]=v[j-1][i-j+3]-sum(k=1,j-1,v1[j-k]*v[k][i-j+2])); v1[i]=v[i][2]-sum(j=1,i-1,v1[j])-1); v1
b(n,m)=my(v, v1, v2, v3); v=f(n); v1=vector(n+1,i,1); v2=v1; v3=vector(n+1,i,0); v3[1]=1; for(i=1,n,for(j=1,n-i+1,v2[j]=v1[j+1]+sum(q=1,j,(m+1)^-(j-q+1)*v[j-q+1]*v1[q])); v1=v2; v3[i+1]=v1[1];); for(i=1,n,v3[i+1]*=(m+1)^i); v3
Null, благодарю за ответ! Буду пробовать и обязательно поделюсь результатом. Функцию

надо просто возводить в степень или же для подсчета ее коэффициентов есть какой-то хитрый ход?