|
%I #15 Sep 08 2022 08:45:52
%S 1,1,1,1,-3,1,1,17,17,1,1,-239,-219,-239,1,1,7169,6933,6933,7169,1,1,
%T -444479,-437307,-437563,-437307,-444479,1,1,56004353,55559877,
%U 55567029,55567029,55559877,56004353,1,1,-14225105663,-14169101307,-14169545803,-14169538395,-14169545803,-14169101307,-14225105663,1
%N Triangle T(n,k) = 1 - A176337(k) - A176337(n-k) + A176337(n) read by rows.
%C Row sums are {1, 2, -1, 36, -695, 28206, -2201133, 334262520, -99297043939, 57953303599938, -66678973493759897, ...}.
%H G. C. Greubel, <a href="/A176339/b176339.txt">Rows n = 0..25 of triangle, flattened</a>
%F T(n,k) = T(n,n-k).
%e Triangle begins as:
%e 1;
%e 1, 1;
%e 1, -3, 1;
%e 1, 17, 17, 1;
%e 1, -239, -219, -239, 1;
%e 1, 7169, 6933, 6933, 7169, 1;
%e 1, -444479, -437307, -437563, -437307, -444479, 1;
%p A176339 := proc(n,m)
%p 1-A176337(m)-A176337(n-m)+A176337(n) ;
%p end proc: # _R. J. Mathar_, May 04 2013
%t b[n_, q_]:= b[n, q]= If[n==0, 0, (1-q^n)*b[n-1, q] +1];
%t T[n_,k_,q_]:= 1 + b[n,q] -b[n-k,q] - b[k,q];
%t Table[T[n,k,2], {n,0,10}, {k,0,n}]//Flatten (* modified by _G. C. Greubel_, Dec 07 2019 *)
%o (PARI) b(n,q) = if(n==0, 0, 1 + (1-q^n)*b(n-1,q) );
%o T(n,k,q) = 1 + b(n,q) - b(n-k,q) - b(k,q);
%o for(n=0,10, for(k=0,n, print1(T(n,k,2), ", "))) \\ _G. C. Greubel_, Dec 07 2019
%o (Magma) function b(n,q)
%o if n eq 0 then return 0;
%o else return 1 - (q^n-1)*b(n-1,q);
%o end if; return b; end function;
%o function T(n,k,q) return 1 + b(n,q) - b(n-k,q) - b(k,q); end function;
%o [ T(n,k,2) : k in [0..n], n in [0..10]]; // _G. C. Greubel_, Dec 07 2019
%o (Sage)
%o @CachedFunction
%o def b(n, q):
%o if (n==0): return 0
%o else: return 1 - (q^n-1)*b(n-1,q)
%o def T(n,k,q): return 1 + b(n,q) - b(n-k,q) - b(k,q)
%o [[T(n,k,2) for k in (0..n)] for n in (0..10)] # _G. C. Greubel_, Dec 07 2019
%o (GAP)
%o b:= function(n,q)
%o if n=0 then return 0;
%o else return 1 - (q^n-1)*b(n-1,q);
%o fi; end;
%o T:= function(n,k,q) return 1 + b(n,q) - b(n-k,q) - b(k,q); end;
%o Flat(List([0..10], n-> List([0..n], k-> T(n,k,2) ))); # _G. C. Greubel_, Dec 07 2019
%Y Cf. A176337, A176338, A176340.
%K sign,tabl
%O 0,5
%A _Roger L. Bagula_, Apr 15 2010
|
