|
MATHEMATICA
|
triamat[li_List] := Block[{len=Sqrt[8Length[li]+1]/2-1/2}, If[IntegerQ[len], Part[li, # ]& /@ Table[If[j>i, j(j-1)/2+i, i(i-1)/2+j], {i, len}, {j, len}], li]]; Table[it=triamat/@ IntegerDigits[Range[0, -1+2^(n(n+1)/2)], 2, n(n+1)/2]; Count[it, (q_)?MatrixQ/; (Det[q]>0)], {n, 5}]
|