Question

我想要（these）大数，并且Matlab完美地做到了这一点。（These是结果。）

因子（NumberOfTypeSym）的结果是另一个符号对象：

>> factor( sym('79228162514264337589248983040') )

ans =

2^32*3*5*17*257*641*65537*6700417

但除此之外我无法做任何事情。

有没有办法访问单个素数和指数？

（这是我无法理解的，为什么他们不会只给我一个nx2矩阵，左边的素数和右边的指数。）

目前我特别需要的是结果中素数的数量，即所有指数的总和。

Answer 1

这个怎么样：

>> x=factor( sym('79228162514264337589248983040') )

x =

2^32*3*5*17*257*641*65537*6700417

>> p=char(x)

p =

2^32*3*5*17*257*641*65537*6700417

>> s=regexp(p,'*','split')

s = 

    '2^32'    '3'    '5'    '17'    '257'    '641'    '65537'    '6700417'

>> exp=regexp(s{1},'\^','split')

exp = 

    '2'    '32'
>> [exp s(2:end)]

ans = 

    '2'    '32'    '3'    '5'    '17'    '257'    '641'    '65537'    '6700417'

现在使用角色，根据需要将它们转换为数字。

Answer 2

很容易使用我在文件交换中找到的vpi工具。这将返回一个直接的整数列表。

x=factor(vpi('79228162514264337589248983040'))
x =
Columns 1 through 7
          2     2     2     2     2     2     2
Columns 8 through 14
          2     2     2     2     2     2     2
Columns 15 through 21
          2     2     2     2     2     2     2
Columns 22 through 28
          2     2     2     2     2     2     2
Columns 29 through 35
          2     2     2     2     3     5    17
Columns 36 through 39
        257   641   65537   6700417

vpi的替换（基本上是准备好发布）是vpij，它花了一半的时间来计算这个数字，但是vpij可以处理比vpi大得多的数字。事实上，最后要写的是一个更加改进的因子版本，但我可能会在完成之前释放vpij。

当然，因子必须也是可变精度整数，因为其中一些因素会非常大。这里有一个有50多位数的数字，有人最近在数学网站上发表了评论。

N = vpij(84)^27 + 1
N =
    9026943488964407632833018690186861978797224381906945

x = factor(N)
x =
   5   17   19  109  367  757  2017  230077  397741265470599434164843152148837

Answer 3

只是P0W提供的答案的延伸：

function [ y ] = symfactor2mat( x )

    Str = char(x) ;

    Parts = regexp(Str,'*','split') ;

    Long = length(Parts) ;

    Mat = sym(ones(Long,2)) ;

    for m=1:Long
        Part = Parts{m} ;
        if isempty(strfind( Part, '^' ))
            Mat(m,1) = sym(Part) ;
        else
            PE = regexp(Part,'\^','split') ;  % prime and exponent
            Mat(m,1) = sym(PE{1}) ;
            Mat(m,2) = sym(PE{2}) ;
        end
    end

    y = Mat ;

end

Answer 4

只是添加另一个选项：

N = sym('79228162514264337589248983040');
F = children(factor(N));
primefactors = F(2:2:end)

  primefactors =

  [ 2, 3, 5, 17, 257, 641, 65537, 6700417]

multiplicities = F(3:2:end)

  multiplicities =

  [ 32, 1, 1, 1, 1, 1, 1, 1]

sum(multiplicities)

  ans =

  39

或者，使用MuPAD documentation,

中的功能

feval(symengine, 'numlib::Omega', N)

  ans =

  39

在Matlab中考虑大数，然后得到结果中的数字

4 个答案: