|
%I #21 Sep 08 2022 08:46:09
%S 2,3,5,7,13,17,43,61,67,73,97,127,131,137,167,241,281,307,353,433,463,
%T 467,541,557,631,641,647,653,661,673,683,821,853,857,907,911,991,1033,
%U 1063,1103,1117,1123,1291,1307,1433,1453,1511,1523,1553,1567,1571,1597,1601,1607
%N Prime numbers P such that 24*P^3-1 is also prime.
%H Pierre CAMI, <a href="/A245640/b245640.txt">Table of n, a(n) for n = 1..10000</a>
%e 24*2^3-1=191 prime so a(1)=2.
%e 24*3^3-1=647 prime so a(2)=3.
%t Select[Prime[Range[300]], PrimeQ[24 #^3 - 1] &] (* _Vincenzo Librandi_, Sep 07 2014 *)
%o (PFGW & SCRIPT)
%o SCRIPT
%o DIM n,0
%o DIMS t
%o OPENFILEOUT myf,a(n).txt
%o LABEL loop1
%o SET n,n+1
%o SETS t,%d,%d\,;n;p(n)
%o PRP 24*p(n)^3-1,t
%o IF ISPRP THEN GOTO a
%o GOTO loop1
%o LABEL a
%o WRITE myf,t
%o GOTO loop1
%o (PARI) select(p->isprime(24*p^3-1), primes(300)) \\ _Colin Barker_, Jul 28 2014
%o (Magma) [p: p in PrimesUpTo(1700)| IsPrime(24*p^3-1)]; // _Vincenzo Librandi_, Sep 07 2014
%Y Cf. A245639, A245641.
%K nonn,easy
%O 1,1
%A _Pierre CAMI_, Jul 28 2014
|
