OFFSET
1,3
MATHEMATICA
Rest[Select[Table[k! k/(Times@@Divisors[k])^2, {k, 30}], IntegerQ]] (* Harvey P. Dale, Oct 17 2020 *)
Select[(# * #! / #^DivisorSigma[0, #]) & /@ Range[1, 30], IntegerQ] (* Amiram Eldar, Jul 16 2022 *)
PROG
(PARI) for(n=1, 30, if((n*(n!))%(n^numdiv(n))==0, print1((n*(n!))/(n^numdiv(n)), ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 08 2002
EXTENSIONS
Name clarified by Amiram Eldar, Jul 16 2022
Missing a(1)=1 inserted by Sean A. Irvine, Apr 06 2024
STATUS
approved