使用埃拉托色尼筛选法(the Sieve of Eratosthenes)在一定范围内求素数及反素数(Emirp)

Programming 1.3 In this problem, you'll be asked to find all the prime numbers from 1 to 1000. Prime numbers are used in all
kinds of circumstances, particularly in fields such as cryptography, hashing among many others. Any method w ill be sufficient for
this problem. The Sieve of Eratosthenes is one algorithm which you can try implementing, but there are plenty of others. A prime is
prime if it is not the product of any lesser natural numbers except for 1 and itself, for example, 1, 2 and 3 should be prime by this
definition, but 4, being the product of 2 and 2, is not.

Programming Challenge 1.1 Prime numbers are fairly useful on their own, but another fun challenge is to find what are known
as emirp primes. Emirp primes are primes which are also primes when reversed. For this challenge, you need to print out the list
of all primes that are also primes (as defined in the challenge above) when reversed. [2]

------------------------------------------------------------------------------------------------------------

实际上素数（prime number）是不包括1的，所以就把1排除掉了。

埃拉托色尼筛选法(the Sieve of Eratosthenes)用于筛选自然数数列中的素数。

反素数(Emirp):一个数是素数且将他反过来也是素数，如13和31。

/* Prime numbers are fairly useful on their ow n, 
but another fun challenge is to find w hat are know n
as emirp primes. Emirp primes are primes w hich are also 
primes w hen reversed. For this challenge, you need to 
print out the list of all primes that are also primes 
(as defined in the challenge above) w hen reversed from 1 - 1000.
*/
#include <stdio.h>
#include <stdlib.h>

#define PRIMECOUNT 1000

void sieveOfEratosthenes(int *n, int number);
void emirp(int *n, int number);
int reverse(int n);

int main(int argc, char **argv){
    /* The primes 1 - 1000, if n is prime,
    primes[n - 1] == 1 and 0 otherwise. */
    int primes[PRIMECOUNT];
    int i;
    for(i = 0; i < PRIMECOUNT; i++){
        primes[i] = 1;
    }
    /* Write a prime checking algorithm here. */
    sieveOfEratosthenes(primes, PRIMECOUNT);
    //printf("%d\n", primes[6]);
    
    emirp(primes, PRIMECOUNT);
    //printf("%d\n", primes[6]);
    /* -------------------------------------- */
    printf("All primes found from 1 - 1000:\n");
    for(i = 0; i < PRIMECOUNT; i++){
        if(primes[i] == 1){
        printf("%d\t",i + 1);
        }
    }
    printf("\n");
    return 0;
}

void sieveOfEratosthenes(int *n, int number) {
    int i = 0, j = 0, k = 2;
    n[0] = 0;
    for (i = 1; i < number; i++) {
        if (n[i] == 0) {
            continue;
        }
        for (j = i, k = 2; (j + 1) * k <= number; k++) {
            n[(j + 1) * k - 1] = 0;
        }
    }
}

void emirp(int *n, int number) {
    int i = 0;
    for (i = 0; i < number; i++) {
        if (n[i] == 1 && n[reverse(i + 1) - 1] == 0) {
            //printf("%d is not!\n", i+1);
            n[i] = 0;
            //n[reverse(i + 1) - 1] = 0;
        }
    }
}

int reverse(int n) {
    int reverseNumber = 0, temp = 0;
    while (n != 0) {
        temp = n % 10;
        n = n / 10;
        reverseNumber = reverseNumber * 10 + temp;
    }
    return reverseNumber;
}

所犯的错误是55行的判断是否是反素数，原先是

void emirp(int *n, int number) {
    int i = 0;
    for (i = 0; i < number; i++) {
        if (!(n[i] == 1 && n[reverse(i + 1) - 1] == 1)) {
            n[i] = 0;
            n[reverse(i + 1) - 1] = 0;
        }
    }
}

这样会使得如2，11，31等被排除反素数，因为20，110，310不是素数。