|
%I #32 Mar 06 2023 05:51:02
%S 1,1,1,2,1,1,1,2,2,1,1,2,1,1,1,2,1,2,1,2,1,1,1,2,2,1,2,2,1,1,1,2,1,1,
%T 1,4,1,1,1,2,1,1,1,2,2,1,1,2,2,2,1,2,1,2,1,2,1,1,1,2,1,1,2,4,1,1,1,2,
%U 1,1,1,4,1,1,2,2,1,1,1,2,2,1,1,2,1,1
%N Multiplicative with a(p^e) = 2^omega(e), where omega = A001221.
%C The number of exponential unitary (or e-unitary) divisors of n and the number of exponential squarefree exponential divisors (or e-squarefree e-divisors) of n. These are divisors of n = Product p(i)^a(i) of the form Product p(i)^b(i) where each b(i) is a unitary divisor of a(i) in the first case, or each b(i) is a squarefree divisor of a(i) in the second case. - _Amiram Eldar_, Dec 29 2018
%H Antti Karttunen, <a href="/A278908/b278908.txt">Table of n, a(n) for n = 1..10000</a>
%H Xiaodong Cao and Wenguang Zahi, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL13/Cao/cao4.html">Some arithmetic functions involving exponential divisors</a>, Journal of Integer Sequences, Vol. 13 (2010), Article 10.3.7, eq (20).
%H Nicusor Minculete and László Tóth, <a href="http://ac.inf.elte.hu/Vol_035_2011/elte_annales_35_jav_vagott.pdf#page=205">Exponential unitary divisors</a>, Annales Univ. Sci. Budapest., Sect. Comp., Vol. 35 (2011), pp. 205-216.
%H László Tóth, <a href="http://ac.inf.elte.hu/Vol_027_2007/155.pdf">On certain arithmetic functions involving exponential divisors, II</a>, Annales Univ. Sci. Budapest., Sect. Comp., Vol. 27 (2007), pp. 155-166; <a href="https://arxiv.org/abs/0708.3557">arXiv preprint</a>, arXiv:0708.3557 [math.NT], 2007-2009.
%H Xiangzhen Zhao, Min Liu, and Yu Huang, <a href="http://fs.unm.edu/ScientiaMagna8no3.pdf#page=116">Mean value for the function t^(e)(n) over square-full numbers</a>, Scientia Magna, Vol. 8, No. 3 (2012), pp. 110-114.
%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>
%F Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = Product_{p prime} (1 + Sum_{k>=2} (2*omega(k) - 2^omega(k-1))/p^k) = 1.5431653193... (Tóth, 2007). - _Amiram Eldar_, Nov 08 2020
%p A278908 := proc(n)
%p local a,p,e;
%p a := 1;
%p if n =1 then
%p ;
%p else
%p for p in ifactors(n)[2] do
%p e := op(2,p) ;
%p a := a*2^A001221(e) ;
%p end do:
%p end if;
%p a ;
%p end proc:
%t Table[Times @@ Apply[Times, FactorInteger[n] /. {p_, e_} /; p > 1 :> 2^PrimeNu[e]], {n, 105}] (* _Michael De Vlieger_, Jul 29 2017 *)
%o (Scheme) (define (A278908 n) (if (= 1 n) n (* (A000079 (A001221 (A067029 n))) (A278908 (A028234 n))))) ;; _Antti Karttunen_, Jul 27 2017
%o (PARI) a(n) = my(f=factor(n)); for (k=1, #f~, f[k,1] = 2^omega(f[k,2]); f[k,2] = 1); factorback(f); \\ _Michel Marcus_, Jul 28 2017
%Y Cf. A001221.
%K nonn,easy,mult
%O 1,4
%A _R. J. Mathar_, Nov 30 2016
|
