Relay-Version: version B 2.10 5/3/83; site utzoo.UUCP
Posting-Version: version B 2.10.1exp 11/4/83; site ihuxl.UUCP
Path: utzoo!linus!decvax!harpo!floyd!clyde!ihnp4!ihuxl!dcn
From: dcn@ihuxl.UUCP (Dave Newkirk)
Newsgroups: net.sources
Subject: Aztec C source - primes.c
Message-ID: <741@ihuxl.UUCP>
Date: Mon, 14-Nov-83 09:31:48 EST
Article-I.D.: ihuxl.741
Posted: Mon Nov 14 09:31:48 1983
Date-Received: Tue, 15-Nov-83 06:37:22 EST
Organization: AT&T Bell Labs, Naperville, Il
Lines: 33


/* Eratosthenes Sieve Prime Number Program */

#define TRUE 1
#define FALSE 0
#define SIZE 8190 

char flags[SIZE+1] ;

main()
{
   int i, prime, k, count, iter ;

   printf("10 iterations\n") ;

   for( iter=1 ; iter <= 10 ; iter++ ) {
	  count = 0 ;
	  for( i=0 ; i <= SIZE ; i++ )
		 flags[i] = TRUE ;
	  for( i=0 ; i<= SIZE ; i++ ) {
		 if( flags[i] ) {
			prime = i + i + 3 ;
/*			printf("\n%d", prime ) ; */
			for( k=i+prime ; k <= SIZE ; k+=prime )
			   flags[k] = FALSE ;
			count++ ;
		 }
	  }
   }

   printf("\n%d primes.\n", count ) ;

} /* end main */