This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A156903 #91 May 08 2025 08:56:16
%S A156903 18,36,72,100,144,162,196,200,288,324,392,400,450,576,648,784,800,882,
%T A156903 900,968,1152,1296,1352,1458,1568,1600,1764,1800,1936,2178,2304,2450,
%U A156903 2500,2592,2704,2916,3042,3136,3200,3528,3600,3872,4050,4356,4608,4624
%N A156903 Abundant numbers (A005101) whose abundance is odd.
%C A156903 Equivalently, abundant numbers with odd sum of divisors.
%C A156903 Complement of A204825 with respect to A005101 (abundant numbers).
%C A156903 Seems to be a proper subset of A083211. - _Robert G. Wilson v_, Mar 30 2010. This sequence is indeed a proper subset of A083211, since the abundance of a number k, A033880(k) = sigma(k) - 2*k, has the same parity as sigma(k). If sigma(k) is odd then the sums of any two complementary subsets of the divisors of k have different parities and thus they cannot be equal. - _Amiram Eldar_, Jun 20 2020
%C A156903 If n is present, so is 2*n. - _Robert G. Wilson v_, Jun 21 2015
%C A156903 If n is in the sequence, so is 100*n (conjectured). - _Sergey Pavlov_, Mar 22 2017. Pavlov's observation trivially follows from the fact that to have odd abundance a number k must be either a square or twice a square. If such a number k is abundant then 100*k = (10^2) * k is abundant as well and has odd abundance. In general, we can say that if k is present, so are t^2*k and 2*t^2*k, for every t>0. - _Giovanni Resta_, Oct 16 2018
%C A156903 Terms are congruent to {0, 2, 4, 8, 9, 14, 16, 18, 20, 26, 28, 32} (mod 36). - _Robert G. Wilson v_, Dec 09 2018
%H A156903 Robert G. Wilson v, <a href="/A156903/b156903.txt">Table of n, a(n) for n = 1..22927</a> (corrected by Michel Marcus)
%H A156903 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/AbundantNumber.html">Abundant Number</a>
%e A156903 k = 18 is in the sequence because its divisors are {1,2,3,6,9,18} which sum to sigma(k) = 39; so its abundance is sigma(k) - 2k = 39 - 36 = 3.
%p A156903 with(numtheory): select(k->sigma(k)>2*k and modp(sigma(k)-2*k,2)=1,[$1..5000]); # _Muniru A Asiru_, Dec 11 2018
%t A156903 abundance[n_] := DivisorSigma[1, n] - 2 n; Select[Range[1000], abundance[#] > 0 && Mod[abundance[#], 2] == 1 &]
%t A156903 abundOddAbundQ[n_] := If[MemberQ[{0, 2, 4, 8, 9, 14, 16, 18, 20, 26, 28, 32}, Mod[n, 36]], a = DivisorSigma[1, n]; OddQ@a && a > 2 n]; Select[ Range@ 5000, abundOddAbundQ@# &] (* _Robert G. Wilson v_, Dec 23 2018 *)
%o A156903 (PARI) is(n)=my(k=sigma(n)-2*n); k>0 && k%2 \\ _Charles R Greathouse IV_, Feb 21 2017
%o A156903 (Python)
%o A156903 from sympy.ntheory import divisor_sigma
%o A156903 def a(n):
%o A156903 return divisor_sigma(n) - 2*n
%o A156903 [n for n in range(18, 5001) if a(n) > 0 and a(n) % 2] # _Indranil Ghosh_, Mar 22 2017
%o A156903 (GAP) Filtered([1..5000],k->Sigma(k)-2*k>0 and (Sigma(k)-2*k) mod 2=1); # _Muniru A Asiru_, Dec 11 2018
%Y A156903 Intersection of A005101 and A028982. - _Amiram Eldar_, Jun 20 2020
%Y A156903 Cf. A000203, A033880, A259231.
%Y A156903 A proper subset of A083211. Setwise difference A083211 \ A171641.
%Y A156903 Positions of odd negative terms in A378600.
%Y A156903 Cf. A204825 (abundant numbers with even sum of divisors), A204826 (deficient numbers with odd sum of divisors), A204827 (deficient numbers with even sum of divisors).
%K A156903 nonn
%O A156903 1,1
%A A156903 _Robert G. Wilson v_, Feb 17 2009
%E A156903 Name edited by _Michel Marcus_ and _Charles R Greathouse IV_, Mar 26 2017
%E A156903 Edited by _N. J. A. Sloane_, Jun 21 2020 at the suggestion of _Amiram Eldar_