http://kenta.blogspot.com/2013/06/bgvtnpai-factorization-of-brents-c158.html The final complete answer is: ? factorint(p252-1) %19 = [2 13] [3 1] [13 1] [23 1] [29 1] [6329 1] [760347109 1] [211898520832851652018708913943317 1] [9409853205696664168149671432955079744397 1] [3035625952640765962086368602839492329385321995258480540549 1] [21014737139156086365509400019795978246705779843330697799891230652189207470396514322958290676782759071 1]