Mathematicaの素因数分解の関数FactorIntegerは、
100%正しく素因数分解出来ていると考えて良いので
しょうか？約20桁の整数の素因数分解を行っています。

ttp://reference.wolfram.com/mathematica/ref/FactorInteger.html

PrimeQの説明の所に、10^16以上の場合は誤動作の可能性が指摘
されているので、少々気になっています。

ttp://reference.wolfram.com/mathematica/note/SomeNotesOnInternalImplementation.html#6849