|
%I #14 Dec 02 2020 03:19:27
%S 16,48,64,80,81,112,144,162,176,192,208,240,256,272,304,320,324,336,
%T 368,400,405,432,448,464,496,512,528,560,567,576,592,624,625,648,656,
%U 688,704,720,729,752,768,784,810,816,832,848,880,891,912,944,960,976,1008
%N Numbers whose prime factorization contains at least one composite exponent.
%C 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.05328066264472198953... (using the method of Shevelev, 2016). - _Amiram Eldar_, Nov 08 2020
%H Alois P. Heinz, <a href="/A322448/b322448.txt">Table of n, a(n) for n = 1..10000</a>
%H Vladimir Shevelev, <a href="https://arxiv.org/abs/1602.04244">A fast computation of density of exponentially S-numbers</a>, arXiv:1602.04244 [math.NT], 2016.
%e 16 = 2^4 is a term because 4 is a composite exponent here.
%p a:= proc(n) option remember; local k; for k from 1+
%p `if`(n=1, 0, a(n-1)) while andmap(i-> i[2]=1 or
%p isprime(i[2]), ifactors(k)[2]) do od; k
%p end:
%p seq(a(n), n=1..70);
%t Select[Range[1000], AnyTrue[FactorInteger[#][[;; , 2]], CompositeQ] &] (* _Amiram Eldar_, Nov 08 2020 *)
%o (PARI) isok(m) = #select(x->((x>1) && !isprime(x)), factor(m)[,2]) > 0; \\ _Michel Marcus_, Dec 02 2020
%Y Cf. A002808, A322449.
%K nonn
%O 1,1
%A _Alois P. Heinz_, Dec 08 2018
|
