git.sur5r.net Git - cc65/blob - samples/sieve.c

   1 /*
   2  * Calculate all primes up to a specific number.
   3  */
   4
   5
   6
   7 #include <stdio.h>
   8 #include <conio.h>
   9
  10
  11
  12 /*****************************************************************************/
  13 /*                                   Data                                    */
  14 /*****************************************************************************/
  15
  16
  17
  18 #define COUNT           8192            /* Up to what number? */
  19 #define SQRT_COUNT      91              /* Sqrt of COUNT */
  20
  21 static unsigned char Sieve[COUNT];
  22
  23
  24
  25 /*****************************************************************************/
  26 /*                                   Code                                    */
  27 /*****************************************************************************/
  28
  29
  30
  31 int main (void)
  32 {
  33     /* This is an example where register variables make sense */
  34     register unsigned char* S;
  35     register unsigned       I;
  36     register unsigned       J;
  37
  38     /* Execute the sieve */
  39     I = 2;
  40     while (I < SQRT_COUNT) {
  41         if (Sieve[I] == 0) {
  42             /* Prime number - mark multiples */
  43             S = &Sieve[J = I*2];
  44             while (J < COUNT) {
  45                 *S = 1;
  46                 S += I;
  47                 J += I;
  48             }
  49         }
  50         ++I;
  51     }
  52
  53     /* Print the result */
  54     for (I = 2; I < COUNT; ++I) {
  55         if (Sieve[I] == 0) {
  56             printf ("%4d\n", I);
  57         }
  58         if (kbhit() && cgetc() == 'q') {
  59             break;
  60         }
  61     }
  62
  63     return 0;
  64 }
  65
  66
  67