login
A skewed version of triangular array A122075.
1

%I #18 May 21 2019 03:48:37

%S 1,0,2,0,1,3,0,0,3,5,0,0,1,7,8,0,0,0,4,15,13,0,0,0,1,12,30,21,0,0,0,0,

%T 5,31,58,34,0,0,0,0,1,18,73,109,55,0,0,0,0,0,6,54,162,201,89,0,0,0,0,

%U 0,1,25,145,344,365,144,0,0,0,0,0,0,7,85,361

%N A skewed version of triangular array A122075.

%C Triangle T(n,k), 0 <= k <= n, read by rows, given by (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (2, -1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

%C Subtriangle of the triangle A122950.

%H Reinhard Zumkeller, <a href="/A236076/b236076.txt">Rows n = 0..125 of triangle, flattened</a>

%H H. Fuks and J.M.G. Soto, <a href="http://arxiv.org/abs/1306.1189">Exponential convergence to equilibrium in cellular automata asymptotically emulating identity</a>, arXiv:1306.1189 [nlin.CG], 2013.

%F G.f.: (1+x*y)/(1 - x*y - x^2*y - x^2*y^2).

%F T(n,k) = T(n-1,k-1) + T(n-2,k-1) + T(n-2,k-2), T(0,0)=1, T(1,0) = 0, T(1,1) = 2, T(n,k) = 0 if k < 0 or if k > n.

%F Sum_{k=0..n} T(n,k) = 2^n = A000079(n).

%F Sum_{n>=k} T(n,k) = A078057(k) = A001333(k+1).

%F T(n,n) = Fibonacci(n+2) = A000045(n+2).

%F T(n+1,n) = A023610(n-1), n >= 1.

%F T(n+2,n) = A129707(n).

%e Triangle begins:

%e 1;

%e 0, 2;

%e 0, 1, 3;

%e 0, 0, 3, 5;

%e 0, 0, 1, 7, 8;

%e 0, 0, 0, 4, 15, 13;

%e 0, 0, 0, 1, 12, 30, 21;

%e 0, 0, 0, 0, 5, 31, 58, 34;

%t T[n_, k_]:= If[k<0 || k>n, 0, If[n==0 && k==0, 1, If[k==0, 0, If[n==1 && k==1, 2, T[n-1, k-1] + T[n-2, k-1] + T[n-2, k-2]]]]]; Table[T[n,k], {n,0,10}, {k,0,n}]//Flatten (* _G. C. Greubel_, May 21 2019 *)

%o (Haskell)

%o a236076 n k = a236076_tabl !! n !! k

%o a236076_row n = a236076_tabl !! n

%o a236076_tabl = [1] : [0, 2] : f [1] [0, 2] where

%o f us vs = ws : f vs ws where

%o ws = [0] ++ zipWith (+) (zipWith (+) ([0] ++ us) (us ++ [0])) vs

%o -- _Reinhard Zumkeller_, Jan 19 2014

%o (PARI)

%o {T(n,k) = if(k<0 || k>n, 0, if(n==0 && k==0, 1, if(k==0, 0, if(n==1 && k==1, 2, T(n-1,k-1) + T(n-2,k-1) + T(n-2,k-2) ))))}; \\ _G. C. Greubel_, May 21 2019

%o (Sage)

%o def T(n, k):

%o if (k<0 or k>n): return 0

%o elif (n==0 and k==0): return 1

%o elif (k==0): return 0

%o elif (n==1 and k==1): return 2

%o else: return T(n-1,k-1) + T(n-2,k-1) + T(n-2,k-2)

%o [[T(n, k) for k in (0..n)] for n in (0..10)] # _G. C. Greubel_, May 21 2019

%Y Cf. variant: A055830, A122075, A122950, A208337.

%Y Cf. A167704 (diagonal sums), A000079 (row sums).

%Y Cf. A111006.

%K easy,nonn,tabl

%O 0,3

%A _Philippe Deléham_, Jan 19 2014