%I#9 2017年10月29日21:02:23
%S 1,3,0,1,0,0,-1,0.0,-1,0,0-1,0,0,0-0,-1,1,0,1.0,0,-1,
%T 0,0,-1,0,0,0-0,0,1,0,0,
%U 0,0,0,-1,0,O,-1,0,0,0,1,0,-1,0,0-0,0,0,0 0,1,0,0,0-0,0,-1,0,0
%g.f.的N个系数A(x),其中-1<=A(N)<=1,对于所有N>1,具有初始项{1,3},使得A(x)^(1/3)完全由整数系数组成。
%C A(z)=0,z=-0.3218604126 5330206348 1666946119 2701743677 6909817084 3826086189 5315539535 3583969883。。。
%H Paul D.Hanna,n表,n=0..1001时的a(n)</a>
%e.G.f.:A(x)=1+3*x+x^3-x^6-x^9-x^12-x^18+x^21-x^24-x^30-x^33+x^39-x^42-x^45+x^48+x^54-x^60-x^66-x^69+x^72-x^75+x^84-x^87-x^90-x^93-x^96-x^102+x^108-x^114-x^126+x^132-x^138+x^141-x^144-x^147-x^153-x^ 159-x^162+x^165-x^168-x^171+x^174-x^180-x^183+x^186-x^189-x ^192+x ^195+。。。
%e其中
%e A(x)^(1/3)=1+x-x ^2+2*x ^3-4*x ^4+9*x ^5-22*x ^6+55*x ^7-142*x ^8+375*x ^9-1009*x ^10+2753*x ^11-7599*x ^12+21178*x ^13-59509*x ^14+168401*x ^15-4777*x ^16+1372536*x ^17-3947678*x ^18+11402376*x ^19-33059314*x ^20+96177750*x^21+…+A196307(n)*x ^n+。。。
%o(PARI){a(n)=局部(a=1+3*x);如果(n=0=0,1,如果(n%3==0,对于(j=1,n,对于(k=-1,1,t=polcoeff((a+k*x^j+x*o(x^j))^(1/3),j);
%o如果(分母(t)==1,A=A+k*x^j;断裂));波尔科夫(A+x*O(x^n),n))}
%Y参考A196307(立方根),A196308(三等分);变体:A106216、A083953。
%K符号
%0、2
%A·保罗·D·汉纳,2011年10月1日
|