cp's OEIS Frontend

A135222 Triangle A049310 + A000012 - I, read by rows.

%I A135222 #9 Mar 27 2022 18:25:49
%S A135222 1,1,1,2,1,1,1,3,1,1,2,1,4,1,1,1,4,1,5,1,1,2,1,7,1,6,1,1,1,5,1,11,1,7,
%T A135222 1,1,2,1,11,1,16,1,8,1,1,1,6,1,21,1,22,1,9,1,1,2,1,16,1,36,1,29,1,10,
%U A135222 1,1,1,7,1,36,1,57,1,37,1,11,1,1,2,1,22,1,71,1,85,1,46,1,12,1,1
%N A135222 Triangle A049310 + A000012 - I, read by rows.
%C A135222 Row sums = A081659: (1, 2 4, 6, 9, 13, 19, 28, ...).
%H A135222 G. C. Greubel, <a href="/A135222/b135222.txt">Rows n = 0..100 of triangle, flattened</a>
%F A135222 T(n,k) = A049310(n,k) + A000012(n,k) - Identity matrix, as infinite lower triangular matrices.
%F A135222 T(n,k) = 1 + abs( ((1+(-1)^(n-k))/2)*binomial((n+k)/2, (n-k)/2)*cos((n-k)*Pi/2) ), with T(n,n) = 1. - _G. C. Greubel_, Nov 20 2019
%e A135222 First few rows of the triangle:
%e A135222   1;
%e A135222   1, 1;
%e A135222   2, 1,  1;
%e A135222   1, 3,  1,  1;
%e A135222   2, 1,  4,  1,  1;
%e A135222   1, 4,  1,  5,  1, 1;
%e A135222   2, 1,  7,  1,  6, 1, 1;
%e A135222   1, 5,  1, 11,  1, 7, 1, 1;
%e A135222   2, 1, 11,  1, 16, 1, 8, 1, 1;
%e A135222 ...
%p A135222 T:= proc(n, k) option remember;
%p A135222       if k=n then 1
%p A135222     else 1 + abs( ((1+(-1)^(n-k))/2)*binomial((n+k)/2, (n-k)/2)*cos((n-k)*Pi/2) )
%p A135222       fi; end:
%p A135222 seq(seq(T(n, k), k=0..n), n=0..15); # _G. C. Greubel_, Nov 20 2019
%t A135222 T[n_, k_]:= T[n, k]= If[k==n, 1, 1 + Abs[Simplify[((1+(-1)^(n-k))/2)* Binomial[(n+k)/2, (n-k)/2]*Cos[(n-k)*Pi/2]]] ]; Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* _G. C. Greubel_, Nov 20 2019 *)
%o A135222 (Sage)
%o A135222 @CachedFunction
%o A135222 def T(n, k):
%o A135222     if (k==n): return 1
%o A135222     else: return 1 + abs( ((1+(-1)^(n-k))/2)*binomial((n+k)/2, (n-k)/2)*cos((n-k)*pi/2) )
%o A135222 [[T(n, k) for k in (0..n)] for n in (0..15)] # _G. C. Greubel_, Nov 20 2019
%Y A135222 Cf. A049310, A081659.
%K A135222 nonn,tabl
%O A135222 0,4
%A A135222 _Gary W. Adamson_, Nov 23 2007
%E A135222 More terms added and offset changed by _G. C. Greubel_, Nov 20 2019