if (i!=d)
            f=f+c;
      }
    }
}

int B(int n)  //计算n的阶乘
{
    if (n==0)
        return (1);
    else
        return (n*B(n-1));
}

void output(B1 *psnew, int fac, int **factor)  //输出psnew指向的节点指向的基拉丁方阵和置换得到的所有其他拉丁方阵
{
    int i, j, k;

    for (i=0; i<N; i++)
    {
      for (j=0; j<N; j++)
          printf("%d ", *(psnew->Pa[i]+j));    //输出基拉丁方阵
          printf("\n");
    }
    printf("\n");

    for (k=1; k<fac; k++)
    {
        for (i=0; i<N; i++)
        {
            for (j=0; j<N; j++)
                printf("%d ", *(factor[k]+(*(psnew->Pa[i]+j)-1)));    //输出置换得到的所有其他拉丁方阵
            printf("\n");
        }
        printf("\n");
    }
}