OFFSET
1,2
COMMENTS
First differs from A339741 in having 1080.
First differs from A382075 in having 18000.
Also numbers whose prime indices can be partitioned into distinct sets.
Differs from A212167, which does not include 18000 = 2^4*3^2*5^3, for example. - R. J. Mathar, Mar 23 2025
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The prime indices of 1080 are {1,1,1,2,2,2,3}, and {{1},{2},{1,2},{1,2,3}} is a partition into a set of sets, so 1080 is in the sequence.
We have 18000 = 2*5*6*10*30, so 18000 is in the sequence.
MAPLE
N:= 1000: # to get all terms <= N
A:= Vector(N):
A[1]:= 1:
for n from 2 to N do
if numtheory:-issqrfree(n) then
S:= [$1..N/n]; T:= n*S; A[T]:= A[T]+A[S]
fi;
od:
remove(t -> A[t]=0, [$1..N]); # Robert Israel, Apr 21 2025
MATHEMATICA
sqfacs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[sqfacs[n/d], Min@@#>d&]], {d, Select[Rest[Divisors[n]], SquareFreeQ]}]];
Select[Range[100], Length[sqfacs[#]]>0&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 21 2025
STATUS
approved
