gap> FactorsPminus1(Factorial(158)+1,100000,1000000);
#I  p-1 for n = 
185327186949373479654360975305107852968236093964410457931173183300920822903860\
684098654886099887970007689751613529715511834491898311282134018439421327637431\
255849363723893479936922146879010853472826978770662654016390805763283535114797\
45865318400000000000000000000000000000000000001
a : 2, Limit1 : 100000, Limit2 : 1000000
#I  First stage
#I  4 - digit factor 5227 was found
#I  4 - digit factor 2879 was found
#I  Second stage
#I  14 - digit factor 18331561438319 was found
#I  13 - digit factor 9561906969931 was found
#I  13 - digit factor 1452486383317 was found
[ [ 2879, 5227, 1452486383317, 9561906969931, 18331561438319, 
      483714299709483760811581110341732950506493218122654853400674921345082310\
906370452295654816571305041217323052879842924826121333143254713674832962773107\
806789945715570386038565256719614524924705165110048148716160964980629081176057\
0095669 ], [  ] ]
gap>