login
Difference between the prime indices of the two factors of the n-th semiprime.
34

%I #24 Dec 21 2020 07:22:07

%S 0,1,0,2,3,1,2,4,0,5,3,6,1,7,4,8,0,5,2,6,9,10,3,7,11,1,12,4,13,8,2,9,

%T 14,5,15,10,6,16,3,0,17,11,12,4,18,13,19,1,7,20,8,21,14,5,22,0,15,23,

%U 16,9,2,24,17,25,6,10,26,3,18,27,11,7,28,19,1,29,12,20,2,21,4,30,8,31,13,22

%N Difference between the prime indices of the two factors of the n-th semiprime.

%C Are there no adjacent equal terms? I have verified this up to n = 10^6. - _Gus Wiseman_, Dec 04 2020

%H Robert Israel, <a href="/A176506/b176506.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A049084(A084127(n)) - A049084(A084126(n)). [corrected by _R. J. Mathar_, Apr 22 2010]

%F a(n) = A338913(n) - A338912(n). - _Gus Wiseman_, Dec 04 2020

%e From _Gus Wiseman_, Dec 04 2020: (Start)

%e The sequence of semiprimes together with the corresponding differences begins:

%e 4: 1 - 1 = 0

%e 6: 2 - 1 = 1

%e 9: 2 - 2 = 0

%e 10: 3 - 1 = 2

%e 14: 4 - 1 = 3

%e 15: 3 - 2 = 1

%e 21: 4 - 2 = 2

%e 22: 5 - 1 = 4

%e 25: 3 - 3 = 0

%e 26: 6 - 1 = 5

%e 33: 5 - 2 = 3

%e (End)

%p isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:

%p A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:

%p A084126 := proc(n) min(op(numtheory[factorset](A001358(n)))) ; end proc:

%p A084127 := proc(n) max(op(numtheory[factorset](A001358(n)))) ; end proc:

%p A176506 := proc(n) numtheory[pi](A084127(n)) - numtheory[pi](A084126(n)) ; end proc: seq(A176506(n),n=1..120) ; # _R. J. Mathar_, Apr 22 2010

%p # Alternative:

%p N:= 500: # to use the first N semiprimes

%p Primes:= select(isprime, [2,seq(i,i=3..N/2,2)]):

%p SP:= NULL:

%p for i from 1 to nops(Primes) do

%p for j from 1 to i do

%p sp:= Primes[i]*Primes[j];

%p if sp > N then break fi;

%p SP:= SP, [sp, i-j]

%p od od:

%p SP:= sort([SP],(s,t) -> s[1]<t[1]):

%p map(t -> t[2], SP); # _Robert Israel_, Jan 17 2019

%t M = 500; (* to use the first M semiprimes *)

%t primes = Select[Join[{2}, Range[3, M/2, 2]], PrimeQ];

%t SP = {};

%t For[i = 1, i <= Length[primes], i++,

%t For[j = 1, j <= i, j++,

%t sp = primes[[i]] primes[[j]];

%t If[sp > M, Break []];

%t AppendTo[SP, {sp, i - j}]

%t ]];

%t SortBy[SP, First][[All, 2]] (* _Jean-François Alcover_, Jul 18 2020, after _Robert Israel_ *)

%t Table[If[!SquareFreeQ[n],0,-Subtract@@PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* _Gus Wiseman_, Dec 04 2020 *)

%o (PARI) lista(nn) = {my(vsp = select(x->(bigomega(x)==2), [1..nn])); vector(#vsp, k, my(f=factor(vsp[k])[,1]); primepi(vecmax(f)) - primepi(vecmin(f)));} \\ _Michel Marcus_, Jul 18 2020

%Y Cf. A109313.

%Y Cf. A001358, A049084, A084126, A084127.

%Y A087794 is product of the same indices.

%Y A176504 is the sum of the same indices.

%Y A115392 lists positions of first appearances.

%Y A128301 lists positions of 0's.

%Y A172348 lists positions of 1's.

%Y A338898 has this sequence as row differences.

%Y A338900 is the squarefree case.

%Y A338912/A338913 give the two prime indices of semiprimes.

%Y A006881 lists squarefree semiprimes.

%Y A024697 is the sum of semiprimes of weight n.

%Y A056239 gives sum of prime indices (Heinz weight).

%Y A087112 groups semiprimes by greater factor.

%Y A270650/A270652/A338899 give the prime indices of squarefree semiprimes.

%Y A338904 groups semiprimes by weight.

%Y A338907/A338906 list semiprimes of odd/even weight.

%Y A339114/A339115 give the least/greatest semiprime of weight n.

%Y Cf. A000040, A001222, A062198, A112141, A112798, A325699, A320655, A339005.

%K nonn,look

%O 1,4

%A _Juri-Stepan Gerasimov_, Apr 19 2010

%E a(51) and a(69) corrected by _R. J. Mathar_, Apr 22 2010