|
%I #29 Oct 31 2020 11:50:52
%S 3,6,12,18,24,36,48,72,96,108
%N Ramanujan's largely composite numbers having 3 as the greatest prime divisor.
%C Theorem. Ramanujan's largely composite numbers (A067128) having the greatest prime divisor p_k = prime(k) do not exceed Product_{2 <= p <= p_k} p^((2*ceiling(log_p(p_(k + 1)) - 1).
%C Proof. Let N be in A067128 with prime power factorization 2^l_1 * 3^l_2 * ... * p_k^l_k.
%C First let us show that l_1 <= 2x_1-1 such that 2^x_1 > p_(k+1).
%C Indeed, consider N_1 = 2^(l_1-x_1)*3^l_2*...*p_k^l_k*p_(k+1).
%C Since 2^x_1 > p_(k+1) then N_1<N.
%C But d(N_1) > d(N) if l_1 >= 2*x_1, so l_1 <= 2x_1-1.
%C Analogously we find l_i <= 2x_i-1 if p_i^x_i > p_(k+1), i <= k.
%C Therefore N <= 2^(2*x_1-1)*3^(2*x_2-1)*...* p_k^(2*x_k-1) and the theorem easily follows.
%C QED
%C The inequality of the theorem gives a way to find the full sequence for every p_k. In particular, in case p_k = 2 we have the sequence {2, 4, 8}. For other cases see A273016, A273018.
%t a = {}; b = {0}; Do[If[# >= Max@ b, AppendTo[a, k] && AppendTo[b, #]] &@ DivisorSigma[0, k], {k, 10^7}]; Select[a, FactorInteger[#][[-1, 1]] == 3 &] (* _Michael De Vlieger_, May 13 2016 *)
%Y Cf. A067128, A065119 (the intersection of these two sequences is the present sequence). Cf. also A003586, A273016, A273018.
%K nonn,fini,full
%O 1,1
%A _Vladimir Shevelev_, May 13 2016
| 