试除法计算最小的N个素数

Trial division

试除法求最小N个素数是一个经典的算法。

可以使用素数分类的性质来加快计算的速度。

除了2以外，素数只能是2n+1，n>=1。

除了2、3和5以外，素数只能是6n+1或6n+5，n>=1。

/*
 * 试除法计算最小的N个素数
 *
 * 素数中除了2以外，均为奇数。
 * 可以用试除法从奇数中挑选出最小的N个素数。
 *
 * 另外，<span style="color: rgb(0, 128, 0);">除了2、3和5之外，</span>素数只可能是6n+1或6n+5,n>=1。
 * 所有有另外一种计算方法，会稍微快一些。
 */

#include <stdio.h>
#include <math.h>

#define N   25

int prime[N] = {2, 3, 5};

void output_result(void);

void primes()
{
    int n = 3;
    int i = 2;
    for(;;) {
        n += 2;
        int j, q =sqrt(n);
        for(j=1; prime[j]<=q; j++)
            if(n % prime[j] == 0)
                break;
        if(prime[j] > q)
            prime[i++] = n;
        if(i == N)
            break;
    }
}

void primes2()
{
    int n1=1, n2 = 5;
    int i = 3;
    for(;;) {
        n1 += 6;
        int j, q =sqrt(n1);
        for(j=1; prime[j]<=q; j++)
            if(n1 % prime[j] == 0)
                break;
        if(prime[j] > q)
            prime[i++] = n1;
        if(i == N)
            break;

        n2 += 6;
        q =sqrt(n2);
        for(j=1; prime[j]<=q; j++)
            if(n2 % prime[j] == 0)
                break;
        if(prime[j] > q)
            prime[i++] = n2;
        if(i == N)
            break;
    }
}

int main(void)
{
    primes();

    output_result();

    primes2();

    output_result();

    return 0;
}

void output_result()
{
    int i;
    for(i=0; i<N; i++)
       printf("%d:%d\n", i+1, prime[i]);
}