怎么判斷一個(gè)數(shù)是素?cái)?shù)？

常規(guī)的方法是枚舉從2開始的數(shù)，看看是否能被整除。

但是，如果要判斷的數(shù)很多的時(shí)候，那么效率會十分低下.。。。

一個(gè)優(yōu)化的方法是不用判斷比這個(gè)數(shù)小的所有數(shù)（到平方根位置），而是判斷比他小的所有素?cái)?shù)。

如果所有小于他的素?cái)?shù)都不能將其整除，那么他就是素?cái)?shù)。

那么如何快速得到小于他的素?cái)?shù)?

有個(gè)Eratosthenes

把最后決定采用的寫法放開頭。。不是素?cái)?shù)的標(biāo)記為1

for(int i=2;i*i<=1000000;i++) { if(!isprime[i])for(int j=i;j*i<=1000000;j++) isprime[i*j]=1; }

根據(jù)劉汝佳的書的代碼如下：

可是當(dāng)n為100萬時(shí)候竟然直接停止工作！

#include<iostream> using namespace std; const int MAXN=1000+10; int main() {int n;cin>>n;int vis[MAXN],primer[MAXN],count=0;memset(vis,0,sizeof(vis));for(int i=2;i<=n;i++){if(!vis[i]) primer[count++]=i; //將素?cái)?shù)存入for(int j=i*i;j<=n;j+=i) //標(biāo)記是否是素?cái)?shù)，如果不存入的話,直接判斷vis[target]是否為0也可以。 vis[j]=1; }for(int i=0;i<count;i++)cout<<primer[i]<<" ";cout<<count<<endl; }

判斷一個(gè)數(shù)是否是兩個(gè)素?cái)?shù)的乘積。

多組數(shù)據(jù)不超過1000組，每組數(shù)據(jù)一行一個(gè)整數(shù)n(2 =< n <=1000000)。

n是否是兩個(gè)素?cái)?shù)的乘積，是輸出Yes，不是輸出No。

下面的寫法是非素?cái)?shù)標(biāo)記為1

聽隊(duì)友的改進(jìn)如下：

#include<cstdio> const int MAXN=1000000+10; bool vis[MAXN]={0};int main() {int count=0;for(int i=2;i*i<=1000000;i++) { for(int j=i;j*i<=1000000;j++) vis[i*j]=1; } int n;while(~scanf("%d",&n)){bool ok=false;for(int i=2;i * i <=n;i+=1){if(vis[i]) continue;int t=n / i;if(t*i==n && !vis[t] ){ok=true;break;}}if(ok)printf("Yes\n");elseprintf("No\n");} }

FZU1563

http://acm.fzu.edu.cn/problem.php?pid=1563

可是在FZU這題直接TLE!搜題解得到如下寫法

#include<cstdio> const int MAXN=1000000+10; const int INF=1000000; bool isprime[MAXN]={0}; int prime[MAXN]; int main() {int pnum=0;for(int i=2;i<=INF;i++) { if(!isprime[i]) prime[pnum++]=i; for(int j=0;j<pnum&&prime[j]*i<=INF;j++) { isprime[prime[j]*i]=1; if(i%prime[j]==0) break; } } int n;scanf("%d",&n);while(n--){int a,b;int cnt=0;scanf("%d%d",&a,&b);for(int i=a;i<=b;i++)if(isprime[i]==0)cnt++;printf("%d\n",cnt);} }

最后改進(jìn)如下：

#include<cstdio> const int MAXN=1000000+10; const int INF=1000000; bool isprime[MAXN]={0}; int prime[MAXN]; int main() {for(int i=2;i*i<=1000000;i++) { if(!isprime[i])for(int j=i;j*i<=1000000;j++) isprime[i*j]=1; } int n;scanf("%d",&n);while(n--){int a,b;int cnt=0;scanf("%d%d",&a,&b);for(int i=a;i<=b;i++)if(isprime[i]==0)cnt++;printf("%d\n",cnt);} }

轉(zhuǎn)載于:https://www.cnblogs.com/murmured/p/5004343.html

總結(jié)

以上是生活随笔為你收集整理的素数表（Eratosthenes）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。