Sieve of Atkin是一种快速的素数筛选算法,算法比较成熟和简单,http://en.wikipedia.org/wiki/Sieve_of_Atkin中的描述已经非常的细致,作者撰写此文的目的在于,对如何把伪代码转为C代码作一个引导,参考如下的示例。
#include <math.h>
#include <stdio.h>
/* limit ← 1000000 */
#define LIMIT (1000000)
#define FALSE (0)
#define TRUE (~FALSE)
int main(int argc, char* argv[])
{
int sieve_list[LIMIT + 1];
int n, r;
int x, y;
int k, i;
/* is_prime(i) ← false, ∀ i ∈ [5, limit] */
for(n = 5; n <= LIMIT; n++)
sieve_list[n] = FALSE;
/* for (x, y) in [1, √limit] × [1, √limit]: */
for(x = 1; x <= (int)sqrt(LIMIT); x++)
for(y = 1; y <= (int)sqrt(LIMIT); y++)
{
/*
n ← 4x²+y²
if (n ≤ limit) and (n mod 12 = 1 or n mod 12 = 5):
is_prime(n) ← ¬is_prime(n)
*/
n = 4 * x * x + y * y;
if(n <= LIMIT && (n % 12 == 1 || n % 12 == 5))
sieve_list[n] =