Как бы вы проверили массив на правильность?

У меня есть массив простых чисел. Длина этого массива составляет 1 000 000. Это означает, что мои последние простые числа элементов [999999] должны иметь 1-миллионное простое число, и на самом деле это должно быть 1 000 001-е простое число с тех пор, как я оставил 2 с начала.

Когда я запускаю свою программу, и она выплевывает последний элемент массива, она выплевывает 999 999-е простое число. Я нахожусь в процессе изучения C, и это небольшой домашний проект для математического класса, который я беру. Я понятия не имею, с чего начать в устранении неполадок этой проблемы. Любая помощь будет принята с благодарностью.

Изменить: Извините, по некоторым причинам я думал, что это может быть простой ответ. Я разместил код ниже. Я распечатал массив простых чисел в варианте 1, и все простые числа верны, а результат правильный. Любой другой вариант, и я получаю ошибочный вывод.

void createArray( int, unsigned long long * );
unsigned long long root(unsigned long long);
void print( unsigned long long *, int );

int main(void){

    unsigned long long *prime = (unsigned long long *)calloc(SIZE6, 
                                    sizeof(unsigned long long));
    unsigned long long *ptrEmpty = &prime[15];
    unsigned long long i;
    unsigned long long j;
    long int sizeComp;
    int flag = 0;
    unsigned long option = 0;
    int length = 0;

    system( "clear" );

    do{ 
        flag = 0;
        option = -1;

        printf( "Welcome to the prime number finder program. Below are a list "
            "of prime numbers\nyou can find as well as their average runni"
            "ng time on the Loki system of UNO.\n\nOption 1: Return the 10"
            "0th prime number. ( Running time: .006s )\nOption 2: Return t"
            "he 1000th prime number. ( Running time: .008s )\nOption 3: Re"
            "turn the 10000th prime number. ( Running time: .089s )\nOptio"
            "n 4: Return the 100000th prime number. ( Running time: 2.188s"
            " )\nOption 5: Return the 1000000th prime number. ( Running ti"
            "me: 57.156s )\nOption 6: Return the 10000000th prime number. "
            "( **CAUTION!** Running time: 35m51.3s )\n\nEnter the number o"
            "f the option you would like the program to\nperform (Enter 0 "
            "to exit):  " );

        scanf( "%lu", &option );

        switch( option ){

            case 0:
                flag = 1;
                break;
            case 1:
                createArray( SIZE1, prime );
                ptrEmpty = &prime[15];
                length = SIZE1;
                break;
            case 2:
                createArray( SIZE2, prime );
                ptrEmpty = &prime[15];
                length = SIZE2;
                break;
            case 3:
                createArray( SIZE3, prime );
                ptrEmpty = &prime[15];
                length = SIZE3;
                break;
            case 4:
                createArray( SIZE4, prime );
                ptrEmpty = &prime[15];
                length = SIZE4;
                break;
            case 5:
                createArray( SIZE5, prime );
                ptrEmpty = &prime[15];
                length = SIZE5;
                break;
            case 6:
                createArray( SIZE6, prime );
                ptrEmpty = &prime[15];
                length = SIZE6;
                break;
            default:
                printf( "Please enter one of the available options ( 1 - 7 ) "
                "or 0 to exit:  " );
            flag= 1;
            break;
        }

        if( flag != 1 ){
            for( i = START, sizeComp = 15; sizeComp <= length; i += 2 ){

                flag = 0;
                for( j = 0; prime[j] < root( i ); j++ ){

                    if( ( i % prime[j] ) == 0 ){
                        flag = 1;
                        break;
                    }
                }

                if( flag == 0 ){
                    *ptrEmpty = i;
                    ++sizeComp;
                    ptrEmpty++;
                }
            }

            printf( "\nThe %dth prime number is %llu.\n\n", length, 
                    prime[ length - 2 ]);

//            print( prime, length );
            printf( "****************************************"
                    "****************************************\n\n" );

    }

    }while( option != 0 );

    free( prime );

    return 0;
}

void print( unsigned long long *prime, int length ){

    int i;
    printf( "\n%4d", 2 );
    for( i = 0; i < length; i++){
        printf( "%4llu ", prime[i] );
    }
    printf( "\n" );
}

void createArray( int length, unsigned long long *prime ){

    unsigned long long *pHolder = NULL;

    pHolder = (unsigned long long *)realloc( prime, length * 
                                        sizeof( unsigned long long ) );

    if(pHolder != NULL){
        prime = pHolder;
    }else{
        free( prime );
        printf( "Error reallocating memory.\n" );
        exit(1);
    }

    prime[0] = 3;
    prime[1] = 5;
    prime[2] = 7;
    prime[3] = 11;
    prime[4] = 13;
    prime[5] = 17;
    prime[6] = 19;
    prime[7] = 23;
    prime[8] = 29;
    prime[9] = 31;
    prime[10] = 37;
    prime[11] = 41;
    prime[12] = 43;
    prime[13] = 47;
    prime[14] = 53;
}

unsigned long long root(unsigned long long a) {

    unsigned long long rem = 0;
    unsigned long long root = 0;
    int i;

    for (i = 0; i < 16; i++) {
        root <<= 1;
        rem <<= 2;
        rem += a >> 30;
        a <<= 2;

        if (root < rem) {
            root++;
            rem -= root;
            root++;
        }
    }

    return (root >> 1);
}

Изменить: РАЗМЕР 1 = 100, РАЗМЕР 2 = 1000, РАЗМЕР 3 = 10000, РАЗМЕР 4 = 100000, РАЗМЕР 5 = 1000000, РАЗМЕР 6 = 10000000

1 ответ

Вы можете заняться отладкой, изменив размер массива на 5 и просмотрев только первые 5 простых чисел, распечатав все 5 из них, и посмотрите, действительно ли это 2 3 5 7 11 (или в вашем случае 3 5 7 11 13) или если у вас есть ошибки там.

Ваша ошибка, вероятно, будет существовать там, и вы можете исправить ее для гораздо меньшего случая, и она должна работать для большего случая.

Другие вопросы по тегам