login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A197680 Numbers whose exponents in their prime power factorization are squares. 17

%I #74 Jan 20 2024 07:26:34

%S 1,2,3,5,6,7,10,11,13,14,15,16,17,19,21,22,23,26,29,30,31,33,34,35,37,

%T 38,39,41,42,43,46,47,48,51,53,55,57,58,59,61,62,65,66,67,69,70,71,73,

%U 74,77,78,79,80,81,82,83,85,86,87,89,91,93,94,95,97,101

%N Numbers whose exponents in their prime power factorization are squares.

%C Numbers whose prime factorization has the form Product_i p_i^e_i where the e_i are all squares.

%C All squarefree numbers (A005117) are in the sequence. - _Vladimir Shevelev_, Nov 16 2015

%C Let h_k be the density of the subsequence of A197680 of numbers whose prime power factorization (PPF) has the form Product_i p_i^e_i where the e_i all squares <= k^2. Then for every k>1 there exists eps_k>0 such that for any x from the interval (h_k-eps_k, h_k) there is no a sequence S of positive integers such that x is the density of numbers whose PPF has the form Product_i p_i^e_i where the e_i are all in S. - For a proof, see [Shevelev], second link. - _Vladimir Shevelev_, Nov 17 2015

%C Numbers with an odd number of exponential divisors (A049419). - _Amiram Eldar_, Nov 05 2021

%H Alois P. Heinz, <a href="/A197680/b197680.txt">Table of n, a(n) for n = 1..10000</a>

%H StackExchange, <a href="http://math.stackexchange.com/questions/73354/two-dirichlets-series-related-to-the-divisor-summatory-function-and-to-the-riem">Question 73354</a>, 2011.

%H Vladimir Shevelev, <a href="http://arxiv.org/abs/1510.05914">Exponentially S-numbers</a>, arXiv:1510.05914 [math.NT], 2015.

%H Vladimir Shevelev, <a href="http://arxiv.org/abs/1511.03860">Set of all densities of exponentially S-numbers</a>, arXiv:1511.03860 [math.NT], 2015.

%H Vladimir Shevelev, <a href="http://dx.doi.org/10.4064/aa8395-5-2016">S-exponential numbers</a>, Acta Arithmetica, Vol. 175 (2016), pp. 385-395.

%F Sum_{i<=x, i is in A197680} 1 = h*x + O(sqrt(x)*log x*e^(c*sqrt(log x)/(log(log x))), where c=4*sqrt(2.4/log 2)=7.44308... and h=Product_{prime p} (1+Sum_{i>=2} (u(i)-u(i-1))/p^i)=0.641115... where u(n) is the characteristic function of sequence A000290. The calculations of h in the formula were done independently by _Juan Arias-de-Reyna_ and _Peter J. C. Moses_. For a proof of the formula, see the first Shevelev link. - _Vladimir Shevelev_, Nov 17 2015

%p a:= proc(n) option remember; local k; for k from 1+

%p `if`(n=1, 0, a(n-1)) while 0=mul(`if`(issqr(

%p i[2]), 1, 0), i=ifactors(k)[2]) do od; k

%p end:

%p seq(a(n), n=1..80); # _Alois P. Heinz_, Jun 30 2016

%t Select[Range[100], Union[IntegerQ /@ Sqrt[Transpose[FactorInteger[#]][[2]]]][[1]] &] (* _T. D. Noe_, Oct 18 2011 *)

%o (PARI) isok(n) = {my(f = factor(n)[,2]); #select(x->issquare(x), f) == #f; } \\ _Michel Marcus_, Oct 23 2015

%Y Cf. A000290, A049419, A005117, A115063, A209061.

%K nonn

%O 1,2

%A _A. Neves_, Oct 17 2011

%E Reformulation of the name by _Vladimir Shevelev_, Oct 14 2015

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 04:23 EDT 2024. Contains 371264 sequences. (Running on oeis4.)