(PARI){a(n)=局部(A1,A2,A3);如果(n<1,5*(n==0),A1=总和(k=0,n,(5*k)!/k!^5*(-x)^k,x*O(x^n));A2=-x*exp);A3=subst(5/A1^2/(1+5^5*x)/(x*A2'/A2)^3,x,serreverse(A2));sumdiv(n,k,moebius(n/k)*polcoff(A3,k))/n^3)}/*迈克尔·索莫斯2004年3月27日*/
(PARI)
cumsum(v)=(i=2,#v,v[i]+=v[i-1]);v;
my(x='x+O('x^(N+1)),h=cumsum(向量(5*N,N,1/N)),
y0=总和(n=0,n,(5*n)/n!^5*x^n),
y1=5*和(n=1,n,(5*n)/n!^5*(h[5*n]-h[n])*x^n),
Qx=x*exp(y1/y0),Xq=serreverse(Qx));
Vec(5*(x*Xq'/Xq)^3/((1-3125*Xq)*sqr(subst(y0,'x,Xq)));
};
序列(N)={
v2=dirmul(矢量(N,N,moebius(N)),矢量(N、N,v1[N+1]));
concat(5,向量(v2,n,v2[n]/n^3));
};