%I #21 May 30 2023 02:21:35
%S 942,945,1572,1575,2202,2205,2832,2835,3462,3465,4092,4095,4722,4725,
%T 5352,5355,5772,5985,6432,6435,6612,6615,6822,6825,7242,7245,7425,
%U 7872,7875,8082,8085,8412,8415,8502,8505,8922,8925,9132,9135,9552,9555,9762,9765
%N Abundant numbers that differ from the next abundant number by 3.
%C Apparently these numbers come up mostly by pairs m, m+3 with m even; the odd terms being a subsequence of A005231. But this is not always the case (e.g., note the term 7425).
%C The numbers of terms not exceeding 10^k, for k = 3, 4, ..., are 2, 43, 393, 3635, 37599, 374092, 3731903, 37338208, 373256850, ... . Apparently, the asymptotic density of this sequence exists and equals 0.00373... . - _Amiram Eldar_, May 30 2023
%H Muniru A Asiru, <a href="/A228382/b228382.txt">Table of n, a(n) for n = 1..10000</a>
%e 942 is abundant, 943 and 944 are deficient, 945 is abundant.
%e 945 is abundant, 946 and 947 are deficient, 948 is abundant.
%p with(numtheory): select(n->sigma(n)>2*n and sigma(n+1)<2*(n+1) and sigma(n+2)<2*(n+2) and sigma(n+3)>2*(n+3),[$1..12000]); # _Muniru A Asiru_, Jun 09 2018
%t With[{abs = Select[Range[10000], DivisorSigma[-1, #] > 2 &]}, abs[[Position[Differences[abs], 3] // Flatten]]] (* _Amiram Eldar_, May 30 2023 *)
%o (PARI) isok(n) = (sigma(n)> 2*n) && (sigma(n+1)< 2*(n+1)) && (sigma(n+2) < 2*(n+2)) && (sigma(n+3) > 2*(n+3)); \\ _Michel Marcus_, Aug 21 2013
%o (GAP) a:=Filtered([1..130000],n->Sigma(n)>2*n and Sigma(n+1)<2*(n+1) and Sigma(n+2)<2*(n+2) and Sigma(n+3)>2*(n+3)); # _Muniru A Asiru_, Jun 09 2018
%Y Subsequence of A005101.
%Y Cf. A005231, A096399.
%K nonn
%O 1,1
%A _Michel Marcus_, Aug 21 2013
|