OFFSET
1,1
MAPLE
with(numtheory); a := []; for i from 2 to 1000 do p := ithprime(i); t1 := factorset(p-1); q := t1[nops(t1)]; if q^4 > p^3 then a := [op(a), p]; fi; od:
MATHEMATICA
Select[Prime[Range[300]], FactorInteger[#-1][[-1, 1]]>Surd[#^3, 4]&] (* Harvey P. Dale, Feb 27 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Aug 23 2002
STATUS
approved