gist-3.c

for (int c0 = 1; c0 < n; c0 += 4)
  for (int c1 = c0 + 1; c1 <= n; c1 += 6)
    s0(c0, c1);