login
Unitary abundant numbers: numbers k such that usigma(k) > 2*k.
42

%I #39 Mar 12 2026 06:56:26

%S 30,42,66,70,78,102,114,138,150,174,186,210,222,246,258,282,294,318,

%T 330,354,366,390,402,420,426,438,462,474,498,510,534,546,570,582,606,

%U 618,630,642,654,660,678,690,714,726,750,762,770,780,786,798,822,834

%N Unitary abundant numbers: numbers k such that usigma(k) > 2*k.

%C If a term n in the sequence ends in neither 0 nor 5, then 10*n is also in the sequence. - _Lekraj Beedassy_, Jun 11 2004

%C The lower asymptotic density of this sequence is larger than 1/18 = 0.0555... which is the density of its subsequence of numbers of the form 6*m where gcd(m, 6) = 1 and m > 1. Numerically, based on counts of terms below 10^n (A302993), it seems that this sequence has an asymptotic density which equals to about 0.070034... - _Amiram Eldar_, Feb 13 2021

%C The asymptotic density of this sequence is in the interval (0.0674, 0.1055) (Wall, 1970). - _Amiram Eldar_, Apr 18 2024

%C All the terms are nonpowerful numbers (A052485). For powerful numbers (A001694) k, usigma(k)/k < 15/Pi^2 = 1.519817... (A082020; the record values are attained at the squares of primorials, A061742). - _Amiram Eldar_, Jul 20 2024

%D C. Sung, Mathematical Buds, "Unitary Divisors", Chap. V, pp. 42-67, Ed. H. D. Ruderman, Mu Alpha Theta OK 1978.

%H Nathaniel Johnston, <a href="/A034683/b034683.txt">Table of n, a(n) for n = 1..10000</a>

%H Charles Robert Wall, <a href="https://trace.tennessee.edu/utk_graddiss/8570/">Topics related to the sum of unitary divisors of an integer</a>, Ph.D. diss., University of Tennessee, 1970.

%p isA034683 := proc(n)

%p is(A034448(n) > 2*n) ;

%p end proc:

%p for n from 1 do

%p if isA034683(n) then

%p print(n);

%p end if;

%p end do: # _R. J. Mathar_, Nov 10 2014

%t usigma[n_] := If[n == 1, 1, Times @@ (1 + Power @@@ FactorInteger[n])];

%t Select[Range[1000], usigma[#] > 2#&] (* _Jean-François Alcover_, Mar 23 2020, after _Giovanni Resta_ in A034448 *)

%o (PARI) is(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + f[i, 1]^f[i, 2]) > 2*n;} \\ _Amiram Eldar_, Apr 18 2024

%Y Subsequence of A005101.

%Y Cf. A034444, A034448, A129485 (odd terms), A129487.

%Y See the comments for the relationships with A001694, A052485, A061742, A082020, A302993,

%K nonn

%O 1,1

%A _Erich Friedman_