(PARI)/*从1+x=2*A(x)^3-A(x)^5开始:*/
{a(n)=局部(a=[1,1]);对于(i=1,n,a=concat(a,0);a[#a]=-Vec(2*Ser(a)^3-Ser(B)^5)[#a]);a[1]}
对于(n=0,25,打印1(a(n),“,”)
(PARI)/*系列逆转:*/
{a(n)=本地(a=1+serreverse(x-4*x^2-8*x^3-5*x^4-x^5+x^2*O(x^n));波尔科夫(a,n)}
对于(n=0,25,打印1(a(n),“,”)