JFIF$        dd7 

Viewing File: /usr/share/doc/perl-Math-BigInt/examples/prime.pl

#!/usr/bin/perl -w

BEGIN { unshift @INC, '../lib'; }       # uncomment to use old, org version

$| = 1;

use Math::BigInt;

# this is a complicated version of the prime number sieve.
# It is not optimized (since we want to benchmark as many features as
# possible).

$amount = Math::BigInt->new( shift || 1000000 );

@primes = (1,1,0); # any not defined number is prime, 0,1 are not, but 2 is
my $prime = Math::BigInt->new (3);      # start

# the loop below is faster in the old version than in the new, since it is
# the worst case for new lib: small numbers and lot's of bstr()/new().
# It also slows down the benchmark too much so we use slightly faster int here
$r = 0; my $a = $amount->numify();
for ($i = 3; $i < $a; $i++)             # int version
  {
  $primes[$i] = $r; $r = 1-$r;
  }

# find primes
OUTER:
while ($prime < $amount)
  {
  # find first unmarked, it is the next prime
  $cur = $prime;
  while ($primes[$cur])
    {
    $cur += 2; last OUTER if $cur >= $amount;   # no more to do
    }
  # $cur is now new prime
  $str = "$cur"; $str =~ s/\+//; # unify output for comapre
  #print "$str $prime $amount\n";
  # now strike out all multiples of $cur
  $add = $cur*2;
  $prime = $cur + 2;                    # next round start two higher
  $cur += $add;
  while ($cur < $amount)
    {
    $primes[$cur] = 1; $cur += $add;
    }
  }

$i = 0;
foreach (@primes)
  {
  push @real_primes, $i if $primes[$i] == 0;
  $i++;
  }
# uncomment to print em:
# foreach (@real_primes) { print "$_\n"; }

print "last prime: $real_primes[-1]\n";

# check against text
open FILE, '1000.txt' or die "Can't read 1000.txt: $!";
my @test;
while (<FILE>)
  {
  next if /^#/;
  next if /^\s*$/;
  $_ =~ s/\s+/ /g;
  $_ =~ s/^\s+//;
  $_ =~ s/\s+$//;
  push @test, split /\s+/,$_;
  }
close FILE;
my $i = 0;
foreach (@real_primes)
  {
  print "oups: $i: $test[$i] != $real_primes[$i]\n"
    if $test[$i] != $real_primes[$i]; $i++;
  last if $i >= 1000;
  }

print "done\n";
Back to Directory  nL+D550H?Mx ,D"v]qv;6*Zqn)ZP0!1 A "#a$2Qr D8 a Ri[f\mIykIw0cuFcRı?lO7к_f˓[C$殷WF<_W ԣsKcëIzyQy/_LKℂ;C",pFA:/]=H  ~,ls/9ć:[=/#f;)x{ٛEQ )~ =𘙲r*2~ a _V=' kumFD}KYYC)({ *g&f`툪ry`=^cJ.I](*`wq1dđ#̩͑0;H]u搂@:~וKL Nsh}OIR*8:2 !lDJVo(3=M(zȰ+i*NAr6KnSl)!JJӁ* %݉?|D}d5:eP0R;{$X'xF@.ÊB {,WJuQɲRI;9QE琯62fT.DUJ;*cP A\ILNj!J۱+O\͔]ޒS߼Jȧc%ANolՎprULZԛerE2=XDXgVQeӓk yP7U*omQIs,K`)6\G3t?pgjrmۛجwluGtfh9uyP0D;Uڽ"OXlif$)&|ML0Zrm1[HXPlPR0'G=i2N+0e2]]9VTPO׮7h(F*癈'=QVZDF,d߬~TX G[`le69CR(!S2!P <0x<!1AQ "Raq02Br#SCTb ?Ζ"]mH5WR7k.ۛ!}Q~+yԏz|@T20S~Kek *zFf^2X*(@8r?CIuI|֓>^ExLgNUY+{.RѪ τV׸YTD I62'8Y27'\TP.6d&˦@Vqi|8-OΕ]ʔ U=TL8=;6c| !qfF3aů&~$l}'NWUs$Uk^SV:U# 6w++s&r+nڐ{@29 gL u"TÙM=6(^"7r}=6YݾlCuhquympǦ GjhsǜNlɻ}o7#S6aw4!OSrD57%|?x>L |/nD6?/8w#[)L7+6〼T ATg!%5MmZ/c-{1_Je"|^$'O&ޱմTrb$w)R$& N1EtdU3Uȉ1pM"N*(DNyd96.(jQ)X 5cQɎMyW?Q*!R>6=7)Xj5`J]e8%t!+'!1Q5 !1 AQaqё#2"0BRb?Gt^## .llQT $v,,m㵜5ubV =sY+@d{N! dnO<.-B;_wJt6;QJd.Qc%p{ 1,sNDdFHI0ГoXшe黅XۢF:)[FGXƹ/w_cMeD,ʡcc.WDtA$j@:) -# u c1<@ۗ9F)KJ-hpP]_x[qBlbpʖw q"LFGdƶ*s+ډ_Zc"?%t[IP 6J]#=ɺVvvCGsGh1 >)6|ey?Lӣm,4GWUi`]uJVoVDG< SB6ϏQ@ TiUlyOU0kfV~~}SZ@*WUUi##; s/[=!7}"WN]'(L! ~y5g9T̅JkbM' +s:S +B)v@Mj e Cf jE 0Y\QnzG1д~Wo{T9?`Rmyhsy3!HAD]mc1~2LSu7xT;j$`}4->L#vzŏILS ֭T{rjGKC;bpU=-`BsK.SFw4Mq]ZdHS0)tLg