2

这是我到目前为止所拥有的：

//project Eular Problem 4: Prime Factors
#include<iostream>
#include<cmath>
typedef unsigned long long int uint_64;
using namespace std;

void storeFactors(uint_64 factors[], uint_64 num)
{
    for(uint_64 i=0;i*i<num;i++
       factors[i]=1;       //assign 1 to all the values

    for(uint_64 j=2; j*j<num;j++){
       if(num%j!=0)
          factors[j]=0;       //assign 0 to non-factors
    } 
}  

//Sieve of Eratosthenes to generate primes
void gen_primes(uint_64 arr[],uint_64 firstElement, uint_64 lastElement, uint_64 size)
{
    for(uint_64 i=0;i<size;i++)   //assigning 1 to all the values
       arr[i]=1;

    for(uint_64 i=2;i*i<=lastElement;i++){   //loop until the square-root of n
       if(arr[i])
          for(uint_64 j=i;j*i<=lastElement;j++)      //eliminate multiples by assigning them 0
            arr[j*i]=0;
       if(firstElement==1)
          arr[firstElement]=0;
    }
}

void arrayComp(uint_64 factors[],uint_64 primeArray[], uint_64 size)
{
    for(uint_64 i=2; i<=size; i++){
        if(factors[i] && primeArray[i]){
                cout<<i<<endl;
        }
    }
}

void processFactors(uint_64 num)
{
    uint_64 size = sqrt(num);
    uint_64 *factors = new uint_64[size];
    uint_64 *primeArray = new uint_64[size];

    storeFactors(factors, num);
    gen_primes(primeArray, 2, num, size);
    arrayComp(factors, primeArray,size);

    delete [] factors;
    delete [] primeArray;
}

int main()
{
    uint_64 number;
    cout<<"Enter a number: "<<endl;
    cin>>number;
    cout<<"The prime factors of "<<number<<" are: "<<endl;
    processFactors(number);

    return 0;
}

我尝试采用筛法来计算因子。所有非因子都分配为 0。如果它们既是输入因子又是素数，则 ArrayComp 将显示这些数字。

我遇到的问题是输出不完整并且程序遇到分段错误。例如，10 的因数是 5 和 2，但它显示 100 的答案相同。

编辑：我不太确定的另一件事是数组的大小。该程序显示 3 是 21 而不是 7 的主要因数，但如果我将大小增加 1，它会显示 3 和 5（不正确）

4

2 回答 2

7

100的质因数与10的质因数相同，只是重数不同。

看来您的表示没有任何方法来存储重复的因素。

于 2011-10-17T15:58:05.100 回答
1 
uint_64 size = sqrt(num);
uint_64 *factors = new uint_64[size];
uint_64 *primeArray = new uint_64[size];

错误的储备大小。

前任。）

数 = 10

大小 = 3

10 的因数 = 2, 5

    if(factors[i] && primeArray[i]){
            cout<<i<<endl;
    }

因素[5] ? 那么数组大小为 3(0,1,2)

于 2011-10-17T23:21:47.237 回答