cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A075706 Safe primes (A005385) (p and (p-1)/2 are primes) such that 8*p+1 (A023228) is also prime.

This page as a plain text file.
%I A075706 #9 Jan 31 2020 12:26:52
%S A075706 5,11,107,179,347,479,1187,1307,1367,1487,1619,2027,2207,2447,2999,
%T A075706 3119,3467,4007,4079,4139,4799,5087,5807,5927,5939,6827,7079,7247,
%U A075706 8699,9587,9839,10607,12107,12539,12659,14207,15299,16139,16187,17027
%N A075706 Safe primes (A005385) (p and (p-1)/2 are primes) such that 8*p+1 (A023228) is also prime.
%H A075706 Harvey P. Dale, <a href="/A075706/b075706.txt">Table of n, a(n) for n = 1..1000</a>
%e A075706 11 is a prime, so is (11-1)/2=5 and also 8*11+1=89; 107 is a prime, (107-1)/2=53 and 8*107+1=857, ...
%p A075706 ts_sg8_var_pras := proc(nmax) local i,tren,atek; tren := 0: for i from 1 to nmax do atek := numtheory[safeprime](i): if (atek > tren) then if (isprime(atek)='true' and isprime(6*atek+1)='true') then tren := atek: fi; fi; od; end: seq(ts_sg8_var_pras(i), i=1..3000);
%t A075706 Select[Prime[Range[2000]],AllTrue[{(#-1)/2,8#+1},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jan 31 2020 *)
%o A075706 (PARI) forprime(p=3,20000,if(isprime((p-1)/2),if(isprime(8*p+1),print1(p","))))
%Y A075706 Cf. A005384, A005385, A059455, A023228.
%K A075706 nonn
%O A075706 1,1
%A A075706 _Jani Melik_, Oct 02 2002
%E A075706 More terms from _Ralf Stephan_, Mar 19 2003