%I #16 Dec 18 2023 01:42:15
%S 1,2,2,2,2,4,2,4,2,4,2,4,2,4,4,2,2,4,2,4,4,4,2,8,2,4,4,4,2,8,2,2,4,4,
%T 4,4,2,4,4,8,2,8,2,4,4,4,2,4,2,4,4,4,2,8,4,8,4,4,2,8,2,4,4,4,4,8,2,4,
%U 4,8,2,8,2,4,4,4,4,8,2,4,2,4,2,8,4,4,4
%N The number of tri-unitary divisors of n.
%C A divisor d of k is tri-unitary if the greatest common bi-unitary divisor of d and k/d is 1.
%C Differs from A037445 at n = 32, 96, 128, 160, 224, ...
%H Amiram Eldar, <a href="/A335385/b335385.txt">Table of n, a(n) for n = 1..10000</a>
%H Graeme L. Cohen, <a href="https://doi.org/10.1090/S0025-5718-1990-0993927-5">On an integer's infinitary divisors</a>, Mathematics of Computation, Vol. 54, No. 189 (1990), pp. 395-411.
%H Pentti Haukkanen, <a href="https://www.fq.math.ca/Scanned/38-5/haukkanen.pdf">On the k-ary convolution of arithmetical functions</a>, The Fibonacci Quarterly, Vol. 38, No. 5 (2000), pp. 440-445.
%F Multiplicative with a(p^e) = 4 if e = 3 or 6, and a(p^e) = 2 otherwise.
%e a(4) = 2 since 4 has 2 tri-unitary divisors, 1 and 4. 2 is not a tri-unitary divisor of 4 since the greatest common bi-unitary divisor of 2 and 4/2 = 2 is 2 and not 1.
%t f[p_, e_] := If[e == 3 || e == 6, 4, 2]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100]
%o (PARI) a(n) = vecprod(apply(x -> if(x == 3 || x == 6, 4, 2), factor(n)[, 2])); \\ _Amiram Eldar_, Dec 18 2023
%Y Cf. A222266, A324706, A324707, A324708, A324709.
%Y Cf. A000005, A034444, A037445, A048105, A049419, A286324, A322483.
%K nonn,easy,mult
%O 1,2
%A _Amiram Eldar_, Jun 04 2020