质数又称素数。一个大于1的自然数,除了1和它自身外,不能被其他自然数整除的数叫做质数;否则称为合数(规定1既不是质数也不是合数)。
普通方法:遍历2~n-1,如果能整除则不是素数,反之则是素数
代码如下
#include<bits/stdc++.h>
using namespace std;
bool ok(int n){
if(n<2){
return 0;
}
for(int i=2;i<n;i++){
if(n%i==0){
return 0;
}
}
return 1;
}
int main()
{
int l,r;
cin>>l>>r;
for(int i=l;i<=r;i++){
if(ok(i)){
cout<<i<<"\n";
}
}
return 0;
}
时间复杂度略优于O(n^2),空间复杂度O(1)。
不难发现,一个数如果能整除i,那它也能整除n/i,这就是优化版的判素数,代码如下(事实上和代码1几乎没有区别,套了一个sqrt)
#include<bits/stdc++.h>
using namespace std;
bool ok(int n){
if(n<2){
return 0;
}
for(int i=2;i<=sqrt(n);i++){
if(n%i==0){
return 0;
}
}
return 1;
}
int main()
{
int l,r;
cin>>l>>r;
for(int i=l;i<=r;i++){
if(ok(i)){
cout<<i<<"\n";
}
}
return 0;
}
时间复杂度优化到了平方根级别,但仍然不够,在10^5左右就会超时,所以就要引入埃氏筛。
埃拉托斯特尼筛法,简称埃氏筛或爱氏筛,是一种由希腊数学家埃拉托斯特尼所提出的一种简单检定素数的算法。要得到自然数n以内的全部素数,必须把不大于根号n的所有素数的倍数剔除,剩下的就是素数。
概括图:
代码如下
#include<bits/stdc++.h>
#define cio ios::sync_with_stdio(0);
using namespace std;
int a[10000005];
void shai(int n){
a[1]=a[0]=1;
for(int i=2;i<=sqrt(n);i++){
for(int j=2;j<=n/i;j++){
a[i*j]=1;
}
}
}
bool ok(int n){
if(!a[n]){
return 1;
}
return 0;
}
int main()
{
int l,r;
cin>>l>>r;
shai(r);
for(int i=l;i<=r;i++){
if(ok(i)){
cout<<i<<"\n";
}
}
return 0;
}
这样可以优化到O(n),但在10^7左右还是会时间超限,还可以做一个小小的优化,每个数从它的平方开始筛,数据越大越有效,10^7左右可以做到原来的33%,代码如下
#include<bits/stdc++.h>
#define cio ios::sync_with_stdio(0);
using namespace std;
int a[10000005];
void shai(int n){
a[1]=a[0]=1;
for(int i=2;i<=sqrt(n);i++){
for(int j=i;j<=n/i;j++){
a[i*j]=1;
}
}
}
bool ok(int n){
if(!a[n]){
return 1;
}
return 0;
}
int main()
{
cio
int l,r;
cin>>l>>r;
shai(r);
for(int i=l;i<=r;i++){
if(ok(i)){
cout<<i<<"\n";
}
}
return 0;
}