打印素数表
开根快速判定法
方法一:算法复杂度 O(sqrt(n))
#include <cmath>
#include <cstdio>
const int maxn = 200;
bool P[maxn]= {0};
bool isPrime(int n) //素数判断
{
int sqr = (int)sqrt(1.0*n); //sqrt的作用是为一个浮点数开根号,乘1.0就是为了使n变为浮点数 需要添加 math.h 头文件
if(n<=1) return false;
else
{
for(int i=2; i<=sqrt(n); i++)
{
if(n%i==0) return false;
}
return true;
}
}
int Prime[maxn],PNum = 0;
void Find_Prime(int n)
{
for(int i=1; i<n; i++)
{
if(isPrime(i)==true)
{
Prime[PNum++]=i;
P[i]=true;
}
}
}
int main()
{
Find_Prime(101);
for(int i=0; i<PNum; i++)
{
printf("%d ",Prime[i]);
}
return 0;
}
埃氏筛法
方法二:埃氏筛法 算法复杂度O(n log logn)
const int maxn=101;
bool P[maxn]={0}; //false代表不需要筛出去,true代表需要被筛出去
int Prime[maxn]={0},Pnum=0; //素数表,素数个数
//唯一需要确定2是素数
void Find_Prime()
{
int m=2;
for(int i=2;i<maxn;i++)
{
if(P[i]==false) //代表是素数不需要被筛出去
{
Prime[Pnum++]=i;
for(int j=i;j<maxn;j=i*m)
{
P[j]=true;
m++;
}
}
}
}
int main()
{
Find_Prime();
for(int i=0;i<Pnum;i++)
{
printf("%d ",Prime[i]);
}
return 0;
}
pat1013 数素数
原题链接
思路:埃氏筛法 难点可能在输出格式吧 找规律 我的比较麻烦
输出优化 参考原贴
注意点:
1.小技巧:空格肉眼看不见,如果报错格式错误,可以把空格换成“#”等其他符号,看看是哪里多了空格
2.考虑到不知道第10^4个素数多大,
不妨将测试上限设的大一点。也可以先用程序测试一下第10^4个素数是多少然后再用这个数作为上限;
#include <iostream>
#include <cmath>
#include <cstdio>
const int maxn=1000000;
bool P[maxn]= {0};
int Prime[maxn]= {0},Pnum=0;
void Find_Prime()
{
for(int i=2; i<=maxn; i++)
{
if(P[i]==false)
{
int m=1;
Prime[Pnum++]=i;
for(int j=i; j<=maxn; j=i*m)
{
P[j]=true;
m++;
}
}
}
}
int main()
{
int m,n;
int cou=0;
scanf("%d %d",&m,&n);
Find_Prime();
// for(int i=m-1; i<n; i++)
// {
// cou++;
// if(n-m+1<10)
// {
// if(i==n-1) printf("%d",Prime[i]);
// else printf("%d ",Prime[i]);
// }
// else
// {
// if(i==n-1) printf("%d",Prime[i]);
// else if(cou%10==0)
// {
// printf("%d",Prime[i]);
// printf("\n");
// }
// else printf("%d ",Prime[i]);
// }
//
// }
for(int i=m-1;i<n;i++)//优化
{
cou++;
printf("%d",Prime[i]);
if(cou%10!=0&&(i!=n-1)) printf(" ");
else printf("\n");
}
return 0;
}