%I #13 Sep 08 2022 08:45:03
%S 23,29,35,39,44,47,53,55,59,62,63,69,71,74,76,79,80,83,87,89,90,95,97,
%T 98,99,103,104,107,109,111,116,118,119,124,125,127,129,131,132,134,
%U 135,139,142,143,146,149,151,152,153,155,159,160,161,164,167,168,169,170
%N Numbers k such that k! is divisible by (k+1)^6.
%H Robert Israel, <a href="/A061754/b061754.txt">Table of n, a(n) for n = 1..10000</a>
%p filter:= proc(n) local q,t,p,i,w,F;
%p F:= ifactors(n+1)[2];
%p for q in F do
%p p:= q[1];
%p t:= 0:
%p for i from 1 do
%p w:= floor(n/p^i);
%p if w = 0 then return false fi;
%p t:= t+w;
%p if t >= 6*q[2] then break fi;
%p od;
%p od;
%p true
%p end proc:
%p select(filter, [$1..200]); # _Robert Israel_, Jul 01 2018
%t Select[Range[200], IntegerQ[ #!/(# + 1)^6] &]
%o (PARI) isok(n) = !(n! % (k+1)^6); \\ _Michel Marcus_, Jul 02 2018
%o (Magma) [n: n in [1..200]| Factorial(n) mod (n+1)^6 eq 0]; // _Vincenzo Librandi_, Jul 02 2018
%K nonn
%O 1,1
%A _Robert G. Wilson v_, Jun 21 2001