|
%I #36 Oct 26 2021 02:14:48
%S 5,11,13,29,43,59,61,67,83,131,139,157,173,211,227,229,277,283,317,
%T 331,347,373,389,419,421,443,461,509,547,563,571,619,643,653,659,661,
%U 691,709,733,787,797,821,853,859,877,907,941,947,997,1019,1021,1069,1091,1093,1109,1123,1163,1181,1213
%N Primes equal to the sum of a pair of consecutive integers which are both squarefree.
%C The associated prime factors will always include 2 and 3.
%C Will every prime number be encountered as a prime factor from the sequence entries?
%C The sequence appears to share many of it terms with A001122.
%C What is the asymptotic behavior?
%C Conjecture: sequence has density A271780/2 = A005597*4/Pi^2 = 0.2675535... in the primes. - _Charles R Greathouse IV_, Jan 24 2018
%C The prime terms of A179017 (except 3). - _Bill McEachen_, Oct 21 2021
%H Charles R Greathouse IV, <a href="/A269844/b269844.txt">Table of n, a(n) for n = 1..10000</a>
%H Bill McEachen, <a href="/A269844/a269844.png">A269844_vs_A001122</a>
%e 277 = 138 + 139 = 2*3*23 + 139 is in the sequence since both terms are squarefree.
%e 281 = 140 + 141 = 2^2*5*7 + 3*47 is not in the sequence since the former term is divisible by 2^2.
%t Select[Prime@ Range[3, 200], PrimeOmega@ # == PrimeNu@ # &[# (# + 1)] &@ Floor[#/2] &] (* _Michael De Vlieger_, Mar 07 2016 *)
%o (PARI)
%o genit(maxx)={for(i5=3,maxx,n=prime(i5);a=factor(floor(n/2.));b=factor(ceil(n/2.));clear=1;for(j5=1,omega(floor(n/2.)),if(a[j5,2]<>1,clear=0));
%o for(j7=1,omega(ceil(n/2.)),if(b[j7,2]<>1,clear=0));if(clear>0,print1(n,",")));}
%o (PARI) is(n)=isprime(n) && issquarefree(n\2) && issquarefree(n\2+1) \\ _Charles R Greathouse IV_, Jan 24 2018
%o (PARI) list(lim)=my(v=List(),t=1); forfactored(k=3,(lim+1)\2, if(vecmax(k[2][,2])>1, t=0, ; if(t && isprime(t=2*k[1]-1), listput(v,t)); t=1)); Vec(v) \\ _Charles R Greathouse IV_, Jan 24 2018
%Y Cf. A001122 (primes with primitive root 2), A179017.
%Y Cf. A005597, A271780.
%K nonn,easy
%O 1,1
%A _Bill McEachen_, Mar 06 2016
|
