两百万以下所有质数的和
传统的算法运行太慢了，好半天才算出来，求效率的算法

#include <stdio.h>
#include <math.h>
int main()
{
    int i,j,sum=0;
    int flag=1;
    for(i=3;i<20000000;i++)
    {
        flag=1;
        for(j=2;j<=sqrt(i);j++)
        {
            if(i%j==0)
            {
                flag=0;
                break;
            }
        }
        if(flag)
            sum+=i;
    }
    printf("%d\n",sum);
    return 0;
}

参考: http://blog.

不好意思 错了
应该是这个页面：
http://blog.

来学习

好大的取值范围。。

基本看不懂 完了

学习了上面那个大牛的博客里的筛选法 写的一个求质数和的算法 真的好快 贴出来 大家看下
看来算法才是王道啊

/*

 *筛选法求大范围内所有质数和 

 * 

 *思路：1. 筛选法 参考：<span style="color: #008000; text-decoration: underline;">http://blog.[/color] 

 *      2. 用bit位标记大范围内所有自然数，0 不是质数 1 是质数 

 * 

 *流程：1. （从质数2开始）将该质数所有的倍数所在位全部清零 

 *      2. 将该质数值加进sum

 *      3. 从该质数所在位向后查找第一个未置零的位(即下一个质数)

 *           若其值在最问题所在范围内 返回1

 *         否则 退出，sum即所求

 * 

 *子方法： 1. BitSearch 从给定位置向后搜索第一个为1的位，返回位置的序号

 *         2. BitClear  对指定位清零 

 */ 

#include <stdio.h>
#include <malloc.h>
#include <string.h>

typedef unsigned long LCNT;
typedef unsigned long CELL;

#define DAT_SIZE 2000000
#define CELL_SIZE (sizeof(CELL)*8)
#define MEM_SIZE (DAT_SIZE/8 + 1)

LCNT BitSearch(CELL bitSet[], LCNT startPos);
void BitClear(CELL bitSet[], LCNT bitPos);

int main()
{
    LCNT i, primSum = 0, primCnt = 2;
    CELL *primSet = (CELL*)malloc(MEM_SIZE);
    
    if(primSet == NULL){
        printf("内存不足\n");
        return 1;
    }
    else{
        memset(primSet, 0xFF, MEM_SIZE);
        primSet[0] ^= 3;  //质数从2开始，自然数[0,1]所在位清零  
    }
    
    while(primCnt < DAT_SIZE){        
        for(i=2; primCnt*i < DAT_SIZE; i++)
            BitClear(primSet, primCnt*i);
            
        primSum += primCnt;
        primCnt = BitSearch(primSet, primCnt+1);    
    }
    
    printf("The Sum is %ld\n", primSum);
    
    return 0;
}

LCNT BitSearch(CELL bitSet[], LCNT startPos)
{
    LCNT arrPos = startPos/CELL_SIZE, bitPos = startPos%CELL_SIZE;
    
    while(arrPos < MEM_SIZE){
        if(bitSet[arrPos] != 0){
            for(bitPos; bitPos < CELL_SIZE; bitPos++)
                if((bitSet[arrPos] >> bitPos)%2 == 1)
                    goto OUT; 
            bitPos = 0;
        }
        arrPos++;
    }

OUT:
    return arrPos*CELL_SIZE + bitPos;
}

void BitClear(CELL bitSet[], LCNT bitPos)
{
    bitSet[bitPos/CELL_SIZE] &= ~((LCNT)1 << bitPos%CELL_SIZE);
}

接分

这么多啊！！！

不好意思 main函数最后忘了释放内存了 加一句 free(primSet); 就没事了
不过在main函数中 也没有影响 程序结束时会自动释放的 要是在其他函数局部域 就内存泄露了
疏忽啊