OFFSET
3,1
MATHEMATICA
<< NumberTheory`NumberTheoryFunctions` a = {}; Do[AppendTo[a, PreviousPrime[Sqrt[(Prime[x - 2])*(Prime[x + 2])]]], {x, 3, 100}]; a
Table[NextPrime[GeometricMean[{Prime[n-2], Prime[n+2]}], -1], {n, 3, 60}] (* Harvey P. Dale, May 14 2015 *)
PROG
(PARI) A126991(n)={ n=sqrtint(prime(n-2)*prime(n+2)); if(0==n%2, n--); while(!isprime(n), n-=2); n } /* then vector(50, n, A126991(n+2)) displays a list of values ; t=3; forprime(p=2, 999, while(A126991(t)<p, t++); if(A126991(t)>p, print(p))) prints primes 2, 17, 41, 59, ... not in this sequence */ \\ M. F. Hasler, Jun 14 2007
CROSSREFS
KEYWORD
nonn,easy,less
AUTHOR
Artur Jasinski, Jan 01 2007
EXTENSIONS
Edited by M. F. Hasler, Jun 14 2007
STATUS
approved