python素數(shù)篩選法淺析
原理:
素數(shù),指在一個大于1的自然數(shù)中,除了1和此整數(shù)自身外,不能被其他自然數(shù)整除的數(shù)。在加密應用中起重要的位置,比如廣為人知的RSA算法中,就是基于大整數(shù)的因式分解難題,尋找兩個超大的素數(shù)然后相乘作為密鑰的。一個比較常見的求素數(shù)的辦法是埃拉托斯特尼篩法(the Sieve of Eratosthenes) ,說簡單一點就是畫表格,然后刪表格,如圖所示:

從2開始依次往后面數(shù),如果當前數(shù)字一個素數(shù),那么就將所有其倍數(shù)的數(shù)從表中刪除或者標記,然后最終得到所有的素數(shù)。
有一個優(yōu)化:
標記2和3的倍數(shù)的時候,6被標記了兩次。所以從i的平方開始標記,減少很多時間。
比如3的倍數(shù)從9開始標記,而不是6,并且每次加6。
除了2以外,所有素數(shù)都是奇數(shù)。奇數(shù)的平方還是奇數(shù),如果再加上奇數(shù)就變成了偶數(shù)一定不會是素數(shù),所以加偶數(shù)(2倍素數(shù))。
預先處理了所有偶數(shù)。
注意:1既不是素數(shù)也不是合數(shù),這里沒有處理1。
#! prime.py
import time
def primes(n):
P = []
f = []
for i in range(n+1):
if i > 2 and i%2 == 0:
f.append(1)
else:
f.append(0)
i = 3
while i*i <= n:
if f[i] == 0:
j = i*i
while j <= n:
f[j] = 1
j += i+i
i += 2
P.append(2)
for i in range(3,n,2):
if f[i] == 0:
P.append(i)
return P
def isPrime(n):
if n > 2 and n%2 == 0:
return 0
i = 3
while i*i <= n:
if n%i == 0:
return 0
i += 2
return 1
def primeCnt(n):
cnt = 0
for i in range(2,n):
if isPrime(i):
cnt += 1
return cnt
if __name__ == '__main__':
start = time.clock()
n = 10000000
P = primes(n);
print("There are %d primes less than %d"%(len(P),n))
#for i in range(10):
# print(P[i])
print("Time: %f"%(time.clock()-start))
#for n in range(2,100000):
# if isPrime(n):
# print("%d is prime"%n)
#print("%d is "%n + ("prime" if isPrime(n) else "not prime"))
start = time.clock()
n = 1000000
print("There are %d primes less than %d"%(primeCnt(n),n))
print("Time: %f"%(time.clock()-start)
用素數(shù)篩選法求1千萬以內的素數(shù)用了5.767s,
普通素數(shù)判斷法求1百萬以內的素數(shù)用了9.642s,
用C++素數(shù)篩選法求1億以內的素數(shù)用了0.948s,
用C++普通素數(shù)判斷法求1千萬以內的素數(shù)用了3.965s,
可見解釋語言確實比編譯語言慢很多。
附C++程序,用了位壓縮優(yōu)化空間
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
#define N 100000001
unsigned f[(N>>5)+5];
int p[5761456],m;
void init()
{
int i,j;
for(i=4;i<N;i+=2)
f[i>>5]|=1<<(i&0x1F);
p[m++]=2;
for(i=3;i*i<N;i+=2)
if(!(f[i>>5]&(1<<(i&0x1F))))
{
p[m++]=i;
for(j=i*i;j<N;j+=i+i)
f[j>>5]|=1<<(j&0x1F);
}
for(;i<N;i+=2)
if(!(f[i>>5]&(1<<(i&0x1F))))
p[m++]=i;
}
int is_prime(int n)
{
int i;
for(i=0;p[i]*p[i]<=n;i++)
if(n%p[i]==0)
return 0;
return 1;
}
int isPrime(int n)
{
if(n>2 && n%2==0)
return 0;
int i=3;
while(i*i<=n)
{
if(n%i==0)
return 0;
i+=2;
}
return 1;
}
int main()
{
int n=0,i;
clock_t st=clock();
init();
/*for(i=2;i<10000000;i++)
if(isPrime(i))
n++;*/
printf("%d %dms\n",m,clock()-st);
/*while(~scanf("%d",&n),n)
{
i=lower_bound(p,p+m,n+1)-p;
printf("%d\n",i);
}*/
return 0;
}
以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。
- 使用Python判斷質數(shù)(素數(shù))的簡單方法講解
- Python編程判斷一個正整數(shù)是否為素數(shù)的方法
- python求素數(shù)示例分享
- Python求出0~100以內的所有素數(shù)
- Python素數(shù)檢測的方法
- python判斷所輸入的任意一個正整數(shù)是否為素數(shù)的兩種方法
- Python 判斷是否為質數(shù)或素數(shù)的實例
- Python實現(xiàn)輸出某區(qū)間范圍內全部素數(shù)的方法
- Python numpy.array()生成相同元素數(shù)組的示例
- Python求區(qū)間正整數(shù)內所有素數(shù)之和的方法實例
相關文章
使用PyCharm配合部署Python的Django框架的配置紀實
這篇文章主要介紹了使用PyCharm配合部署Python的Django框架的配置紀實,PyCharm是一款強大的Python的IDE,需要的朋友可以參考下2015-11-11
Python實現(xiàn)批量讀取HDF多波段柵格數(shù)據(jù)并繪制像元直方圖
這篇文章主要為大家詳細介紹了如何基于Python語言gdal模塊,實現(xiàn)多波段HDF柵格圖像文件的讀取、處理與像元值可視化(直方圖繪制)等操作,需要的可以參考一下2023-03-03

