OFFSET
1,1
COMMENTS
Includes the first 10 primes.
If p is prime, the next four primes are p+2, p+6, p+12 and p+14, and p^2+6*p+48 is prime, then the twin primes p and p+2 are in the sequence, with the same value of A341937.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(5) = 11 is a term because (p,q,r,s)=(11,13,17,19) are consecutive primes with 2*p*q-2*q*r+r*s = 167, which is prime.
MAPLE
P:= select(isprime, [2, seq(i, i=3..10000, 2)]):
B:= select(i -> isprime(P[i+2]*P[i+3]-2*P[i+1]*(P[i+2]-P[i])), [$1..nops(P)-3]):
P[B];
PROG
(Python)
from sympy import isprime, nextprime
p, q, r, s, A341934_list=2, 3, 5, 7, []
while p < 10**6:
if isprime(2*q*(p-r)+r*s):
A341934_list.append(p)
p, q, r, s = q, r, s, nextprime(s) # Chai Wah Wu, Feb 24 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 23 2021
STATUS
approved