A128874 Numbers n such that 16 + the product of the first n lesser twin primes (A097489) is prime.

1, 2, 3, 469

Offset

1,2

Comments

It is interesting from the example that the first three primes 19,31 and 181 greater than twinl#(n) + 16 are all greater twin primes. The next PRP is the 1824 digit number twinl#(469) + 16 and is not part of a twin prime pair.

Links

Table of n, a(n) for n=1..4.

Formula

Define twinl#(n)as the product of the first n lesser twin primes. Then if twinl#(n)+ 16 is prime, list n.

Example

Twinl#(3) + 16 = 3*5*11+16 = 181 is prime so 3 is in the sequence.

Prog

(PARI)
twinl(n) = { local(c, x); c=0; x=1; while(c<n, if(isprime(prime(x)+2), c++); x++; ); return(prime(x-1)) } \\ The n-th lower twin prime.
twiprimesl2(n, a) = { local(pr, x, y, j); for(a=0, n, for(j=1, n, pr=1; for(x=1, j, pr*=twinl(x); ); y=pr+2^a; if(ispseudoprime(y), print1(y", "); break ) ) ) } \\ Cino Hilliard, May 08 2007
(PARI) lista(nn) = { for(n=1, nn, if (ispseudoprime(16+prod(i=1, n, twinl(i))), print1(n, ", "))); } \\ Michel Marcus, Feb 10 2014

Crossrefs

Sequence in context: A109855 A306370 A355644 * A196070 A173342 A090510

Adjacent sequences: A128871 A128872 A128873 * A128875 A128876 A128877

Keyword

more,nonn

Author

Cino Hilliard, May 08 2007

Extensions

Partially edited by Michel Marcus, Feb 10 2014

Status

approved