(PARI)Lp(n)=我的(x='x+O('x^(n+4)));波尔科夫((1+x)*(1-6*x^2+7*x^4-(1-3*x^2)*sqrt(1-6x^2+5*x^4))/(2*x^4*(1-x)),n)\\A039660型
M(n)=我的(A);如果(n<1,0,n--;A=O(x);对于(k=0,n\2,A=1/(1-x-x^2/(1+x-x^2*A));极系数(A,n))\\A055879美元
N(N)=波尔科夫((1-x-sqrt(1-6*x+5*x^2+x^2*O(x^N)))/2,N+1)\\A002212年
b(n)=n(n+3)-9*n(n+2)+25*n(n+1)-21*n(n)+(M(n+3.)-M(n+2.)-3*M(n+1)+3*M(n)+Lp(n))/2;