|
%I #89 Sep 23 2020 15:06:57
%S 1,3,4,5,8,21,7,12,35,96,9,16,49,144,410,11,20,63,192,574,1680,13,24,
%T 77,240,738,2240,6692,15,28,91,288,902,2800,8604,26112,17,32,105,336,
%U 1066,3360,10516,32640,100296,19,36,119,384,1230,3920,12428,39168,122584,380480
%N Triangle read by rows: T(n,k) = 2*n+1 for k = 0 and otherwise T(n,k) = Sum_{i=n-k..n, j=0..i-n+k, i<>n or j<>k} T(i,j).
%H Oboifeng Dira, <a href="/A335436/b335436.txt">Table of n, a(n) for n = 0..54</a>
%F T(n,0) = 2*n+1 for k=0;
%F T(n,k) = ((2+sqrt(2))^(k+1)-(2-sqrt(2))^(k+1))*(2*n-k+1)/(4*sqrt(2)) for 1<=k<=n.
%e Triangle begins:
%e 1;
%e 3, 4;
%e 5, 8, 21;
%e 7, 12, 35, 96;
%e 9, 16, 49, 144, 410;
%e ...
%e T(3,2) = ((2+sqrt(2))^3-(2-sqrt(2))^3)*(6-2+1)/(4*sqrt(2)) = (28*sqrt(2))*(5)/(4*sqrt(2) = 35.
%p T := proc (n, k) if k = 0 and 0 <= n then 2*n+1 elif 1 <= k and k <= n then round((((2+sqrt(2))^(k+1)-(2-sqrt(2))^(k+1))*(2*n-k+1)/(4*sqrt(2)))) else 0 end if end proc:seq(print(seq(T(n, k), k=0..n)), n=0..9);
%o (PARI) T(n,k) = if (k==0, 2*n+1, if (k<=n, sum(i=n-k, n, sum(j=0, i-n+k, if ((i==n) && (j==k), 0, T(i,j)), 0))));
%o matrix(10, 10, n, k, T(n-1, k-1)) \\ _Michel Marcus_, Sep 08 2020
%o (PARI) T(n, k) = if (k==0, 2*n+1, if (k>n, 0, my(w=quadgen(8, 'w)); ((2+w)^(k+1)-(2-w)^(k+1))*(2*n-k+1)/(4*w)));
%o matrix(10, 10, n, k, T(n-1, k-1)) \\ _Michel Marcus_, Sep 10 2020
%Y Cf. A053199, A294285, A064339.
%K nonn,tabl
%O 0,2
%A _Oboifeng Dira_, Jul 14 2020
|
