|
%I #24 Oct 20 2023 08:41:08
%S 5,12,18,42,44,50,58,67,74,83,90,95,126,142,194,210,211,275,367,382,
%T 479,655,734,818,839,926,939,1047
%N Indices of semiprimes in A000073.
%C a(29) >= 1091. A000073(1091) is a 288-digit composite number with unknown factorization. Other possible terms are 1390, 1630, 2034, 2131, 2339, ... - _Tyler Busby_, Feb 09 2023
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Semiprime.html">Semiprime</a>
%F A000073(a(n)) = A101758(n).
%t Flatten[Position[LinearRecurrence[{1,1,1},{1,1,2},100],_?(PrimeOmega[ #]==2&)]]+1 (* This generates the first 12 terms of the sequence. To generate more, increase 100 to a higher number, but it will take increasingly long to run. *) (* _Harvey P. Dale_, May 31 2013 *)
%Y Cf. A000073, A001358, A101758.
%Y Equivalent sequence for primes: A303263.
%K nonn,more,less
%O 1,1
%A _Jonathan Vos Post_ and _Ray Chandler_, Dec 26 2004
%E More terms from _David Wasserman_, Apr 03 2008
|
