triple factorial primes

[Maximilian Hasler]

Neil and dear seqfans,
for

A037083  	 	 Numbers n such that n!!! + 1 is prime.  	 	 +0
3
	1, 2, 4, 5, 6, 7, 9, 10, 11, 17, 24, 29, 39, 40, 57, 58, 59, 91, 155,
175, 977, 1077, 1327, 2076

I find several more terms between 175 and 977. Can anyone confirm
this, or is this a bug in my PARI program or its isprime() function ?
(IIRC, in PARI isprime() is NOT ispeudoprime() as in Maple.)
(The corresponding link in the page
http://www.prothsearch.net/multifac.html cited in A037083 is broken.)
In case there is no objection, I would like to ask Neil to add the
missing terms, PARI code (esp. for A007661) and the b-files (well, for
b037083 this is maybe not very useful, the "table" link will do).

(PARI)
A007661(n) = prod(i=1,(n-1)\3,n-=3,n+!n)
c=0;for(n=1,999,if(isprime(A007661(n)+1),print1("b037083.txt",c++," "n)))
c=0;for(n=1,999,if(isprime(t=A007661(n)+1),write("b037082.txt",c++," "t)))

Regards,
Maximilian Hasler

1 1
2 2
3 4
4 5
5 6
6 7
7 9
8 10
9 11
10 17
11 24
12 29
13 39
14 40
15 57
16 58
17 59
18 91
19 155
20 175
21 245
22 359
23 372
24 597
25 864
26 977
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: b037083.txt
URL: <http://list.seqfan.eu/pipermail/seqfan/attachments/20071126/6e28c873/attachment.txt>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: b037082.txt
URL: <http://list.seqfan.eu/pipermail/seqfan/attachments/20071126/6e28c873/attachment-0001.txt>