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

A346880 Sum of the divisors, except the smallest and the largest, of the n-th positive even number.

Original entry on oeis.org

0, 2, 5, 6, 7, 15, 9, 14, 20, 21, 13, 35, 15, 27, 41, 30, 19, 54, 21, 49, 53, 39, 25, 75, 42, 45, 65, 63, 31, 107, 33, 62, 77, 57, 73, 122, 39, 63, 89, 105, 43, 139, 45, 91, 143, 75, 49, 155, 72, 116, 113, 105, 55, 171, 105, 135, 125, 93, 61, 239, 63, 99, 185, 126, 121
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

a(n) has a symmetric representation.

Examples

			For n = 5 the 5th even number is 10 and the divisors of 10 are [1, 2, 5, 10] and the sum of the divisors of 10 except the smaller and the largest is 2 + 5 = 7, so a(5) = 7.

Crossrefs

Bisection of A048050.
Partial sums give A346870.
Cf. A000203, A005843, A062731, A237593, A245092, A244049, A326124, A346879.

Programs

  • Mathematica
    a[n_] := DivisorSigma[1, 2*n] - 2*n - 1; Array[a, 100] (* Amiram Eldar, Aug 19 2021 *)
  • Python
    from sympy import divisors
def a(n): return sum(divisors(2*n)[1:-1])
print([a(n) for n in range(1, 66)]) # Michael S. Branicky, Aug 19 2021

Formula

a(n) = A048050(2*n).
Sum_{k=1..n} a(k) = (5*Pi^2/24 - 1) * n^2 + O(n*log(n)). - Amiram Eldar, Mar 21 2024

A346869 Sum of all divisors, except the smallest and the largest of every number, of the first n odd numbers.

Original entry on oeis.org

0, 0, 0, 0, 3, 3, 3, 11, 11, 11, 21, 21, 26, 38, 38, 38, 52, 64, 64, 80, 80, 80, 112, 112, 119, 139, 139, 155, 177, 177, 177, 217, 235, 235, 261, 261, 261, 309, 327, 327, 366, 366, 388, 420, 420, 440, 474, 498, 498, 554, 554, 554, 640, 640, 640, 680, 680, 708, 772, 796
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

Partial sums of the odd-indexed terms of Chowla's function A048050.
a(n) has a symmetric representation.

Crossrefs

Partial sums of A346879.
Cf. A000203, A005408, A008438, A048050, A237593, A245092, A244049, A326123, A346870.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 0,
      a(n-1)+numtheory[sigma](2*n-1)-2*n)
    end:
seq(a(n), n=1..60);  # Alois P. Heinz, Aug 19 2021
  • Mathematica
    s[1] = 0; s[n_] := DivisorSigma[1, 2*n - 1] - 2*n; Accumulate @ Array[s, 50] (* Amiram Eldar, Aug 19 2021 *)
Accumulate[Join[{0},Table[DivisorSigma[1,n]-n-1,{n,3,151,2}]]] (* Harvey P. Dale, Jul 29 2023 *)
  • Python
    from sympy import divisors
from itertools import accumulate
def A346879(n): return sum(divisors(2*n-1)[1:-1])
def aupton(nn): return list(accumulate(A346879(n) for n in range(1, nn+1)))
print(aupton(60)) # Michael S. Branicky, Aug 19 2021

A346877 Sum of the divisors, except for the largest, of the n-th odd number.

Original entry on oeis.org

0, 1, 1, 1, 4, 1, 1, 9, 1, 1, 11, 1, 6, 13, 1, 1, 15, 13, 1, 17, 1, 1, 33, 1, 8, 21, 1, 17, 23, 1, 1, 41, 19, 1, 27, 1, 1, 49, 19, 1, 40, 1, 23, 33, 1, 21, 35, 25, 1, 57, 1, 1, 87, 1, 1, 41, 1, 29, 65, 25, 12, 45, 31, 1, 47, 1, 27, 105, 1, 1, 51, 25, 35, 81, 1, 1, 81, 37
Offset: 1

Views

Author

Omar E. Pol, Aug 20 2021

Keywords

Comments

Sum of aliquot divisors (or aliquot parts) of the n-th odd number.
a(n) has a symmetric representation.

Examples

			For n = 5 the 5th odd number is 9 and the divisors of 9 are [1, 3, 9] and the sum of the divisors of 9 except for the largest is 1 + 3 = 4, so a(5) = 4.

Crossrefs

Bisection of A001065.
Partial sums give A347153.
Cf. A000203, A005408, A008438, A057427, A153485, A237593, A245092, A244049, A326123, A346869, A346878, A346879, A347154.

Programs

  • Mathematica
    a[n_] := DivisorSigma[1, 2*n - 1] - 2*n + 1; Array[a, 100] (* Amiram Eldar, Aug 20 2021 *)
Total[Most[Divisors[#]]]&/@Range[1,161,2] (* Harvey P. Dale, Sep 29 2024 *)
  • PARI
    a(n) = sigma(2*n-1) - (2*n-1); \\ Michel Marcus, Aug 20 2021
  • Python
    from sympy import divisors
def a(n): return sum(divisors(2*n-1)[:-1])
print([a(n) for n in range(1, 79)]) # Michael S. Branicky, Aug 20 2021

Formula

a(n) = A001065(2*n-1).
a(n) = A057427(n-1) + A346879(n).
G.f.: Sum_{k>=0} (2*k + 1) * x^(3*k + 2) / (1 - x^(2*k + 1)). - Ilya Gutkovskiy, Aug 20 2021
Sum_{k=1..n} a(k) = (Pi^2/8 - 1)*n^2 + O(n*log(n)). - Amiram Eldar, Mar 17 2024

A347153 Sum of all divisors, except the largest of every number, of the first n odd numbers.

Original entry on oeis.org

0, 1, 2, 3, 7, 8, 9, 18, 19, 20, 31, 32, 38, 51, 52, 53, 68, 81, 82, 99, 100, 101, 134, 135, 143, 164, 165, 182, 205, 206, 207, 248, 267, 268, 295, 296, 297, 346, 365, 366, 406, 407, 430, 463, 464, 485, 520, 545, 546, 603, 604, 605, 692, 693, 694, 735, 736, 765, 830, 855
Offset: 1

Views

Author

Omar E. Pol, Aug 20 2021

Keywords

Comments

Sum of all aliquot divisors (or aliquot parts) of the first n odd numbers.
Partial sums of the odd-indexed terms of A001065.
a(n) has a symmetric representation.

Crossrefs

Partial sums of A346877.
Cf. A000203, A001065, A001477, A005408, A008438, A048050, A153485, A237593, A245092, A244049, A326123, A346869, A346878, A346879, A347154.

Programs

  • Mathematica
    s[n_] := DivisorSigma[1, 2*n - 1] - 2*n + 1; Accumulate @ Array[s, 100] (* Amiram Eldar, Aug 20 2021 *)
  • PARI
    a(n) = sum(k=1, n, k = 2*k-1; sigma(k)-k); \\ Michel Marcus, Aug 20 2021
  • Python
    from sympy import divisors
from itertools import accumulate
def A346877(n): return sum(divisors(2*n-1)[:-1])
def aupton(nn): return list(accumulate(A346877(n) for n in range(1, nn+1)))
print(aupton(60)) # Michael S. Branicky, Aug 20 2021

Formula

a(n) = A001477(n-1) + A346869(n).
G.f.: (1/(1 - x)) * Sum_{k>=0} (2*k + 1) * x^(3*k + 2) / (1 - x^(2*k + 1)). - Ilya Gutkovskiy, Aug 20 2021
a(n) = (Pi^2/8 - 1)*n^2 + O(n*log(n)). - Amiram Eldar, Mar 21 2024
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.