上三角蛇形矩阵

void matrix(int n)
{
 int i;
 for (i = 0; i < n; i++)
 {
  int j;
  for (j = 0; j < n-i; j++)
  {
   ++i;
   j++;
   printf("%-3d",i*(i+1)/2+j*(j+1)/2+i*(j-1)-i-j+1);
   --i;
   --j;
  }
  puts("");
 }
}