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 A113645 #23 Sep 17 2020 09:01:22 %S A113645 4,8,12,16,18,24,27,32,36,48,54,64,72,80,81,96,108,120,128,144,160, %T A113645 162,180,192,200,216,240,243,256,270,288,300,320,324,360,384,400,405, %U A113645 432,448,450,480,486,500,512,540,576,600,640,648,672,675,720,729,750,768 %N A113645 Numbers k such that sum of exponents in prime factorization of k (i.e., A001222(k)) is >= each prime divisor of k. %H A113645 Amiram Eldar, <a href="/A113645/b113645.txt">Table of n, a(n) for n = 1..10000</a> %F A113645 A number k is included if A006530(k) <= A001222(k). %e A113645 12 = 2^2 *3^1. Since the sum of the prime factorization exponents, 2+1 = 3, is >= the largest prime dividing 12, which is 3, then 12 is included in the sequence. %t A113645 fQ[n_] := Block[{f = FactorInteger@n}, Plus @@ Last /@ f >= Max[First /@ f]]; Select[ Range[2, 800], fQ@ # &] (* _Robert G. Wilson v_, Jan 16 2006 *) %t A113645 qu[n_]:=n>1&&Block[{f=Transpose@FactorInteger@n, s}, s=Plus@@f[[2]];s>=Max@f[[1]]]; L ={};Do[If[qu[n], Print[n];AppendTo[L, n]], {n, 1000}];L (* _Giovanni Resta_, Jan 16 2006 *) %o A113645 (PARI) isok(m) = {my(f=factor(m)); #select(x->(x>bigomega(f)), f[,1]~) == 0;} \\ _Michel Marcus_, Sep 17 2020 %Y A113645 A proper subset of A068936. %Y A113645 Cf. A001222, A006530. %K A113645 nonn %O A113645 1,1 %A A113645 _Leroy Quet_, Jan 15 2006 %E A113645 More terms from _Robert G. Wilson v_ and _Giovanni Resta_, Jan 16 2006