login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A131108
T(n,k) = 2*A007318(n,k) - A097806(n,k).
4
1, 1, 1, 2, 3, 1, 2, 6, 5, 1, 2, 8, 12, 7, 1, 2, 10, 20, 20, 9, 1, 2, 12, 30, 40, 30, 11, 1, 2, 14, 42, 70, 70, 42, 13, 1, 2, 16, 56, 112, 140, 112, 56, 15, 1, 2, 18, 72, 168, 252, 252, 168, 72, 17, 1, 2, 20, 90, 240, 420, 504, 420, 240, 90, 19, 1
OFFSET
0,4
COMMENTS
Row sums give A095121.
Triangle T(n,k), 0 <= k <= n, read by rows given by [1, 1, -2, 1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1,0,0,1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 18 2007
FORMULA
Twice Pascal's triangle minus A097806, the pairwise operator.
G.f.: (1-x*y+x^2+x^2*y)/((-1+x+x*y)*(x*y-1)). - R. J. Mathar, Aug 11 2015
EXAMPLE
First few rows of the triangle are:
1;
1, 1;
2, 3, 1;
2, 6, 5, 1;
2, 8, 12, 7, 1;
2, 10, 20, 20, 9, 1;
...
MAPLE
seq(seq( `if`(k=n-1, 2*n-1, `if`(k=n, 1, 2*binomial(n, k))), k=0..n), n=0..12); # G. C. Greubel, Nov 18 2019
MATHEMATICA
Table[If[k==n-1, 2*n-1, If[k==n, 1, 2*Binomial[n, k]]], {n, 0, 12}, {k, 0, n}]//Flatten (* G. C. Greubel, Nov 18 2019 *)
PROG
(PARI) T(n, k) = if(k==n-1, 2*n-1, if(k==n, 1, 2*binomial(n, k))); \\ G. C. Greubel, Nov 18 2019
(Magma)
function T(n, k)
if k eq n-1 then return 2*n-1;
elif k eq n then return 1;
else return 2*Binomial(n, k);
end if;
return T;
end function;
[T(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Nov 18 2019
(Sage)
@CachedFunction
def T(n, k):
if (k==n-1): return 2*n-1
elif (k==n): return 1
else: return 2*binomial(n, k)
[[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 18 2019
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Gary W. Adamson, Jun 15 2007
EXTENSIONS
Corrected by Philippe Deléham, Dec 17 2007
More terms added and data corrected by G. C. Greubel, Nov 18 2019
STATUS
approved