1

我试图通过尝试在 Project Eular 页面上做示例来了解更多关于 C 代码的信息：http ://projecteuler.net/problems 目前我试图从数字 600851475143 中计算出最高素数。这是我的代码：

/* The prime factors of 13195 are 5, 7, 13 and 29.

What is the largest prime factor of the number 600851475143 ?*/
#include <stdio.h>

unsigned long long GetPrimeFactor(unsigned long long number)
{
    /* printf("The number you typed was %d\n", number);*/
    unsigned long long prime = 2;
    unsigned long long LargestPrime = 0;

    while (prime != 0)
    {
        if(number%prime == 0)
        {
            if(number == prime)
            {
                if(LargestPrime < prime)
                    LargestPrime = prime;

                printf("End of prime factors \nLargest Prime: %d\nEnding Number: %d\n", LargestPrime, number);
                break;
            }
            else if(LargestPrime < prime)
                LargestPrime = prime;

            //divide number by prime

            number = number/prime;
            printf("Prime: %d\n", prime);
        }

        else
        {
            //get new prime
            prime = nextPrimeNumber(prime);
        }
    }
}

unsigned long long nextPrimeNumber(unsigned long long input)
{
    unsigned long long nextPrimeNumber;

    while(input !=0)
    {
        if( input < 0)
        {
            printf("The number you have entered is negative\n");
        }
        else if (input > 0)
        {
            nextPrimeNumber = input + 1;

            if(nextPrimeNumber%2 == 0 && nextPrimeNumber != 2)
            {
                nextPrimeNumber += 1;
            }

            while(!isPrime(nextPrimeNumber))
            {
                nextPrimeNumber += 2;
            }

            return nextPrimeNumber;
        }
    }
}

int isPrime(int number)
{
    int i;
    int prime = 1; //true

    if(number == 2)
        prime = 0; //false

    if(number%2 == 0 || number <= 1)
        prime = 0;
    else
    {
        for(i=3; i<sqrt(number) && prime == 1; i+=2)
            if(number%i == 0)
                prime = 0;
    }

    return prime;
}

int main(void)
{
    unsigned long long number;

    printf("Enter a Number: \n");
    scanf("%d", &number);
    GetPrimeFactor(number);

    return 0;
}

此代码适用于示例中的数字 13195，但是当我尝试更大的数字并且我不确定出了什么问题时会变坏。

输出大数 600851475143：素数：3 素数：29 素数：5108231 等。这不可能是正确的，因为该数字无论如何都不能被 3 整除。

无法弄清楚出了什么问题，所以任何帮助都会很棒

4

2 回答 2

2

签名ints最高 2147483647，unsigned ints最高 4294967295，您可以使用unsigned long long最高 18446744073709551615，如果您需要更大的东西，请使用任意精度/bignum 库，如 GnuMP（http://gmplib.org/）

更新您编辑的问题：

你使用的 scanf 的格式说明符甚至没有读完整个unsigned long long，我读到了你必须使用的地方%llu，但是在我的 gcc 上这不起作用，我必须使用：

printf("Enter a Number: \n");
scanf("%I64d", &number);
printf("%I64d", number);  // <- make sure that this prints the number you 
                 // gave it before even calling GetPrimeFactor()
于 2012-02-16T13:05:12.910 回答
0

你有没有改变scanf？

scanf("%L", &number)
于 2012-02-16T13:29:16.350 回答