OFFSET
1,2
COMMENTS
If n is in the sequence and m < n is squarefree and coprime to n, then m*n is in the sequence. - Robert Israel, Aug 12 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local t;
andmap(t -> t[2]=1 and t[1]^2<=n, ifactors(n)[2])
end proc:
select(filter, [$1..1000]); # Robert Israel, Aug 12 2019
MATHEMATICA
selQ[n_] := SquareFreeQ[n] && AllTrue[FactorInteger[n][[All, 1]], #^2<=n&];
Select[Range[1000], selQ] (* Jean-François Alcover, Jul 30 2020 *)
PROG
(PARI) isok(n) = (n==1) || (issquarefree(n) && (vecmax(factor(n)[, 1])^2 < n)); \\ Michel Marcus, Aug 12 2019
KEYWORD
nonn
AUTHOR
Max Alekseyev, Nov 14 2017
STATUS
approved
