1for (int c0 = 2; c0 <= n + m; c0 += 1)
2  for (int c1 = max(-m + c0, 1); c1 <= min(n, c0 - 1); c1 += 1)
3    S1(c1, c0 - c1);
4