%I #11 Jun 16 2021 09:27:04
%S 1,1,2,3,4,2,5,6,7,4,8,6,9,5,8,10,11,7,12,12,10,8,13,12,14,9,15,15,16,
%T 8,17,18,16,11,20,21,19,12,18,24,20,10,21,24,28,13,22,20,23,14,22,27,
%U 24,15,32,30,24,16,25,24,26,17,35,27,36,16,28,33,26,20
%N a(n) is the product of indices of unitary prime power divisors of n.
%C Equivalently: replace each prime power p^e in the prime factorisation of n by its index in A246655. - _M. F. Hasler_, Jun 16 2021
%H Robert Israel, <a href="/A333235/b333235.txt">Table of n, a(n) for n = 1..10000</a>
%F If n = Product (p_j^k_j) then a(n) = Product (A025528(p_j^k_j)).
%F a(prime(n)) = A027883(n).
%F a(2^n) = A182908(n).
%F a(A246655(n)) = n.
%e a(600) = a(2^3 * 3 * 5^2) = a(A246655(6) * A246655(2) * A246655(14)) = 6 * 2 * 14 = 168.
%p N:= 1000: # for a(1)..a(N)
%p R:= NULL: p:= 2:
%p while p < N do
%p R:= R, seq(p^k,k=1..ilog[p](N));
%p p:= nextprime(p);
%p od:
%p L:= sort([R]):
%p f:= proc(n) local F, t;
%p F:= ifactors(n)[2];
%p mul(ListTools:-BinarySearch(L,t[1]^t[2]),t=F)
%p end proc:
%p map(f, [$1..N]); # _Robert Israel_, Feb 11 2021
%t PrimePowerPi[n_] := Sum[Boole[PrimePowerQ[k]], {k, 1, n}]; a[1] = 1; a[n_] := Times @@ (PrimePowerPi[#[[1]]^#[[2]]] & /@ FactorInteger[n]); Table[a[n], {n, 1, 70}]
%o (PARI) apply( {A333235(n)=vecprod([A322981(f[1]^f[2])|f<-factor(n)~])}, [1..99]) \\ _M. F. Hasler_, Jun 16 2021
%Y Cf. A003963, A025528, A027883, A141128, A141809, A156061, A182908, A246655.
%Y Cf. A322981 (the index of n = p^e in A246655).
%K nonn,mult,look
%O 1,3
%A _Ilya Gutkovskiy_, Mar 12 2020
|