(PARI){a(n)=局部(a=1+x);对于(i=1,n,a=和(m=0,n,x^m/prod(k=1,m,1-k*x*a+x*O(x^n));polcoff(a,n)}
对于(n=0,30,打印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))^(m-k));波尔科夫(a,n)}
对于(n=0,20,打印1(a(n),“,”)