matlab nextpow,大数的素性检测(用MATLAB仿真)

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function f = factor(n)

%FACTOR Prime factors.

% FACTOR(N) returns a vector containing the prime factors of N.%% This function uses the simple sieve approach. It may require large

% memory allocation if the number given is too big. Technically it is% possible to improve this algorithm, allocating less memory for most% cases and resulting in a faster execution time. However, it will still

% have problems in the worst case.

% % Class support for input N:

% float: double, single

% integer: uint8, int8, uint16, int16, uint32, int32, uint64, int64

%% See also PRIMES, ISPRIME.

% Copyright 1984-2012 The MathWorks, Inc.

% $Revision: 1.14.4.5 $ $Date: 2012/10/15 20:09:54 $

if ~isscalar(n)

rror(message('MATLAB:factor:NonScalarInput'));

end

if ~isreal(n) || (n < 0) || (floor(n) ~= n)

error(message('MATLAB:factor:InputNotPosInt'));

end

if (isfloat(n) && n > flintmax(class(n)))

error(message('MATLAB:factor:InputOutOfRange'));

end

if n < 4

f=n;

return

else

f = [];

end

if (isa(n,'uint64') || isa(n,'int64')) && n > flintmax

p = primes(2.^(nextpow2(n)/2));

else

p = primes(cast(sqrt(double(n)),class(n)));

end

while n>1,

d = find(rem(n,p)==0);

if isempty(d)

f = [f n];

break;

end

p = p(d);

f = [f p];

n = n/prod(p);

end

f = sort(f);