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

A062731 Sum of divisors of 2*n.

Original entry on oeis.org

3, 7, 12, 15, 18, 28, 24, 31, 39, 42, 36, 60, 42, 56, 72, 63, 54, 91, 60, 90, 96, 84, 72, 124, 93, 98, 120, 120, 90, 168, 96, 127, 144, 126, 144, 195, 114, 140, 168, 186, 126, 224, 132, 180, 234, 168, 144, 252, 171, 217, 216, 210, 162, 280, 216, 248, 240, 210
Offset: 1

Views

Author

Jason Earls, Jul 11 2001

Keywords

Comments

a(n) is also the total number of parts in all partitions of 2*n into equal parts. - Omar E. Pol, Feb 14 2021

Crossrefs

Sigma(k*n): A000203 (k=1), A144613 (k=3), A193553 (k=4, even bisection), A283118 (k=5), A224613 (k=6), A283078 (k=7), A283122 (k=8), A283123 (k=9).
Cf. A008438, A074400, A182818, A239052 (odd bisection), A326124 (partial sums), A054784, A215947, A336923, A346870, A346878, A346880, A355750.
Row 2 of A319526. Column & Row 2 of A216626. Row 1 of A355927.
Shallow diagonal (2n,n) of A265652. See also A244658.

Programs

Formula

a(n) = A000203(2*n). - R. J. Mathar, Apr 06 2011
a(n) = A000203(n) + A054785(n). - R. J. Mathar, May 19 2020
From Vaclav Kotesovec, Aug 07 2022: (Start)
Dirichlet g.f.: zeta(s) * zeta(s-1) * (3 - 2^(1-s)).
Sum_{k=1..n} a(k) ~ 5 * Pi^2 * n^2 / 24. (End)
From Miles Wilson, Sep 30 2024: (Start)
G.f.: Sum_{k>=1} k*x^(k/gcd(k, 2))/(1 - x^(k/gcd(k, 2))).
G.f.: Sum_{k>=1} k*x^(2*k/(3 + (-1)^k))/(1 - x^(2*k/(3 + (-1)^k))). (End)

Extensions

Zero removed and offset corrected by Omar E. Pol, Jul 17 2009

A346878 Sum of the divisors, except for the largest, of the n-th positive even number.

Original entry on oeis.org

1, 3, 6, 7, 8, 16, 10, 15, 21, 22, 14, 36, 16, 28, 42, 31, 20, 55, 22, 50, 54, 40, 26, 76, 43, 46, 66, 64, 32, 108, 34, 63, 78, 58, 74, 123, 40, 64, 90, 106, 44, 140, 46, 92, 144, 76, 50, 156, 73, 117, 114, 106, 56, 172, 106, 136, 126, 94, 62, 240, 64, 100, 186, 127
Offset: 1

Views

Author

Omar E. Pol, Aug 20 2021

Keywords

Comments

Sum of aliquot divisors (or aliquot parts) of the n-th positive even number.
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 for the largest is 1 + 2 + 5 = 8, so a(5) = 8.
		

Crossrefs

Programs

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

Formula

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

A346870 Sum of all divisors, except the smallest and the largest of every number, of the first n positive even numbers.

Original entry on oeis.org

0, 2, 7, 13, 20, 35, 44, 58, 78, 99, 112, 147, 162, 189, 230, 260, 279, 333, 354, 403, 456, 495, 520, 595, 637, 682, 747, 810, 841, 948, 981, 1043, 1120, 1177, 1250, 1372, 1411, 1474, 1563, 1668, 1711, 1850, 1895, 1986, 2129, 2204, 2253, 2408, 2480, 2596, 2709, 2814
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

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

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 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[n_] := DivisorSigma[1, 2*n] - 2*n - 1; Accumulate @ Array[s, 50] (* Amiram Eldar, Aug 19 2021 *)
  • Python
    from sympy import divisors
    from itertools import accumulate
    def A346880(n): return sum(divisors(2*n)[1:-1])
    def aupton(nn): return list(accumulate(A346880(n) for n in range(1, nn+1)))
    print(aupton(52)) # Michael S. Branicky, Aug 19 2021
    
  • Python
    from math import isqrt
    def A346870(n): return (t:=isqrt(m:=n>>1))**2*(t+1) - sum((q:=m//k)*((k<<1)+q+1) for k in range(1,t+1))-3*((s:=isqrt(n))**2*(s+1) - sum((q:=n//k)*((k<<1)+q+1) for k in range(1,s+1))>>1)-n*(n+2) # Chai Wah Wu, Nov 02 2023

Formula

a(n) = (5*Pi^2/24 - 1) * n^2 + O(n*log(n)). - Amiram Eldar, May 15 2023

A346879 Sum of the divisors, except the smallest and the largest, of the n-th odd number.

Original entry on oeis.org

0, 0, 0, 0, 3, 0, 0, 8, 0, 0, 10, 0, 5, 12, 0, 0, 14, 12, 0, 16, 0, 0, 32, 0, 7, 20, 0, 16, 22, 0, 0, 40, 18, 0, 26, 0, 0, 48, 18, 0, 39, 0, 22, 32, 0, 20, 34, 24, 0, 56, 0, 0, 86, 0, 0, 40, 0, 28, 64, 24, 11, 44, 30, 0, 46, 0, 26, 104, 0, 0, 50, 24, 34, 80, 0, 0, 80, 36
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

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 the smaller and the largest is 3, so a(5) = 3.
For n = 6 the 6th odd number is 11 and the divisors of 11 are [1, 11] and the sum of the divisors of 11 except the smaller and the largest is 0, so a(6) = 0.
		

Crossrefs

Bisection of A048050.
Partial sums give A346869.

Programs

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

Formula

a(n) = A048050(2*n-1).

A347154 Sum of all divisors, except the largest of every number, of the first n positive even numbers.

Original entry on oeis.org

1, 4, 10, 17, 25, 41, 51, 66, 87, 109, 123, 159, 175, 203, 245, 276, 296, 351, 373, 423, 477, 517, 543, 619, 662, 708, 774, 838, 870, 978, 1012, 1075, 1153, 1211, 1285, 1408, 1448, 1512, 1602, 1708, 1752, 1892, 1938, 2030, 2174, 2250, 2300, 2456, 2529, 2646, 2760
Offset: 1

Views

Author

Omar E. Pol, Aug 20 2021

Keywords

Comments

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

Crossrefs

Programs

  • Mathematica
    s[n_] := DivisorSigma[1, 2*n] - 2*n; Accumulate @ Array[s, 100] (* Amiram Eldar, Aug 20 2021 *)
  • PARI
    a(n) = sum(k=1, n, k*=2; sigma(k)-k); \\ Michel Marcus, Aug 20 2021
    
  • Python
    from sympy import divisors
    from itertools import accumulate
    def A346878(n): return sum(divisors(2*n)[:-1])
    def aupton(nn): return list(accumulate(A346878(n) for n in range(1, nn+1)))
    print(aupton(51)) # Michael S. Branicky, Aug 20 2021
    
  • Python
    from math import isqrt
    def A347154(n): return (t:=isqrt(m:=n>>1))**2*(t+1) - sum((q:=m//k)*((k<<1)+q+1) for k in range(1,t+1))-3*((s:=isqrt(n))**2*(s+1) - sum((q:=n//k)*((k<<1)+q+1) for k in range(1,s+1))>>1)-n*(n+1) # Chai Wah Wu, Nov 02 2023

Formula

a(n) = n + A346870(n).
a(n) = (5*Pi^2/24 - 1) * n^2 + O(n*log(n)). - Amiram Eldar, May 15 2023

A348554 Irregular triangle read by rows: row n gives the divisors d of 2*n with 1 < d < 2*n, for n >= 2.

Original entry on oeis.org

2, 2, 3, 2, 4, 2, 5, 2, 3, 4, 6, 2, 7, 2, 4, 8, 2, 3, 6, 9, 2, 4, 5, 10, 2, 11, 2, 3, 4, 6, 8, 12, 2, 13, 2, 4, 7, 14, 2, 3, 5, 6, 10, 15, 2, 4, 8, 16, 2, 17, 2, 3, 4, 6, 9, 12, 18, 2, 19, 2, 4, 5, 8, 10, 20, 2, 3, 6, 7, 14, 21, 2, 4, 11, 22, 2, 23, 2, 3, 4, 6, 8, 12, 16, 24, 2, 5, 10, 25
Offset: 2

Views

Author

Wolfdieter Lang, Oct 22 2021

Keywords

Comments

This gives the rows 2*n of A137510, for n >= 2.
The length of row n is A069930(n) = tau(2*n) - 2 = A099777(n) - 2.
The sum of row n is A346880(n) = A062731(n) - (2*n + 1).

Examples

			The irregular triangle T(n, k) begins:
n, 2*n / k 1  2  3  4  5  6  7 ...
----------------------------------
2,   4:    2
3,   6:    2  3
4,   8:    2  4
5,  10:    2  5
6   12:    2  3  4 6
7,  14:    2  7
8,  16:    2  4  8
9,  18:    2  3  6  9
10, 20:    2  4  5 10
11, 22:    2 11
12, 24:    2  3  4  6  8 12
13, 26:    2 13
14, 28:    2  4  7 14
15, 30:    2  3  5  6 10 15
16, 32:    2  4  8  1
17, 34:    2 17
18, 36:    2  3  4  6  9 12 18
19, 38:    2 19
20, 40:    2  4  5  8 10 20
...
		

Crossrefs

Programs

  • Mathematica
    Flatten@Table[Select[Divisors[2n],1<#<2n&],{n,2,25}] (* Giorgos Kalogeropoulos, Oct 22 2021 *)
  • PARI
    row(n) = select(x->((x>1) && (x<2*n)), divisors(2*n)); \\ Michel Marcus, Oct 23 2021

Formula

T(n, k) = A137510(2*n, k), for n >= 2 and k = 1, 2, ..., A069930(n).
Showing 1-6 of 6 results.