cp's OEIS Frontend

A377591 Powerful numbers k that are not prime powers such that there exist no numbers m such that rad(m) | k and Omega(m) > Omega(k), where rad = A007947 and Omega = A001222.

%I A377591 #6 Nov 13 2024 17:17:22
%S A377591 225,675,1225,2025,3025,5929,6075,6125,8281,8575,14161,15125,18225,
%T A377591 20449,30625,34969,41503,42875,43681,48841,54675,57967,60025,61009,
%U A377591 64009,65219,75625,89401,99127,101761,104329,107653,116281,142129,152881,153125,162409,164025
%N A377591 Powerful numbers k that are not prime powers such that there exist no numbers m such that rad(m) | k and Omega(m) > Omega(k), where rad = A007947 and Omega = A001222.
%C A377591 Terms are odd; proper subset of A363217, which is a proper subset of A286708, itself contained in A001694.
%C A377591 Proper subset of A377590.
%H A377591 Michael De Vlieger, <a href="/A377591/b377591.txt">Table of n, a(n) for n = 1..10000</a>
%H A377591 Michael De Vlieger, <a href="https://doi.org/10.13140/RG.2.2.18645.64480">Numbers k for which floor(log k / log lpf(k)) <= bigomega(k)</a>, 2024.
%e A377591 36 is not in the sequence since 2^5 < 36, Omega(32) = 5, but Omega(36) = 4.
%e A377591 72 is not in the sequence since 2^6 < 72, but Omega(72) = 5.
%e A377591 225 is in the sequence since 3^4 < 225, Omega(81) = Omega(225) = 4.
%e A377591 441 is not in the sequence since 3^5 < 441, Omega(243) = 5, but Omega(441) = 4, etc.
%t A377591 Select[With[{nn = 200000},
%t A377591   Rest@ Select[
%t A377591     Union@ Flatten@ Table[a^2*b^3, {b, Surd[nn, 3]}, {a, Sqrt[nn/b^3]}],
%t A377591     Not@*PrimePowerQ] ],
%t A377591   Function[{n, k},
%t A377591     NoneTrue[FactorInteger[n][[All, 1]],
%t A377591       Floor@ Log[#, n] > k &]] @@ {#, PrimeOmega[#]} &]
%Y A377591 Cf. A001222, A001694, A007947, A286708, A363217, A376846, A377590.
%K A377591 nonn,easy
%O A377591 1,1
%A A377591 _Michael De Vlieger_, Nov 02 2024