ehi*_*ime 1 php algorithm bash primes
嘿伙计们,我正试图从我的发电机中挤出更多的质数,但是在第60次运行中已达到+ - 14,000,000左右的极限值,我希望尽可能推高可能达到25-30万.关于如何完成这一壮举的任何想法?
这是我的PHP代码
<?php
$i = 2;
$primes = array();
while(true)
{
$prime=true;
$sqrt=floor(sqrt($i));
foreach($primes as $num)
{
if($i%$num==0)
{
$prime=false;
break;
}
if($num>$sqrt) break;
}
if($prime) echo "$i\n";
$i++;
}
Run Code Online (Sandbox Code Playgroud)
这是运行它的小bash脚本
#!/bin/bash
outfile="$1.out"
`php $1 > $outfile &`
sleep 60
killall php
Run Code Online (Sandbox Code Playgroud)
编辑
这是昨晚的pcntl_forked版本只是为了看看会更快.出于一些奇怪的原因,我做得越少越好,如果是,我发现瓶颈来自gmp_strval(gmp_nextprime($start))哪个非常慢.
<?php
$primeCount = 0;
for ($i = 1; $i <= 1; ++$i) {
$start = $i;
$pid = pcntl_fork();
if (!$pid) {
while(true) echo $start = gmp_strval(gmp_nextprime($start)) . "\n";
exit($i);
}
}
while (pcntl_waitpid(0, $status) != -1) {
$status = pcntl_wexitstatus($status);
#echo "Child $status , $primeCount completed\n";
}
die;
Run Code Online (Sandbox Code Playgroud)
小智 8
除了2和3之外,所有素数都是高于6的倍数或者低于6的倍数.这应该将您的搜索空间减少到大约1/3而不是每次递增1.用2,3筛选筛子并使用缩小的搜索空间来加快时间.
如果您需要一个简单的解释:
1 + n*6 (This is one above, and I state it's got primes)
2 + n*6 (multiple of 2)
3 + n*6 (multiple of 3)
4 + n*6 (multiple of 2)
5 + n*6 (This is one below, and I state it's got primes)
6 + n*6 (multiple of 2,3,6)
Run Code Online (Sandbox Code Playgroud)