|
MATHEMATICA
|
aa = {}; Do[If[(Sqrt[n! - 1] /. Sqrt[_] -> 1) > 1, Print[n]; AppendTo[aa, n], {n, 1, 1000}]; aa (*Artur Jasinski*)
nfdsQ[n_]:=AnyTrue[Rest[Divisors[n!-1]], IntegerQ[Sqrt[#]]&]; Select[Range[ 2, 610], nfdsQ] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, May 24 2020 *)
|