|
%I #7 May 05 2021 17:01:11
%S 5,5,5,5,7,7,11,11,11,11,23,23,23,23,23,23,23,23,23,23,23,23,47,47,47,
%T 47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47,59,59,
%U 59,59,59,59,59,59,59,59,59,59,83,83,83,83,83,83,83,83,83
%N Smallest safe prime ((p-1)/2 is also prime) > n.
%H Michael S. Branicky, <a href="/A033299/b033299.txt">Table of n, a(n) for n = 1..10000</a>
%p with(numtheory); safeprime(n);
%o (Python)
%o from sympy import isprime, nextprime
%o def a(n):
%o p = nextprime(n)
%o while not isprime((p-1)//2): p = nextprime(p)
%o return p
%o print([a(n) for n in range(1, 68)]) # _Michael S. Branicky_, May 05 2021
%K nonn
%O 1,1
%A _N. J. A. Sloane_.
%E a(57) and beyond from _Michael S. Branicky_, May 05 2021
|
