OFFSET
0,5
REFERENCES
Carl Erik Froeberg, Numerical Mathematics, Benjamin/Cummings Pu.Co. 1985, ISBN 0-8053-2530-1, Chapter 17.2.
Charles Jordan, Calculus of Finite Differences, Chelsea 1965, p. 513.
LINKS
Alois P. Heinz, Rows n = 0..100, flattened
W. M. Johnson, On Cotesian numbers: their history, computation and values to n=20, Quart. J. Pure Appl. Math., 46 (1914), 52-65. [Annotated scanned copy]
EXAMPLE
0, 1/2, 1/2, 1/6, 2/3, 1/6, 1/8, 3/8, 3/8, 1/8, 7/90, 16/45, 2/15, 16/45, 7/90, 19/288, 25/96, 25/144, 25/144, 25/96, 19/288, 41/840, 9/35, 9/280, 34/105, 9/280, 9/35, 41/840, ... = A100640/A100641 = A100642/A002176 (the latter is not in lowest terms)
Triangle begins
0;
1, 1;
1, 4, 1;
1, 3, 3, 1;
7, 32, 12, 32, 7;
MAPLE
# (This defines the Cotesian numbers C(n, i))
with(combinat); C:=proc(n, i) if i=0 or i=n then RETURN( (1/n!)*add(n^a*stirling1(n, a)/(a+1), a=1..n+1) ); fi; (1/n!)*binomial(n, i)* add( add( n^(a+b)*stirling1(i, a)*stirling1(n-i, b)/((b+1)*binomial(a+b+1, b+1)), b=1..n-i+1), a=1..i+1); end;
den:=proc(n) local t1, i; t1:=1; for i from 0 to n do t1:=ilcm(t1, denom(C(n, i))); od: t1; end;
# Then den(n)*C(n, k) gives the current sequence
seq(seq(den(n, k)*C(n, k), k=0..n), n=0..10);
MATHEMATICA
c[n_, i_] /; i == 0 || i == n = (1/n!)*Sum[n^a*StirlingS1[n, a]/(a+1), {a, 1, n+1}]; c[n_, i_] = (1/n!)*Binomial[n, i]*Sum[n^(a + b)*StirlingS1[i, a]*StirlingS1[n-i, b]/((b+1)*Binomial[a+b+1, b+1]), {b, 1, n}, {a, 1, i+1}]; den[n_] := (For[t1 = 1; i = 0, i <= n, i++, t1 = LCM[t1, c[n, i] // Denominator]]; t1); Table[den[n]*c[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Apr 11 2013, after Maple *)
CROSSREFS
KEYWORD
AUTHOR
N. J. A. Sloane, Dec 04 2004
STATUS
approved