火烧赤壁

Editing: primes.awk

#!/usr/bin/mawk -f

# primes.awk
#
#  mawk -f primes.awk  [START]  STOP
# find all primes    between 2 and STOP
#       or START and STOP
#

function usage()
{ ustr = sprintf("usage: %s  [start] stop", ARGV[0])
  system( "echo " ustr) 
  exit 1
}

BEGIN { if (ARGC == 1 || ARGC > 3 ) usage()
        if ( ARGC == 2 )  { start = 2  ; stop = ARGV[1]+0 }
	else
	if ( ARGC == 3 )  { start = ARGV[1]+0 ; stop = ARGV[2]+0 }

if ( start < 2 ) start = 2
   if ( stop < start ) stop = start

prime[ p_cnt = 1 ] =  3  # keep primes in prime[]

# keep track of integer part of square root by adding
# odd integers 
   odd = test = 5
   root = 2
   squares = 9

while ( test <= stop )
{
   if ( test >= squares )
   { root++
     odd += 2
     squares += odd 
   }

flag = 1
   for ( i = 1 ; prime[i] <= root ; i++ )
   	if ( test % prime[i] == 0 )  #  not prime
	{ flag = 0 ; break }

if ( flag )  prime[ ++p_cnt ] = test

test += 2
}

prime[0] = 2

for( i = 0 ; prime[i] < start ; i++)  ;

for (  ;  i <= p_cnt ; i++ )  print prime[i]

}

Back