此题根据洛谷上的素数方阵的题目给出的代码,稍有改动,以下是链接。
https://www.luogu.com.cn/problem/B3738
样例输入:4
样例输出:
2 3 5 7
37 41 43 11
31 53 47 13
29 23 19 17
实现的原理看下面手绘图,字不怎么样,请见谅。
下面是代码,里面有注释。
#include<stdio.h>
int a[200][200],p[50000],b=1;//这里规定n最大为200。
int prime(int n)//用于找寻素数
{
if(n<=1) return 0;
for(int i=2;i*i<=n;i++){
if(n%i==0) return 0;
}
return 1;
}
int main()
{
int n,k=1,count=0;
scanf("%d",&n);
for(int i=2;count!=n*n&&i<1000000000;i++){
if(prime(i))//提前将需要用到的素数储存到p数组里面,n*n个素数就刚刚好。
{
p[b++]=i;
count++;
}
}
int S=n,N=1,W=1,E=n; //枚举上下左右四个方向。
int f=1;
for(int i=1;i<=n/2+1;i++){ //n为偶数时,循环n/2次,n为奇数时,循环(n+1)/2,直接用n/2+1次即可。
if(f%4==1)//用%操作,完成四个方向的循环。
{
for(int i=W;i<=E;i++)//开始时,模拟成从方阵的上部,从左往右铺设素数,当到达方阵的边界时,for循环会停止。
a[W][i]=p[k++];//保证铺设的素数从小到大。
f++;//进行下一个方向,在方阵右部,从上往下铺设素数。
N++;//此时为了避免重复铺设,上边界由原来的x,变为x+1。
}
if(f%4==2)
{
for(int i=N;i<=S;i++)//。。。
a[i][S]=p[k++];
f++;
E--;//此时方阵的右边界收缩,右边界由原来的x,变为x-1。
}
if(f%4==3)//以下两个方向大致同理
{
for(int i=E;i>=W;i--)
a[S][i]=p[k++];
f++;
S--;
}
if(f%4==0)
{
for(int i=S;i>=N;i--)
a[i][W]=p[k++];
W++;
f++;
}
}
for(int i=1;i<=n;i++){//最后打印出来,大功告成
for(int j=1;j<=n;j++)
printf("%d ",a[i][j]);
puts("");
}
return 0;
}