#include<stdio.h> #include<math.h> int main() { int N,m; do{ scanf("%d",&N); }while(N<=1); printf("2 "); for(m=3;m<=N;m++){ int i; int tmp=(int)sqrt(m); for(i=2;i<=m;i++) if(m%i==0) break; if(i>tmp) printf("%d ",m); } return 0; } 优化该代码
时间: 2023-12-05 13:03:39 浏览: 132
上述代码是朴素的暴力枚举法,每次判断一个数是否为素数时都需要从2到该数的平方根进行循环,时间复杂度为O(n*sqrt(n)),效率较低。下面介绍两种优化方法:
1. 埃氏筛法优化
我们可以使用埃氏筛法来优化。具体做法是:首先定义一个长度为N+1的数组is_prime[],其中is_prime[i]表示i是否是素数,初始时全部初始化为true。然后从2开始枚举每个素数i,将i的倍数j(j≥2且i*j≤n)标记为合数,即令is_prime[i*j]=false。最后输出所有is_prime[i]为true的i,即为小于等于n的所有素数。时间复杂度为O(n*loglog(n)),效率更高。
下面是优化后的代码:
```c
#include <stdio.h>
#include <stdlib.h>
#include <stdbool.h>
#include <math.h>
#define MAXN 1000000
int main()
{
int n;
bool is_prime[MAXN + 1] = {false}; // 初始化为false
scanf("%d", &n);
printf("2 ");
for (int i = 3; i <= n; i += 2) // 从3开始,每次加2
{
if (!is_prime[i])
{
printf("%d ", i);
for (int j = i * i; j <= n; j += 2 * i) // 从i的平方开始标记,每次加2*i
{
is_prime[j] = true;
}
}
}
return 0;
}
```
2. 埃氏筛法进一步优化
在埃氏筛法的基础上,我们还可以进一步优化,减少时间和空间的使用。具体做法是:首先定义一个长度为N+1的数组is_prime[],其中is_prime[i]表示i是否是素数,初始时全部初始化为true。然后从2开始枚举每个素数i,将i的倍数j(j≥2且i*j≤n)标记为合数,即令is_prime[i*j]=false。最后输出所有is_prime[i]为true的i,即为小于等于n的所有素数。
优化的关键在于:对于一个素数i,它的倍数2i、3i、4i、……已经被之前的素数标记过了,因此,在枚举素数i时,我们可以从i的平方开始标记,而不是从2i开始标记。这样可以减少一半的标记次数。
下面是优化后的代码:
```c
#include <stdio.h>
#include <stdlib.h>
#include <stdbool.h>
#include <math.h>
#define MAXN 1000000
int main()
{
int n, cnt = 0;
bool is_prime[MAXN + 1] = {false}; // 初始化为false
scanf("%d", &n);
printf("2 ");
for (int i = 3; i <= n; i += 2) // 从3开始,每次加2
{
if (!is_prime[i])
{
printf("%d ", i);
cnt++;
if (cnt % 10 == 0) // 每10个数换行
{
printf("\n");
}
for (int j = i * i; j <= n; j += 2 * i) // 从i的平方开始标记,每次加2*i
{
is_prime[j] = true;
}
}
}
return 0;
}
```
阅读全文
相关推荐
![-](https://img-home.csdnimg.cn/images/20241231044955.png)
![-](https://img-home.csdnimg.cn/images/20241231044901.png)
![-](https://img-home.csdnimg.cn/images/20241231044930.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)