容易就是使用我的
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