OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local q, t, p, i, w, F;
F:= ifactors(n+1)[2];
for q in F do
p:= q[1];
t:= 0:
for i from 1 do
w:= floor(n/p^i);
if w = 0 then return false fi;
t:= t+w;
if t >= 6*q[2] then break fi;
od;
od;
true
end proc:
select(filter, [$1..200]); # Robert Israel, Jul 01 2018
MATHEMATICA
Select[Range[200], IntegerQ[ #!/(# + 1)^6] &]
PROG
(PARI) isok(n) = !(n! % (k+1)^6); \\ Michel Marcus, Jul 02 2018
(Magma) [n: n in [1..200]| Factorial(n) mod (n+1)^6 eq 0]; // Vincenzo Librandi, Jul 02 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jun 21 2001
STATUS
approved