cp's OEIS Frontend

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.

Showing 1-3 of 3 results.

A294898 Deficiency minus binary weight: a(n) = A033879(n) - A000120(n) = A005187(n) - A000203(n).

Original entry on oeis.org

0, 0, 0, 0, 2, -2, 3, 0, 3, 0, 7, -6, 9, 1, 2, 0, 14, -5, 15, -4, 7, 5, 18, -14, 16, 7, 10, -3, 24, -16, 25, 0, 16, 12, 19, -21, 33, 13, 18, -12, 37, -15, 38, 1, 8, 16, 41, -30, 38, 4, 26, 3, 48, -16, 33, -11, 30, 22, 53, -52, 55, 23, 16, 0, 44, -14, 63, 8, 39, -7, 66, -53, 69, 31, 22, 9, 54, -16, 73, -28, 38, 35, 78, -59, 58
Offset: 1

Views

Author

Antti Karttunen, Nov 25 2017

Keywords

Comments

"Least deficient numbers" or "almost perfect numbers" are those k for which A033879(k) = 1, or equally, for which a(k) = -A048881(k-1). The only known solutions are powers of 2 (A000079), all present also in A295296. See also A235796 and A378988. - Antti Karttunen, Dec 16 2024

Crossrefs

Cf. A000120, A000203, A001065, A005187, A011371, A013661, A033879, A048881, A235796, A294896, A294899, A297114 (Möbius transform), A317844 (difference from a(n)), A326133, A326138, A324348 (a(n) applied to Doudna sequence), A379008 (a(n) applied to prime shift array), A378988.
Cf. A295296 (positions of zeros), A295297 (parity of a(n)).

Programs

Formula

a(n) = A005187(n) - A000203(n).
a(n) = A011371(n) - A001065(n).
a(n) = A033879(n) - A000120(n).
Sum_{k=1..n} a(k) ~ c * n^2, where c = 1 - zeta(2)/2 = 0.177532... . - Amiram Eldar, Feb 22 2024

Extensions

Name edited by Antti Karttunen, Dec 16 2024

A326131 Positive numbers n for which A000120(n) = k*A294898(n), with k < 0; numbers for which A326130(n) = sigma(n) - A005187(n).

Original entry on oeis.org

6, 28, 110, 496, 884, 8128, 18632, 85936, 116624, 15370304, 33550336, 73995392, 815634435, 3915380170, 5556840416, 6800695312, 8589869056
Offset: 1

Views

Author

Antti Karttunen, Jun 09 2019

Keywords

Comments

No further terms below 2^31.
See also comments in A326133.
The quotients A000120(k)/(sigma(k)-A005187(k)) for these terms are: 1, 1, 5, 1, 3, 1, 5, 9, 2, 2, 1, 2, 2. Ones occur at the positions of perfect numbers.
a(18) > 10^11. - Amiram Eldar, Jan 03 2021

Examples

			110 is "1101110" in binary, thus A000120(110) = 5. Sigma(110) = 216, while A005187(110) = 215, thus as 5 = 5*(216-215), 110 is included in this sequence.
		

Crossrefs

Intersection of A326132 and A326133, also of A326132 and A326138.
Cf. also A325981, A326141.

Programs

  • Mathematica
    q[n_] := Module[{bw = DigitCount[n, 2, 1], ab = DivisorSigma[1, n] - 2*n, sum}, (sum = ab + bw) > 0 && Divisible[bw, sum]]; Select[Range[10^5], q] (* Amiram Eldar, Jan 03 2021 *)
  • PARI
    A005187(n) = { my(s=n); while(n>>=1, s+=n); s; };
    isA326131(n) = { my(t=sigma(n)-A005187(n)); (gcd(hammingweight(n), t) == t); };

Extensions

a(14)-a(17) from Amiram Eldar, Jan 03 2021

A326133 Numbers n for which sigma(n) > A005187(n).

Original entry on oeis.org

6, 12, 18, 20, 24, 28, 30, 36, 40, 42, 48, 54, 56, 60, 66, 70, 72, 78, 80, 84, 88, 90, 96, 100, 102, 104, 108, 110, 112, 114, 120, 126, 132, 138, 140, 144, 150, 156, 160, 162, 168, 174, 176, 180, 186, 192, 196, 198, 200, 204, 208, 210, 216, 220, 222, 224, 228, 234, 240, 246, 252, 258, 260, 264, 270, 272, 276, 280, 282, 288, 294, 300
Offset: 1

Views

Author

Antti Karttunen, Jun 11 2019

Keywords

Comments

Differs from A023196 for the first time at the 28th term, which here is 110, which is not included in A023196.
Note that as there is at least one odd number (815634435) in A326138, it means that A005231 does not contain all odd terms of this sequence.

Crossrefs

Positions of negative terms in A294898.
Cf. A000396, A005231, A083207, A111592, A326131, A326138 (subsequences).

Programs

  • Mathematica
    Select[Range[300], DivisorSigma[1, #] > 2*# - DigitCount[2*#, 2, 1] &] (* Amiram Eldar, Aug 06 2023 *)
  • PARI
    A005187(n) = { my(s=n); while(n>>=1, s+=n); s; };
    isA326133(n) = (sigma(n)>A005187(n));
Showing 1-3 of 3 results.