| S3(2, 1); | |
| S1(3, 1); | |
| for (int c0 = 4; c0 <= M + 1; c0 += 1) { | |
| S1(c0, 1); | |
| for (int c1 = 2; c1 < (c0 + 1) / 2; c1 += 1) | |
| S2(c0, c1); | |
| if (c0 % 2 == 0) | |
| S4(c0, c0 / 2); | |
| } | |
| for (int c0 = M + 2; c0 <= 2 * M; c0 += 1) { | |
| for (int c1 = -M + c0; c1 < (c0 + 1) / 2; c1 += 1) | |
| S2(c0, c1); | |
| if (c0 % 2 == 0) | |
| S4(c0, c0 / 2); | |
| } |