gist-5.c

for (int c0 = 1; c0 <= n; c0 += 12)
  for (int c1 = c0; c1 <= n; c1 += 8)
    s0(c0, c1);