OFFSET
1,1
EXAMPLE
37^3 - 1 = 2^2*3^3*7*67, 37^3 + 1 = 2*19*31*43 and 37*3 > 43 and 37*3 > 67
MATHEMATICA
Select[Prime[Range[PrimePi[3000]]], 3 # > FactorInteger[#^3 - 1][[-1, 1]] && 3 # > FactorInteger[#^3 + 1][[-1, 1]] &] (* T. D. Noe, Nov 07 2012 *)
PROG
(PARI) forprime(h=3, 5001, for(n=3, 3, a=h^n; b=a-1; c=a+1; d=vecmax(factor(b)[, 1]~); e=vecmax(factor(c)[, 1]~); g=h*n; if(g>d && g>e, print1(h, ", "))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Robin Garcia, Nov 07 2012
STATUS
approved