%I #9 Nov 04 2012 00:54:19
%S 1,1,1,2,1,1,2,2,1,1,1,0,2,2,2,2,1,1,1,2,1,1,0,0,2,2,2,0,2,2,2,2,1,1,
%T 1,2,1,1,2,2,1,1,1,0,0,0,0,0,2,2,2,0,2,2,0,0,2,2,2,0,2,2,2,2,1,1,1,2,
%U 1,1,2,2,1,1,1,0,2,2,2,2,1,1,1,2,1,1,0,0,0,0,0,0,0,0,0,0
%N a(n) = A183161(n) (mod 4), n>=0.
%C Conjecture: a(n) never equals 3.
%C A183161(n) is defined by the convolution:
%C Sum_{k=0..n} A183161(n-k)*A183161(k) = Sum_{k=0..n} C(n+k,n-k)*C(2*n-k,k).
%C The g.f. F(x) of A183161 satisfies: F(x) = 1/sqrt(1 - 2*x*G(x)^2 - 3*x^2*G(x)^4), where G(x) = 1 + x*G(x)^3 = g.f. of A001764.
%H Paul D. Hanna, <a href="/A218622/b218622.txt">Table of n, a(n) for n = 0..2048</a>
%e Formatting the terms into groups of 8 reveals complex binary patterns:
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 0,0,0,0,0,0,0,0,
%e 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 0,0,0,0,0,0,0,0,
%e 2,2,2,0,2,2,0,0, 2,2,2,0,0,0,0,0, 2,2,2,0,2,2,0,0, 2,2,2,0,2,2,2,2,
%e 1,1,1,2,1,1,2,2, 1,1,1,0,2,2,2,2, 1,1,1,2,1,1,0,0, 2,2,2,0,2,2,2,2, ...
%o (PARI) {a(n)=local(A2=sum(m=0, n, sum(k=0, m, binomial(m+k, m-k)*binomial(2*m-k, k))*x^m+x*O(x^n))); polcoeff(A2^(1/2), n)%4}
%o (PARI) {a(n)=local(G=1); for(i=0, n, G=1+x*G^3+O(x^(n+1))); polcoeff(1/sqrt(1-2*x*G^2-3*x^2*G^4), n)%4}
%o (PARI) /* Using Central Trinomial Coefficients A002426: */
%o {A002426(n)=sum(k=0, n\2, binomial(n, 2*k)*binomial(2*k, k))}
%o {a(n)=if(n==0, 1, sum(k=0, n, A002426(k)*binomial(3*n-k, n-k)*2*k/(3*n-k)))%4}
%o /* Format Print of a(n) into 4 columns of 8 terms each: */
%o for(n=0,1024,if(n>0,if(n%32==0,print(""),if(n%8==0,print1(" "))));print1(a(n),","))
%Y Cf. A183161.
%K nonn
%O 0,4
%A _Paul D. Hanna_, Nov 03 2012
|