cp's OEIS Frontend

A327827 Sum of multinomials M(n; lambda), where lambda ranges over all partitions of n into parts incorporating 1.

%I A327827 #12 Dec 09 2020 15:06:11
%S A327827 0,1,2,9,40,235,1476,11214,91848,859527,8710300,97675138,1179954612,
%T A327827 15490520786,217602374458,3280028076615,52571985879600,
%U A327827 895913825750191,16140560853800556,307048409240931810,6143666813617775100,129096480664676773542,2840750997343361802150
%N A327827 Sum of multinomials M(n; lambda), where lambda ranges over all partitions of n into parts incorporating 1.
%H A327827 Alois P. Heinz, <a href="/A327827/b327827.txt">Table of n, a(n) for n = 0..450</a>
%H A327827 Wikipedia, <a href="https://en.wikipedia.org/wiki/Multinomial_theorem#Multinomial_coefficients">Multinomial coefficients</a>
%H A327827 Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_(number_theory)">Partition (number theory)</a>
%F A327827 a(n) ~ c * n!, where c = A247551 = 2.5294774720791526481801161542539542411787... - _Vaclav Kotesovec_, Sep 28 2019
%p A327827 b:= proc(n, i, k) option remember; `if`(n=0, 1,
%p A327827      `if`(i>n, 0, b(n, i+1, `if`(i=k, 0, k))+
%p A327827      `if`(i=k, 0, b(n-i, i, k)*binomial(n, i))))
%p A327827     end:
%p A327827 a:= n-> b(n, 1, 0)-b(n, 1$2):
%p A327827 seq(a(n), n=0..23);
%t A327827 b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 2, 0, b[n, i - 1, If[i == k, 0, k]]] + If[i == k, 0, b[n - i, Min[n - i, i], k]/i!]];
%t A327827 T[n_, k_] := n! (b[n, n, 0] - If[k == 0, 0, b[n, n, k]]);
%t A327827 a[n_] := T[n, 1];
%t A327827 a /@ Range[0, 23] (* _Jean-François Alcover_, Dec 09 2020, after _Alois P. Heinz_ *)
%Y A327827 Column k=1 of A327801.
%K A327827 nonn
%O A327827 0,3
%A A327827 _Alois P. Heinz_, Sep 26 2019