{{{id=1|
# replace /home/tom/maths/free_sigma/ with the path to cutjoin.sage on your computer
load "/home/tom/maths/free_sigma/cutjoin.sage"
///
}}}

{{{id=2|
#check [2]*[2] then  [2]*[2]*[2]
n=7
print multSnAlg([2],[2],n)
print multSnAlgArr(multSnAlg([2],[2],n),[[1,[2]]],n)
///
[[21, [1, 1, 1, 1, 1, 1, 1]], [3, [3, 1, 1, 1, 1]], [2, [2, 2, 1, 1, 1]]]
[[71, [2, 1, 1, 1, 1, 1]], [16, [4, 1, 1, 1]], [9, [3, 2, 1, 1]], [6, [2, 2, 2, 1]]]
}}}

{{{id=3|
#check [2]*[2]*[2]
for n in range(6,15):
    print binomial(n,2)^3,(3*n^2+n-12)/2*binomial(n,2)+96*binomial(n,4) + 180*binomial(n,5) + 90*binomial(n,6)
///
3375 3375
9261 9261
21952 21952
46656 46656
91125 91125
166375 166375
287496 287496
474552 474552
753571 753571
}}}