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.

Previous Showing 21-30 of 71 results. Next

A035188 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 6.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 0, 1, 1, 2, 0, 1, 0, 0, 2, 1, 0, 1, 2, 2, 0, 0, 2, 1, 3, 0, 1, 0, 2, 2, 0, 1, 0, 0, 0, 1, 0, 2, 0, 2, 0, 0, 2, 0, 2, 2, 2, 1, 1, 3, 0, 0, 2, 1, 0, 0, 2, 2, 0, 2, 0, 0, 0, 1, 0, 0, 2, 0, 2, 0, 2, 1, 2, 0, 3, 2, 0, 0, 0, 2, 1
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 24. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[6, #] &]]; Table[ a[n], {n, 1, 100}] (* G. C. Greubel, Apr 27 2018 *)
  • PARI
    my(m=6); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(6, d)); \\ Amiram Eldar, Nov 20 2023

Formula

From Amiram Eldar, Oct 17 2022: (Start)
a(n) = Sum_{d|n} Kronecker(6, d).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = log(5+2*sqrt(6)) / sqrt(6) = 0.935881... . (End)
Multiplicative with a(p^e) = 1 if Kronecker(6, p) = 0 (p = 2 or 3), a(p^e) = (1+(-1)^e)/2 if Kronecker(6, p) = -1 (p is in A038877), and a(p^e) = e+1 if Kronecker(6, p) = 1 (p is in A097934). - Amiram Eldar, Nov 20 2023

A035195 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 13.

Original entry on oeis.org

1, 0, 2, 1, 0, 0, 0, 0, 3, 0, 0, 2, 1, 0, 0, 1, 2, 0, 0, 0, 0, 0, 2, 0, 1, 0, 4, 0, 2, 0, 0, 0, 0, 0, 0, 3, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 2, 1, 0, 4, 1, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 2, 4, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 5
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 13. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[13, #] &]]; Table[ a[n], {n, 1, 100}] (* G. C. Greubel, Apr 27 2018 *)
  • PARI
    my(m=13); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(13, d)); \\ Amiram Eldar, Nov 18 2023

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*log((3+sqrt(13))/2)/sqrt(13) = 0.662735... . - Amiram Eldar, Oct 11 2022
From Amiram Eldar, Nov 18 2023: (Start)
a(n) = Sum_{d|n} Kronecker(13, d).
Multiplicative with a(13^e) = 1, a(p^e) = (1+(-1)^e)/2 if Kronecker(13, p) = -1 (p is in A038884), and a(p^e) = e+1 if Kronecker(13, p) = 1 (p is in A038883 \ {13}). (End)

A035199 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 17.

Original entry on oeis.org

1, 2, 0, 3, 0, 0, 0, 4, 1, 0, 0, 0, 2, 0, 0, 5, 1, 2, 2, 0, 0, 0, 0, 0, 1, 4, 0, 0, 0, 0, 0, 6, 0, 2, 0, 3, 0, 4, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 1, 2, 0, 6, 2, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 7, 0, 0, 2, 3, 0, 0, 0, 4, 0, 0, 0, 6, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 17. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[17, #] &]]; Table[ a[n], {n, 1, 100}] (* G. C. Greubel, Apr 27 2018 *)
  • PARI
    my(m=17); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(17, d)); \\ Amiram Eldar, Nov 18 2023

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*log(4+sqrt(17))/sqrt(17) = 1.016084... . - Amiram Eldar, Oct 11 2022
From Amiram Eldar, Nov 18 2023: (Start)
a(n) = Sum_{d|n} Kronecker(17, d).
Multiplicative with a(17^e) = 1, a(p^e) = (1+(-1)^e)/2 if Kronecker(17, p) = -1 (p is in A038890), and a(p^e) = e+1 if Kronecker(17, p) = 1 (p is in A038889 \ {17}). (End)

A112605 Number of representations of n as a sum of a square and six times a triangular number.

Original entry on oeis.org

1, 2, 0, 0, 2, 0, 1, 2, 0, 2, 2, 0, 0, 0, 0, 2, 2, 0, 1, 2, 0, 0, 4, 0, 0, 2, 0, 2, 0, 0, 0, 2, 0, 0, 2, 0, 3, 2, 0, 0, 2, 0, 2, 2, 0, 2, 0, 0, 0, 2, 0, 0, 2, 0, 2, 2, 0, 0, 0, 0, 1, 4, 0, 0, 4, 0, 0, 2, 0, 2, 2, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 2, 2, 0, 0, 4, 0, 2, 0, 0, 1, 2, 0, 0, 2, 0, 2, 0, 0, 4, 4, 0, 0, 0, 0
Offset: 0

Views

Author

James Sellers, Dec 21 2005

Keywords

Comments

Ramanujan theta functions: f(q) := Prod_{k>=1} (1-(-q)^k) (see A121373), phi(q) := theta_3(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), psi(q) := Sum_{k=0..oo} q^(k*(k+1)/2) (A010054), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).

Examples

			a(22) = 4 since we can write 22 = 4^2 + 6*1 = (-4)^2 + 6*1 = 2^2 + 6*3 = (-2)^2 + 6*3.
G.f. = 1 + 2*x + 2*x^4 + x^6 + 2*x^7 + 2*x^9 + 2*x^10 + 2*x^15 + 2*x^16 + ... - _Michael Somos_, Aug 11 2009
G.f. = q^3 + 2*q^7 + 2*q^19 + q^27 + 2*q^31 + 2*q^39 + 2*q^43 + 2*q^63 + ... - _Michael Somos_, Aug 11 2009
		

Crossrefs

A112608(n) = a(2*n). 2 * A112609(n) = a(2*n + 1). A112604(n) = a(3*n). 2 * A121361(n) = a(3*n + 1). A112606(n) = a(6*n). 2 * A131962(n) = a(6*n + 1). 2 * A112607(n) = a(6*n + 3). 2 * A131964(n) = a(6*n + 4). - Michael Somos, Aug 11 2009

Programs

  • Mathematica
    a[n_] := DivisorSum[4n+3, KroneckerSymbol[-3, #]&]; Table[a[n], {n, 0, 104}] (* Jean-François Alcover, Dec 04 2015, adapted from PARI *)
  • PARI
    {a(n) = if(n<0, 0, sumdiv(4*n+3, d, kronecker(-3, d)))}; /* Michael Somos, May 20 2006 */
    
  • PARI
    {a(n) = my(A); if(n<0, 0, A = x*O(x^n); polcoeff( eta(x^2+A)^5*eta(x^12+A)^2 / eta(x+A)^2 / eta(x^4+A)^2 / eta(x^6+A), n))}; /* Michael Somos, May 20 2006 */

Formula

a(n) = d_{1, 3}(4n+3) - d_{2, 3}(4n+3) where d_{a, m}(n) equals the number of divisors of n which are congruent to a mod m.
Expansion of q^(-3/4)*eta(q^2)^5*eta(q^12)^2/(eta(q)^2*eta(q^4)^2*eta(q^6)) in powers of q. - Michael Somos, May 20 2006
Euler transform of period 12 sequence [ 2, -3, 2, -1, 2, -2, 2, -1, 2, -3, 2, -2, ...]. - Michael Somos, May 20 2006
a(n)=A002324(4n+3). - Michael Somos, May 20 2006
Expansion of phi(q)*psi(q^6) in powers of q where phi(),psi() are Ramanujan theta functions. - Michael Somos, May 20 2006, Sep 29 2006
G.f. is a period 1 Fourier series which satisfies f(-1 / (48 t)) = 3^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is g.f. for A164273. - Michael Somos, Aug 11 2009
a(3*n + 2) = 0. - Michael Somos, Aug 11 2009

A035192 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 10.

Original entry on oeis.org

1, 1, 2, 1, 1, 2, 0, 1, 3, 1, 0, 2, 2, 0, 2, 1, 0, 3, 0, 1, 0, 0, 0, 2, 1, 2, 4, 0, 0, 2, 2, 1, 0, 0, 0, 3, 2, 0, 4, 1, 2, 0, 2, 0, 3, 0, 0, 2, 1, 1, 0, 2, 2, 4, 0, 0, 0, 0, 0, 2, 0, 2, 0, 1, 2, 0, 2, 0, 0, 0, 2, 3, 0, 2, 2, 0, 0, 4, 2, 1, 5
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 40. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[10, #] &]]; Table[ a[n], {n, 1, 100}] (* G. C. Greubel, Apr 27 2018 *)
  • PARI
    my(m=10); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(10, d)); \\ Amiram Eldar, Nov 18 2023

Formula

From Amiram Eldar, Nov 18 2023: (Start)
a(n) = Sum_{d|n} Kronecker(10, d).
Multiplicative with a(p^e) = 1 if Kronecker(10, p) = 0 (p = 2 or 5), a(p^e) = (1+(-1)^e)/2 if Kronecker(10, p) = -1 (p is in A038880), and a(p^e) = e+1 if Kronecker(10, p) = 1 (p is in A097955).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*log(sqrt(10)+3)/sqrt(10) = 1.1500865228... . (End)

A035194 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 12.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 2, 1, 2, 0, 0, 1, 0, 1, 0, 0, 0, 2, 2, 1, 1, 2, 1, 0, 0, 0, 0, 1, 2, 0, 0, 1, 2, 0, 2, 0, 0, 0, 0, 2, 0, 2, 2, 1, 1, 1, 0, 2, 0, 1, 0, 0, 0, 0, 2, 0, 2, 0, 0, 1, 0, 2, 0, 0, 2, 0, 2, 1, 2, 2, 1, 0, 0, 2, 0, 0, 1
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 12. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[12, #] &]]; Table[ a[n], {n, 1, 100}] (* G. C. Greubel, Apr 27 2018 *)
  • PARI
    my(m=12); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(12, d)); \\ Amiram Eldar, Nov 18 2023

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = log(2+sqrt(3))/sqrt(3) = 0.760345... (A196530). - Amiram Eldar, Oct 11 2022
From Amiram Eldar, Nov 18 2023: (Start)
a(n) = Sum_{d|n} Kronecker(12, d).
Multiplicative with a(p^e) = 1 if Kronecker(12, p) = 0 (p = 2 or 3), a(p^e) = (1+(-1)^e)/2 if Kronecker(12, p) = -1 (p is in A003630), and a(p^e) = e+1 if Kronecker(12, p) = 1 (p is in A097933). (End)

A035210 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 28.

Original entry on oeis.org

1, 1, 2, 1, 0, 2, 1, 1, 3, 0, 0, 2, 0, 1, 0, 1, 0, 3, 2, 0, 2, 0, 0, 2, 1, 0, 4, 1, 2, 0, 2, 1, 0, 0, 0, 3, 2, 2, 0, 0, 0, 2, 0, 0, 0, 0, 2, 2, 1, 1, 0, 0, 2, 4, 0, 1, 4, 2, 2, 0, 0, 2, 3, 1, 0, 0, 0, 0, 0, 0, 0, 3, 0, 2, 2, 2, 0, 0, 0, 0, 5
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 28. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := DivisorSum[n, KroneckerSymbol[28, #] &]; Array[a, 100] (* Amiram Eldar, Nov 19 2023 *)
  • PARI
    my(m = 28); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(28, d)); \\ Amiram Eldar, Nov 19 2023

Formula

From Amiram Eldar, Nov 19 2023: (Start)
a(n) = Sum_{d|n} Kronecker(28, d).
Multiplicative with a(p^e) = 1 if Kronecker(28, p) = 0 (p = 2 or 7), a(p^e) = (1+(-1)^e)/2 if Kronecker(28, p) = -1 (p is in A003632), and a(p^e) = e+1 if Kronecker(28, p) = 1 (p is in A296934).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = log(3*sqrt(7)+8)/sqrt(7) = 1.046454884756... . (End)

A035211 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 29.

Original entry on oeis.org

1, 0, 0, 1, 2, 0, 2, 0, 1, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 2, 0, 0, 2, 0, 3, 0, 0, 2, 1, 0, 0, 0, 0, 0, 4, 1, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 3, 0, 0, 2, 2, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 1, 4, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 29. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := DivisorSum[n, KroneckerSymbol[29, #] &]; Array[a, 100] (* Amiram Eldar, Nov 19 2023 *)
  • PARI
    my(m = 29); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(29, d)); \\ Amiram Eldar, Nov 19 2023

Formula

From Amiram Eldar, Nov 19 2023: (Start)
a(n) = Sum_{d|n} Kronecker(29, d).
Multiplicative with a(29^e) = 1, a(p^e) = (1+(-1)^e)/2 if Kronecker(29, p) = -1 (p is in A038902), and a(p^e) = e+1 if Kronecker(29, p) = 1 (p is in A191022).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*log((sqrt(29)+5)/2)/sqrt(29) = 0.611766289562... . (End)

A035219 Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 37.

Original entry on oeis.org

1, 0, 2, 1, 0, 0, 2, 0, 3, 0, 2, 2, 0, 0, 0, 1, 0, 0, 0, 0, 4, 0, 0, 0, 1, 0, 4, 2, 0, 0, 0, 0, 4, 0, 0, 3, 1, 0, 0, 0, 2, 0, 0, 2, 0, 0, 2, 2, 3, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 1, 0, 0, 2, 0, 0, 0, 2, 0, 2, 0, 2, 0, 4, 0, 0, 0, 5
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Dedekind zeta function for the quadratic number field of discriminant 37. See A002324 for formula and Maple code. - N. J. A. Sloane, Mar 22 2022

Crossrefs

Dedekind zeta functions for imaginary quadratic number fields of discriminants -3, -4, -7, -8, -11, -15, -19, -20 are A002324, A002654, A035182, A002325, A035179, A035175, A035171, A035170, respectively.
Dedekind zeta functions for real quadratic number fields of discriminants 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 40 are A035187, A035185, A035194, A035195, A035199, A035203, A035188, A035210, A035211, A035215, A035219, A035192, respectively.

Programs

  • Mathematica
    a[n_] := DivisorSum[n, KroneckerSymbol[37, #] &]; Array[a, 100] (* Amiram Eldar, Nov 20 2023 *)
  • PARI
    my(m = 37); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
    
  • PARI
    a(n) = sumdiv(n, d, kronecker(37, d)); \\ Amiram Eldar, Nov 20 2023

Formula

From Amiram Eldar, Nov 20 2023: (Start)
a(n) = Sum_{d|n} Kronecker(37, d).
Multiplicative with a(37^e) = 1, a(p^e) = (1+(-1)^e)/2 if Kronecker(37, p) = -1 (p is in A038914), and a(p^e) = e+1 if Kronecker(37, p) = 1 (p is in A191027).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*log(sqrt(37)+6)/sqrt(37) = 0.819292168725... . (End)

A046913 Sum of divisors of n not congruent to 0 mod 3.

Original entry on oeis.org

1, 3, 1, 7, 6, 3, 8, 15, 1, 18, 12, 7, 14, 24, 6, 31, 18, 3, 20, 42, 8, 36, 24, 15, 31, 42, 1, 56, 30, 18, 32, 63, 12, 54, 48, 7, 38, 60, 14, 90, 42, 24, 44, 84, 6, 72, 48, 31, 57, 93, 18, 98, 54, 3, 72, 120, 20, 90, 60, 42, 62, 96, 8, 127, 84, 36, 68, 126
Offset: 1

Views

Author

Keywords

Examples

			Divisors of 12 are 1 2 3 4 6 12 and discarding 3 6 and 12 we get a(12) = 1 + 2 + 4 = 7.
x + 3*x^2 + x^3 + 7*x^4 + 6*x^5 + 3*x^6 + 8*x^7 + 15*x^8 + x^9 + 18*x^10 + ...
		

Crossrefs

Programs

  • Magma
    [SumOfDivisors(3*k)-3*SumOfDivisors(k):k in [1..70]]; // Marius A. Burtea, Jun 01 2019
  • Mathematica
    Table[DivisorSigma[1, 3*w]-3*DivisorSigma[1, w], {w, 1, 256}]
    DivisorSum[#1, # &, Mod[#, 3] != 0 &] & /@ Range[68] (* Jayanta Basu, Jun 30 2013 *)
    f[p_, e_] := If[p == 3, 1, (p^(e+1)-1)/(p-1)]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 17 2020 *)
  • PARI
    {a(n) = if( n<1, 0, sigma(3*n) - 3 * sigma(n))} /* Michael Somos, Jul 19 2004 */
    
  • PARI
    a(n) = sigma(n \ 3^valuation(n, 3)) \\ David A. Corneth, Jun 01 2019
    

Formula

Multiplicative with a(3^e) = 1, a(p^e) = (p^(e+1)-1)/(p-1) for p<>3. - Vladeta Jovovic, Sep 11 2002
G.f.: Sum_{k>0} x^k*(1+2*x^k+2*x^(3*k)+x^(4*k))/(1-x^(3*k))^2. - Vladeta Jovovic, Dec 18 2002
a(n) = A000203(3n)-3*A000203(n). - Labos Elemer, Aug 14 2003
Inverse Mobius transform of A091684. - Gary W. Adamson, Jul 03 2008
Dirichlet g.f.: zeta(s)*zeta(s-1)*(1-1/3^(s-1)). - R. J. Mathar, Feb 10 2011
G.f. A(x) satisfies: 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w)= u^2 + 9 * v^2 + 16 * w^2 - 6 * u*v + 4 * u*w - 24 * v*w - v + w. - Michael Somos, Jul 19 2004
L.g.f.: log(Product_{k>=1} (1 - x^(3*k))/(1 - x^k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Mar 14 2018
a(n) = A002324(n) + 3*Sum_{j=1, n-1} A002324(j)*A002324(n-j). See Farkas and Guerzhoy links. - Michel Marcus, Jun 01 2019
a(3*n) = a(n). - David A. Corneth, Jun 01 2019
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 / 18. - Vaclav Kotesovec, Sep 17 2020
Previous Showing 21-30 of 71 results. Next