筛选一万亿内的素数.docVIP

筛选一万亿内任意一个区段的素数 函数参数：as_1? ? ? ? ? ? ? decimal? value? ? ? 筛选范围区段起点? ? ? ? ? as_2? ? ? ? ? ? ? decimal? value? ? ? 筛选范围区段终点 ? ? ? ? as_return[]? ? ? ? decimal? reference? 存放种子素数? ? ? ? ? as_return_value[]? decimal? reference? 存放筛选结果 调用详解：1.必须定义的变量? dec as_return[],as_return_value[] ? ? ? ? 2.as_return,as_return_value切勿赋值。 ? ? ? ? 3.筛选区段一般不大于100万为佳，即as_2－as_1=1000000, ? ? ? ? ? 函数允许区段超过100万。 ? ? ? ? 4.调用本函数f_prime_number(as_1,as_2,as_return,as_return_value) ? ? ? ? 5.as_return_value[]就是生成的某区段内的素数 函数返回：=0 所找到素数的个数 ? ? ? ? ? 0 error 实际调用：求一万亿最后一个一万区段的素数 ? ? ? ? f_prime_number(999999990001.0,1000000000000.0,as_return,as_return_value) 要点及原理： ? ? 一.筛100万内素数： ? ? ? 1.筛法原理，一个素数不能被它平方根以内的所有素数的倍数筛掉，它就是素数。 ? ? ? 2.100万的平方根是1000，所有只用1000以内的素数就可以筛出100万内所有素数 ? ? ? 3.根据上面所说，程序分两步，第一步1000以内进行筛选，第二部把1000以上的筛选 ? ? ? ? 剩下的素数挑选出来 ? ? ? 4.第一步的具体思路： ? ? ? ? a.定义数组：p[1000001]，数组下标表示数本身，初值为0，筛选时把素数的整数倍置为-1， ? ? ? ? 剩下为0的元素即为素数 ? ? ? ? b.2是素数，其倍数是偶数，由程序滤掉，不再筛 ? ? ? ? c.筛选从种子的平方开始，因为其平方以内的的数以被前面的素数筛过了 ? ? ? ? d.筛选的步长用种子的2倍，因为用种子作步长有一半是偶数，筛过了 ? 二.筛一万亿内的素数 ? ? ? 1.一万亿的内的素数不宜一下子全筛出来，可以一次以100万为单位，一个区段一个区段筛。 ? ? ? 2.算法要点如下： ? ? ? ? a.一万亿的平方根为100万，用100万以内的素数可以筛出一万亿内的所有素数，所以函数 ? ? ? ? 首先自行计算一次100万以内的素数存入数组as_return。 ? ? ? ? b.定义数组p[1000001],其下标与区段内的自然序号相对应，元素值0为素数，1不是 ? ? ? ? c.根据as_return内的素数筛选至种子的平方大于筛选区间的上限 ? ? ? ? d.计算出每个种子在指定区间内筛掉的第一个数，然后按种子长度筛下去。 ? ? ? ? 程序中语句为“j=as_return -mod(k1,as_return)”，如i=12，as_return＝37 ? ? ? ? 则j＝21，所以令p[21]＝1，即将k1+21筛掉（k1为区间起始值）再按步长37筛下去 ? ? ? ? e.偶数仍由程序筛掉，所以素数2不用 ? ? ? ? f.数组p元素加上区间起始值k1就是实际要求的素数 ? 三.求更大的素数，如求1000万平方(100万亿)内的素数 ? ? ? 方法一，把本函数涉及到100万的地方扩大10倍，即ls_rangp,这样 ? ? ? ? ? ? ? 就可以求出1000万平方内的所有素数，但会有一定的时间和资源开销，具体要看 ? ? 实际硬件配置。 ? ? ? 方法二，使用本函数计算出1000万内的素数作为种子存放在数组或文件中，再使用 ? ? ? ? ? ? ? 本函数的后半部分（内部函数2）筛选1万亿以上部分的素数。 ? ? ? 方法二的实际调用示例： ? ? ? //求1万亿后第一个一万区段的素数 ? ? ? dec as_return[],as_return_value[],temp[],k ? ? ? f_prime_number(1as_return,as_return_value) ? ? ? as_return=as_return_valu