编写一个程序，将自然数1～n² 按“蛇形”填入n×n矩阵中。例如，当n=5时，构造的方阵如下。

      

    （1）编程思路1。

        分析数的填法，是按“从右上到左下”的”蛇形”或“从左下到右上”的”蛇形”，沿平行于副对角线的各条对角线上，将自然数从小到大填写。

        设用（i,j）表示矩阵元素的坐标，初始时，i=0，j=0。

        先从右上到左下填数，当从右上到左下时，坐标i增加（i++），坐标j减小（j--），当j减到小于0（副对角线及以上部分从右上到左下填写时），或i加到大于n-1（副对角线以下部分从右上到左下填写时）时结束。然后调整坐标（i,j），调整方法为：当在副对角线及以上时（此时j<0且i<n），j=0；当过副对角线后（此时i>n-1），i=n-1，j=j+2。调整坐标后，再开始从左下向右上填数。

        当从左下到右上时，坐标i减小（i--），坐标j增加（j++），当i减到小于0（副对角线及以上部分从左下到右上填写时），或j加到大于n-1（副对角线以下部分从左下到右上填写时）时结束。然后调整坐标（i,j），调整方法为：当在副对角线及以上时（此时i<0且j<n），i=0；当过副对角线后（此时j>n-1），i=i+2，j=n-1。调整坐标后，又开始从右上向左下填数。

        如此循环，直到n*n个数填完为止。

        （2）源程序1。

#include <stdio.h>
#define  N  19
int main()
{
    int a[N][N];
    int n,i=0,j=0,k=1;            // i,j是矩阵元素的下标,k是要填入的自然数
    scanf("%d",&n);
    while(i<n && j<n)
    {
        while(i<n && j>-1)       // 从右上向左下填数
        { a[i][j]=k++; i++ ;j--; }
        if ((j<0)&&(i<n)) j=0;  // 副对角线及以上部分的新i,j坐标
        else {j=j+2; i=n-1;}    // 副对角线以下的新的i,j坐标.
        while(i>-1 && j<n)      // 从左下向右上
        { a[i][j]=k++; i--; j++; }
        if(i<0 && j<n) i=0;
        else{i=i+2; j=n-1;}
    }
    for (i=0;i<n;i++)
    {
        for (j=0;j<n;j++)
            printf("%4d",a[i][j]);
        printf("\n");
    }
    return 0;
}

        （3）编程思路2。 

        n×n矩阵有2*n-1条平行于副对角线的对角线（下面简称为斜线）。例如，n=5时，有9条斜线。可将这9条斜线按从上到下用编号k依次标记为1~9。K=5时为副对角线，k<5时，斜线位于副对角线上面。

        由题目给出的5×5方阵结果知，该方阵有9条斜线。第1条斜线上填写1个数，第2条斜线上填写2个数，…，第5条斜线上填写5个数，第6条斜线上填写4个数，…，第9条斜线上填写1个数。

        设第k条斜线上需要填写的数字个数为q。显然有

        当k<=n时，q=k；  当k>n时，q=2*n-k。

        设用（i,j）表示矩阵元素的坐标，1≤i≤n,1≤j≤n。

        第k条斜线上的坐标（i,j）一定满足等式： i+j=k+1

        例如，第3条斜线上的3个元素的坐标从右上到左下依次为 （1,3）、（2,2）、（3,1），满足i+j=k+1=4。同时，可以看出，第3条斜线上的元素坐标i从1到3。

        再如，第4条斜线上的4个元素的坐标从左下到右上依次为 （4,1）、（3,2）、（2,3）、（1,4），满足i+j=k+1=5。同时，可以看出，第4条斜线上的元素坐标j从1到4。

        由此，可找出第k条斜线上填写的q个元素中的第p个元素的坐标（i,j）与k、p、q的关系。

        当k为奇数时，斜线上的数字从右上到左下填写，i=p，j=q+1-p；

        当k为偶数时，斜线上的数字从左下到右上填写，i=q+1-j，j=p。

        但当k>n时，此时斜线在副对角线的下方，需要对坐标再修正一下。

        例如，第6条斜线上的4个元素的坐标从左下到右上依次为 （5,2）、（4,3）、（3,4）、（2,5），同样满足i+j=k+1=7。但是，第6条斜线上的元素坐标j从2到5，i从5递减到2。而按上面的关系计算出的第6条斜线上的元素坐标j从1到4，i从4递减到1。此时，只需要将计算出的坐标i和j的值均加上1（也即n-q）即可。

        同样，考察第7条斜线上的3个元素的坐标从右上到左下依次为 （3,5）、（4,4）、（5,3），同样满足i+j=k+1=8。但是，第7条斜线上的元素坐标i从3到5，j从5递减到3。而按上面的关系计算出的第7条斜线上的元素坐标i从1到3，j从3递减到1。此时，只需要将计算出的坐标i和j的值均加上2（也即n-q）即可。

        计算出了每条斜线上应该填写的每个元素的坐标，在该坐标位置填上相应的数字。采用循环将这2*n-1条斜线上均填写上数字，则“蛇形”方阵就可以构造完成。

        （4）源程序2。

#include <stdio.h>
#define  N  20
int main()
{
    int a[N][N];
    int n,i,j,k=1,m,p,q;
    scanf("%d",&n);
    m=1;                  // m是要填入的自然数
    for (k=1;k<2*n; k++)  // 第k条平行与副对角线的对角线
    {
       if (k<n) q=k;      // q是第k条对角线上填写数字的个数
       else  q=2*n-k;
       for (p=1;p<=q;p++)
       {
           if (k%2) { i=p; j=q+1-p;}
           else  { i=q+1-p; j=p; }
           if (k>n) { i=i+n-q;  j=j+n-q; }
           a[i][j]=m++;
       }
    }
    for (i=1;i<=n;i++)
    {
        for (j=1;j<=n;j++)
            printf("%4d",a[i][j]);
        printf("\n");
    }
    return 0;
}

         （5）编程思路3。 

        观察图1所示的蛇形阵可知，方阵在逐个填数构造的过程中，是沿两种斜线进行的，一种是“右上到左下”（斜向下），一种是“左下到右上”（斜向上）。     

      

                                                        图1 蛇形方阵的构造示意图

        设当前已填入数字的位置坐标为（i,j），i和j均在0~n-1之间。若按斜向下填写，则下一位置为row++、col-- ；若按斜向上填写，则下一位置为row--、col++。由于下一位置可能超出方阵的边界，因此有时需要调整。调整有4种情况，下面分别进行说明。

        斜向上填写时，会出现两种情况：

        （1）超过了首行的位置（即i==-1），如图1（a）所示，3填写好后，下一个数4的计算位置越界了，此时进行调整，方法为i++。

        （2）超过最右列的位置（即j==n），如图1（b）所示，19填写好后，下一个数20的计算位置越界了，此时进行调整，方法为i+=2、j--。

        一种特例，在如图1（c）所示的4阶方阵中，10填写好后，下一个数11的计算位置，行和列都越界了，其处理方法同列越界，因此在程序中应先处理j==n的情况，再处理i==-1的情况。这样对于这种特例，由于处理了j==n后，i加了2，不会越界，因此不会再处理i==-1的情况。

        斜向下填写时，也会出现两种情况：

        （1）超过最左列的位置（即j==-1），如图1（d）所示，6填写好后，下一个数7的计算位置越界了，此时进行调整，方法为j++。

        （2）超过了底行的位置（即 i==n），如图1（e）所示，22填写好后，下一个数23的计算位置越界了，此时进行调整，方法为 i--,j+=2。

        一种特例，15填写好后，下一个数16的计算位置，行和列都越界了，但处理方法同行越界，因此在程序中应先处理i==n的情况，再处理j==-1的情况。这样对于这种特例，由于处理了i==n后，j加了2，不会越界，因此不会再处理j==-1的情况。

        每次进行越界调整填数后，填数的方向也会发生变化。因此，可设置一个方向变量up，当up=1时，表示斜向下填数；up=0时，表示斜向上填数。

        初始化时，令up=1、i=0、j=0；在当前位置填上1（即a[i][j]=1），之后从 m=2开始进行循环，直到m==n*n，全部n*n个数填写完毕。循环中，总是先按up的方向，确定下一个位置，然后填上相应的数。

        （6）源程序3。

#include <stdio.h>
#define  N  20
int main()
{
    int a[N][N];
    int n,i,j,m;
    scanf("%d",&n);
    int up=1;
    i=0;     j=0;
    a[i][j]=1;
    m=2;      // m是要填入的自然数
    while (m<=n*n)
    {
        if (up) { i++;  j--; }         // 从右上到左下
        else    { i--;  j++; }         // 从左下到右上
        if (j==n)                      // 超过最右列的位置
        {    i=i+2; j--;   up=!up;  }
        if (i==-1)                    // 超过首行的位置
        {    i++;   up=!up;    }
        if (i==n)                     // 超过底行的位置
        {   i--;   j+=2; up=!up; }
        if (j==-1)                    // 超过最左列的位置
        {    j++;    up=!up;   }
        a[i][j]=m++;
    }
    for (i=0;i<n;i++)
    {
        for (j=0;j<n;j++)
            printf("%4d",a[i][j]);
        printf("\n");
    }
    return 0;
}