j'utilise le facteur (n) dans matlab mais cela ne fonctionne pas pour les nombres comme ((10^11-1)/11), pour le même nombre que puis-je faire?

Question

J'utilise factor(n) dans matlab mais cela ne fonctionne pas pour les nombres comme ((10^11-1)/11) - que puis-je faire? Ma source est ici.

m=input('Enter your number: '); 
i=0; 

while (i<m) 
    if(isprime(i)) 
    % sum=factor((10^(i-1)-1)); 
     sum=factor((10^(i-1)-1)/i); 
     disp(sum); 
    end 
    i =i+1; 

end 

mais pour un grand n, il renvoie des erreurs !!

>> FactorGen 

Enter your number: 45 

3  3 

3  3 11 

3  3 11 101 

3  3  3  7 11 13 37 

     3   3   11   41   271  9091 

     3   3   3   7   11   13   37   101  9901 

Error using factor (line 26) When n is single or double, its maximum 
allowed value is FLINTMAX. 

Error in FactorGen (line 7) sum=factor((10^(i-1)-1)); 

Je veux la fonction factor((10^(i-1)-1)) de travailler pour autant. Comment puis-je résoudre mon problème?

Answer 1

Je pense que cela peut être partiellement atténué en convertissant votre grand nombre en format uint64. Pour R2014b entier maximum qui peut être manipulé est:

n = intmax('uint64') 
n = 1.8447e+19 

Alors que le double maximum qui peut être manipulé est:

n = flintmax('double') 
n = 9.0072e+15 

Cela peut être vérifié par exemple simple. Utilisons factor sur le nombre plus grand que flintmax. Tout d'abord, essayez double:

factor(10^16) 
Error using factor (line 26) 
When n is single or double, its maximum allowed value is FLINTMAX. 

Maintenant, nous essayons uint64:

factor(uint64(10^16)) 
ans = 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 

Ainsi, en utilisant factor(uint64(x)) augmente votre limite de plusieurs ordres. En effet, cela n'a aucun sens d'utiliser le format double pour la fonction factor.

Par ailleurs, puisque l'utilisation factor sur un grand nombre peut geler votre système, je suggère d'utiliser cette fonction:

function f = Factorize(n) 
i = 0; 
while mod(n,2) == 0 
    i = i + 1; 
    n = n/2; 
    f(i) = 2; 
    disp([2,n]) 
end 
q = round(sqrt(double(n))); 
q = q + 1 - mod(q,2); 
for j = 3:2:q 
    while mod(n,j) == 0 
     i = i + 1; 
     f(i) = j; 
     n = n/j; 
    end 
end 
if n > 2; 
    i = i + 1; 
    f(i) = n; 
end 

Il est beaucoup plus rapide pour un grand nombre et ne surcharge pas le système dans son ensemble n espoir cela aide