OFFSET
1,1
COMMENTS
The notation # refers to the primorials A002110, the partial products of primes.
Roughly 75% of the entries are nonzero.
For roughly 50% of the solutions (that is roughly 1/3 of all entries if zeros are included) the q are smaller than prime(n)/5.
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..750
EXAMPLE
2*2 + 1 = 5 (prime) with q=2, q#=2, prime(n)# = 2 so a(1)=2.
2*2*3 + 1 = 13 (prime) with q=2, q#=2, prime(2)# = 2*3 so a(2)=2.
2*2*3*5 + 1 = 61 (prime) with q=2, q#=2, prime(3)# = 2*3*5 so a(3)=2.
2*2*3*5*7 + 1 = 421 (prime) with q=2, q#=2, prime(4)# = 2*3*5*7 so a(4)=2.
MAPLE
A002110 := proc(n) option remember; mul(ithprime(i), i=1..n) ; end proc:
A190617 := proc(n) local psharp ; psharp := A002110(n) ; for i from 1 to n do if isprime(1+psharp*A002110(i)) then return ithprime(i) ; end if; end do: return 0 ; end proc:
seq(A190617(n), n=1..80) ; # R. J. Mathar, Jun 02 2011
PROG
PFGW from Primeformgroup for prime search and certification
pfgw64 -f in.txt , results in pfgw-prime.log and pfgw.log
in.txt scriptyfile
SCRIPT
DIM nn, 0
DIM kk
DIM mm
DIM jj
DIMS tt
LABEL loopn
SET nn, nn+1
IF nn>750 THEN END
SET kk, p(nn)
SET mm, 0
LABEL loopm
SET mm, mm+1
IF mm>nn THEN GOTO loopn
SET jj, p(mm)
SETS tt, %d, %d\,; kk; jj
PRP kk#*(jj#)+1, tt
IF ISPRIME THEN GOTO loopn
IF ISPRP THEN GOTO loopn
goto loopm
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, May 14 2011
STATUS
approved