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”).
%I #19 May 10 2023 06:15:09
%S 1,2,12,56,297,1632,9531,58634,379371,2574254,18276457,135463074,
%T 1046041114,8399533370,70013963418,604840440328,5407301690915,
%U 49958478263502,476403955991034,4683463406478004,47414166201239781,493803423334040824,5285548108715948453
%N Number of partitions of {1...n} containing 5 strings of 3 consecutive integers, where each string is counted within a block and a string of more than 3 consecutive integers are counted three at a time.
%H Augustine O. Munagi, <a href="https://doi.org/10.1155/IJMMS.2005.215">Set Partitions with Successions and Separations</a>, Int. J. Math and Math. Sc., 2005:3 (2005), 451-463.
%F a(n) = Sum_{k=1..n} c(n, k, 5), where c(n, k, 5) is the case r=5 of c(n, k, r) given by c(n, k, r)=c(n-1, k-1, r)+(k-1)c(n-1, k, r)+c(n-2, k-1, r)+(k-1)c(n-2, k, r)+c(n-1, k, r-1)-c(n-2, k-1, r-1)-(k-1)c(n-2, k, r-1), r=0, 1, .., n-k-1, k=1, 2, .., n-2r, c(n, k, 0) = Sum_{j= 0..floor(n/2)} binomial(n-j, j)*S2(n-j-1, k-1).
%e a(8) = 2 because the partitions of {1,...,8} with 5 strings of 3 consecutive integers are 1234567/8, 1/2345678.
%p c := proc(n,k,r) option remember ; local j ; if r =0 then add(binomial(n-j,j)*combinat[stirling2](n-j-1,k-1),j=0..floor(n/2)) ; else if r <0 or r > n-k-1 then RETURN(0) fi ; if n <1 then RETURN(0) fi ; if k <1 then RETURN(0) fi ; RETURN( c(n-1,k-1,r)+(k-1)*c(n-1,k,r)+c(n-2,k-1,r)+(k-1)*c(n-2,k,r) +c(n-1,k,r-1)-c(n-2,k-1,r-1)-(k-1)*c(n-2,k,r-1) ) ; fi ; end: A105487 := proc(n) local k ; add(c(n,k,5),k=1..n) ; end: for n from 7 to 30 do printf("%d, ",A105487(n)) ; od ; # _R. J. Mathar_, Feb 20 2007
%t S2[_, -1] = 0;
%t S2[n_, k_] = StirlingS2[n, k];
%t c[n_, k_, r_] := c[n, k, r] = Which[
%t r == 0, Sum[Binomial[n - j, j]*S2[n - j - 1, k - 1],
%t {j, 0, Floor[n/2]}],
%t r < 0 || r > n - k - 1, 0,
%t n < 1, 0,
%t k < 1, 0,
%t True, c[n - 1, k - 1, r] +
%t (k - 1)*c[n - 1, k, r] +
%t c[n - 2, k - 1, r] +
%t (k - 1)*c[n - 2, k, r] +
%t c[n - 1, k, r - 1] -
%t c[n - 2, k - 1, r - 1] -
%t (k - 1)*c[n - 2, k, r - 1]];
%t A105487[n_] := Sum[c[n, k, 5], {k, 1, n}];
%t Table[A105487[n], {n, 7, 30}] (* _Jean-François Alcover_, May 10 2023, after _R. J. Mathar_ *)
%Y Cf. A105486, A105482.
%K nonn
%O 7,2
%A _Augustine O. Munagi_, Apr 10 2005
%E More terms from _R. J. Mathar_, Feb 20 2007