OFFSET
1,1
COMMENTS
LINKS
EXAMPLE
MATHEMATICA
fct = Array[#! &, 10]; prev = {}; jp = fct; While[jp != prev, prev = jp; jp = Select[Union @@ Outer[Times, jp, fct], # <= fct[[-1]] &]]; fctdiv[n_] := Module[{m = 1, k = 1}, While[Divisible[n, m], k++; m *= k]; m /= k; n/m]; Select[jp, FixedPoint[fctdiv, #] != 1 &] (* Amiram Eldar, May 22 2021 *)
PROG
(PARI)
search_up_to = 2^22;
A076934(n) = for(k=2, oo , if(n%k, return(n), n /= k));
A001013list(lim, mx=lim)=if(lim<2, return([1])); my(v=[1], t=1); for(n=2, mx, t*=n; if(t>lim, break); v=concat(v, t*A001013list(lim\t, t))); Set(v) \\ From A001013
v001013 = A001013list(search_up_to);
A001013(n) = v001013[n];
isA344179(n) = if(v001013[#v001013]<n, -(1/0), ((1!=A093411(n))&&vecsearch(v001013, n)));
for(n=1, search_up_to, if(isA344179(n), print1(n, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 18 2021
STATUS
approved