The OEIS is supported by the many generous donors to the OEIS Foundation.
%I #15 Aug 02 2024 08:57:35
%S 16,32,48,64,80,81,96,112,128,144,160,162,176,192,208,224,240,243,256,
%T 272,288,304,320,324,336,352,368,384,400,405,416,432,448,464,480,486,
%U 496,512,528,544,560,567,576,592,608,624,625,640,648,656,672,688,704,720
%N Numbers whose prime factorization has exactly one exponent that is larger than 3.
%C Subsequence of A046101 and first differs from it at n = 98: A046101(98) = 1296 = 2^4 * 3^4 is not a term of this sequence.
%C Numbers k such that the powerful part of k, A057521(k), is a prime power whose exponent is larger than 3 (A246550).
%C The asymptotic density of this sequence is (1/zeta(4)) * Sum_{p prime} 1/(p^4-1) = 0.075131780079404733755... .
%H Amiram Eldar, <a href="/A375146/b375146.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>.
%e 16 = 2^4 is a term since its prime factorization has exactly one exponent, 4, that is larger than 3.
%t q[n_] := Count[FactorInteger[n][[;;, 2]], _?(# > 3 &)] == 1; Select[Range[1000], q]
%o (PARI) is(k) = #select(x -> x > 3, factor(k)[, 2]) == 1;
%Y Subsequence of A046101.
%Y Cf. A013662, A057521, A190641, A246550, A375145.
%K nonn,easy
%O 1,1
%A _Amiram Eldar_, Aug 01 2024