什么是求素数
素数指的是因子只有1和本身的数(1不是素数),求解素数在数学上应用非常广泛,而求解n以内的素数也是我们编程时常遇到的问题,在这个问题上,筛选法求解素数运行得非常快。
i在2到n-1之间任取一个数,如果n能被整除则不是素数,否则就是素数
称筛法
筛选法又称筛法,是求不超过自然数N(N>1)的所有质数的一种方法。据说是古希腊的埃拉托斯特尼(Eratosthenes,约公元前274~194年)发明的,又称埃拉托斯特尼筛子。
具体做法是:
先把N个自然数按次序排列起来。1不是质数,也不是合数,要划去。第二个数2是质数留下来,而把2后面所有能被2整除的数都划去。2后面第一个没划去的数是3,把3留下,再把3后面所有能被3整除的数都划去。3后面第一个没划去的数是5,把5留下,再把5后面所有能被5整除的数都划去。这样一直做下去,就会把不超过N的全部合数都筛掉,留下的就是不超过N的全部质数。因为希腊人是把数写在涂腊的板上,每要划去一个数,就在上面记以小点,寻求质数的工作完毕后,这许多小点就像一个筛子,所以就把埃拉托斯特尼的方法叫做“埃拉托斯特尼筛”,简称“筛法”。(另一种解释是当时的数写在纸草上,每要划去一个数,就把这个数挖去,寻求质数的工作完毕后,这许多小洞就像一个筛子。)
普通枚举法:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
#include <iostream>
#include <string>
#include <cmath>
#include <cstring>
using namespace std;
bool isPlain( int x){
if (x<2) return false ;
else {
for ( int i=2;i<x;i++)
{
if (!(x%i))
return false ;
}
}
return true ;
}
int main()
{
int n;
cin>>n;
int cot=0;
for ( int j=0;j<n;j++){
if (isPlain(j)){
cout<<j<<((++cot%7==0)? "\n" : "\t" );
}
}
}
|
筛选法:
原始版本:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
#include <iostream>
#include <string>
#include <cmath>
#include <cstring>
using namespace std;
int main()
{
int n;
cin>>n;
bool * ans= new bool [n];
memset (ans, true , sizeof ( bool )*n); //
ans[0]= false ;
ans[1]= false ;
for ( int i=2;i<n;i++){
if (ans[i]){
for ( int j=i*2;j<n;j+=i){ //倍数取整
ans[j]= false ;
}
}
}
int col = 0;
for ( int i=0;i<n;i++){
if (ans[i]){
cout<<i<< " " ;
}
}
return 0;
}
|
改进版本
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
#include <iostream>
#include <string>
#include <cmath>
#include <cstring>
#include <bitset>
using namespace std;
int main()
{
int n;
cin>>n;
bitset<100000> ans;
ans.set(0);
ans.set(1);
for ( int j=2; j<= sqrt (n); j++)
{
for ( int i=2*j; i < n; i+=j)
{
ans.set(i);
}
}
int cot=0;
for ( int i=0; i<n; i++)
{
if (ans[i]!=1)
{
cout<<i<<((++cot%7==0)? "\n" : "\t" );
}
}
}
|
总结
以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,如果有疑问大家可以留言交流,谢谢大家对服务器之家的支持。
原文链接:http://www.cnblogs.com/dgwblog/p/8035784.html