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.
%I A349181 #7 Nov 10 2021 07:07:40 %S A349181 100,300,700,1100,1225,1300,1700,1900,2100,2300,2450,2900,3100,3300, %T A349181 3675,3700,3900,4100,4225,4300,4700,5100,5300,5700,5900,6100,6700, %U A349181 6900,7100,7300,7350,7700,7900,8300,8450,8700,8900,9100,9300,9700,10100,10300,10700,10900 %N A349181 Powerful harmonic numbers: numbers k such that the set of powerful divisors of k that are larger than 1 has more than one element and that the harmonic mean of this set is an integer. %C A349181 Numbers with a single powerful divisor > 1 are A060687 and trivially have an integer harmonic mean. %C A349181 The least term that is not divisible by 5 (or 25) is a(5446) = 1413721. %H A349181 Amiram Eldar, <a href="/A349181/b349181.txt">Table of n, a(n) for n = 1..10000</a> %e A349181 100 is a term since its powerful divisors > 1 are 4, 25 and 100 and their harmonic mean, 10, is an integer. %t A349181 powQ[n_] := Min[FactorInteger[n][[;; , 2]]] > 1; powHarmQ[n_] := Module[{d = Select[Divisors[n], powQ]}, Length[d] > 1 && IntegerQ[HarmonicMean[d]]]; Select[Range[10^4], powHarmQ] %Y A349181 Cf. A001694, A060687. %Y A349181 Similar sequences: A001599, A006086, A063947, A286325, A319745, A348964. %K A349181 nonn %O A349181 1,1 %A A349181 _Amiram Eldar_, Nov 09 2021