A120880
G.f. satisfies: A(x) = A(x^3)*(1 + 2*x + x^2); thus a(n) = 2^A062756(n), where A062756(n) is the number of 1's in the ternary expansion of n.
Original entry on oeis.org
1, 2, 1, 2, 4, 2, 1, 2, 1, 2, 4, 2, 4, 8, 4, 2, 4, 2, 1, 2, 1, 2, 4, 2, 1, 2, 1, 2, 4, 2, 4, 8, 4, 2, 4, 2, 4, 8, 4, 8, 16, 8, 4, 8, 4, 2, 4, 2, 4, 8, 4, 2, 4, 2, 1, 2, 1, 2, 4, 2, 1, 2, 1, 2, 4, 2, 4, 8, 4, 2, 4, 2, 1, 2, 1, 2, 4, 2, 1, 2, 1, 2, 4, 2, 4, 8, 4, 2, 4, 2, 4, 8, 4, 8, 16, 8, 4, 8, 4, 2, 4, 2, 4, 8, 4
Offset: 0
Records are 2^n at positions: 0,1,4,13,40,121,...,(3^n-1)/2,... (n>=0).
A(x) = 1 + 2*x + x^2 + 2*x^3 + 4*x^4 + 2*x^5 + x^6 + 2*x^7 + x^8 +...
-
a120880 n = sum $ map (a039966 . (n -)) $ takeWhile (<= n) a005836_list
-- Reinhard Zumkeller, Jan 28 2012
-
Nest[ Join[#, 2 #, #] &, {1}, 5] (* Robert G. Wilson v, Jul 27 2014 *)
-
a(n)=if(n==0,1,a(n\3)*2^((n%3)%2))
A117939
Triangle related to powers of 3 partitions of n.
Original entry on oeis.org
1, 2, 1, 1, -2, 1, 2, 0, 0, 1, 4, 2, 0, 2, 1, 2, -4, 2, 1, -2, 1, 1, 0, 0, -2, 0, 0, 1, 2, 1, 0, -4, -2, 0, 2, 1, 1, -2, 1, -2, 4, -2, 1, -2, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 4, 2, 0, 0, 0, 0, 0, 0, 0, 2, 1, 2, -4, 2, 0, 0, 0, 0, 0, 0, 1, -2, 1, 4, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 0, 1, 8, 4, 0, 4, 2, 0, 0, 0, 0, 4, 2, 0, 2, 1
Offset: 0
Triangle begins
1;
2, 1;
1, -2, 1;
2, 0, 0, 1;
4, 2, 0, 2, 1;
2, -4, 2, 1, -2, 1;
1, 0, 0, -2, 0, 0, 1;
2, 1, 0, -4, -2, 0, 2, 1;
1, -2, 1, -2, 4, -2, 1, -2, 1;
-
T[n_, k_]:= Sum[JacobiSymbol[Binomial[n, j], 3]*JacobiSymbol[Binomial[n-j, k], 3], {j, 0, n}]; Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 29 2021 *)
-
T(n,k)=(matrix(n+1,n+1,r,c,(binomial(r-1,c-1)+1)%3-1)^2)[n+1,k+1] \\ Paul D. Hanna, Jul 08 2006
-
def A117939(n, k): return sum(jacobi_symbol(binomial(n, j), 3)*jacobi_symbol(binomial(n-j, k), 3) for j in (0..n))
flatten([[A117939(n, k) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Oct 29 2021
A117963
Antidiagonal sums of a Legendre-binomial triangle for p = 3.
Original entry on oeis.org
1, 1, 2, 0, 2, 2, 1, 3, 4, -2, 2, 0, 2, 2, 4, 0, 4, 4, -1, 3, 2, 2, 4, 6, 1, 7, 8, -6, 2, -4, 4, 0, 4, -2, 2, 0, 2, 2, 4, 0, 4, 4, 2, 6, 8, -4, 4, 0, 4, 4, 8, 0, 8, 8, -5, 3, -2, 4, 2, 6, -1, 5, 4, 0, 4, 4, 2, 6, 8, 2, 10, 12, -5, 7, 2, 6, 8, 14, 1, 15, 16, -14, 2, -12, 8, -4, 4, -6, -2, -8, 8, 0, 8, -4, 4, 0, 4, 4, 8, -6, 2
Offset: 0
The triangle mentioned in the name starts:
{1},
{1, 1},
{1,-1, 1},
{1, 0, 0, 1},
{1, 1, 0, 1, 1},
{1,-1, 1, 1,-1, 1},
{1, 0, 0,-1, 0, 0, 1}.
-
a[n_] := Sum[JacobiSymbol[Binomial[n - k, k], 3], {k, 0, n/2}];
a /@ Range[0, 100] (* Jean-François Alcover, Oct 14 2019 *)
-
{a(n)=local(A=1+x+x*O(x^n)); for(i=1,#binary(n), A=subst(A,x,x^3+x*O(x^n)) *(1-4*x^3-x^6)/(1-x-x^2+x*O(x^n))); polcoeff(A,n,x)} \\ Paul D. Hanna, Jul 11 2006
-
A117963list(upto_n) = { my(A=1+x+x*O(x^upto_n)); for(i=1, #binary(upto_n), print1(i,", "); A=subst(A, x, x^3+x*O(x^upto_n)) *(1-4*x^3-x^6)/(1-x-x^2+x*O(x^upto_n))); print(); vector(upto_n,n,polcoeff(A,n-1)); }; \\ Antti Karttunen, Jan 01 2023, after Paul D. Hanna's program above.
-
A117963(n) = sum(k=0,n\2,kronecker(binomial(n-k,k),3)); \\ Antti Karttunen, Jan 01 2023
A120854
Matrix log of A117939, read by rows, consisting only of 0's, 3's and signed 2's.
Original entry on oeis.org
0, 2, 0, 3, -2, 0, 2, 0, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 3, -2, 0, 3, 0, 0, -2, 0, 0, 0, 0, 3, 0, 0, -2, 0, 2, 0, 0, 0, 3, 0, 0, -2, 3, -2, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 3, -2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 2, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 3, -2, 0, 0, 0, 0, 0
Offset: 0
Triangle begins:
0;
2, 0;
3,-2, 0;
2, 0, 0, 0;
0, 2, 0, 2, 0;
0, 0, 2, 3,-2, 0;
3, 0, 0,-2, 0, 0, 0;
0, 3, 0, 0,-2, 0, 2, 0;
0, 0, 3, 0, 0,-2, 3,-2, 0;
2, 0, 0, 0, 0, 0, 0, 0, 0, 0; ...
Matrix exponentiation gives A117939:
1;
2, 1;
1,-2, 1;
2, 0, 0, 1;
4, 2, 0, 2, 1;
2,-4, 2, 1,-2, 1;
1, 0, 0,-2, 0, 0, 1;
2, 1, 0,-4,-2, 0, 2, 1;
1,-2, 1,-2, 4,-2, 1,-2, 1; ...
and A117939 is the matrix square of A117947:
1;
1, 1;
1,-1, 1;
1, 0, 0, 1;
1, 1, 0, 1, 1;
1,-1, 1, 1,-1, 1;
1, 0, 0,-1, 0, 0, 1;
1, 1, 0,-1,-1, 0, 1, 1;
1,-1, 1,-1, 1,-1, 1,-1, 1; ...
-
/* Generated as the Matrix LOG of A117939: */ T(n,k)=local(M=matrix(n+1,n+1,r,c,(binomial(r-1,c-1)+1)%3-1)^2, L=sum(i=1,#M,-(M^0-M)^i/i));return(L[n+1,k+1])
-
/* Generated as the Ternary Fractal: */ T(n,k)=local(r=n,c=k,s=floor(log(n+1)/log(3))+1,u=vector(s),v=vector(s),d,e); if(n<=k,0,if(n<3&k<3,[0,0,0;2,0,0;3,-2,0][n+1,k+1], for(i=1,#u,u[i]=r%3;r=r\3);for(i=1,#v,v[i]=c%3;c=c\3); d=0;for(i=1,#v,if(u[i]!=v[i],d+=1;e=i));if(d==1,T(u[e],v[e]),0)))
A120855
Row sums of triangle A120854, which is the matrix log of triangle A117939.
Original entry on oeis.org
0, 2, 1, 2, 4, 3, 1, 3, 2, 2, 4, 3, 4, 6, 5, 3, 5, 4, 1, 3, 2, 3, 5, 4, 2, 4, 3, 2, 4, 3, 4, 6, 5, 3, 5, 4, 4, 6, 5, 6, 8, 7, 5, 7, 6, 3, 5, 4, 5, 7, 6, 4, 6, 5, 1, 3, 2, 3, 5, 4, 2, 4, 3, 3, 5, 4, 5, 7, 6, 4, 6, 5, 2, 4, 3, 4, 6, 5, 3, 5, 4, 2, 4, 3, 4, 6, 5, 3, 5, 4, 4, 6, 5, 6, 8, 7, 5, 7, 6, 3, 5, 4, 5, 7, 6
Offset: 0
-
f[n_] := DigitCount[n, 3] /. {a_, b_, c_} -> 2a + b + 0c; Array[f, 105, 0] (* Robert G. Wilson v, Jul 31 2012 *)
-
{a(n)=local(M=matrix(n+1,n+1,r,c,(binomial(r-1,c-1)+1)%3-1)^2, L=sum(i=1,#M,-(M^0-M)^i/i));return(sum(k=0,n,L[n+1,k+1]))}
Showing 1-5 of 5 results.
Comments