A364997 Numbers k neither squarefree nor prime power such that rad(k)*A119288(k) > k but rad(k)*A053669(k) < k.
40, 45, 56, 63, 88, 99, 104, 117, 136, 152, 153, 171, 175, 176, 184, 207, 208, 232, 248, 261, 272, 275, 279, 280, 296, 297, 304, 315, 325, 328, 333, 344, 351, 368, 369, 376, 387, 423, 424, 425, 440, 459, 464, 472, 475, 477, 488, 495, 496, 513, 520, 531, 536, 539
Offset: 1
Keywords
Examples
Let b(n) = A126706(n), S = A360767, and T = A360765. b(1) = 12 is not in the sequence since p*r = 3*6 = 18 and q*r = 5*6 = 30; both exceed 12, thus 12 is in S but not in T. b(2) = 18 is not in the sequence since p*r = 3*6 = 18 and q*r = 5*6 = 30. Indeed, neither is less than 18, hence 18 is not in S but is in T. b(6) = 36 is not in the sequence since p*r = 3*6 = 18 and q*r = 5*6, and both do not exceed 36, therefore 36 is not in S but is in T. b(7) = a(1) = 40 since p*r = 5*10 = 50 and q*r = 3*10 = 30. We have both 50 > 40 and 30 < 40, thus 40 is in both S and T, etc.
Links
- Michael De Vlieger, Table of n, a(n) for n = 1..10000
- Michael De Vlieger, Annotated plot of b(n) = A126706(n), n = 20*(y-1) + x at (x, -y), for x = 1..20 and y = 1..20, thus for n = 1..400. Terms in this sequence are colored black, those in A364999 in blue, in A364998 in gold, and in A361098 in red.
- Michael De Vlieger, Plot of b(n), n = 120*(y-1) + x at (x, -y), for x = 1..120 and y = 1..120, thus for n = 1..14400 using the same color scheme as immediately above.
- Michael De Vlieger, Plot of b(n), with n = 1016*(y-1) + x at (x, -y), for x = 1..1016 and y = 1..1016, thus showing 1032256 terms. Terms b(n) in this sequence are colored black, else white.
Crossrefs
Programs
-
Mathematica
Select[Select[Range[500], Nor[PrimePowerQ[#], SquareFreeQ[#]] &], Function[{k, f}, Function[{p, q, r}, And[p r > k, q r < k]] @@ {f[[2, 1]], SelectFirst[Prime@ Range[PrimePi[f[[-1, 1]]] + 1], ! Divisible[k, #] &], Times @@ f[[All, 1]]}] @@ {#, FactorInteger[#]} &]
Comments