cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A106221 Self-convolution 4th power equals A106220, which consists entirely of digits {0,1,2,3} after the initial terms {1,4}.

This page as a plain text file.
%I A106221 #7 Mar 14 2015 09:57:56
%S A106221 1,1,-1,2,-4,10,-26,71,-199,569,-1652,4855,-14413,43153,-130143,
%T A106221 394967,-1205268,3695771,-11381215,35183209,-109138163,339599993,
%U A106221 -1059702401,3315256789,-10396158911,32671424776,-102879610571,324557399534,-1025643986057,3246330348415,-10290418283163
%N A106221 Self-convolution 4th power equals A106220, which consists entirely of digits {0,1,2,3} after the initial terms {1,4}.
%F A106221 Limit a(n+1)/a(n) = -3.30697774878897620974321728382452592372871...
%e A106221 A(x) = 1 + x - x^2 + 2*x^3 - 4*x^4 + 10*x^5 - 26*x^6 + 71*x^7 -+...
%e A106221 A(x)^4 = 1 + 4*x + 2*x^2 + 3*x^4 + 2*x^6 + x^8 + 2*x^14 +...
%e A106221 A106220 = {1,4,2,0,3,0,2,0,1,0,0,0,0,0,2,0,0,0,2,...}.
%o A106221 (PARI) {a(n)=local(A=1+4*x);if(n==0,1, for(j=1,n, for(k=0,3,t=polcoeff((A+k*x^j+x*O(x^j))^(1/4),j); if(denominator(t)==1,A=A+k*x^j;break))); return(polcoeff((A+x*O(x^n))^(1/4),n)))}
%Y A106221 Cf. A106220, A106219, A106223, A106225.
%K A106221 sign,base
%O A106221 0,4
%A A106221 _Paul D. Hanna_, May 01 2005