Return to primes.p CVS log

Up to [CSRG BSD Unix] / 43BSDTahoe / ucb / pascal / tstpx / src

Annotation of 43BSDTahoe/ucb/pascal/tstpx/src/primes.p, revision 1.1.1.1

1.1       root        1: program primes(output);
                      2: const n = 1000; n1 = 33;   (*n1 = sqrt(n)*)
                      3: var i,k,x,inc,lim,square,l: integer;
                      4:     prim: boolean;
                      5:     p,v: array[1..n1] of integer;
                      6: begin writeln(wallclock);
                      7:    write(2:6, 3:6); l := 2;
                      8:    x := 1; inc := 4; lim := 1; square := 9;
                      9:    for i := 3 to n do
                     10:    begin (*find next prime*)
                     11:       repeat x := x + inc; inc := 6-inc;
                     12:          if square <= x then
                     13:             begin lim := lim+1;
                     14:                v[lim] := square; square := sqr(p[lim+1])
                     15:             end ;
                     16:          k := 2; prim := true;
                     17:          while prim and (k<lim) do
                     18:          begin k := k+1;
                     19:             if v[k] < x then v[k] := v[k] + 2*p[k];
                     20:             prim := x <> v[k]
                     21:          end
                     22:       until prim;
                     23:       if i <= n1 then p[i] := x;
                     24:       write(x:6); l := l+1;
                     25:       if l = 10 then
                     26:          begin writeln; l := 0
                     27:          end
                     28:    end ;
                     29:    writeln; writeln(wallclock)
                     30: end .

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.