%I #54 Sep 14 2021 12:11:57
%S 9,14,21,26,403,12367,41303,66893,68297,73147,111607,116813,118003,
%T 130133,146873,222757,260497,418307,429491,439097,478061,559003,
%U 628241,729007,822397,1116707,1239869,1595683,1887239,2148589,2225669,2481463,2502977,2539553
%N The lesser of two semiprime brothers.
%C Semiprime brothers are two consecutive semiprimes (A001358) whose prime factors are consecutive primes (A000040).
%C The first several examples of semiprime brothers are {9, 10}, {14, 15}, {21, 22}, {26, 33} & {403, 407}.
%C The only square term is 9 and the only even terms are 14 and 26.
%C Obviously the difference between the primepi of the factors of the two consecutive semiprimes is either {-1, 1} or {1, -1}.
%C Number of terms < 10^n: 1, 4, 5, 5, 10, 25, 62, 143, 319, 761, 2010, 5275, etc.
%C Only the first three terms have as the next semiprime the next integer making them twins. - _Robert G. Wilson v_, Jun 21 2018
%H Jonathan Vos Post, Robert G. Wilson v, and Giovanni Resta, <a href="/A289548/b289548.txt">Table of n, a(n) for n = 1..5275</a> (terms < 10^12, terms > 10^10 from G. Resta)
%e 26 is in the sequence because 26 = 2*13 and the next semiprime is 33 = 3*11 with 2 & 3 consecutive primes and 11 & 13 consecutive primes.
%e 403 is in the sequence because 403 = 13*31 and the next semiprime is 407 = 11*37 with 11 & 13 and 31 & 37 being consecutive primes.
%t p = q = 4; fp = fq = {1, 1}; lst = {}; While[p < 26000000, While[fq = Flatten[ Table[#1, {#2}] & @@@ FactorInteger@ q]; Length@ fq != 2, q++]; d = Sort[{fp, fq}]; If[ NextPrime[ d[[1, 1]]] == d[[2, 1]] && NextPrime[ d[[2, 2]]] == d[[1, 2]], AppendTo[lst, p]]; p = q; fp = fq; q++]; lst
%o (PARI) isok(p, q) = (nextprime(p+1) == q) || (nextprime(q+1) == p);
%o pairp(n) = if (issquare(n), vector(2, k, sqrtint(n)), (factor(n)[,1])~);
%o lista(nn) = {na = 2; while (na < nn, if (bigomega(na) != 2, na++, nb = na + 1; while (bigomega(nb) != 2, nb++); fpa = pairp(na); fpb = pairp(nb); if (isok(fpa[1], fpb[1]) && isok(fpa[2], fpb[2]), print1(na, ", ")); na = nb;););} \\ _Michel Marcus_, Jul 11 2017
%o (Python)
%o from sympy import factorint, nextprime
%o def is_semiprime(n):
%o return sum(e for e in factorint(n).values()) == 2
%o def next_semiprime(n):
%o nxt = n + 1
%o while not is_semiprime(nxt): nxt += 1
%o return nxt
%o def are_consecutive(p, q):
%o return max(p, q) == nextprime(min(p, q))
%o def ok(n):
%o if not is_semiprime(n): return False
%o nextsp = next_semiprime(n)
%o fn, fm = factorint(n, multiple=True), factorint(nextsp, multiple=True)
%o return are_consecutive(fn[0], fm[0]) and are_consecutive(fn[1], fm[1])
%o print(list(filter(ok, range(150000)))) # _Michael S. Branicky_, Sep 14 2021
%Y Cf. A000040, A001358.
%K nonn
%O 1,1
%A _Jonathan Vos Post_ and _Robert G. Wilson v_, Jul 07 2017
|