proposed
approved
proposed
approved
editing
proposed
b[args_List] := b[args] = Module[{nargs = Length[args]}, If[args[[1]] == 0, If[nargs == 2, 1, b[args // Rest]], Sum[If[args[[j]] - Last[args] < 0, 0, b[Append[Sort[Flatten[Table[args[[i]] - If[i == j, Last[args], 0], {i, 1, nargs-1}]]], Last[args]-1]]], {j, 1, nargs-1}]]];
a[n_] := If[Mod[#, 3] == 0, b[{#/3, #/3, #/3, n}]/6, 0]&[n(n+1)/2];
Array[a, 42] (* Jean-François Alcover, Oct 30 2020, after Alois P. Heinz *)
approved
editing
Alois P. Heinz, <a href="/A112972/b112972_1.txt">Table of n, a(n) for n = 1..125</a>
editing
approved
Alois P. Heinz, <a href="/A112972/b112972_1.txt">Table of n, a(n) for n = 1..100125</a>
approved
editing
editing
approved
editing
approved
b:= proc() option remember; local i, j, t; `if`(args[1]=0, `if`(nargs=2, 1, b(args[t] $t=2..nargs)), add(`if`(args[j] -args[nargs] <0, 0, b(sort([seq(args[i] -`if`(i=j, args[nargs], 0), i=1..nargs-1)])[], args[nargs]-1)), j=1..nargs-1)) end: a:= proc(n) local m; m:= n*(n+1)/2; `if`(irem(m, 3)=0, b((m/3)$3, n)/6, 0) end: seq(a(n), n=1..42); # Alois P. Heinz, Sep 03 2009
b:= proc() option remember; local i, j, t; `if`(args[1]=0,
`if`(nargs=2, 1, b(args[t] $t=2..nargs)), add(
`if`(args[j] -args[nargs]<0, 0, b(sort([seq(args[i]-
`if`(i=j, args[nargs], 0), i=1..nargs-1)])[],
args[nargs]-1)), j=1..nargs-1))
end:
a:= n-> (m-> `if`(irem(m, 3)=0, b((m/3)$3, n)/6, 0))(n*(n+1)/2):
seq(a(n), n=1..42); # Alois P. Heinz, Sep 03 2009
approved
editing