Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #14 Feb 14 2020 03:21:44
%S 1,5,12,50,60,84,300,420,450,756,900,1950,3780,7800,9900,33150,49140,
%T 54600,100800,132600,265200,491400,928200,1856400,8353800,8884200,
%U 16707600,52211250,65995776,78566400,182739375,183783600,208845000,280348992,293046000,329978880
%N Numbers k that are unitary harmonic in Gaussian integers: k * A332476(k) is divisible by A332472(k) + i*A332473(k) (where i is the imaginary unit).
%C Analogous to unitary harmonic numbers (A006086), with the number and sum of unitary divisors functions generalized for Gaussian integers (A332476, A332472 + i * A332473) instead of the number and sum of unitary divisors functions (A034444, A034448).
%H Amiram Eldar, <a href="/A332477/b332477.txt">Table of n, a(n) for n = 1..75</a>
%e 5 is a term since 5 * A332476(5)/(A332472(5) + i*A332473(5)) = 5 * 4/(4 + 8*i) = 1 - 2*i is a Gaussian integer.
%t sigma[p_, e_] := If[Abs[p] == 1, 1, (p^e + 1)]; tau[p_, e_] := If[Abs[p] == 1, 1, 2]; unitaryHarmoincQ[n_] := Divisible[n * Times @@ tau @@@ (f = FactorInteger[n, GaussianIntegers -> True]), Times @@ sigma @@@ f]; Select[Range[10^6], unitaryHarmoincQ]
%Y Cf. A034444, A034448, A006086, A332317, A332472, A332473, A332476.
%K nonn
%O 1,2
%A _Amiram Eldar_, Feb 13 2020