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.

A372962 a(n) = Sum_{1 <= x_1, x_2, x_3 <= n} ( n/gcd(x_1, x_2, x_3, n) )^2.

Original entry on oeis.org

1, 29, 235, 925, 3101, 6815, 16759, 29597, 57097, 89929, 160931, 217375, 371125, 486011, 728735, 947101, 1419569, 1655813, 2475739, 2868425, 3938365, 4666999, 6435815, 6955295, 9690601, 10762625, 13874563, 15502075, 20510309, 21133315, 28628191, 30307229, 37818785
Offset: 1

Views

Author

Seiichi Manyama, May 18 2024

Keywords

Links

Crossrefs

Cf. A068963, A084218, A372963.
Cf. A001160, A008683.
Cf. A013662, A013664, A088246.
Cf. A350156, A372964.
Cf. A371492, A372952.

Programs

  • Mathematica
    f[p_, e_] := (p^(5*e+5) - p^(5*e+2) + p^2 - 1)/(p^5-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, May 21 2024 *)
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*(n/d)^2*sigma(d, 5));

Formula

a(n) = Sum_{d|n} mu(n/d) * (n/d)^2 * sigma_5(d).
From Amiram Eldar, May 21 2024: (Start)
Multiplicative with a(p^e) = (p^(5*e+5) - p^(5*e+2) + p^2 - 1)/(p^5-1).
Dirichlet g.f.: zeta(s)*zeta(s-5)/zeta(s-2).
Sum_{k=1..n} a(k) ~ c * n^6 / 6, where c = zeta(6)/zeta(4) = 2*Pi^2/21 = 0.939962323... (1/A088246). (End)
a(n) = Sum_{d|n} phi(n/d) * (n/d)^4 * sigma_4(d^2)/sigma_2(d^2). - Seiichi Manyama, May 24 2024
a(n) = Sum_{1 <= x_1, x_2, x_3 <= n} ( gcd(x_1, n)/gcd(x_1, x_2, x_3, n) )^3. - Seiichi Manyama, May 25 2024

A350156 Inverse Moebius transform of A000056.

Original entry on oeis.org

1, 7, 25, 55, 121, 175, 337, 439, 673, 847, 1321, 1375, 2185, 2359, 3025, 3511, 4897, 4711, 6841, 6655, 8425, 9247, 12145, 10975, 15121, 15295, 18169, 18535, 24361, 21175, 29761, 28087, 33025, 34279, 40777, 37015, 50617, 47887, 54625, 53119, 68881, 58975, 79465, 72655, 81433, 85015
Offset: 1

Views

Author

Werner Schulte, Jan 19 2022

Keywords

Comments

Let f be an arbitrary arithmetic function. Define the sequence a(f; n) by a(f; n) = Sum_{i=1..n, k=1..n} f(n / gcd(gcd(i,k),n)) for n > 0. Then a(f; n) equals inverse Moebius transform of f(n) * A007434(n) for n > 0; if f is multiplicative then a(f; n) is multiplicative; this sequence uses f(n) = n (see formula section).

Crossrefs

Column k=2 of A372968.
Cf. A000010, A000056, A000578, A001158, A007434, A055615, A057660.
Cf. A372962, A372964.

Programs

  • Mathematica
    f[p_, e_] := p^(3*e) - (p - 1)*(p^(3*e) - 1)/(p^3 - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Jan 19 2022 *)
  • Python
    from math import prod
from sympy import factorint
def A350156(n): return prod((q:=p**(3*e))-(p-1)*(q-1)//(p**3-1) for p,e in factorint(n).items()) # Chai Wah Wu, Mar 04 2025

Formula

Multiplicative with a(p^e) = p^(3*e) - (p-1) * (p^(3*e) - 1) / (p^3 - 1) for prime p and e >= 0.
Dirichlet g.f.: Sum_{n>0} a(n) / n^s = zeta(s-3) * zeta(s) / zeta(s-1).
a(n) = Sum_{i=1..n, k=1..n} n / gcd(gcd(i,k),n) for n > 0.
Dirichlet convolution with A000010 equals A000578.
Dirichlet convolution of A001158 and A055615.
Sum_{k=1..n} a(k) ~ c * n^4, where c = Pi^4/(360*zeta(3)) = 0.225098... . - Amiram Eldar, Oct 16 2022
a(n) = Sum_{d|n} phi(n/d) * (n/d)^2 * sigma_2(d^2)/sigma(d^2). - Seiichi Manyama, May 24 2024
a(n) = Sum_{1 <= x_1, x_2 <= n} ( gcd(x_1, n)/gcd(x_1, x_2, n) )^2. - Seiichi Manyama, May 25 2024

A371491 a(n) = Sum_{1 <= x_1, x_2, x_3, x_4, x_5 <= n} ( n/gcd(x_1, x_2, x_3, x_4, x_5, n) )^3.

Original entry on oeis.org

1, 249, 6535, 63737, 390501, 1627215, 5764459, 16316665, 42876109, 97234749, 214357551, 416521295, 815728525, 1435350291, 2551924035, 4177066233, 6975752529, 10676151141, 16983556183, 24889362237, 37670739565, 53375030199, 78310973115, 106629405775
Offset: 1

Views

Author

Seiichi Manyama, May 24 2024

Keywords

Links

Crossrefs

Cf. A068970, A084220, A372950, A372964.
Cf. A372952, A372963.
Cf. A372966.
Cf. A013664, A013667.

Programs

  • Mathematica
    f[p_, e_] := (p^(8*e + 8) - p^(8*e + 3) + p^3 - 1)/(p^8 - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 24] (* Amiram Eldar, May 24 2024 *)
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*(n/d)^3*sigma(d, 8));
  • PARI
    a(n) = sumdiv(n,d, eulerphi(n/d)*(n/d)^3*sigma(d^2, 8)/sigma(d^2, 4));

Formula

a(n) = Sum_{d|n} mu(n/d) * (n/d)^3 * sigma_8(d).
a(n) = Sum_{d|n} phi(n/d) * (n/d)^3 * sigma_8(d^2)/sigma_4(d^2).
From Amiram Eldar, May 24 2024: (Start)
Multiplicative with a(p^e) = (p^(8*e+8) - p^(8*e+3) + p^3 - 1)/(p^8-1).
Dirichlet g.f.: zeta(s)*zeta(s-8)/zeta(s-3).
Sum_{k=1..n} a(k) ~ c * n^9 / 9, where c = zeta(9)/zeta(6) = 0.984926747... . (End)
a(n) = Sum_{1 <= x_1, x_2, x_3, x_4, x_5 <= n} ( gcd(x_1, x_2, n)/gcd(x_1, x_2, x_3, x_4, x_5, n) )^5. - Seiichi Manyama, May 25 2024

A373103 a(n) = Sum_{1 <= x_1, x_2, x_3, x_4, x_5 <= n} ( n/gcd(x_1, x_2, x_3, x_4, x_5, n) )^4.

Original entry on oeis.org

1, 497, 19603, 254449, 1952501, 9742691, 40351207, 130277873, 385845769, 970392997, 2357933051, 4987963747, 10604470813, 20054549879, 38274877103, 66702270961, 118587792977, 191765347193, 322687567459, 496811926949, 791004710821, 1171892726347, 1801152381623
Offset: 1

Views

Author

Seiichi Manyama, May 25 2024

Keywords

Links

Crossrefs

Cf. A371491, A371878, A373007, A373105.
Cf. A350156, A372962, A372964.
Cf. A013664, A013668.

Programs

  • Mathematica
    f[p_, e_] :=  (p^(9*e+9) - p^(9*e+4) + p^4 - 1)/(p^9-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 24] (* Amiram Eldar, May 25 2024 *)
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*(n/d)^4*sigma(d, 9));
  • PARI
    a(n) = sumdiv(n, d, eulerphi(n/d)*(n/d)^8*sigma(d^2, 8)/sigma(d^2, 4));

Formula

a(n) = Sum_{1 <= x_1, x_2, x_3, x_4, x_5 <= n} ( gcd(x_1, n)/gcd(x_1, x_2, x_3, x_4, x_5, n) )^5.
a(n) = Sum_{d|n} mu(n/d) * (n/d)^4 * sigma_9(d).
a(n) = Sum_{d|n} phi(n/d) * (n/d)^8 * sigma_8(d^2)/sigma_4(d^2).
From Amiram Eldar, May 25 2024: (Start)
Multiplicative with a(p^e) = (p^(9*e+9) - p^(9*e+4) + p^4 - 1)/(p^9-1).
Dirichlet g.f.: zeta(s)*zeta(s-9)/zeta(s-4).
Sum_{k=1..n} a(k) ~ c * n^10 / 10, where c = zeta(10)/zeta(6) = Pi^4/99 = 0.983930212464... . (End)

A372950 a(n) = Sum_{1 <= x_1, x_2 <= n} ( n/gcd(x_1, x_2, n) )^3.

Original entry on oeis.org

1, 25, 217, 793, 3001, 5425, 16465, 25369, 52705, 75025, 159721, 172081, 369097, 411625, 651217, 811801, 1414945, 1317625, 2469241, 2379793, 3572905, 3993025, 6424177, 5505073, 9378001, 9227425, 12807289, 13056745, 20486761, 16280425, 28599361, 25977625, 34659457
Offset: 1

Views

Author

Seiichi Manyama, May 18 2024

Keywords

Links

Crossrefs

Cf. A084218, A350156.
Cf. A068970, A084220, A372964.
Cf. A001160, A008683.
Cf. A002117, A013664, A347328.

Programs

  • Mathematica
    f[p_, e_] := (p^(5*e+5) - p^(5*e+3) + p^3 - 1)/(p^5-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, May 21 2024 *)
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*(n/d)^3*sigma(d, 5));

Formula

a(n) = Sum_{d|n} mu(n/d) * (n/d)^3 * sigma_5(d).
From Amiram Eldar, May 21 2024: (Start)
Multiplicative with a(p^e) = (p^(5*e+5) - p^(5*e+3) + p^3 - 1)/(p^5-1).
Dirichlet g.f.: zeta(s)*zeta(s-5)/zeta(s-3).
Sum_{k=1..n} a(k) ~ c * n^6 / 6, where c = zeta(6)/zeta(3) = 0.846335... (A347328). (End)
Dirichlet convolution of A334659 and A001160. - R. J. Mathar, Jul 14 2025

A371628 a(n) = Sum_{1 <= x_1, x_2, x_3, x_4 <= n} ( gcd(x_1, n)/gcd(x_1, x_2, x_3, x_4, n) )^3.

Original entry on oeis.org

1, 65, 757, 4225, 16001, 49205, 119365, 271489, 554797, 1040065, 1783541, 3198325, 4850977, 7758725, 12112757, 17392769, 24211265, 36061805, 47162485, 67604225, 90359305, 115930165, 148291397, 205517173, 250266001, 315313505, 404686153, 504317125, 595481825
Offset: 1

Views

Author

Seiichi Manyama, May 24 2024

Keywords

Links

Crossrefs

Cf. A373059, A371492.
Cf. A372963, A372964.
Cf. A084220.
Cf. A002117, A013662, A013665.

Programs

  • Mathematica
    f[p_, e_] := (p^(6*e+2)*(p^4+p^3+2*p^2+p+1) - p^(4*e+2)*(p^2-p+1) + p^2+p+1)/((p+1)^2*(p^2+1)*(p^2-p+1)); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 30] (* Amiram Eldar, May 24 2024 *)
  • PARI
    a(n) = sumdiv(n, d, eulerphi(n/d)*(n/d)^3*sigma(d^2, 6)/sigma(d^2, 3));

Formula

a(n) = Sum_{d|n} phi(n/d) * (n/d)^3 * sigma_6(d^2)/sigma_3(d^2).
From Amiram Eldar, May 24 2024: (Start)
Multiplicative with a(p^e) = (p^(6*e+2)*(p^4+p^3+2*p^2+p+1) - p^(4*e+2)*(p^2-p+1) + p^2+p+1)/((p+1)^2*(p^2+1)*(p^2-p+1)).
Dirichlet g.f.: zeta(s)*zeta(s-4)*zeta(s-6)/zeta(s-3)^2.
Sum_{k=1..n} a(k) ~ c * n^7 / 7, where c = zeta(3)*zeta(7)/zeta(4)^2 = 1.034718122... . (End)
Showing 1-6 of 6 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.