第一部分:素性判定
素数【质数】——就是约数只有1和它本身的数,所以从2到n-1,如果存在某个数m,使得n%m==0 那n就不是质数了。。。
另外,如果i是n的约数,那么n/i 也一定会是n的约数,所以只需要在2到根号(n)中去找是否存在能整除n的数,如果存在,那么n不素数,否则就是
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
using namespace std;
bool is_prime(int x)
{
if(x<=1)
return false;
int k=(floor)(sqrt((double)x)+0.5); 刘汝佳大神教的小技巧,能避免系统误差。。。
for(int i=2;i<=k;i++)
if(x%i==0) //如果能整除,那么i是x的约数
return false; //找到了1和本身之外的约数,直接返回,肯定不是素数了。。。
return true; //能到达这一步,说明上面的测试全部通过,没有其他约数
}
int main()
{
int x;
while(scanf("%d",&x)!=EOF)
{
if(is_prime(x))
printf("Yes\n");
else
printf("No\n");
}
return 0;
}
第二部分:埃氏筛法
定义一个数组[0,n)—— bool is_peime[n] ; 表示i是否为素数
初始化:is_prime[0]和is_peime[1]为false,表示0,1不是素数。。。其他全为true
从2开始,如果这个数时素数,那么把n范围内,2的倍数全部标记为false,一定部位素数。。。
3也为素数,然后把3的倍数全部标记为false。。。以此类推。。。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
using namespace std;
#define MAX_N 100000
bool is_prime [MAX_N];
int main()
{
for(int i=0;i<MAX_N;i++)
is_prime[i]=true; //初始化全为素数
is_prime[0]=is_prime[1]=false; //0和1肯定不是素数,标记为false
int N; //求N范围内的素数
while(scanf("%d",&N)!=EOF)
{
for(int i=2;i<=N;i++) //从2开始,0和1果断不是素数。。。
if(is_prime[i]) //如果是素数
for(int j=2;i*j<=N;j++)
is_prime[i*j]=false; //它的倍数全部标记为false
int countn=0;
for(int i=0;i<=N;i++)
{
if(is_prime[i]) //如果是质数
{
printf("%5d",i);
countn++;
if(countn%10==0) //每行打印10个
printf("\n");
}
}
printf("\n");
}
return 0;
}