cp's OEIS Frontend

A224610 Smallest j such that j*2*prime(n)^3-1 and j*2*prime(n)*q^2-1 are prime.

%I A224610 #22 Apr 22 2013 12:08:03
%S A224610 2,2,5,7,59,142,264,25,8,21,124,33,60,87,9,231,5,6,82,155,7,66,72,21,
%T A224610 42,105,15,48,250,68,222,54,47,195,255,360,205,6,83,26,5,1,50,220,173,
%U A224610 1,976,30,228,130,30,129,46,1106,65,62,15,109,24,41,922,15,132,89
%N A224610 Smallest j such that j*2*prime(n)^3-1 and j*2*prime(n)*q^2-1 are prime.
%H A224610 Pierre CAMI, <a href="/A224610/b224610.txt">Table of n, a(n) for n = 1..10000</a>
%t A224610 a[n_] := For[p = Prime[n]; j = 1, j < 10^6, j++, If[PrimeQ[q = j*2*p^3 - 1] && PrimeQ[j*p*2*q^2 - 1], Return[j]]]; Table[a[n], {n, 1, 75}] (* _Jean-François Alcover_, Apr 22 2013 *)
%o A224610 (PFGW & SCRIPTIFY)
%o A224610 SCRIPT
%o A224610 DIM n,1
%o A224610 DIM j
%o A224610 DIM q
%o A224610 DIMS t
%o A224610 OPENFILEOUT myf,a(n).txt
%o A224610 LABEL a
%o A224610 SET n,n+1
%o A224610 SET j,0
%o A224610 LABEL b
%o A224610 SET j,j+1
%o A224610 SET q,j*2*p(n)^3-1
%o A224610 SETS t,%d\,;p(n)
%o A224610 PRP q,t
%o A224610 IF ISPRP THEN GOTO c
%o A224610 GOTO b
%o A224610 LABEL c
%o A224610 SET q,j*2*p(n)*q^2-1
%o A224610 PRP q,t
%o A224610 IF ISPRP THEN GOTO d
%o A224610 GOTO b
%o A224610 LABEL d
%o A224610 WRITE myf,t
%o A224610 GOTO a
%Y A224610 Cf. A224490, A224609, A224611, A224612.
%K A224610 nonn
%O A224610 1,1
%A A224610 _Pierre CAMI_, Apr 12 2013