c中的溢出错误以查找十二位数字的质因子

问题描述

我想通过以下程序解决此问题：https://projecteuler.net/problem=3，但是我不确定使用long long int的正确方法。

#include <stdio.h>
#include <stdlib.h>

int main(int argc,char *argv[]) {
    long long int result = factors(6051475143);
    printf("%lld",result); 
    return 0;
}

void factors(int number) {
    int factor[100000];
    int index = 0,i = 1;
    for (; i < number; i++) {
        if (number % i == 0) {
            factor[index] = i;
            index += 1;
        }
    }
    findprime(factor,100000);
}

int findprime(int prime[],int size) {
    int i = 0,j;
    long long int latestprime;
    for(; i < size; i++) {
        if (prime[i] == 0)
            break;
        int is_prime = 1;
        for (j = 2; j < prime[i]; j++) {
            if (prime[i] % j == 0 && prime[i] != j) {   
                is_prime = 0;
                break;
            }
        }
        if (is_prime == 1)
            latestprime = prime[i];
    } 
    return latestprime;
}

如果我尝试使用10位数字，那么它会起作用，但是当我尝试使用12位数字时，它将返回零。

解决方法

要解决此问题，您确实需要一个足以表示数字的类型。 unsigned long long被指定为至少具有64个值位，即18446744073709551615。

您的程序具有未定义的行为，即使是很小的数字：result = factors(6051475143)也不能可靠地存储任何有用的东西，因为factors被定义为void函数，并且确实甚至没有return语句。之所以会偶然，是因为最后一条语句findprime(factor,100000);将findprime的返回值留在了寄存器中，其中main检索了它期望的{{1 }}，这是编译器根据缺少原型推断出的返回类型。

要找到最大的素因数，只要找到一个将其均分的因子，就应尝试并减少其数目。

这是修改后的版本：

int

,

#include <stdio.h>
#include <stdlib.h>


int main() {

long long int result =factors(6051475143);
printf("%lld",result); 
return 0;
}

factors(long long int number)
{
    long int factor[1000000];
    long long int x;
    long int index=0,i=1;
    for(;i<800000; i++)
    {
        if (number%i==0)
        {
            factor[index] = i;
            index+=1;
        
        }
    }
    x = findprime(factor,index);
    return x;
}

int findprime(long int prime[],long int size)
{
    long int i = 0,j;
    long long int latestprime;
    for(;i<size;i++)
    {
        if (prime[i] == 0)
            break;
        int is_prime=1;
        for(j=2;j<prime[i];j++)
        {
            if (prime[i]%j==0 && prime[i]!=j)
            {   
                is_prime=0;
                break;
            }
        }
        if (is_prime==1)
            latestprime = prime[i];
    } 
    return latestprime;
}