OFFSET
1,2
COMMENTS
Aebi and Cairns proved that if (p, p+2) are twin primes, then 2^(p+2) == 3p+8 (mod p(p+2)). This sequence contains the quotients of this congruence.
LINKS
Christian Aebi and Grant Cairns, Catalan numbers, primes, and twin primes, Elemente der Mathematik, Vol. 63, No. 4 (2008), pp. 153-164.
EXAMPLE
a(2) = 3, because A001359(2) = 5 and Q(5) = (2^(5+2)-3*5-8)/(5*(5+2)) = 3.
MATHEMATICA
s={}; Do[If[PrimeQ[n] && PrimeQ[n+2], q = (2^(n+2)-3n-8)/(n(n+2)); AppendTo[s, q]], {n, 1, 1000}]; s
PROG
(PARI) lista(nn) = {forprime(p=2, nn, if (isprime(p+2), print1((2^(p+2)-3*p-8)/(p*(p+2)), ", "); ); ); } \\ Michel Marcus, Aug 04 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Aug 04 2019
STATUS
approved