素数之恋

　　 http://primes.utm.edu/bios/page.php?id=1498 http://primes.utm.edu/ 以上网页说的是目前找的最大素数是2 ^ 43112609 - 1 下面是在Unix-Center.net的机器上验证程序，可以验证1-E9内的质数，如果改为E10的话程序会出现Segmentation Fault. 希望有人能告诉我再往上该怎么验证。

　　

-bash-3.2$ more prime.c
#include <stdio.h>
#include <time.h>
#include <math.h>
#include <string.h>
#include <stdlib.h>

#define  E9        1072000000 
#define  E10       10000000000
#define  E11       100000000000
#define  E12       1000000000000

char table[E9];

int num = 0;
//int prime[E8];

void cal_table()
{
        int j, k;
        time_t t1, t2;

        t1 = time(NULL);

        memset(table, 1, E9);
        for (j = 2 ; j < E9; j++) {
                if ( table[j] ) {
                        num++;
                        for ( k = j + j; k < E9; k += j )
                        {
                                table[k] = 0;
                        }
                }
        }
        t2 = time(NULL);

        printf("Totaly %d primes until E9, cost %d time_t.\n", num, t2 - t1);
}

/*
void init_prime()
{
        int i;
        num = 0;
        printf("Init prime start...");
        for (i = 2; i < 2*E9; i++) {
            if (table[i])
                prime[num++] = i;
        }
        printf("Init prime finished...");
}
*/


/*
int e9_to_e12(int min, int max)
{
    int i, j, k;
    time_t t1, t2;

    t1 = time(NULL);

    for (i = 0; i < 1000; i++) {
        memset(table, 1, E9);
        for (j = 2 ; j < E9; j++) {
                if ( table[j] ) {
                        num++;
                        for ( k = j + j; k < E9; k += j )
                        {
                                table[k] = 0;
                        }
                }
        }
    }

    t2 = time(NULL);

    printf("Totaly %d primes until 1 000 000 000, cost %d time_t.\n", num, t2 - t1);
}
*/


/*
void e9_to_e10()
{
    int i;
    for (i = E9 + 1; i < E10; i++) {
    }
}
*/


/*
int isprime(int i)
{
        int j, k;
        k = (int)sqrt(i);
        for (j = 0;  prime[j] <= k; j++) {
                if (i % prime[j] == 0)
                        return 0;
        }
        return 1;
}
*/

int main()
{
    int i;

    cal_table();
    //init_prime();

    return 0;
}

-bash-3.2$ gcc prime.c
-bash-3.2$ ./a.out
Totaly 54316419 primes until E9, cost 179 time_t.
-bash-3.2$

使人感兴趣的是发现黎曼那相当令人悲哀的和有点让人可怜的特性。但是这些可能只是这个人可见的外表和举止。在羞怯、孤独外表的内部，是才华卓越和惊人勇敢的头脑。无论他在偶然见到他的人面前显得怎样胆怯和倦怠，黎曼的数学却具有拿破仑的一场战役那样无畏的冲击力和能量。当然，他的数学方面的朋友和同事都知道这一点，并且尊敬他。