1 #include <stdio.h>
 2  #define MAX 1000
 3  #define null1 0
 4  #define NEXT(x)  x=next[x]
 5  #define REMOVE(x) {   previous[next[x]]=previous[x];   \
 6                        next[previous[x]]=next[x];       \
 7                    }
 8  
 9  #define INITIAL(n)  { unsigned long i;                    \
10                        for(i=2;i<=n;i++)                   \
11                            previous[i]=i-1,next[i]=i+1;    \
12                        previous[2]=next[n]=null1;           \
13                      }
14  
15  int main()
16  {
17      unsigned long previous[MAX+1]={0};
18      unsigned long next[MAX+1]={0};
19      unsigned long prime,fact,i,mult;
20      unsigned long n;
21      unsigned long count=0;
22      
23      scanf("%lu",&n);
24  
25      INITIAL(n); //initial the array
26  
27      for(prime=2;prime*prime<=n;NEXT(prime))
28      {
29          for(fact=prime;prime*fact<=n;NEXT(fact)) 
30          {
31              for(mult=prime*fact;mult<=n;mult*=prime) 
32                  REMOVE(mult);
33          }
34      }
35      for(i=2;i!=null1;NEXT(i))
36          printf("%lu ",i),count++;
37      printf("\nThe sum of the prime numbers is %lu\n",count);
38  }