cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A335387 Tri-unitary harmonic numbers: numbers k such that the harmonic mean of the tri-unitary divisors of k is an integer.

This page as a plain text file.
%I A335387 #9 Jul 22 2020 02:20:50
%S A335387 1,6,45,60,90,270,420,630,2970,5460,8190,9100,15925,27300,36720,40950,
%T A335387 46494,47520,54600,81900,95550,136500,163800,172900,204750,232470,
%U A335387 245700,257040,332640,409500,464940,491400,646425,716625,790398,791700,819000,900900,929880
%N A335387 Tri-unitary harmonic numbers: numbers k such that the harmonic mean of the tri-unitary divisors of k is an integer.
%C A335387 Equivalently, numbers k such that A324706(k) | (k * A335385(k)).
%C A335387 Differs from A063947 from n >= 18.
%H A335387 Amiram Eldar, <a href="/A335387/b335387.txt">Table of n, a(n) for n = 1..200</a>
%e A335387 45 is a term since its tri-unitary divisors are {1, 5, 9, 45} and their harmonic mean, 3, in an integer.
%t A335387 f1[p_, e_] := If[e == 3 || e == 6, 4, 2]; f2[p_, e_] := If[e == 3, (p^4 - 1)/(p - 1), If[e == 6, (p^8 - 1)/(p^2 - 1), p^e + 1]]; f[p_, e_] := p^e * f1[p, e]/f2[p, e]; tuhQ[1] = True; tuhQ[n_] := IntegerQ[Times @@ (f @@@ FactorInteger[n])]; Select[Range[10^4], tuhQ]
%Y A335387 A324707 is a subsequence.
%Y A335387 Analogous sequences: A001599 (harmonic numbers), A006086 (unitary), A063947 (infinitary), A286325 (bi-unitary), A319745 (nonunitary).
%Y A335387 Cf. A324706, A335385.
%K A335387 nonn
%O A335387 1,2
%A A335387 _Amiram Eldar_, Jun 04 2020