A164949 - OEIS

%I #13 Jun 08 2018 08:06:13

%S 0,0,0,0,0,0,1,3,9,23,67,203,693,2584,9929,37480,137067,522854,

%T 2052657,8199728,33456333,137831268,574295984,2392149818,9950364020,

%U 41860671346,177512155194,757447761138,3254519322231,14049972380612,60960849334377,265354255338637

%N Number of different ways to select 4 disjoint subsets from {1..n} with equal element sum.

%e a(7) = 1, because {1,6}, {2,5}, {3,4}, {7} are disjoint subsets of {1..7} with element sum 7.

%e a(8) = 3: {1,6}, {2,5}, {3,4}, {7} have element sum 7, {1,7}, {2,6}, {3,5}, {8} have element sum 8, and {1,8}, {2,7}, {3,6}, {4,5} have element sum 9.

%p b:= proc() option remember; local i, j; `if`(args[1]=0 and args[2]=0 and args[3]=0 and args[4]=0, 1, `if`(add(args[j], j=1..4)> args[5] *(args[5]-1)/2, 0, b(args[j]$j=1..4, args[5]-1)) +add(`if`(args[j] -args[5]<0, 0, b(sort([seq(args[i] -`if`(i=j, args[5], 0), i=1..4)])[], args[5]-1)), j=1..4)) end: a:= n-> add(b(k$4, n), k=7..floor(n*(n+1)/8)) /24: seq(a(n), n=1..20);

%t b[l_, n_, k_] := b[l, n, k] = Module[{i, j}, If[l == Array[0&, k], 1, If[ Total[l] > n(n-1)/2, 0, b[l, n-1, k]] + Sum[If[l[[j]]-n < 0, 0, b[Sort[ Table[l[[i]] - If[i==j, n, 0], {i, 1, k}]], n-1, k]], {j, 1, k}]]];

%t T[n_, k_] := Sum[b[Array[t&, k], n, k], {t, 2k-1, Floor[n(n+1)/(2k)]}]/k!;

%t a[n_] := T[n, 4];

%t Array[a, 20] (* _Jean-François Alcover_, Jun 08 2018, after _Alois P. Heinz_'s Maple code in A196231 *)

%Y Column k=4 of A196231.

%Y Cf. A161943, A164934.

%K nonn

%O 1,8

%A _Alois P. Heinz_, Sep 01 2009