A121682 -id:A121682 - cp's OEIS frontend

A337001 a(n) = n! * Sum_{k=0..n} k^3 / k!.

0, 1, 10, 57, 292, 1585, 9726, 68425, 547912, 4931937, 49320370, 542525401, 6510306540, 84633987217, 1184875823782, 17773137360105, 284370197765776, 4834293362023105, 87017280516421722, 1653328329812019577, 33066566596240399540, 694397898521048399601

Offset: 0

Ilya Gutkovskiy, Aug 10 2020

nonn

Exponential convolution of cubes (A000578) and factorial numbers (A000142).

Seiichi Manyama, Table of n, a(n) for n = 0..449

Cf. A000142, A000522, A000578, A007526, A030297, A256016, A337002.

Row sums of A121682.

Table[n! Sum[k^3/k!, {k, 0, n}], {n, 0, 21}]
nmax = 21; CoefficientList[Series[x (1 + 3 x + x^2) Exp[x]/(1 - x), {x, 0, nmax}], x] Range[0, nmax]!
a[0] = 0; a[n_] := a[n] = n (n^2 + a[n - 1]); Table[a[n], {n, 0, 21}]

a(n) = n! * sum(k=0, n, k^3/k!); \\ Michel Marcus, Aug 12 2020

E.g.f.: x * (1 + 3*x + x^2) * exp(x) / (1 - x).
a(0) = 0; a(n) = n * (n^2 + a(n-1)).
a(n) ~ 5*exp(1)*n!. - Vaclav Kotesovec, Jan 13 2024

A121662 Triangle read by rows: T(i,j) for the recurrence T(i,j) = (T(i-1,j) + 1)*i.

Original entry on oeis.org

1, 4, 2, 15, 9, 3, 64, 40, 16, 4, 325, 205, 85, 25, 5, 1956, 1236, 516, 156, 36, 6, 13699, 8659, 3619, 1099, 259, 49, 7, 109600, 69280, 28960, 8800, 2080, 400, 64, 8, 986409, 623529, 260649, 79209, 18729, 3609, 585, 81, 9, 9864100, 6235300, 2606500, 792100, 187300, 36100, 5860, 820, 100, 10

Offset: 1

Thomas Wieder, Aug 15 2006

The first column is A007526 = "the number of (nonnull) "variations" of n distinct objects, namely the number of permutations of nonempty subsets of {1,...,n}." E.g. for n=3 there are 15 subsets: {a}, {b}, {c}, {ab}, {ba}, {ac}, {ca}, {bc}, {cb}, {abc}, {acb}, {bac}, {bca}, {cab}, {cba}. These are subsets with a number of elements l=1,...,n. The second column excludes all subsets with l=n elements. For n=3 one has therefore only the 9 subsets {a}, {b}, {c}, {ab}, {ba}, {ac}, {ca}, {bc}, {cb}. The third column excludes all subsets with l>=n-1 elements. For n=3 one has therefore only the 3 subsets {a}, {b},{c}. See also A121684. The second column is A038156 = n!*Sum(1/k!, k=1..n-1). The first lower diagonal are the squares A000290 = n^2. The second lower diagonal (15, 40, 85...) is A053698 = n^3 + n^2 + n + 1. The row sum is A030297 = a(n) = n*(n+a(n-1)).

T(i, j) is the total number of ordered sets of size 1 to i-j+1 that can be created from i distinct items. - Manfred Boergens, Jun 22 2022

			Triangle T(i,j) begins:
       1
       4     2
      15     9     3
      64    40    16     4
     325   205    85    25    5
    1956  1236   516   156   36   6
   13699  8659  3619  1099  259  49  7
   ...

Cf. A007526, A030297, A053698, A038156, A000290, A121682.

Mirror of triangle A285268.

T:= proc(i, j) option remember;
      `if`(j<1 or j>i, 0, T(i-1, j)*i+i)
    end:
seq(seq(T(n, k), k=1..n), n=1..10);  # Alois P. Heinz, Jun 22 2022

Table[Sum[m!/(m - i)!, {i, n}], {m, 9}, {n, m, 1, -1}] // Flatten (* Michael De Vlieger, Apr 22 2017 *)
(* Sum-free code *)
b[j_] = If[j == 0, 0, Floor[j! E - 1]];
T[i_, j_] = b[i] - i! b[j - 1]/(j - 1)!;
Table[T[i, j], {i, 24}, {j, i}] // Flatten
(* Manfred Boergens, Jun 22 2022 *)

From Manfred Boergens, Jun 22 2022: (Start)
T(i, j) = Sum_{k=1..i-j+1} i!/(i-k)! = Sum_{k=j-1..i-1} i!/k!.
Sum-free formula: T(i, j) = b(i) - i!*b(j-1)/(j-1)! where b(0)=0, b(j)=floor(j!*e-1) for j>0.
(End)

Edited by N. J. A. Sloane, Sep 15 2006

Formula in name corrected by Alois P. Heinz, Jun 22 2022

cp's OEIS Frontend

A337001 a(n) = n! * Sum_{k=0..n} k^3 / k!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula