OFFSET
1,2
COMMENTS
The sequence is formed by 1 and the numbers whose prime power factorization contains only odd exponents.
The density of the sequence is the constant given by A065463.
Except for the first term the same as A002035. - R. J. Mathar, Feb 07 2016
Also numbers k all of whose divisors are bi-unitary divisors (i.e., A286324(k) = A000005(k)). - Amiram Eldar, Dec 19 2018
The term "exponentially odd integers" was apparently coined by Cohen (1960). These numbers were also called "unitarily 2-free", or "2-skew", by Cohen (1961). - Amiram Eldar, Jan 22 2024
LINKS
Peter J. C. Moses, Table of n, a(n) for n = 1..2000
Eckford Cohen, Arithmetical functions associated with the unitary divisors of an integer, Mathematische Zeitschrift, Vol. 74 (1960), pp. 66-80.
Eckford Cohen, Some sets of integers related to the k-free integers, Acta Sci. Math. (Szeged), Vol. 22, No. 3-4 (1961), pp. 223-233.
Vladimir Shevelev, Exponentially S-numbers, arXiv:1510.05914 [math.NT], 2015.
Vladimir Shevelev, Set of all densities of exponentially S-numbers, arXiv preprint arXiv:1511.03860 [math.NT], 2015.
Vladimir Shevelev, S-exponential numbers, Acta Arithmetica, Vol. 175(2016), 385-395.
D. Suryanarayana and R. Sita Rama Chandra Rao, Distribution of unitarily k-free integers, Journal of the Australian Mathematical Society, Vol. 20 , No. 2 (1975), pp. 129-141.
FORMULA
Sum_{a(n)<=x} 1 = C*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 C = Product_{prime p} (1 - 1/p*(p + 1)) = 0.7044422009991... (A065463).
Sum_{n>=1} 1/a(n)^s = zeta(2*s) * Product_{p prime} (1 + 1/p^s - 1/p^(2*s)), s>1. - Amiram Eldar, Sep 26 2023
MATHEMATICA
Select[Range@ 100, AllTrue[Last /@ FactorInteger@ #, OddQ] &] (* Version 10, or *)
Select[Range@ 100, Times @@ Boole[OddQ /@ Last /@ FactorInteger@ #] == 1 &] (* Michael De Vlieger, Feb 02 2016 *)
PROG
(PARI) isok(n)=my(f = factor(n)); for (k=1, #f~, if (!(f[k, 2] % 2), return (0))); 1; \\ Michel Marcus, Feb 02 2016
(Python)
from itertools import count, islice
from sympy import factorint
def A268335_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n:all(e&1 for e in factorint(n).values()), count(max(startvalue, 1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Feb 01 2016
STATUS
approved