|
MATHEMATICA
|
q[n_, k_] := n + k + 1; (* A094727 *)
r[0] = 1; r[k_] := Sum[q[k - 1, i] r[k - 1 - i], {i, 0, k - 1}]
p[n_, k_] := n!/(k! (n - k)!); (* Pascal's triangle *)
v[n_] := Sum[p[n, k] r[n - k], {k, 0, n}]
Table[v[n], {n, 0, 18}] (* A193659 *)
TableForm[Table[q[i, k], {i, 0, 4}, {k, 0, i}]]
Table[r[k], {k, 0, 8}] (* A193668 *)
TableForm[Table[p[n, k], {n, 0, 4}, {k, 0, n}]]
|