改進的篩素數法

2021-08-21 14:42:23 字數 2259 閱讀 1394

最簡單的篩素數法方法就是從2開始,將所以2的倍數去掉,然後從3開始,將3的倍數去掉。根據這樣很容易寫出**,下面**就是是篩素數法得到100以內的素數並儲存到primes陣列中。

const int maxn = 100;

bool flag[maxn];

int primes[maxn / 3], pi;

void getprime_1()

}

可以看出這種會有很多重複訪問,如在訪問flag[2]和flag[5]時會各訪問flag[10]一次。因此最好想方法來減少這種重複訪問,讓flag陣列的每個元素只被訪問一次。可以這樣考慮——簡單的篩素數法是利用乙個素數的倍數必須不是素數,同樣任何乙個數與其它所有素數的乘積必然也不是素數(這是因為每個合數必有乙個最小素因子)。

為了試驗這種想法,先用2到10之間的數來驗證下。

2,3,4,5,6,7,8,9,10      初始時所以flag都是無標記的。

第一步 訪問2,flag[2]無標記所以將2加入素數表中,然後將2與素數表中的所有數相乘得到的數必定不是素數,2*2=4因此標記flag[4]。

2,3,4,5,6,7,8,9,10

第二步 訪問3,flag[3]無標記所以將3加入素數表中,將3與素數表中的所有數相乘得到的數必定不是素數,3*2=6,3*3=9因此標記flag[6]和flag[9]。

2,3,4,5,6,7,8,9,10

第三步 訪問4,flag[4]有標記所以4不加入素數表中,將4與素數表中的所有數相乘得到的數必定不是素數, 4*2=8,4*3=12因此標記flag[8]。

2,3,4,5,6,7,8,9,10

第四步 訪問5,flag[5]無標記所以將5加入素數表中,將5與素數表中的所有數相乘得到的數必定不是素數,5*2=10,5*3=15因此標記flag[10]。

2,3,4,5,6,7,8,9,10

第五步 訪問6,flag[6]有標記所以6不加入素數表中,將6與素數表中的所有數相乘得到的數必定不是素數, 6*2=12,6*3=18,6*5=30。

2,3,4,5,6,7,8,9,10

後面幾步類似,**不難寫出:

const int maxn = 100;

bool flag[maxn];

int primes[maxn / 3], pi;

void getprime_2()

}

這份**對不對了?仔細回顧下分析過程,可以發現有些資料還是被訪問多次了,這當然不是我們希望的結果,我們的要求是讓每個合數僅被它的最小素因子篩去一次。比如12,它的最小素因子是2,所以就只應該被在計算6*2時去訪問,而且不應該在計算4*3時去訪問,同理18也只應該被在計算9*2時去訪問,而且不應該在計算6*3時去訪問。

找到原因後,再來思考如何解決。6*3不行而9*2可以了,是因為6是2的倍數,所以在計算6*2之後就不能再將6與比2大的素數相乘,這些相乘的結果必定會導致重複計算。因此對於任何數來說,如果它如果是該素數的倍數那麼它就不能再與素數表中該素數之後的素數相乘了,如9是3的倍數,所以在9*3之後就不能再去用計算9*5了。因此在**中再增加一行判斷語句:

const int maxn = 100;

bool flag[maxn];

int primes[maxn / 3], pi;

void getprime_2()

}}

想知道這二種篩素數法方法的區別嗎?現在對求2到1億之間的素數進行測試,看看區別到底會有多大,測試**如下:

// 普通的篩素數方法與改進之後的效率對比

// by morewindows( )

#include #include #include #include const int maxn = 100000000;

bool flag[maxn];

int primes[maxn / 3], pi;

// 利用對每個素數的倍數必定不是素數來篩選

void getprime_1()

}// 利用了每個合數必有乙個最小素因子來篩選

void getprime_2() }}

int main()

測試結果如圖所示:

可以看出,效率有4倍之差。改進還是比較可觀。還有空間壓縮技巧來改進後的篩素數法方進行空間壓縮。(位運算)

1.普通的篩素數的原理是乙個素數的倍數必須不是素數。

2.改進的篩素數的原理是每個合數必有乙個最小素因子,根據每個最小素因子去訪問合數就能防止合數被重複訪問。

改進的篩素數法

最簡單的篩素數法方法就是從2開始,將所以2的倍數去掉,然後從3開始,將3的倍數去掉。根據這樣很容易寫出 下面 就是是篩素數法得到100以內的素數並儲存到primes陣列中。cpp view plain copy by morewindows const intmaxn 100 bool flag m...

篩素數法(改進版)

非常感謝作者提供的學習機會,這裡只是copy過來方便自己學習。最簡單的篩素數法方法就是從2開始,將所以2的倍數去掉,然後從3開始,將3的倍數去掉。根據這樣很容易寫出 下面 就是是篩素數法得到100以內的素數並儲存到primes陣列中。cpp view plain copy by morewindow...

改進的素數篩

最簡單的篩素數法方法就是從2開始,將所以2的倍數去掉,然後從3開始,將3的倍數去掉。根據這樣很容易寫出 下面 就是是篩素數法得到100以內的素數並儲存到primes陣列中。cpp view plain copy by morewindows const intmaxn 100 bool flag m...