OFFSET
1,2
COMMENTS
Subsequence of A209061, A386803 and A388972, and first differs from them at n = 62: A209061(62) = A386803(62) = A388972(62) = 64 = 2^6 is not a term in this sequence.
First differs from A115063 and A369939 at n = 122: a(122) = 128 = 2^7 is neither a term in A115063 nor a term in A369939. Also, A115063(243) = A369939(243) = 256 = 2^8 is the least term in A115063 and A369939 that is not a term in this sequence.
Except for a(1) = 1, a subsequence of A140823 and differs from it by not having the terms 48, 64, 80, 112, 144, ... .
The squarefree numbers (A005117) and numbers whose prime factorization exponents are all primes (A056166) are all terms in this sequence. Each term in this sequence has a unique representation as the product of two coprime numbers, one is squarefree and the other is a term in A056166.
The asymptotic density of this sequence is Product_{p prime} (1 - 1/p^4 + Sum_{q prime >= 5} (1/p^q - 1/p^(q+1))) = 0.94671933735527801046... . [corrected by Jason Yuen, Nov 11 2025]
Numbers for which their sets of exponential divisors (A322791), exponential unitary divisors (A361255), and exponential semiproper divisors (defined in A323308 and A323309) coincide. - Amiram Eldar, Dec 05 2025
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[n_] := AllTrue[FactorInteger[n][[;; , 2]], ! CompositeQ[#] &]; Select[Range[100], q]
PROG
(PARI) isok(k) = vecsum(apply(x -> if(x == 1 || isprime(x), 0, 1), factor(k)[, 2])) == 0;
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Nov 05 2025
STATUS
approved
