A251726 Numbers n > 1 for which gpf(n) < lpf(n)^2, where lpf and gpf (least and greatest prime factor of n) are given by A020639(n) and A006530(n).
2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 15, 16, 17, 18, 19, 21, 23, 24, 25, 27, 29, 31, 32, 35, 36, 37, 41, 43, 45, 47, 48, 49, 53, 54, 55, 59, 61, 63, 64, 65, 67, 71, 72, 73, 75, 77, 79, 81, 83, 85, 89, 91, 95, 96, 97, 101, 103, 105, 107, 108, 109, 113, 115, 119, 121, 125, 127, 128, 131, 133, 135, 137, 139, 143, 144
Offset: 1
Keywords
Examples
For 35 = 5*7, 7 is less than 5^2, thus 35 is included. For 90 = 2*3*3*5, 5 is not less than 2^2, thus 90 is NOT included. For 105 = 3*5*7, 7 is less than 3^2, thus 105 is included.
Links
- Antti Karttunen, Table of n, a(n) for n = 1..10000
- Ratio A251726(n)/A251727(n) (plotted with OEIS-server's Plot2-utility)
Crossrefs
Programs
-
Mathematica
pfQ[n_]:=Module[{f=FactorInteger[n]},f[[-1,1]]
Harvey P. Dale, May 01 2015 *) -
PARI
for(n=2, 150, if(vecmax(factor(n)[,1]) < vecmin(factor(n)[,1])^2, print1(n,", "))) \\ Indranil Ghosh, Mar 24 2017
-
Python
from sympy import primefactors print([n for n in range(2, 150) if max(primefactors(n))
Indranil Ghosh, Mar 24 2017
Formula
Extensions
A new simpler definition found Jan 01 2015 and the original definition moved to the Comments section.
Comments