OFFSET
1,2
COMMENTS
A006530(a(n))^2 <= a(n). - Reinhard Zumkeller, Oct 12 2011
This set (say S) has density d(S) = 1-log(2) and multiplicative density m(S) = 1-exp(-Gamma). Multiplicative density: let A be a set of numbers, A(x) = { k in A | gpf(k) <=x } where gpf(k) denotes the greatest prime factor of k and let m(x)(A) = Product_{p<=x} (1 - 1/p)*Sum_{k in A(x)} 1/k. If lim_{x->infinity} m(x)(A) exists = m(A), this limit is called "multiplicative density" of A (Erdős and Davenport, 1951). - Benoit Cloitre, Jun 12 2002
LINKS
T. D. Noe and William A. Tedeschi, Table of n, a(n) for n = 1..10000 (first 1000 terms computed by T. D. Noe)
H. Davenport and P. Erdős, On sequences of positive integers, J. Indian Math. Soc. 15 (1951), pp. 19-24.
Eric Weisstein's World of Mathematics, Greatest Prime Factor
Eric Weisstein's World of Mathematics, Round Number
MATHEMATICA
gpf[n_] := FactorInteger[n][[-1, 1]]; A048098 = {}; For[n = 1, n <= 200, n++, If[ gpf[n] <= Sqrt[n], AppendTo[ A048098, n] ] ]; A048098 (* Jean-François Alcover, Jan 26 2012 *)
PROG
(PARI)
print1(1, ", "); for(n=2, 1000, if(vecmax(factor(n)[, 1])<=sqrt(n), print1(n, ", ")))
(Haskell)
a048098 n = a048098_list !! (n-1)
a048098_list = [x | x <- [1..], a006530 x ^ 2 <= x]
-- Reinhard Zumkeller, Oct 12 2011
(Python)
from sympy import factorint
def ok(n):
return n == 1 if n < 2 else max(factorint(n))**2 <= n
print([k for k in range(196) if ok(k)]) # Michael S. Branicky, Dec 22 2021
(Python)
from math import isqrt
from sympy import primepi
def A048098(n):
def f(x): return int(n+sum(primepi(x//i)-primepi(i) for i in range(1, isqrt(x)+1)))
def bisection(f, kmin=0, kmax=1):
while f(kmax) > kmax: kmax <<= 1
while kmax-kmin > 1:
kmid = kmax+kmin>>1
if f(kmid) <= kmid:
kmax = kmid
else:
kmin = kmid
return kmax
return bisection(f) # Chai Wah Wu, Sep 01 2024
CROSSREFS
KEYWORD
easy,nonn,nice
AUTHOR
EXTENSIONS
More terms from James A. Sellers, Apr 22 2000
Edited by Charles R Greathouse IV, Nov 08 2010
STATUS
approved