|
%I #25 Oct 27 2015 18:05:57
%S 10,1060,1795,1885,2965,3245,3335,4065,4325,5015,5875,6985,7605,7905,
%T 9785,11315,12045,12360,14390,14970,15285,15500,15885,17195,18220,
%U 20670,20695,22160,24915,25645,25955,26025,29410,29910,32925,35530,36280
%N Numbers n such that each of p=6*n+1, q=6*p+1, r=6*q+1 and s=6*r+1 is prime.
%C Each p is a starting prime of a complete generalized Cunningham chain p(k)=6*p(k-1)+1.
%C All terms are multiples of 5. Hence t = 6s+1 = 1555+7776n are always composite, and the chains are indeed "complete."
%C Subsequence of A263310 (and as such of A263309 and of A024899).
%H Zak Seidov, <a href="/A263311/b263311.txt">Table of n, a(n) for n = 1..20000</a>
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Cunningham_chain">Cunningham chain</a>
%p isA263311 := proc(n)
%p return isprime(6*n+1) and isprime(36*n+7) and isprime(216*n+43) and isprime(1296*n+259) ;
%p end proc:
%p for n from 1 to 30000 do
%p if isA263311(n) then
%p printf("%d,",n);
%p end if;
%p end do; # _R. J. Mathar_, Oct 17 2015
%t Select[Range[10,100000,5],PrimeQ[p=6*#+1]&&PrimeQ[q=6*p+1]&&PrimeQ[r=6*q+1]&&PrimeQ[s=6*r+1]&]
%o (PARI) for(n=1, 1e5, if(isprime(p=6*n+1) && isprime(q=6*p+1) && isprime(r=6*q+1) && isprime(6*r+1), print1(n", "))) \\ _Altug Alkan_, Oct 17 2015
%Y Cf. A024899, A263309, A263310.
%K nonn
%O 1,1
%A _Zak Seidov_, Oct 13 2015
|
