|
%I #33 May 02 2023 02:28:22
%S 5,13,61,181,349,421,661,829,1453,1621,1669,2029,2269,3469,3853,4021,
%T 4093,4261,4789,6781,6829,6949,7549,8221,8293,8821,9421,10069,10093,
%U 10141,10501,10861,12253,12613,13933,14389,14629,14869,16069,16189,16981,17389,17749
%N Greater of the pairs of twin primes in A001122.
%C Primes p such that both p - 2 and p are both in A001122.
%C Apart from the first term, all terms are congruent to 13 mod 24, since terms in A006512 are congruent to 1 mod 6 apart from the first one, and terms in A001122 are congruent to 3 or 5 mod 8.
%C Note that "there are infinitely many pairs of twin primes" and "there are infinitely many primes with primitive root 2" are two famous and unsolved problems, so a stronger conjecture implying both of them is that this sequence is infinite.
%C Also note that a pair of cousin primes can't both appear in A001122, while a pair of sexy primes can.
%H Amiram Eldar, <a href="/A319249/b319249.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..406 from Jianing Song)
%F a(n) = A319248(n) + 2.
%F For n >= 2, a(n) = 24*A319250(n-1) + 13.
%e 11 and 13 is a pair of twin primes both having 2 as a primitive root, so 13 is a term.
%e 59 and 61 is a pair of twin primes both having 2 as a primitive root, so 61 is a term.
%e Although 137 and 139 is a pair of twin primes, 139 has 2 as a primitive root while 137 doesn't, so 139 is not a term.
%t Select[Prime[Range[2^11]], PrimeQ[# - 2] && PrimitiveRoot[# - 2] == 2 && PrimitiveRoot[#] == 2 &] (* _Amiram Eldar_, May 02 2023 *)
%o (PARI) forprime(p=3, 10000, if(znorder(Mod(2,p))==p-1 && znorder(Mod(2,p+2))==p+1, print1(p+2, ", ")))
%o (Python)
%o from itertools import islice
%o from sympy import isprime, nextprime, is_primitive_root
%o def A319249_gen(): # generator of terms
%o p = 2
%o while (p:=nextprime(p)):
%o if isprime(p+2) and is_primitive_root(2,p) and is_primitive_root(2,p+2):
%o yield p+2
%o A319249_list = list(islice(A319249_gen(),30)) # _Chai Wah Wu_, Feb 13 2023
%Y Cf. A001122, A001359, A006512.
%Y A319248 gives p-2, A319250 gives (p-13)/24.
%K nonn
%O 1,1
%A _Jianing Song_, Sep 15 2018
|
