經常有初學者詢問求解n內所有素數(質數)的問題,對此,網上的解答也很多,但很多要麼不夠專業,要麼只有程式沒有演算法解析,所以三藏大廈對此問題做個小結,**一下求解素數的常見演算法,同時給出相應的c語言程式及其解析。為了方便初學者理解,本文將從易到難闡述不同演算法,高手可以直接看後面的高效演算法
質數的定義
乙個數,如果只有1和它本身兩個因數,這樣的數叫做質數,又稱素數。
試除判斷法
演算法描述
:從上述定義可知,素數不能被1和它本身之外的數整除,所以,判斷乙個數x是否素數只要看它是否能被2~sqrt(x)間的數整除即可;而求n內所有素數則是迴圈重複上述過程。
c語言實現
:複製內容到剪貼簿
**:#include
#include
#define n 100000
//簡單試除判斷法
ver1
int s i m p l edivisionv1(int n)
// i
不能被2~sqrt(i)
間的數整除
,素數也
if (j > sqrt(i))}//
因輸出費時,且和演算法核心相關不大,故略
//釋放記憶體,別忘了傳說中的記憶體洩漏
free(primes);
return count;
}void main()
計算結果:
[100000]以內素數個數:9592, 計算用時:468毫秒
[1000000]以內素數個數:78498, 計算用時:10859毫秒
[5000000]以內素數個數:348513, 計算用時:103560毫秒
噢噢,算算十萬還行,百萬就10秒多了,而且時間增長很快,這不行,得優化一下!
優化分析
:仔細研究一下s i m p l edivisionv1我們可以發現以下幾個問題:
1.在迴圈條件中重複呼叫sqrt(i)顯然是比較浪費時間的
2.判斷素數,真的需要拿2~sqrt(i)間的所有整數去除嗎?我們知道,合數都可以分解成若干質數,所以只要2~sqrt(i)間的質數不能整除i即可
根據上面兩點,我們可將s i m p l edivisionv1公升級為s
i m p l edivisionv2,如下
複製內容到剪貼簿
**://
簡單試除判斷法
ver2
int s i m p l edivisionv2(int n)
// i
不能被2~sqrt(i)
間的素數整除,自然也不能被其他數整除
,素數也
if (j == stop)}//
因輸出費時,且和演算法核心相關不大,故略
//釋放記憶體,別忘了傳說中的記憶體洩漏
free(primes);
return count;
}然後將main中呼叫的函式替換為s
i m p l edivisionv2,在看一下執行結果:
[100000]以內素數個數:9592, 計算用時:46毫秒
[1000000]以內素數個數:78498, 計算用時:546毫秒
[5000000]以內素數個數:348513, 計算用時:3515毫秒
[10000000]以內素數個數:664579, 計算用時:8000毫秒
很開心的看到,經過優化,速度提高了幾十倍,尤其是時間增長曲線的坡度變小了,n值越大,v2函式比v1的效率就越高
對於試除判斷這種質數演算法來說,三藏認為s i m p l edivisionv2基本已經接近極限,不大可能有量級上的突破了,有興趣的朋友可以自己進一步優化。初學者除了參看上述例子外,可以嘗試做各種修改及細節優化,也可以將除法變乘法,多加練習是學習程式設計的好方法。
雖然,上例中v2已經比v1快了很多了,但隨著n的增大,耗時還是不少,那麼我們還有更好的方法嗎?
計算素數演算法的一些優化(程式設計珠璣閱讀感想)
這篇文章同樣是關於我讀程式設計珠璣的一些體會與感想,這是關於計算素數演算法優化的一些感想。在學c語言的時候,用程式列印出某個範圍內的所有素數的練習題大家都不陌生,我當時想到的演算法類似下面這樣的。int prime int n 但是作者還是覺得迴圈是累贅,有提供了下面的演算法,思想大概是這樣,對能被...
計算一定範圍內素數個數的演算法
問題 給定乙個大整數n,計算開區間 1,n 的素數有多少?include include include include include undef true define true 1 undef false define false 0 typedef int bool32 a不能整除小於等於根...
Java計算素數
素數又稱質數,指的是,除了1和它本身,沒有第三個數能夠整除它。例如 注意 1本身只能夠被1整除,不算質數。根據質數 素數 的定義不難得出,要計算乙個數是不是質數,需要明確是不是除了1和本身以外,還有其他除數。由此可以有乙個計算思路 大致流程圖如下 根據上面流程圖分析,下圖為計算 判斷乙個數是不是質數...