(PARI){a(n)=局部(a=1+x);对于(i=1,n,a=和(m=0,n,x^m*a^m/prod(k=1,m,1-k*x+x*O(x^n));polcoeff(a,n)}
对于(n=0,20,打印1(a(n),“,”)
(PARI){斯特林2(n,k)=n!*polceoff(((exp(x+x*O(x^n))-1)^k)/k!,n)}
{a(n)=局部(a=1+x);对于(i=0,n,a=和(m=0,n,x^m*和(k=0,m,斯特林2(m,k)*(a+x*O(x^n)));波尔科夫(a,n)}
对于(n=0,20,打印1(a(n),“,”)