login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A127492 Indices m of primes such that Sum_{k=0..2, k<j<=3} prime(m+k)*prime(m+j)*prime(m+j+1) is twice a prime. 2
2, 10, 17, 49, 71, 72, 75, 145, 161, 167, 170, 184, 244, 250, 257, 266, 267, 282, 286, 301, 307, 325, 343, 391, 405, 429, 450, 537, 556, 561, 584, 685, 710, 743, 790, 835, 861, 904, 928, 953 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Let p_0 .. p_4 be five consecutive primes, starting with the m-th prime. The index m is in the sequence if the absolute value [x^0] of the polynomial (x-p_0)*[(x-p_1)*(x-p_2) + (x-p_2)*(x-p_3) + (x-p_3)*(x-p_4)] + (x-p_1)*[(x-p_2)*(x-p_3) + (x-p_3)*(x-p_4)] + (x-p_2)*(x-p_3)*(x-p_4) is two times a prime. The correspondence with A127491: the coefficient [x^2] of the polynomial (x-p_0)*(x-p_1)*..*(x-p_4) is the sum of 10 products of a set of 3 out of the 5 primes. Here the sum is restricted to the 6 products where the two largest of the 3 primes are consecutive. - R. J. Mathar, Apr 23 2023
LINKS
MAPLE
isA127492 := proc(k)
local x, j ;
(x-ithprime(k))* mul( x-ithprime(k+j), j=1..2)
+(x-ithprime(k))* mul( x-ithprime(k+j), j=2..3)
+(x-ithprime(k))* mul( x-ithprime(k+j), j=3..4)
+(x-ithprime(k+1))* mul( x-ithprime(k+j), j=2..3)
+(x-ithprime(k+1))* mul( x-ithprime(k+j), j=3..4)
+(x-ithprime(k+2))* mul( x-ithprime(k+j), j=3..4) ;
p := abs(coeff(expand(%/2), x, 0)) ;
if type(p, 'integer') then
isprime(p) ;
else
false ;
end if ;
end proc:
for k from 1 to 900 do
if isA127492(k) then
printf("%a, ", k) ;
end if ;
end do: # R. J. Mathar, Apr 23 2023
MATHEMATICA
a = {}; Do[If[PrimeQ[(Prime[x] Prime[x + 1]Prime[x + 2] + Prime[x] Prime[x + 2]Prime[x + 3] + Prime[x] Prime[x + 3] Prime[x + 4] + Prime[x + 1] Prime[x + 2]Prime[x + 3] + Prime[x + 1] Prime[x + 3]Prime[x + 4] + Prime[x + 2] Prime[x + 3] Prime[x + 4])/2], AppendTo[a, x]], {x, 1, 1000}]; a
prQ[{a_, b_, c_, d_, e_}]:=PrimeQ[(a b c+a c d+a d e+b c d+b d e+c d e)/2]; PrimePi/@Select[ Partition[ Prime[Range[1000]], 5, 1], prQ][[;; , 1]] (* Harvey P. Dale, Apr 21 2023 *)
CROSSREFS
Sequence in context: A304809 A304805 A214086 * A258974 A366508 A077247
KEYWORD
nonn,uned,obsc
AUTHOR
Artur Jasinski, Jan 16 2007
EXTENSIONS
Definition simplified by R. J. Mathar, Apr 23 2023
Edited by Jon E. Schoenfield, Jul 23 2023
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 06:14 EDT 2024. Contains 371964 sequences. (Running on oeis4.)