login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A295084 Number of sqrt(n)-smooth numbers <= n. 3

%I

%S 1,1,1,3,3,3,3,4,7,7,7,8,8,8,8,9,9,10,10,10,10,10,10,11,16,16,17,17,

%T 17,18,18,19,19,19,19,20,20,20,20,21,21,21,21,21,22,22,22,23,30,31,31,

%U 31,31,32,32,33,33,33,33,34,34,34,35,36,36,36,36,36,36,37,37,38,38,38,39,39,39,39,39,40

%N Number of sqrt(n)-smooth numbers <= n.

%C a(n) = number of positive integers m<=n such that A006530(m) <= sqrt(n).

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Smooth_number">Smooth number</a>

%F a(n) = n - A241419(n).

%F If n is in A063539, then a(n)=a(n-1)+1; if n is in A001248, i.e., n=p^2 for prime p, then a(n)=a(n-1)+p; otherwise a(n)=a(n-1).

%F a(n) = (1 - log(2))*n + O(n/log(n)) as n -> infinity. - _Robert Israel_, Nov 14 2017

%p N:= 100: # to get a(1)..a(N)

%p G:= [0,seq(max(numtheory:-factorset(n)),n=2..N)]:

%p seq(nops(select(t -> t^2 <= n, G[1..n])),n=1..N); # _Robert Israel_, Nov 14 2017

%o (PARI) A295084(n) = my(r=n); forprime(p=sqrtint(n)+1,n, r-=n\p); r;

%Y Cf. A048098 (indices of records), A063539, A241419.

%K nonn,look

%O 1,4

%A _Max Alekseyev_, Nov 13 2017

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 27 17:58 EST 2020. Contains 331296 sequences. (Running on oeis4.)