cp's OEIS Frontend

A295920 Number of twice-factorizations of n of type (P,R,R).

%I A295920 #20 Dec 06 2018 16:33:00
%S A295920 1,1,1,3,1,1,1,3,3,1,1,1,1,1,1,8,1,1,1,1,1,1,1,1,3,1,3,1,1,1,1,3,1,1,
%T A295920 1,3,1,1,1,1,1,1,1,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,17,1,1,1,1,
%U A295920 1,1,1,1,1,1,1,1,1,1,1,1,8,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,1,1,1,1,1
%N A295920 Number of twice-factorizations of n of type (P,R,R).
%C A295920 a(n) is also the number of ways to choose a perfect divisor d|n and then a sequence of log_d(n) perfect divisors of d.
%H A295920 Antti Karttunen, <a href="/A295920/b295920.txt">Table of n, a(n) for n = 1..16384</a>
%H A295920 <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>
%F A295920 a(n) = Sum_{d|A052409(n)} A000005(A052409(n^(1/d)))^d. - _Antti Karttunen_, Dec 06 2018, after Mathematica-code
%e A295920 The a(64) = 17 twice-factorizations are:
%e A295920 (2)*(2)*(2)*(2)*(2)*(2)  (2*2)*(2*2)*(2*2)  (2*2*2)*(2*2*2)  (2*2*2*2*2*2)
%e A295920 (2*2)*(2*2)*(4)          (2*2)*(4)*(2*2)    (4)*(2*2)*(2*2)
%e A295920 (2*2)*(4)*(4)            (4)*(2*2)*(4)      (4)*(4)*(2*2)
%e A295920 (2*2*2)*(8)              (8)*(2*2*2)
%e A295920 (4)*(4)*(4)              (4*4*4)
%e A295920 (8)*(8)                  (8*8)
%e A295920 (64)
%t A295920 Table[Sum[Length[Divisors[GCD@@FactorInteger[n^(1/d)][[All,2]]]]^d,{d,Divisors[GCD@@FactorInteger[n][[All,2]]]}],{n,100}]
%o A295920 (PARI)
%o A295920 A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
%o A295920 A295920(n) = if(1==n,n,my(r); sumdiv(A052409(n), d, if(!ispower(n,d,&r),(1/0),numdiv(A052409(r))^d))); \\ _Antti Karttunen_, Dec 06 2018, after Mathematica-code
%Y A295920 Cf. A000005, A001055, A052409, A052410, A089723, A279789, A281113, A295923, A295924, A295931, A295935.
%K A295920 nonn
%O A295920 1,4
%A A295920 _Gus Wiseman_, Nov 30 2017
%E A295920 More terms from _Antti Karttunen_, Dec 06 2018