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-2 of 2 results.

A374941 a(n) = (sum of proper prime factors of n) + Sum_{d = proper composite factor of n} a(d).

Original entry on oeis.org

0, 0, 0, 2, 0, 5, 0, 4, 3, 7, 0, 12, 0, 9, 8, 8, 0, 13, 0, 16, 10, 13, 0, 28, 5, 15, 6, 20, 0, 30, 0, 16, 14, 19, 12, 40, 0, 21, 16, 36, 0, 36, 0, 28, 19, 25, 0, 64, 7, 19, 20, 32, 0, 32, 16, 44, 22, 31, 0, 90, 0, 33, 23, 32, 18, 48, 0, 40, 26, 42, 0, 112, 0, 39
Offset: 1

Views

Author

Tanmaya Mohanty, Jul 24 2024

Keywords

Examples

			For n=24, the tree of recurrences "a(d)" is
         24
  /  /  /  \  \  \
  2  3  4  6   8  12
       /  / \  / \  \ \ \ \
      2  2  3  2  4  2 3 4 6
                 /      /  / \
                2      2   2  3
a(24) = 2 + 3 + a(4) + a(6) + a(8) + a(12)
= 5 + 2 + 2 + 3 + 2 + a(4) + 2 + 3 + a(4) + a(6)
= 14 + 2 + 5 + 2 + 2 + 3
= 28
		

Crossrefs

Cf. A371075 (fixed points).

Programs

  • Python
    from sympy import divisors, isprime
    def a(n): return sum(di if isprime(di) else a(di) for di in divisors(n)[1:-1])
    print([a(n) for n in range(1, 75)]) # Michael S. Branicky, Jul 24 2024

Formula

a(p) = 0, iff p = 1 or a prime number.
a(p^k) == 2^(k-2)*p for p prime, k > 1. - Michael S. Branicky, Aug 01 2024

Extensions

More terms from Michael S. Branicky, Jul 24 2024

A376445 Numbers k such that A374941(k) > k.

Original entry on oeis.org

24, 36, 48, 60, 72, 84, 90, 96, 108, 120, 132, 144, 160, 168, 180, 192, 210, 216, 240, 252, 264, 270, 280, 288, 300, 312, 320, 336, 360, 384, 396, 408, 420, 432, 456, 468, 480, 504, 528, 540, 552, 560, 576, 600, 612, 624, 630, 640, 648, 660, 672, 684, 696, 720
Offset: 1

Views

Author

Tanmaya Mohanty, Sep 22 2024

Keywords

Examples

			24 is a term because A374941(24) = 28 > 24.
		

Crossrefs

Programs

  • Python
    from sympy import divisors, isprime
    from functools import cache
    @cache
    def f(n): return sum(di if isprime(di) else f(di) for di in divisors(n)[1:-1])
    def ok(n): return f(n) > n
    print([k for k in range(1, 801) if ok(k)]) # Michael S. Branicky, Sep 23 2024
Showing 1-2 of 2 results.