OFFSET
1,1
COMMENTS
Any term of A391088 is of the form k*m where k is a term of this sequence and m is a squarefree number coprime to k. Therefore, A391088 can be generated from this sequence by multiplying terms with coprime squarefree numbers, and the asymptotic density of A391088 can be evaluated from the terms of this sequence (see the Comments section of A391088).
The least odd term is a(1458) = A321147(1) = 225450225 = (3 * 5 * 7 * 11 * 13)^2, which is the least odd exponential abundant number.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10127 (terms below 10^10)
FORMULA
36 = 2^2 * 3^2 is a term since it is a powerful number, and its exponential unitary divisors, {6, 12, 18, 36}, can be partitioned into 2 disjoint sets whose sum is equal: 6 + 12 + 18 = 36.
MATHEMATICA
pows[max_] := Union[Flatten[Table[i^2*j^3, {j, 1, Surd[max, 3]}, {i, 1, Sqrt[max/j^3]}]]];
seq[max_] := Select[pows[max], euzQ]; seq[200000] (* using the function "euzQ" from A391088 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 28 2025
STATUS
approved
