有关n阶拉丁方阵的问题

#include<stdio.h>
#include<conio.h>

int as[20][20]; /*拉丁方阵在此二维数组上规划*/
int a1[20]; /*每行可出现的数都在里面*/
int yy=0;


int do1(int x,int y,int n)
{
int i=1;
for(;i<x;i++)
{
if(as[y][i]==n)
{
 return 0;
}
 
}

for(i=1;i<y;i++)
{
if(as[i][x]==n)
{
return 0;
}

}

return 1;
}



void print(int a[20][20],int n)
{
int i=1,j=1;
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
if(j==1)
{
printf("\n");
}
printf("%d ",a[i][j]);
}
}

}




void fun(int x,int y,int n) /*x y 为当前处理的坐标*/
{
int i=1;
for(;i<=n;i++)
{
if(do1(x,y,a1[i])) /*do1()是判断此处可否为a1[i]*/
{
as[y][x]=a1[i];

if(x==n&&y==n)
{
yy++;
print(as,n); 
printf("\n"); 
  return;//返回 
}

if(x==n)
{
   fun(1,y+1,n);
}
else
{
  fun(x+1,y,n);
}

}
}
}




int main()
{ int n=0,i=1;
printf("需要求的次方:");
scanf("%d",&n);

for(;i<=n;i++)
{
a1[i]=i;
}

fun(1,1,n);

printf("\n\n\n It is %d.",yy);
return 0;
}