2 for (int c1 = 1; c1 <= M; c1 += 1) {
4 for (int c2 = c1 + 1; c2 <= M; c2 += 1)
7 for (int c0 = 1; c0 < 3 * M - 1; c0 += 3) {
9 for (int c1 = (c0 + 5) / 3; c1 <= M; c1 += 1) {
11 for (int c4 = (c0 + 5) / 3; c4 < c1; c4 += 1)
12 S5(c4, c1, (c0 + 2) / 3);
14 for (int c1 = (c0 + 5) / 3; c1 <= M; c1 += 1)