|
%I #20 Jun 28 2017 20:37:28
%S 1,2,4,5,8,10,12,14,16,21,22,25,27,29,33,35,38,41,45,47,50,52,55,58,
%T 62,64,67,70,73,76,82,84,87,88,93,96,100,104,107,111,113,115,120,121,
%U 124,126,133,137,141,142,143,147,149,155,158,162,168,169,174,176,178,183
%N First appearance of n-th prime as prime factor in list of semiprimes.
%C The positions of even semiprimes in A001358. - _Juri-Stepan Gerasimov_, Apr 11 2010
%H Zak Seidov, <a href="/A115392/b115392.txt">Table of n, a(n) for n = 1..1000</a>
%F A001358(a(n)) = 2*A000040(n) - _Zak Seidov_, Jun 27 2017
%e a(5)=8 because 5th prime, 11, first appeared as a prime factor in 8th semiprime 22=2*11;
%e a(6)=10 because 6th prime, 13, first appeared as a prime factor in 10th semiprime 26=2*13.
%t sp = Select[Range[4,20000],2 == PrimeOmega[#]&]; Table[Position[Mod[sp, Prime[k]],0][[1,1]],{k,1000}] (* For first 1000 terms. - _Zak Seidov_, Jun 27 2017 *)
%o (PARI) first(n) = my(l = List([4,1]), u = 2*prime(n), res = vector(n), t=0); forprime(p = 2, t++; sqrt(2*prime(n)), forprime(q = p+1, u\p, listput(l, [p*q, t])); listsort(l); for(i=1, #l, if(res[l[i][3]]==0, res[l[i][3]] = i)) \\ _David A. Corneth_, Jun 28 2017
%Y Cf. A000040, A001358.
%K nonn
%O 1,2
%A _Zak Seidov_, Mar 08 2006
%E Edited by _Zak Seidov_, Jun 27 2017
|
