OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
3660 is in the sequence because it has 4 distinct prime factors (2, 3, 5 and 61) and 61 > sqrt(3660).
MAPLE
with(numtheory): a:=proc(n) if nops(factorset(n))=4 and factorset(n)[4]^2>=n then n else fi end: seq(a(n), n=1..4500);
MATHEMATICA
pf4Q[n_]:=Module[{f=FactorInteger[n]}, Length[f]==4 && f[[-1, 1]] >= Sqrt[ n]]; Select[Range[5000], pf4Q] (* Harvey P. Dale, Sep 13 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Emeric Deutsch, Feb 02 2006
STATUS
approved