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 24 results. Next

A210648 Continued fraction of A137245 = sum( 1/(p log p), p prime ) = 1.63661632335126...

Original entry on oeis.org

1, 1, 1, 1, 3, 32, 2, 3, 3, 1, 3, 2, 3, 6, 1, 3, 1, 3, 6, 8, 1, 4, 17, 1, 2, 3, 3, 1, 2, 2, 3, 5, 15, 1, 1, 5, 1, 1, 3, 2, 1, 1, 1, 3, 1, 2, 2, 2, 1, 1, 1, 25, 1, 8, 10, 1, 1, 11, 2, 2, 25, 2, 16, 1, 2, 1, 2, 1, 1, 1, 2, 11, 4, 35, 3, 1, 1, 9, 59, 1, 2, 2, 1
Offset: 0

Views

Author

M. F. Hasler, Mar 27 2012

Keywords

Programs

  • PARI
    s=1.6366163233512608685696580039218636711815970761312; e=.1^default(realprecision); a=[contfrac(s+e),contfrac(s-e)]; for(n=1,min(#a[1],#a[2]), a[1][n] == a[2][n] & print1(a[1][n]","))

Extensions

More terms from Vaclav Kotesovec, Jun 12 2022

A000040 The prime numbers.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271
Offset: 1

Views

Author

Keywords

Comments

See A065091 for comments, formulas etc. concerning only odd primes. For all information concerning prime powers, see A000961. For contributions concerning "almost primes" see A002808.
A number p is prime if (and only if) it is greater than 1 and has no positive divisors except 1 and p.
A natural number is prime if and only if it has exactly two (positive) divisors.
A prime has exactly one proper positive divisor, 1.
The paper by Kaoru Motose starts as follows: "Let q be a prime divisor of a Mersenne number 2^p-1 where p is prime. Then p is the order of 2 (mod q). Thus p is a divisor of q - 1 and q > p. This shows that there exist infinitely many prime numbers." - Pieter Moree, Oct 14 2004
1 is not a prime, for if the primes included 1, then the factorization of a natural number n into a product of primes would not be unique, since n = n*1.
Prime(n) and pi(n) are inverse functions: A000720(a(n)) = n and a(n) is the least number m such that a(A000720(m)) = a(n). a(A000720(n)) = n if (and only if) n is prime.
Second sequence ever computed by electronic computer, on EDSAC, May 09 1949 (see Renwick link). - Russ Cox, Apr 20 2006
Every prime p > 3 is a linear combination of previous primes prime(n) with nonzero coefficients c(n) and |c(n)| < prime(n). - Amarnath Murthy, Franklin T. Adams-Watters and Joshua Zucker, May 17 2006; clarified by Chayim Lowen, Jul 17 2015
The Greek transliteration of 'Prime Number' is 'Protos Arithmos'. - Daniel Forgues, May 08 2009 [Edited by Petros Hadjicostas, Nov 18 2019]
A number n is prime if and only if it is different from zero and different from a unit and each multiple of n decomposes into factors such that n divides at least one of the factors. This applies equally to the integers (where a prime has exactly four divisors (the definition of divisors is relaxed such that they can be negative)) and the positive integers (where a prime has exactly two distinct divisors). - Peter Luschny, Oct 09 2012
Motivated by his conjecture on representations of integers by alternating sums of consecutive primes, for any positive integer n, Zhi-Wei Sun conjectured that the polynomial P_n(x) = Sum_{k=0..n} a(k+1)*x^k is irreducible over the field of rational numbers with the Galois group S_n, and moreover P_n(x) is irreducible mod a(m) for some m <= n(n+1)/2. It seems that no known criterion on irreducibility of polynomials implies this conjecture. - Zhi-Wei Sun, Mar 23 2013
Questions on a(2n) and Ramanujan primes are in A233739. - Jonathan Sondow, Dec 16 2013
From Hieronymus Fischer, Apr 02 2014: (Start)
Natural numbers such that there is exactly one base b such that the base-b alternate digital sum is 0 (see A239707).
Equivalently: Numbers p > 1 such that b = p-1 is the only base >= 1 for which the base-b alternate digital sum is 0.
Equivalently: Numbers p > 1 such that the base-b alternate digital sum is <> 0 for all bases 1 <= b < p-1. (End)
An integer n > 1 is a prime if and only if it is not the sum of positive integers in arithmetic progression with common difference 2. - Jean-Christophe Hervé, Jun 01 2014
Conjecture: Numbers having prime factors <= prime(n+1) are {k|k^f(n) mod primorial(n)=1}, where f(n) = lcm(prime(i)-1, i=1..n) = A058254(n) and primorial(n) = A002110(n). For example, numbers with no prime divisor <= prime(7) = 17 are {k|k^60 mod 30030=1}. - Gary Detlefs, Jun 07 2014
Cramer conjecture prime(n+1) - prime(n) < C log^2 prime(n) is equivalent to the inequality (log prime(n+1)/log prime(n))^n < e^C, as n tend to infinity, where C is an absolute constant. - Thomas Ordowski, Oct 06 2014
I conjecture that for any positive rational number r there are finitely many primes q_1,...,q_k such that r = Sum_{j=1..k} 1/(q_j-1). For example, 2 = 1/(2-1) + 1/(3-1) + 1/(5-1) + 1/(7-1) + 1/(13-1) with 2, 3, 5, 7 and 13 all prime, 1/7 = 1/(13-1) + 1/(29-1) + 1/(43-1) with 13, 29 and 43 all prime, and 5/7 = 1/(3-1) + 1/(7-1) + 1/(31-1) + 1/(71-1) with 3, 7, 31 and 71 all prime. - Zhi-Wei Sun, Sep 09 2015
I also conjecture that for any positive rational number r there are finitely many primes p_1,...,p_k such that r = Sum_{j=1..k} 1/(p_j+1). For example, 1 = 1/(2+1) + 1/(3+1) + 1/(5+1) + 1/(7+1) + 1/(11+1) + 1/(23+1) with 2, 3, 5, 7, 11 and 23 all prime, and 10/11 = 1/(2+1) + 1/(3+1) + 1/(5+1) + 1/(7+1) + 1/(43+1) + 1/(131+1) + 1/(263+1) with 2, 3, 5, 7, 43, 131 and 263 all prime. - Zhi-Wei Sun, Sep 13 2015
Numbers k such that ((k-2)!!)^2 == +-1 (mod k). - Thomas Ordowski, Aug 27 2016
Does not satisfy Benford's law [Diaconis, 1977; Cohen-Katz, 1984; Berger-Hill, 2017]. - N. J. A. Sloane, Feb 07 2017
Prime numbers are the integer roots of 1 - sin(Pi*Gamma(s)/s)/sin(Pi/s). - Peter Luschny, Feb 23 2018
Conjecture: log log a(n+1) - log log a(n) < 1/n. - Thomas Ordowski, Feb 17 2023

Examples

			From _David A. Corneth_, Oct 22 2024: (Start)
7 is a prime number as it has exactly two divisors, 1 and 7.
8 is not a prime number as it does not have exactly two divisors (it has 1, 2, 4 and 8 as divisors though it is sufficient to find one other divisor than 1 and 8)
55 is not a prime number as it does not have exactly two divisors. One other divisor than 1 and 55 is 5.
59 is a prime number as it has exactly two divisors; 1 and 59. (End)
		

References

  • M. Aigner and G. M. Ziegler, Proofs from The Book, Springer-Verlag, Berlin, 2nd. ed., 2001; see p. 3.
  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 2.
  • E. Bach and Jeffrey Shallit, Algorithmic Number Theory, I, Chaps. 8, 9.
  • D. M. Bressoud, Factorization and Primality Testing, Springer-Verlag NY 1989.
  • M. Cipolla, "La determinazione asintotica dell'n-mo numero primo.", Rend. d. R. Acc. di sc. fis. e mat. di Napoli, s. 3, VIII (1902), pp. 132-166.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 127-149.
  • R. Crandall and C. Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see p. 1.
  • Martin Davis, "Algorithms, Equations, and Logic", pp. 4-15 of S. Barry Cooper and Andrew Hodges, Eds., "The Once and Future Turing: Computing the World", Cambridge 2016.
  • J.-P. Delahaye, Merveilleux nombres premiers, Pour la Science-Belin Paris, 2000.
  • J.-P. Delahaye, Savoir si un nombre est premier: facile, Pour La Science, 303(1) 2003, pp. 98-102.
  • M. Dietzfelbinger, Primality Testing in Polynomial Time, Springer NY 2004.
  • M. du Sautoy, The Music of the Primes, Fourth Estate / HarperCollins, 2003; see p. 5.
  • J. Elie, "L'algorithme AKS", in 'Quadrature', No. 60, pp. 22-32, 2006 EDP-sciences, Les Ulis (France);
  • W. & F. Ellison, Prime Numbers, Hermann Paris 1985
  • T. Estermann, Introduction to Modern Prime Number Theory, Camb. Univ. Press, 1969.
  • J. M. Gandhi, Formulae for the nth prime. Proc. Washington State Univ. Conf. on Number Theory, 96-106. Wash. St. Univ., Pullman, Wash., 1971.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §3.2 Prime Numbers, pp. 77-78.
  • R. K. Guy, Unsolved Problems Number Theory, Section A.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 2.
  • Peter Hilton and Jean Pedersen, A Mathematical Tapestry: Demonstrating the Beautiful Unity of Mathematics, Cambridge University Press, 2010, pp. (260-264).
  • H. D. Huskey, Derrick Henry Lehmer [1905-1991]. IEEE Ann. Hist. Comput. 17 (1995), no. 2, 64-68. Math. Rev. 96b:01035, cf. http://www.ams.org/mathscinet-getitem?mr=1336709
  • M. N. Huxley, The Distribution of Prime Numbers, Oxford Univ. Press, 1972.
  • D. S. Jandu, Prime Numbers And Factorization, Infinite Bandwidth Publishing, N. Hollywood CA 2007.
  • E. Landau, Handbuch der Lehre von der Verteilung der Primzahlen, Chelsea, NY, 1974.
  • D. H. Lehmer, The sieve problem for all-purpose computers. Math. Tables and Other Aids to Computation, Math. Tables and Other Aids to Computation, 7, (1953). 6-14. Math. Rev. 14:691e
  • D. N. Lehmer, "List of Prime Numbers from 1 to 10,006,721", Carnegie Institute, Washington, D.C. 1909.
  • W. J. LeVeque, Topics in Number Theory. Addison-Wesley, Reading, MA, 2 vols., 1956, Vol. 1, Chap. 6.
  • H. Lifchitz, Table des nombres premiers de 0 à 20 millions (Tomes I & II), Albert Blanchard, Paris 1971.
  • R. F. Lukes, C. D. Patterson and H. C. Williams, Numerical sieving devices: their history and some applications. Nieuw Arch. Wisk. (4) 13 (1995), no. 1, 113-139. Math. Rev. 96m:11082, cf http://www.ams.org/mathscinet-getitem?mr=96m:11082
  • P. Ribenboim, The New Book of Prime Number Records, Springer-Verlag NY 1995.
  • P. Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004.
  • H. Riesel, Prime Numbers and Computer Methods for Factorization, Birkhäuser Boston, Cambridge MA 1994.
  • B. Rittaud, "31415879. Ce nombre est-il premier?" ['Is this number prime?'], La Recherche, Vol. 361, pp. 70-73, Feb 15 2003, Paris.
  • D. Shanks, Solved and Unsolved Problems in Number Theory, 2nd. ed., Chelsea, 1978, Chap. 1.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 107-119.
  • D. Wells, Prime Numbers: The Most Mysterious Figures In Math, J. Wiley NY 2005.
  • H. C. Williams and Jeffrey Shallit, Factoring integers before computers. Mathematics of Computation 1943-1993: a half-century of computational mathematics (Vancouver, BC, 1993), 481-531, Proc. Sympos. Appl. Math., 48, AMS, Providence, RI, 1994. Math. Rev. 95m:11143

Crossrefs

For is_prime and next_prime, see A010051 and A151800.
Cf. A000720 ("pi"), A001223 (differences between primes), A002476, A002808, A003627, A006879, A006880, A008578, A080339, A233588.
Cf. primes in lexicographic order: A210757, A210758, A210759, A210760, A210761.
Cf. A003558, A179480 (relating to the Quasi-order theorem of Hilton and Pedersen).
Boustrophedon transforms: A000747, A000732, A230953.
a(2n) = A104272(n) - A233739(n).
Related sequences:
Primes (p) and composites (c): A002808, A000720, A065855.
Primes between p(n) and 2*p(n): A063124, A070046; between c(n) and 2*c(n): A376761; between n and 2*n: A035250, A060715, A077463, A108954.
Composites between p(n) and 2*p(n): A246514; between c(n) and 2*c(n): A376760; between n and 2*n: A075084, A307912, A307989, A376759.

Programs

  • GAP
    A000040:=Filtered([1..10^5],IsPrime); # Muniru A Asiru, Sep 04 2017
    
  • Haskell
    -- See also Haskell Wiki Link.
    import Data.List (genericIndex)
    a000040 n = genericIndex a000040_list (n - 1)
    a000040_list = base ++ larger where
    base = [2,3,5,7,11,13,17]
    larger = p : filter prime more
    prime n = all ((> 0) . mod n) $ takeWhile (\x -> x*x <= n) larger
    _ : p : more = roll $ makeWheels base
    roll (Wheel n rs) = [n * k + r | k <- [0..], r <- rs]
    makeWheels = foldl nextSize (Wheel 1 [1])
    nextSize (Wheel size bs) p = Wheel (size * p)
    [r | k <- [0..p-1], b <- bs, let r = size*k+b, mod r p > 0]
    data Wheel = Wheel Integer [Integer]
    -- Reinhard Zumkeller, Apr 07 2014
    
  • Magma
    [n : n in [2..500] | IsPrime(n)];
    
  • Magma
    a := func< n | NthPrime(n) >;
    
  • Maple
    A000040 := n->ithprime(n); [ seq(ithprime(i),i=1..100) ];
    # For illustration purposes only:
    isPrime := s -> is(1 = sin(Pi*GAMMA(s)/s)/sin(Pi/s)):
    select(isPrime, [$2..100]); # Peter Luschny, Feb 23 2018
  • Mathematica
    Prime[Range[60]]
  • Maxima
    A000040(n) := block(
    if n = 1 then return(2),
    return( next_prime(A000040(n-1)))
    )$ /* recursive, to be replaced if possible - R. J. Mathar, Feb 27 2012 */
    
  • PARI
    {a(n) = if( n<1, 0, prime(n))};
    
  • PARI
    /* The following functions provide asymptotic approximations, one based on the asymptotic formula cited above (slight overestimate for n > 10^8), the other one based on pi(x) ~ li(x) = Ei(log(x)) (slight underestimate): */
    prime1(n)=n*(log(n)+log(log(n))-1+(log(log(n))-2)/log(n)-((log(log(n))-6)*log(log(n))+11)/log(n)^2/2)
    prime2(n)=solve(X=n*log(n)/2,2*n*log(n),real(eint1(-log(X)))+n)
    \\ M. F. Hasler, Oct 21 2013
    
  • PARI
    forprime(p=2, 10^3, print1(p, ", ")) \\ Felix Fröhlich, Jun 30 2014
    
  • PARI
    primes(10^5) \\ Altug Alkan, Mar 26 2018
    
  • Python
    from sympy import primerange
    print(list(primerange(2, 272))) # Michael S. Branicky, Apr 30 2022
  • Sage
    a = sloane.A000040
    a.list(58)  # Jaap Spies, 2007
    
  • Sage
    prime_range(1, 300)  # Zerinvary Lajos, May 27 2009
    

Formula

The prime number theorem is the statement that a(n) ~ n * log n as n -> infinity (Hardy and Wright, page 10).
For n >= 2, n*(log n + log log n - 3/2) < a(n); for n >= 20, a(n) < n*(log n + log log n - 1/2). [Rosser and Schoenfeld]
For all n, a(n) > n log n. [Rosser]
n log(n) + n (log log n - 1) < a(n) < n log n + n log log n for n >= 6. [Dusart, quoted in the Wikipedia article]
a(n) = n log n + n log log n + (n/log n)*(log log n - log n - 2) + O( n (log log n)^2/ (log n)^2). [Cipolla, see also Cesàro or the "Prime number theorem" Wikipedia article for more terms in the expansion]
a(n) = 2 + Sum_{k = 2..floor(2n*log(n)+2)} (1-floor(pi(k)/n)), for n > 1, where the formula for pi(k) is given in A000720 (Ruiz and Sondow 2002). - Jonathan Sondow, Mar 06 2004
I conjecture that Sum_{i>=1} (1/(prime(i)*log(prime(i)))) = Pi/2 = 1.570796327...; Sum_{i=1..100000} (1/(prime(i)*log(prime(i)))) = 1.565585514... It converges very slowly. - Miklos Kristof, Feb 12 2007
The last conjecture has been discussed by the math.research newsgroup recently. The sum, which is greater than Pi/2, is shown in sequence A137245. - T. D. Noe, Jan 13 2009
A000005(a(n)) = 2; A002033(a(n+1)) = 1. - Juri-Stepan Gerasimov, Oct 17 2009
A001222(a(n)) = 1. - Juri-Stepan Gerasimov, Nov 10 2009
From Gary Detlefs, Sep 10 2010: (Start)
Conjecture:
a(n) = {n| n! mod n^2 = n(n-1)}, n <> 4.
a(n) = {n| n!*h(n) mod n = n-1}, n <> 4, where h(n) = Sum_{k=1..n} 1/k. (End)
For n = 1..15, a(n) = p + abs(p-3/2) + 1/2, where p = m + int((m-3)/2), and m = n + int((n-2)/8) + int((n-4)/8). - Timothy Hopper, Oct 23 2010
a(2n) <= A104272(n) - 2 for n > 1, and a(2n) ~ A104272(n) as n -> infinity. - Jonathan Sondow, Dec 16 2013
Conjecture: Sequence = {5 and n <> 5| ( Fibonacci(n) mod n = 1 or Fibonacci(n) mod n = n - 1) and 2^(n-1) mod n = 1}. - Gary Detlefs, May 25 2014
Conjecture: Sequence = {5 and n <> 5| ( Fibonacci(n) mod n = 1 or Fibonacci(n) mod n = n - 1) and 2^(3*n) mod 3*n = 8}. - Gary Detlefs, May 28 2014
Satisfies a(n) = 2*n + Sum_{k=1..(a(n)-1)} cot(k*Pi/a(n))*sin(2*k*n^a(n)*Pi/a(n)). - Ilya Gutkovskiy, Jun 29 2016
Sum_{n>=1} 1/a(n)^s = P(s), where P(s) is the prime zeta function. - Eric W. Weisstein, Nov 08 2016
a(n) = floor(1 - log(-1/2 + Sum_{ d | A002110(n-1) } mu(d)/(2^d-1))/log(2)) where mu(d) = A008683(d) [Ghandi, 1971] (see Ribenboim). Golomb gave a proof in 1974: Give each positive integer a probability of W(n) = 1/2^n, then the probability M(d) of the integer multiple of number d equals 1/(2^d-1). Suppose Q = a(1)*a(2)*...*a(n-1) = A002110(n-1), then the probability of random integers that are mutually prime with Q is Sum_{ d | Q } mu(d)*M(d) = Sum_{ d | Q } mu(d)/(2^d-1) = Sum_{ gcd(m, Q) = 1 } W(m) = 1/2 + 1/2^a(n) + 1/2^a(n+1) + 1/2^a(n+2) + ... So ((Sum_{ d | Q } mu(d)/(2^d-1)) - 1/2)*2^a(n) = 1 + x(n), which means that a(n) is the only integer so that 1 < ((Sum_{ d | Q } mu(d)/(2^d-1)) - 1/2)*2^a(n) < 2. - Jinyuan Wang, Apr 08 2019
Conjecture: n * (log(n)+log(log(n))-1+((log(log(n))-A)/log(n))) is asymptotic to a(n) if and only if A=2. - Alain Rocchelli, Feb 12 2025
From Stefano Spezia, Apr 13 2025: (Start)
a(n) = 1 + Sum_{m=1..2^n} floor(floor(n/Sum_{j=1..m} A080339(j))^(1/n)) [Willans, 1964].
a(n) = 1 + Sum_{m=1..2^n} floor(floor(n/(1 + A000720(m)))^(1/n)) [Willans, 1964]. (End)

A349394 a(p^e) = p^(e-1) for prime powers, a(n) = 0 for all other n; Dirichlet convolution of A003415 (arithmetic derivative of n) with A055615 (Dirichlet inverse of n).

Original entry on oeis.org

0, 1, 1, 2, 1, 0, 1, 4, 3, 0, 1, 0, 1, 0, 0, 8, 1, 0, 1, 0, 0, 0, 1, 0, 5, 0, 9, 0, 1, 0, 1, 16, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 7, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 32, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 27, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2021

Keywords

Comments

Dirichlet convolution of this sequence with Euler phi (A000010) is A300251.
Convolving this sequence with sigma (A000203) produces A319684.
With a(1) = 1 instead of 0, this would be the Dirichlet convolution of A129283 (A003415(n)+n) with A055615. Thus when we subtract A063524 from that convolution, we get this sequence. (See also A349434). Compare also to the convolution of A069359 (sequence agreeing with A003415 on squarefree numbers) with A055615, which is the characteristic function of primes, A010051. - Antti Karttunen, Nov 20 2021

Crossrefs

Programs

  • Haskell
    import Math.NumberTheory.Primes
    a n = case factorise n of
        [(p,e)] -> unPrime p^(e-1) :: Int
         -> 0 -- _Sebastian Karlsson, Nov 19 2021
  • Mathematica
    f[p_, e_] := e/p; d[1] = 0; d[n_] := n * Plus @@ f @@@ FactorInteger[n]; a[n_] := DivisorSum[n, # * MoebiusMu[#] * d[n/#] &]; Array[a, 100] (* Amiram Eldar, Nov 19 2021 *)
  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
    A055615(n) = (n*moebius(n));
    A349394(n) = sumdiv(n,d,A003415(n/d)*A055615(d));
    
  • PARI
    A349394(n) = { my(p=0,e); if((e=isprimepower(n,&p)),p^(e-1),0); }; \\ (After Sebastian Karlsson's new formula) - Antti Karttunen, Nov 20 2021
    

Formula

a(n) = Sum_{d|n} A003415(n/d) * A055615(d).
a(n) = 0 unless n is a prime power (A246655), in which case a(p^e) = p^(e-1). - Sebastian Karlsson, Nov 19 2021
a(n) = A003557(n) * A069513(n). [From above] - Antti Karttunen, Nov 20 2021
Dirichlet g.f.: Sum_{p prime} 1/(p^s-p) [Follows from the D.g.f. of A003415 proved by Haukkanen et al.]. - Sebastian Karlsson, Nov 25 2021
Sum_{k=1..n} a(k) has an average value c*n, where c = A137245 = Sum_{primes p} 1/(p*log(p)) = 1.63661632335... - Vaclav Kotesovec, Mar 03 2023

Extensions

Added Sebastian Karlsson's formula as the new primary definition - Antti Karttunen, Nov 20 2021

A115563 Decimal expansion of Sum_{n>=2} 1/(n*log(n)^2).

Original entry on oeis.org

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

Views

Author

Pierre CAMI, Mar 11 2006

Keywords

Comments

Sum_{n>1} 1/(n*log(n)^2) is a tiny bit greater than (zeta(2))^(3/2) = (Pi^2 / 6)^(3/2) = 2.109709908063657.... - Daniel Forgues, Mar 30 2012
From Bernard Schott, Oct 03 2021: (Start)
Theorem: Bertrand series Sum_{n>=2} 1/(n*log(n)^q) is convergent iff q > 1 (for q = 3, 4, 5 see respectively A145419, A145420, A145421).
As H(n) ~ log(n), compare with A347145. (End)

Examples

			2.10974280123689197447925719761655132638553198439474202264991560319281...
		

Crossrefs

Cf. A137245, A257812. A097906 is a similar sum.

Programs

  • Mathematica
    digits = 150; NSum[1/(n*Log[n]^2), {n, 2, Infinity}, NSumTerms -> 200000, WorkingPrecision -> digits + 5, Method -> {"EulerMaclaurin", Method -> {"NIntegrate", "MaxRecursion" -> 20}}] (* Vaclav Kotesovec, Mar 01 2016, after Jean-François Alcover *)
    maxiter = 20; nn = 10000; alfa = 2; bas = Sum[1/(k*Log[k]^alfa), {k, 2, nn}] + 1/((alfa - 1)*Log[nn + 1/2]^(alfa - 1)); sub = 0; Do[sub = sub + 1/4^s/(2*s + 1)! * NSum[(D[1/(x*Log[x]^alfa), {x, 2 s}]) /. x -> k, {k, nn + 1, Infinity}, WorkingPrecision -> 120, NSumTerms -> 100000, PrecisionGoal -> 120, Method -> {"NIntegrate", "MaxRecursion" -> 100}]; Print[bas - sub], {s, 1, maxiter}] (* Vaclav Kotesovec, Jun 12 2022 *)

Extensions

Removed incorrect speculations about relations to A097906 - R. J. Mathar, Oct 14 2010
More terms from Robert G. Wilson v, Dec 12 2012
Corrected a(55) and beyond, Vaclav Kotesovec, Mar 01 2016

A221711 Decimal expansion of sum 1/(p^2 * log p) over the primes p=2,3,5,7,11,...

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Jan 26 2013

Keywords

Examples

			0.50778218785919931877437510379470557...
		

References

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

Crossrefs

Cf. A137245.

Programs

  • Mathematica
    digits = 106; precision = digits + 15;
    tmax = 400; (* integrand considered negligible beyond tmax *)
    kmax = 400; (* f(k) considered negligible beyond kmax *)
    InLogZeta[k_] := NIntegrate[Log[Zeta[t]], {t, k, tmax},
       WorkingPrecision -> precision, MaxRecursion -> 20,
       AccuracyGoal -> precision];
    f[k_] := With[{mu = MoebiusMu[k]}, If[mu == 0, 0, (mu/k^2)*InLogZeta[2k]]];
    s = 0; Do[s = s + f[k]; Print[k, " ", s], {k, 1, kmax}];
    RealDigits[s][[1]][[1 ;; digits]] (* Jean-François Alcover, Feb 06 2021, updated Jun 23 2022 *)
  • PARI
    \\ See Belabas, Cohen link. Run as SumEulerlog(2) after setting the required precision.
    
  • PARI
    default(realprecision, 200); s=0; for(k=1, 300, s = s + moebius(k)/k^2 * intnum(x=2*k,[[1], 1], log(zeta(x))); print(s)); \\ Vaclav Kotesovec, Jun 12 2022

Extensions

More terms from Hugo Pfoertner, Feb 01 2020
More digits from Vaclav Kotesovec, Jun 12 2022

A319231 Decimal expansion of Sum_{p = prime} 1/(p*log(p)^2).

Original entry on oeis.org

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

Views

Author

R. J. Mathar, Sep 14 2018

Keywords

Comments

Computed by expanding the formalism of arXiv:0811.4739 to double integrals over the Riemann zeta function.

Examples

			1/(2*A253191) + 1/(3*A175478) +1/(5*2.59029...) +1/(7*3.7865)+ ... = 1.52097043...
		

Crossrefs

Programs

  • Mathematica
    digits = 105; precision = digits + 10;
    tmax = 500; (* integrand considered negligible beyond tmax *)
    kmax = 500; (* f(k) considered negligible beyond kmax *)
    InLogZeta[k_] := NIntegrate[(t - k) Log[Zeta[t]], {t, k, tmax}, WorkingPrecision -> precision, MaxRecursion -> 20, AccuracyGoal -> precision];
    f[k_] := With[{mu = MoebiusMu[k]}, If[mu == 0, 0, (mu/k^3)*InLogZeta[k]]];
    s = 0;
    Do[s = s + f[k]; Print[k, " ", s], {k, 1, kmax}];
    RealDigits[s][[1]][[1 ;; digits]] (* Jean-François Alcover, Jun 21 2022, after Vaclav Kotesovec *)
  • PARI
    default(realprecision, 200); s=0; for(k=1, 500, s=s+moebius(k)/k^3 * intnum(x=k,[[1], 1],(x-k)*log(zeta(x))); print(s)); \\ Vaclav Kotesovec, Jun 12 2022

Extensions

More digits from Vaclav Kotesovec, Jun 12 2022

A319232 Decimal expansion of Sum_{p = prime} 1/(p*log p)^2.

Original entry on oeis.org

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

Views

Author

R. J. Mathar, Sep 14 2018

Keywords

Comments

Obtained by expanding the formalism of arXiv:0811.4739 to double integrals over the Riemann zeta function.

Examples

			1/A016627^2 + 1/A016650^2 + 1/8.047189^2 + ... = 0.637056184074676....
		

Crossrefs

Programs

  • Mathematica
    digits = 106; precision = digits + 10;
    tmax = 500; (* integrand considered negligible beyond tmax *)
    kmax = 300; (* f(k) considered negligible beyond kmax *)
    InLogZeta[k_] := NIntegrate[(t - 2k) Log[Zeta[t]], {t, 2k, tmax}, WorkingPrecision -> precision, MaxRecursion -> 20, AccuracyGoal -> precision];
    f[k_] := With[{mu = MoebiusMu[k]}, If[mu == 0, 0, (mu/k^3)*InLogZeta[k]]];
    s = 0;
    Do[s = s + f[k]; Print[k, " ", s], {k, 1, kmax}];
    RealDigits[s][[1]][[1 ;; digits]] (* Jean-François Alcover, Jun 21 2022, after Vaclav Kotesovec *)
  • PARI
    default(realprecision, 200); s=0; for(k=1, 300, s = s + moebius(k)/k^3 * intnum(x=2*k,[[1], 1], (x-2*k)*log(zeta(x))); print(s)); \\ Vaclav Kotesovec, Jun 12 2022

Extensions

More terms from Vaclav Kotesovec, Jun 12 2022

A354917 Decimal expansion of Sum_{p = primes} 1 / (p * log(p)^3).

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, Jun 12 2022

Keywords

Examples

			1.8461474193664495...
		

Crossrefs

Programs

  • Mathematica
    digits = 105; precision = digits + 15;
    tmax = 500; (* integrand considered negligible beyond tmax *)
    kmax = 500; (* f(k) considered negligible beyond kmax *)
    InLogZeta[k_] := NIntegrate[(t - k)^2 Log[Zeta[t]], {t, k, tmax}, WorkingPrecision -> precision, MaxRecursion -> 20, AccuracyGoal -> precision];
    f[k_] := With[{mu = MoebiusMu[k]}, If[mu == 0, 0, (mu/(2 k^4))*InLogZeta[k]]];
    s = 0;
    Do[s = s + f[k]; Print[k, " ", s], {k, 1, kmax}];
    RealDigits[s][[1]][[1 ;; digits]] (* Jean-François Alcover, Jun 21 2022, after Vaclav Kotesovec *)
  • PARI
    default(realprecision, 200); s=0; for(k=1, 500, s = s + moebius(k)/(2*k^4) * intnum(x=k,[[1], 1], (x-k)^2 * log(zeta(x))); print(s));

Extensions

Last digit corrected by Jean-François Alcover and confirmed by Vaclav Kotesovec, Jun 22 2022

A370093 Decimal expansion of Lichtman constant f(N*(2)).

Original entry on oeis.org

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

Views

Author

Artur Jasinski, Feb 09 2024

Keywords

Comments

Definition:
f(N*(k)) = Integral_{s>=1} P_k*(s), where P_k*(s) = Sum_{n>1 and (big) Omega(n)=k} mu(n)^2/n^s, where mu is Möbius (or Moebius) Mu function see A008683, and (big) Omega is number of prime divisors of n counted with multiplicity see A001222.
Lichtman constant f(N*(1)) see A137245.
Lichtman constant f(N*(2)) this sequence.
Lichtman constant f(N*(3)) see A370112.
Lichtman constant f(N*(4)) see A370113.
Limit_{k->oo} f(N*(k)) = 6/Pi^2 = 0.607927101854... see A059956.
Value computed and communicated by Bill Allombert.

Examples

			0.890925479476318332...
		

Crossrefs

Programs

  • PARI
    pz(x)= sum(n=1,max(2,bitprecision(x)/x),my(a=moebius(n));if(a!=0,a*log(zeta(n*x))/n));
    Lichtman(n)=intnum(s=1,[oo,log(2)],exp(-sum(i=1,n,pz(i*s)*x^i/i)+O(x^(n+1)))-1)
    Lichtman(20)
    \\ Bill Allombert, Feb 14 2014 [via Artur Jasinski]

A354887 Decimal expansion of Sum_{primes p} log(log(p)) / (p*log(p)).

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, Jun 10 2022

Keywords

Examples

			0.6410802156599846604833518891513999518913451...
		

Crossrefs

Showing 1-10 of 24 results. Next