29 lines
681 B
Plaintext
29 lines
681 B
Plaintext
|
C Eratosthenes Sieve from BYTE magazine
|
||
|
|
program sieve
|
||
|
|
logical flags( 8191 )
|
||
|
|
integer*2 i, prime, k, count
|
||
|
|
integer*2 iter
|
||
|
|
|
||
|
|
write( *, 50 )
|
||
|
|
50 format( ' 10 iterations' )
|
||
|
|
do 92 iter = 1, 10
|
||
|
|
count = 0
|
||
|
|
do 10 i = 0, 8190
|
||
|
|
10 flags( i ) = .true.
|
||
|
|
do 91 i = 0, 8190
|
||
|
|
if ( .not. flags( i ) ) go to 91
|
||
|
|
prime = i + i + 3
|
||
|
|
k = i + prime
|
||
|
|
20 if ( k .gt. 8190 ) go to 90
|
||
|
|
flags( k ) = .false.
|
||
|
|
k = k + prime
|
||
|
|
go to 20
|
||
|
|
90 count = count + 1
|
||
|
|
91 continue
|
||
|
|
92 continue
|
||
|
|
write( *, 200 ) count
|
||
|
|
200 format( 1X, I6, ' primes' )
|
||
|
|
stop
|
||
|
|
100 format( 1X, I6 )
|
||
|
|
end
|