%I #37 Oct 11 2023 10:17:48
%S 30,60,70,84,90,105,120,126,132,140,154,165,168,180,182,195,198,210,
%T 220,231,234,252,260,264,273,280,286,306,308,312,315,330,336,340,357,
%U 360,364,374,380,385,390,396,399,408,418,420,429,440,442,455,456,462
%N Numbers k such that if one writes k = Product p_i^e_i (p_i primes) and P = max p_i^e_i, then k/P > P.
%C Numbers k such that A346418(k) > 1 (conjectured). - _Amiram Eldar_, Jul 16 2021
%C Numbers k for which A354933(k) > A034699(k), i.e., where A354995(k) > 0 (conjectured).- _Antti Karttunen_, Jun 19 2022
%H Reinhard Zumkeller, <a href="/A051283/b051283.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A080170(n) + 1 (conjectured). - _Ralf Stephan_, Feb 20 2004
%e 120 = 2^3*3^1*5^1, P = 2^3 = 8. 120 is included because 120/8 = 15 > 8.
%t ok[n_] := n > Max[Power @@@ FactorInteger[n]]^2; Select[Range[465], ok] (* _Jean-François Alcover_, Apr 11 2011 *)
%o (Haskell)
%o a051283 n = a051283_list !! (n-1)
%o a051283_list = filter (\x -> (a034699 x) ^ 2 < x) [1..]
%o -- _Reinhard Zumkeller_, May 30 2013
%o (PARI) lista(nn) = for(k=2, nn, f=factor(k); if(k>vecmax(vector(#f[, 1], i, f[i, 1]^f[i, 2]))^2, print1(k, ", "))); \\ _Jinyuan Wang_, Feb 28 2020
%Y Cf. A034699, A080170, A110612, A346418, A354933, A354995, A354996 (characteristic function).
%K nonn,easy,nice
%O 1,1
%A _Leroy Quet_
%E More terms from _James A. Sellers_, Dec 11 1999
|