OFFSET
1,1
COMMENTS
First differs from its subsequence A381315 at n = 40: a(40) = 432 = 2^4 * 3^3 is not a term of A381315.
Numbers k such that A295883(k) = 1.
The asymptotic density of this sequence is Product_{p prime} (1 - 1/p^3 + 1/p^4) * Sum_{p prime} (p-1)/(p^4 - p + 1) = 0.092831691827595439609... (Elma and Martin, 2024).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Ertan Elma and Greg Martin, Distribution of the number of prime factors with a given multiplicity, Canadian Mathematical Bulletin, Vol. 67, No. 4 (2024), pp. 1107-1122; arXiv preprint, arXiv:2406.04574 [math.NT], 2024.
MATHEMATICA
f[p_, e_] := If[e == 3, 1, 0]; s[1] = 0; s[n_] := Plus @@ f @@@ FactorInteger[n]; Select[Range[300], s[#] == 1 &]
PROG
(PARI) isok(k) = vecsum(apply(x -> if(x == 3, 1, 0), factor(k)[, 2])) == 1;
CROSSREFS
A381315 is subsequence.
Cf. A295883.
Numbers that have exactly one exponent in their prime factorization that is equal to k: A119251 (k=1), A386796 (k=2), this sequence (k=3), A386804 (k=4), A386808 (k=5).
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Aug 03 2025
STATUS
approved
