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-10 of 41 results. Next

A085548 Decimal expansion of the prime zeta function at 2: Sum_{p prime} 1/p^2.

Original entry on oeis.org

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

Views

Author

Cino Hilliard, Jul 03 2003

Keywords

Comments

Mathar's Table 1 (cited below) lists expansions of the prime zeta function at integers s in 10..39. - Jason Kimberley, Jan 05 2017

Examples

			0.4522474200410654985065... = 1/2^2 + 1/3^2 + 1/5^2 +1/7^2 + 1/11^2 + 1/13^2 + ...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, pp. 94-98.

Crossrefs

Decimal expansion of the prime zeta function: this sequence (at 2), A085541 (at 3), A085964 (at 4) to A085969 (at 9).
Cf. A136271 (derivative), A117543 (semiprimes), A222056, A209329, A124012.

Programs

  • Magma
    R := RealField(106);
    PrimeZeta := func;
    Reverse(IntegerToSequence(Floor(PrimeZeta(2,173)*10^105)));
    // Jason Kimberley, Dec 30 2016
  • Mathematica
    RealDigits[PrimeZetaP[2], 10, 105][[1]]  (* Jean-François Alcover, Jun 24 2011, updated May 06 2021 *)
  • PARI
    recip2(n) = { v=0; p=1; forprime(y=2,n, v=v+1./y^2; ); print(v) }
    
  • PARI
    eps()=my(p=default(realprecision)); precision(2.>>(32*ceil(p*38539962/371253907)),9)
    lm=lambertw(log(4)/eps())\log(4);
    sum(k=1,lm, moebius(k)/k*log(abs(zeta(2*k)))) \\ Charles R Greathouse IV, Jul 19 2013
    
  • PARI
    sumeulerrat(1/p,2) \\ Hugo Pfoertner, Feb 03 2020
    

Formula

P(2) = Sum_{p prime} 1/p^2 = Sum_{n>=1} mobius(n)*log(zeta(2*n))/n. - Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003
Equals A085991 + A086032 + 1/4. - R. J. Mathar, Jul 22 2010
Equals Sum_{k>=1} 1/A001248(k). - Amiram Eldar, Jul 27 2020
Equals Sum_{k>=2} pi(k)*(2*k+1)/(k^2*(k+1)^2), where pi(k) = A000720(k) (Shamos, 2011, p. 9). - Amiram Eldar, Mar 12 2024

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003
Offset corrected by R. J. Mathar, Feb 05 2009

A005596 Decimal expansion of Artin's constant Product_{p=prime} (1-1/(p^2-p)).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

On Simon Plouffe's web page (and in the book freely available at Gutenberg project) the value is given with an error of +1e-31, as "...651641..." instead of "...641641...". In the reference [Wrench, 1961] cited there, these digits are correct. They are also correct on the Plouffe's Inverter page, as computed by Oliveira e Silva, who comments it took 1 hour at 200 MHz with Mathematica. Using Amiram Eldar's PARI program, the same 500 digits are computed instantly (less than 0.1 sec). - M. F. Hasler, Apr 20 2021
Named after the Austrian mathematician Emil Artin (1898-1962). - Amiram Eldar, Jun 20 2021

Examples

			0.37395581361920228805472805434641641511162924860615...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 169.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    a = Exp[-NSum[ (LucasL[n] - 1)/n PrimeZetaP[n], {n, 2, Infinity}, PrecisionGoal -> 500, WorkingPrecision -> 500, NSumTerms -> 100000]]; RealDigits[a, 10, 111][[1]] (* Robert G. Wilson v, Sep 03 2014 taken from Mathematica's Help file on PrimeZetaP *)
  • PARI
    prodinf(n=2,1/zeta(n)^(sumdiv(n, d, moebius(n/d)*(fibonacci(d-1)+fibonacci(d+1)))/n)) \\ Charles R Greathouse IV, Aug 27 2014
    
  • PARI
    prodeulerrat(1-1/(p^2-p)) \\ Amiram Eldar, Mar 12 2021

Formula

Equals Product_{j>=2} 1/Zeta(j)^A006206(j), where Zeta = A013661, A002117 etc. is Riemann's zeta function. - R. J. Mathar, Feb 14 2009
Equals Sum_{k>=1} mu(k)/(k*phi(k)), where mu is the Moebius function (A008683) and phi is the Euler totient function (A000010). - Amiram Eldar, Mar 11 2020
Equals 1/A065488. - Vaclav Kotesovec, Jul 17 2021

Extensions

More terms from Tomás Oliveira e Silva (http://www.ieeta.pt/~tos)

A136141 Decimal expansion of Sum_{p prime} 1/(p*(p-1)).

Original entry on oeis.org

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

Views

Author

R. J. Mathar, Mar 09 2008

Keywords

Comments

Excess of prime factors with multiplicity over distinct prime factors for random (large) integers. - Charles R Greathouse IV, Sep 06 2011
Sum of reciprocals of (proper) prime powers. The sum of reciprocals of all proper powers is A072102. - Charles R Greathouse IV, Apr 24 2012
Decimal expansion of the infinite sum of the reciprocals of the prime powers which are not prime (A246547). - Robert G. Wilson v, May 13 2019
See the second 'Applications' example under the Mathematica help file for the function PrimePowerQ. - Robert G. Wilson v, May 13 2019
It easy to prove that this constant < 1 (Sum_{p prime} 1/(p*(p-1)) < Sum_{k>=2} 1/(k*(k-1)) = 1). Luthar (1969) asks for a better upper bound. The solution shows that this constant is < 3/2 - log(2) = 0.80685... . - Amiram Eldar, Feb 14 2025

Examples

			Equals 1/2 + 1/(3*2) + 1/(5*4) + 1/(7*6) + ...
= 0.7731566690497951278643674598559423956187413360831860483110060673567...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • Steven R. Finch, Mathematical Constants, Cambridge Univ. Press, 2003, Meissel-Mertens constants, p. 94.

Crossrefs

Cf. A152447 (over the semiprimes), A000040, A000720, A001248, A046660 (excess, see first comment), A072102, A077761, A083342, A179119, A246547.
Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), A085964 (at 4) to A085969 (at 9).

Programs

  • Magma
    R := RealField(105);
    c := &+[R|(EulerPhi(n)-MoebiusMu(n))/n*Log(ZetaFunction(R,n)):n in[2..360]];
    Reverse(IntegerToSequence(Floor(c*10^103))); // Jason Kimberley, Jan 12 2017
  • Mathematica
    digits = 103; sp = NSum[PrimeZetaP[n], {n, 2, Infinity}, WorkingPrecision -> digits + 10, NSumTerms -> 2*digits]; RealDigits[sp, 10, digits] // First (* Jean-François Alcover, Sep 02 2015 *)
  • PARI
    W(x)=solve(y=log(x)/2,max(1,log(x)),y*exp(y)-x)
    eps()=2. >> (32*ceil(default(realprecision)/9.63))
    primezeta(s)=my(t=s*log(2),iter=W(t/eps())\t);sum(k=1,iter, moebius(k)/k*log(abs(zeta(k*s))))
    a(lim,e)={ \\ choose parameters to maximize speed and precision
        my(x,y=exp(W(lim)-.5));
        x=lim^e*(e*log(y))^e*(y*log(y))^-e*incgam(-e,e*log(y));
        forprime(p=2,lim,x+=1/((p*1.)^e*(p-1)));
        x+sum(n=2,e,primezeta(n))
    }; \\ Charles R Greathouse IV, Sep 07 2011
    
  • PARI
    sumeulerrat(1/(p*(p-1))) \\ Amiram Eldar, Mar 18 2021
    

Formula

Equals Sum_{n>=1} 1/A036689(n).
Equals Sum_{s>=2} P(s), where P is the prime zeta function. - Charles R Greathouse IV, Sep 06 2011
Equals A083342 - A077761, that is, Sum_{n>=2} ((EulerPhi(n) - MoebiusMu(n))/n) * log(zeta(n)). - Jean-François Alcover, Sep 02 2015
Equals 2 * Sum_{k>=2} pi(k)/(k^3-k), where pi(k) = A000720(k) (Shamos, 2011, p. 8). - Amiram Eldar, Mar 12 2024

Extensions

More terms from D. S. McNeil, Sep 06 2011
More digits from Jean-François Alcover, Sep 02 2015

A085541 Decimal expansion of the prime zeta function at 3.

Original entry on oeis.org

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

Views

Author

Cino Hilliard, Jul 02 2003

Keywords

Comments

Mathar's Table 1 (cited below) lists expansions of the prime zeta function at integers s in 10..39. - Jason Kimberley, Jan 05 2017

Examples

			0.1747626392994435364231...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • J. W. L. Glaisher, On the Sums of Inverse Powers of the Prime Numbers, Quart. J. Math. 25, 347-362, 1891.

Crossrefs

Decimal expansion of the prime zeta function: A085548 (at 2), this sequence (at 3), A085964 (at 4) to A085969 (at 9).

Programs

  • Magma
    R := RealField(106);
    PrimeZeta := func;
    Reverse(IntegerToSequence(Floor(PrimeZeta(3,117)*10^105)));
    // Jason Kimberley, Dec 30 2016
  • Mathematica
    (* If Mathematica version >= 7.0 then RealDigits[PrimeZetaP[3]//N[#,105]&][[1]] else : *) m = 200; $MaxExtraPrecision = 200; PrimeZetaP[s_] := NSum[MoebiusMu[k]*Log[Zeta[k*s]]/k, {k, 1, m}, AccuracyGoal -> m, NSumTerms -> m, PrecisionGoal -> m, WorkingPrecision -> m]; RealDigits[PrimeZetaP[3]][[1]][[1 ;; 105]] (* Jean-François Alcover, Jun 24 2011 *)
  • PARI
    recip3(n) = { v=0; p=1; forprime(y=2,n, v=v+1./y^3; ); print(v) }
    
  • PARI
    sumeulerrat(1/p,3) \\ Hugo Pfoertner, Feb 03 2020
    

Formula

P(3) = Sum_{p prime} 1/p^3 = Sum_{n>=1} mobius(n)*log(zeta(3*n))/n. - Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003
Equals A086033 + A085992 + 1/8. - R. J. Mathar, Jul 22 2010
Equals Sum_{k>=1} 1/A030078(k). - Amiram Eldar, Jul 27 2020

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003

A065474 Decimal expansion of Product_{p prime} (1 - 2/p^2).

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Nov 19 2001

Keywords

Comments

Density of A007674, squarefree n such that n + 1 is squarefree. - Charles R Greathouse IV, Aug 10 2011
Product_{k>=1} (1 - 2/k^2) = sin(sqrt(2)*Pi) / (sqrt(2)*Pi). - Vaclav Kotesovec, May 23 2020
The asymptotic probability that, for two integers k and m, 0 < k <= m, we have gcd(k*(k+1), m) = 1 (when k and m are chosen at random in the range 1..n and n->oo) (Tóth and Sándor, 1989). - Amiram Eldar, Apr 29 2023

Examples

			0.322634098939244670579531692548...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.

Crossrefs

Programs

  • Mathematica
    $MaxExtraPrecision = 800; digits = 98; terms = 800; P[n_] := PrimeZetaP[n]; LR = Join[{0, 0}, LinearRecurrence[{0, 2}, {-4, 0}, terms + 10]]; r[n_Integer] := LR[[n]]; Exp[NSum[r[n]*P[n - 1]/(n - 1), {n, 3, terms}, NSumTerms -> terms, WorkingPrecision -> digits + 10]] // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Apr 18 2016 *)
  • PARI
    prodeulerrat(1 - 2/p^2) \\ Amiram Eldar, Mar 16 2021

Extensions

Edited by Dean Hickerson, Sep 10 2002
More digits from Vaclav Kotesovec, Dec 18 2019

A085964 Decimal expansion of the prime zeta function at 4.

Original entry on oeis.org

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

Views

Author

Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003

Keywords

Comments

Mathar's Table 1 (cited below) lists expansions of the prime zeta function at integers s in 10..39. - Jason Kimberley, Jan 05 2017

Examples

			0.0769931397642468449426...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • J. W. L. Glaisher, On the Sums of Inverse Powers of the Prime Numbers, Quart. J. Math. 25, 347-362, 1891.

Crossrefs

Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), this sequence (at 4), A085965 (at 5) to A085969 (at 9).

Programs

  • Magma
    R := RealField(106);
    PrimeZeta := func;
    [0]cat Reverse(IntegerToSequence(Floor(PrimeZeta(4,87)*10^105)));
    // Jason Kimberley, Dec 30 2016
    
  • Mathematica
    s[n_] := s[n] = Sum[ MoebiusMu[k]*Log[Zeta[4*k]]/k, {k, 1, n}] // RealDigits[#, 10, 104]& // First // Prepend[#, 0]&; s[100]; s[n = 200]; While[s[n] != s[n - 100], n = n + 100]; s[n] (* Jean-François Alcover, Feb 14 2013 *)
    RealDigits[ PrimeZetaP[ 4], 10, 111][[1]] (* Robert G. Wilson v, Sep 03 2014 *)
  • PARI
    sumeulerrat(1/p,4) \\ Hugo Pfoertner, Feb 03 2020

Formula

P(4) = Sum_{p prime} 1/p^4 = Sum_{n>=1} mobius(n)*log(zeta(4*n))/n
Equals A086034 + A085993 + 1/16. - R. J. Mathar, Jul 22 2010
Equals Sum_{k>=1} 1/A030514(k). - Amiram Eldar, Jul 27 2020

A005597 Decimal expansion of the twin prime constant C_2 = Product_{ p prime >= 3 } (1-1/(p-1)^2).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

C_2 = Product_{ p prime > 2} (p * (p-2) / (p-1)^2) is the 2-tuple case of the Hardy-Littlewood prime k-tuple constant (part of First H-L Conjecture): C_k = Product_{ p prime > k} (p^(k-1) * (p-k) / (p-1)^k).
Although C_2 is commonly called the twin prime constant, it is actually the prime 2-tuple constant (prime pair constant) which is relevant to prime pairs (p, p+2m), m >= 1.
The Hardy-Littlewood asymptotic conjecture for Pi_2m(n), the number of prime pairs (p, p+2m), m >= 1, with p <= n, claims that Pi_2m(n) ~ C_2(2m) * Li_2(n), where Li_2(n) = Integral_{2, n} (dx/log^2(x)) and C_2(2m) = 2 * C_2 * Product_{p prime > 2, p | m} (p-1)/(p-2), which gives: C_2(2) = 2 * C_2 as the prime pair (p, p+2) constant, C_2(4) = 2 * C_2 as the prime pair (p, p+4) constant, C_2(6) = 2* (2/1) * C_2 as the prime pair (p, p+6) constant, C_2(8) = 2 * C_2 as the prime pair (p, p+8) constant, C_2(10) = 2 * (4/3) * C_2 as the prime pair (p, p+10) constant, C_2(12) = 2 * (2/1) * C_2 as the prime pair (p, p+12) constant, C_2(14) = 2 * (6/5) * C_2 as the prime pair (p, p+14) constant, C_2(16) = 2 * C_2 as the prime pair (p, p+16) constant, ... and, for i >= 1, C_2(2^i) = 2 * C_2 as the prime pair (p, p+2^i) constant.
C_2 also occurs as part of other Hardy-Littlewood conjectures related to prime pairs, e.g., the Hardy-Littlewood conjecture concerning the distribution of the Sophie Germain primes (A156874) on primes p such that 2p+1 is also prime.
Another constant related to the twin primes is Viggo Brun's constant B (sometimes also called the twin primes Viggo Brun's constant B_2) A065421, where B_2 = Sum (1/p + 1/q) as (p,q) runs through the twin primes.
Reciprocal of the Selberg-Delange constant A167864. See A167864 for additional comments and references. - Jonathan Sondow, Nov 18 2009
C_2 = Product_{prime p>2} (p-2)p/(p-1)^2 is an analog for primes of Wallis' product 2/Pi = Product_{n=1 to oo} (2n-1)(2n+1)/(2n)^2. - Jonathan Sondow, Nov 18 2009
One can compute a cubic variant, product_{primes >2} (1-1/(p-1)^3) = 0.855392... = (2/3) * 0.6601618...* 1.943596... by multiplying this constant with 2/3 and A082695. - R. J. Mathar, Apr 03 2011
Cohen (1998, p. 7) referred to this number as the "twin prime and Goldbach constant" and noted that, conjecturally, the number of twin prime pairs (p,p+2) with p <= X tends to 2*C_2*X/log(X)^2 as X tends to infinity. - Artur Jasinski, Feb 01 2021

Examples

			0.6601618158468695739278121100145557784326233602847334133194484233354056423...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • Richard Crandall and Carl Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see p. 11.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, Vol. 94, Cambridge University Press, 2003, pp. 84-93, 133.
  • R. K. Guy, Unsolved Problems in Number Theory, Section A8.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, ch. 22.20.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 194, 263-264.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A065645 (continued fraction), A065646 (denominators of convergents to twin prime constant), A065647 (numerators of convergents to twin prime constant), A062270, A062271, A114907, A065418 (C_3), A167864, A000010, A008683.

Programs

  • Mathematica
    s[n_] := (1/n)*N[ Sum[ MoebiusMu[d]*2^(n/d), {d, Divisors[n]}], 160]; C2 = (175/256)*Product[ (Zeta[n]*(1 - 2^(-n))*(1 - 3^(-n))*(1 - 5^(-n))*(1 - 7^(-n)))^(-s[n]), {n, 2, 160}]; RealDigits[C2][[1]][[1 ;; 105]] (* Jean-François Alcover, Oct 15 2012, after PARI *)
    digits = 105; f[n_] := -2*(2^n-1)/(n+1); C2 = Exp[NSum[f[n]*(PrimeZetaP[n+1] - 1/2^(n+1)), {n, 1, Infinity}, NSumTerms -> 5 digits, WorkingPrecision -> 5 digits]]; RealDigits[C2, 10, digits][[1]] (* Jean-François Alcover, Apr 16 2016, updated Apr 24 2018 *)
  • PARI
    \p1000; 175/256*prod(k=2,500,(zeta(k)*(1-1/2^k)*(1-1/3^k)*(1-1/5^k)*(1-1/7^k))^(-sumdiv(k,d,moebius(d)*2^(k/d))/k))
    
  • PARI
    prodeulerrat(1-1/(p-1)^2, 1, 3) \\ Amiram Eldar, Mar 12 2021

Formula

Equals Product_{k>=2} (zeta(k)*(1-1/2^k))^(-Sum_{d|k} mu(d)*2^(k/d)/k). - Benoit Cloitre, Aug 06 2003
Equals 1/A167864. - Jonathan Sondow, Nov 18 2009
Equals Sum_{k>=1} mu(2*k-1)/phi(2*k-1)^2, where mu is the Möbius function (A008683) and phi is the Euler totient function (A000010) (Bruckman, 2001). - Amiram Eldar, Jan 14 2022

Extensions

More terms from Vladeta Jovovic, Nov 08 2001
Commented and edited by Daniel Forgues, Jul 28 2009, Aug 04 2009, Aug 12 2009
PARI code removed by D. S. McNeil, Dec 26 2010

A086242 Decimal expansion of the sum of 1/(p-1)^2 over all primes p.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Jul 13 2003

Keywords

Examples

			1.37506499474863528791725313052243969917959996017...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, Vol. 94, Cambridge University Press, 2003, pp. 94-98.

Crossrefs

Programs

  • Mathematica
    digits = 116; Np = NSum[(n-1)*PrimeZetaP[n], {n, 2, Infinity}, NSumTerms -> 3*digits, WorkingPrecision -> digits+10]; RealDigits[Np, 10, digits] // First (* Jean-François Alcover, Sep 02 2015 *)
  • PARI
    default(realprecision,256);
    (f(k)=return(sum(n=1,1024,moebius(n)/n*log(zeta(k*n)))));
    sum(k=2,1024,(k-1)*f(k)) /* Robert Gerbicz, Sep 12 2012 */
    
  • PARI
    sumeulerrat(1/(p-1)^2) \\ Amiram Eldar, Mar 19 2021

Formula

Equals Sum_{k>=2} (k-1)*primezeta(k). - Robert Gerbicz, Sep 12 2012
Equals lim_{n -> oo} A119686(n)/A334746(n). - Petros Hadjicostas, May 11 2020
Equals Sum_{k>=2} (J_2(k)-phi(k)) * log(zeta(k)) / k, where J_2 = A007434 and phi = A000010 (Jakimczuk, 2017). - Amiram Eldar, Mar 18 2024

Extensions

More digits copied from Cohen's paper by R. J. Mathar, Dec 05 2008
More terms from Robert Gerbicz, Sep 12 2012

A085966 Decimal expansion of the prime zeta function at 6.

Original entry on oeis.org

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

Views

Author

Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003

Keywords

Comments

Mathar's Table 1 (cited below) lists expansions of the prime zeta function at integers s in 10..39. - Jason Kimberley, Jan 07 2017

Examples

			0.0170700868506365129541...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • J. W. L. Glaisher, On the Sums of Inverse Powers of the Prime Numbers, Quart. J. Math. 25, 347-362, 1891.

Crossrefs

Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), A085964 (at 4), A085965 (at 5), this sequence (at 6), A085967 (at 7) to A085969 (at 9).

Programs

  • Magma
    R := RealField(106);
    PrimeZeta := func;
    [0]cat Reverse(IntegerToSequence(Floor(PrimeZeta(6,57)*10^105)));
    // Jason Kimberley, Dec 30 2016
    
  • Mathematica
    s[n_] := s[n] = Sum[ MoebiusMu[k]*Log[Zeta[6*k]]/k, {k, 1, n}] // RealDigits[#, 10, 104]& // First // Prepend[#, 0]&; s[100]; s[n = 200]; While[s[n] != s[n - 100], n = n + 100]; s[n] (* Jean-François Alcover, Feb 14 2013 *)
    RealDigits[ PrimeZetaP[ 6], 10, 111][[1]] (* Robert G. Wilson v, Sep 03 2014 *)
  • PARI
    sumeulerrat(1/p,6) \\ Hugo Pfoertner, Feb 03 2020

Formula

P(6) = Sum_{p prime} 1/p^6 = Sum_{n>=1} mobius(n)*log(zeta(6*n))/n
Equals 1/2^6 + A085995 + A086036. - R. J. Mathar, Jul 14 2012
Equals Sum_{k>=1} 1/A030516(k). - Amiram Eldar, Jul 27 2020

A085965 Decimal expansion of the prime zeta function at 5.

Original entry on oeis.org

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

Views

Author

Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jul 06 2003

Keywords

Comments

Mathar's Table 1 (cited below) lists expansions of the prime zeta function at integers s in 10..39. - Jason Kimberley, Jan 05 2017

Examples

			0.0357550174839242571328...
		

References

  • Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
  • J. W. L. Glaisher, On the Sums of Inverse Powers of the Prime Numbers, Quart. J. Math. 25, 347-362, 1891.

Crossrefs

Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), A085964 (at 4), this sequence (at 5), A085966 (at 6) to A085969 (at 9).

Programs

  • Magma
    R := RealField(106);
    PrimeZeta := func;
    [0] cat Reverse(IntegerToSequence(Floor(PrimeZeta(5,69)*10^105)));
    // Jason Kimberley, Dec 30 2016
    
  • Mathematica
    s[n_] := s[n] = Sum[ MoebiusMu[k]*Log[Zeta[5*k]]/k, {k, 1, n}] // RealDigits[#, 10, 104]& // First // Prepend[#, 0]&; s[100]; s[n=200]; While[s[n] != s[n-100], n = n+100]; s[n] (* Jean-François Alcover, Feb 14 2013, from 1st formula *)
    RealDigits[ PrimeZetaP[ 5], 10, 111][[1]] (* Robert G. Wilson v, Sep 03 2014 *)
  • PARI
    sumeulerrat(1/p,5) \\ Hugo Pfoertner, Feb 03 2020

Formula

P(5) = Sum_{p prime} 1/p^5 = Sum_{n>=1} mobius(n)*log(zeta(5*n))/n.
Equals 1/2^5 + A085994 + A086035. - R. J. Mathar, Jul 14 2012
Equals Sum_{k>=1} 1/A050997(k). - Amiram Eldar, Jul 27 2020
Showing 1-10 of 41 results. Next