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

A030078 Cubes of primes.

Original entry on oeis.org

8, 27, 125, 343, 1331, 2197, 4913, 6859, 12167, 24389, 29791, 50653, 68921, 79507, 103823, 148877, 205379, 226981, 300763, 357911, 389017, 493039, 571787, 704969, 912673, 1030301, 1092727, 1225043, 1295029, 1442897, 2048383, 2248091, 2571353, 2685619, 3307949
Offset: 1

Views

Author

Keywords

Comments

Numbers with exactly three factorizations: A001055(a(n)) = 3 (e.g., a(4) = 1*343 = 7*49 = 7*7*7). - Reinhard Zumkeller, Dec 29 2001
Intersection of A014612 and A000578. Intersection of A014612 and A030513. - Wesley Ivan Hurt, Sep 10 2013
Let r(n) = (a(n)-1)/(a(n)+1) if a(n) mod 4 = 1, (a(n)+1)/(a(n)-1) otherwise; then Product_{n>=1} r(n) = (9/7) * (28/26) * (124/126) * (344/342) * (1332/1330) * ... = 48/35. - Dimitris Valianatos, Mar 06 2020
There exist 5 groups of order p^3, when p prime, so this is a subsequence of A054397. Three of them are abelian: C_p^3, C_p^2 X C_p and C_p X C_p X C_p = (C_p)^3. For 8 = 2^3, the 2 nonabelian groups are D_8 and Q_8; for odd prime p, the 2 nonabelian groups are (C_p x C_p) : C_p, and C_p^2 : C_p (remark, for p = 2, these two semi-direct products are isomorphic to D_8). Here C, D, Q mean Cyclic, Dihedral, Quaternion groups of the stated order; the symbols X and : mean direct and semidirect products respectively. - Bernard Schott, Dec 11 2021

Examples

			a(3) = 125; since the 3rd prime is 5, a(3) = 5^3 = 125.
		

References

  • Edmund Landau, Elementary Number Theory, translation by Jacob E. Goodman of Elementare Zahlentheorie (Vol. I_1 (1927) of Vorlesungen über Zahlentheorie), by Edmund Landau, with added exercises by Paul T. Bateman and E. E. Kohlbecker, Chelsea Publishing Co., New York, 1958, pp. 31-32.

Crossrefs

Other sequences that are k-th powers of primes are: A000040 (k=1), A001248 (k=2), this sequence (k=3), A030514 (k=4), A050997 (k=5), A030516 (k=6), A092759 (k=7), A179645 (k=8), A179665 (k=9), A030629 (k=10), A079395 (k=11), A030631 (k=12), A138031 (k=13), A030635 (k=16), A138032 (k=17), A030637 (k=18).
Cf. A060800, A131991, A000578, subsequence of A046099.
Subsequence of A007422 and of A054397.

Programs

Formula

n such that A062799(n) = 3. - Benoit Cloitre, Apr 06 2002
a(n) = A000040(n)^3. - Omar E. Pol, Jul 27 2009
A064380(a(n)) = A000010(a(n)). - Vladimir Shevelev, Apr 19 2010
A003415(a(n)) = A079705(n). - Reinhard Zumkeller, Jun 26 2011
A056595(a(n)) = 2. - Reinhard Zumkeller, Aug 15 2011
A000005(a(n)) = 4. - Wesley Ivan Hurt, Sep 10 2013
a(n) = A119959(n) * A008864(n) -1.- R. J. Mathar, Aug 13 2019
Sum_{n>=1} 1/a(n) = P(3) = 0.1747626392... (A085541). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(3)/zeta(6) (A157289).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(3) (A088453). (End)

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

A003627 Primes of the form 3n-1.

Original entry on oeis.org

2, 5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, 101, 107, 113, 131, 137, 149, 167, 173, 179, 191, 197, 227, 233, 239, 251, 257, 263, 269, 281, 293, 311, 317, 347, 353, 359, 383, 389, 401, 419, 431, 443, 449, 461, 467, 479, 491, 503, 509, 521, 557, 563, 569, 587
Offset: 1

Views

Author

Keywords

Comments

Inert rational primes in the field Q(sqrt(-3)). - N. J. A. Sloane, Dec 25 2017
Primes p such that 1+x+x^2 is irreducible over GF(p). - Joerg Arndt, Aug 10 2011
Primes p dividing sum(k=0,p,C(2k,k)) -1 = A006134(p)-1. - Benoit Cloitre, Feb 08 2003
A039701(A049084(a(n))) = 2; A134323(A049084(a(n))) = -1. - Reinhard Zumkeller, Oct 21 2007
The set of primes of the form 3n - 1 is a superset of the set of lesser of twin primes larger than three (A001359). - Paul Muljadi, Jun 05 2008
Primes of this form do not occur in or as divisors of {n^2+n+1}. See A002383 (n^2+n+1 = prime), A162471 (prime divisors of n^2+n+1 not in A002383), and A002061 (numbers of the form n^2-n+1). - Daniel Tisdale, Jul 04 2009
Or, primes not in A007645. A003627 UNION A007645 = A000040. Also, primes of the form 6*k-5/2-+3/2. - Juri-Stepan Gerasimov, Jan 28 2010
Except for first term "2", all these prime numbers are of the form: 6*n-1. - Vladimir Joseph Stephan Orlovsky, Jul 13 2011
A088534(a(n)) = 0. - Reinhard Zumkeller, Oct 30 2011
For n>1: Numbers k such that (k-4)! mod k =(-1)^(floor(k/3)+1)*floor((k+1)/6), k>4. - Gary Detlefs, Jan 02 2012
Binomial(a(n),3)/a(n)= (3*A024893(n)^2+A024893(n))/2, n>1. - Gary Detlefs, May 06 2012
For every prime p in this sequence, 3 is a 9th power mod p. See Williams link. - Michel Marcus, Nov 12 2017
2 adjoined to A007528. - David A. Corneth, Nov 12 2017
For n >= 2 there exists a polygonal number P_s(3) = 3s - 3 = a(n) + 1. These are the only primes p with P_s(k) = p + 1, s >= 3, k >= 3, since P_s(k) - 1 is composite for k > 3. - Ralf Steiner, May 17 2018

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Primes of form 3n+1 give A002476.
These are the primes arising in A024893, A087370, A088879. A091177 gives prime index.
Subsequence of A034020.

Programs

  • Haskell
    a003627 n = a003627_list !! (n-1)
    a003627_list = filter ((== 2) . (`mod` 3)) a000040_list
    -- Reinhard Zumkeller, Oct 30 2011
    
  • Magma
    [n: n in PrimesUpTo(720) | n mod 3 eq 2]; // Bruno Berselli, Apr 05 2011
    
  • Maple
    t1 := {}; for n from 0 to 500 do if isprime(3*n+2) then t1 := {op(t1),3*n+2}; fi; od: A003627 := convert(t1,list);
  • Mathematica
    Select[Range[-1, 600, 3], PrimeQ[#] &] (* Vincenzo Librandi, Jun 17 2015 *)
    Select[Prime[Range[200]],Mod[#,3]==2&] (* Harvey P. Dale, Jan 31 2023 *)
  • PARI
    is(n)=n%3==2 && isprime(n) \\ Charles R Greathouse IV, Mar 20 2013

Formula

From R. J. Mathar, Apr 03 2011: (Start)
Sum_{n>=1} 1/a(n)^2 = 0.30792... = A085548 - 1/9 - A175644.
Sum_{n>=1} 1/a(n)^3 = 0.134125... = A085541 - 1/27 - A175645. (End)

A077761 Decimal expansion of Mertens's constant, which is the limit of (Sum_{i=1..k} 1/prime(i)) - log(log(prime(k))) as k goes to infinity, where prime(i) is the i-th prime number.

Original entry on oeis.org

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

Views

Author

T. D. Noe, Nov 14 2002

Keywords

Comments

Graham, Knuth & Patashnik incorrectly give this constant as 0.261972128. - Robert G. Wilson v, Dec 02 2005 [This was corrected in the second edition (1994). - T. D. Noe, Mar 11 2017]
Also the average deviation of the number of distinct prime factors: sum_{n < x} omega(n) = x log log x + B_1 x + O(x) where B_1 is this constant, see (e.g.) Hardy & Wright. - Charles R Greathouse IV, Mar 05 2021
Named after the Polish mathematician Franz Mertens (1840-1927). Sometimes called Meissel-Mertens constant, after Mertens and the German astronomer Ernst Meissel (1826-1895). - Amiram Eldar, Jun 16 2021

Examples

			0.26149721284764278375542683860869585905156664826119920619206421392...
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2004, pp. 94-98
  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, A Foundation For Computer Science, Addison-Wesley, Reading, MA, 1989, p. 23.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 4th ed. (1975). Oxford, England: Oxford University Press. See 22.10, "The number of prime factors of n".
  • József Sándor, Dragoslav S. Mitrinovic, Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Section VII.28, p. 257.

Crossrefs

Cf. A001620.

Programs

  • Mathematica
    $MaxExtraPrecision = 400; RealDigits[ N[EulerGamma + NSum[(MoebiusMu[m]/m)*Log[N[Zeta[m], 120]], {m, 2, 1000}, Method -> "EulerMaclaurin", AccuracyGoal -> 120, NSumTerms -> 1000, PrecisionGoal -> 120, WorkingPrecision -> 120] , 120]][[1, 1 ;; 105]]
    (* or, from version 7 up: *) digits = 105; M = EulerGamma - NSum[ PrimeZetaP[n] / n, {n, 2, Infinity}, WorkingPrecision -> digits+10, NSumTerms -> 3*digits]; RealDigits[M, 10, digits] // First (* Jean-François Alcover, Mar 16 2011, updated Sep 01 2015 *)

Formula

Equals A001620 - Sum_{n>=2} zeta_prime(n)/n where the zeta prime sequence is A085548, A085541, A085964, A085965, A085966 etc. [Sebah and Gourdon] - R. J. Mathar, Apr 29 2006
Equals gamma + Sum_{p prime} (log(1-1/p) + 1/p), where gamma is Euler's constant (A001620). - Amiram Eldar, Dec 25 2021
Equals lim_{k->oo} -k + Sum_{p prime} 1/(p*log(p)^(1/k)) conjectured by Meissel in 1866 and proven by Peter Lindqvist and Jaak Peetre in 1997 see links - Artur Jasinski, Mar 11 2025

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

A179119 Decimal expansion of Sum_{p prime} 1/(p*(p+1)).

Original entry on oeis.org

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

Views

Author

R. J. Mathar, Jan 21 2013

Keywords

Examples

			0.33022992626420324101.. = 1/(2*3) +1/(3*4) +1/(5*6) + 1/(7*8) +... = sum_{n>=1} 1/ (A000040(n)*A008864(n)).
		

Crossrefs

Cf. A136141 for 1/(p(p-1)), A085548 for 1/p^2.
Decimal expansion of the prime zeta function: A085548 (at 2), A085541 (at 3), A085964 (at 4) to A085969 (at 9).
Cf. A307379.

Programs

  • Magma
    R:=RealField(103);
    ExhaustSum :=
      function(
        k_min, term
      : IZ := func)
        c:=R!0; k:=k_min;
        repeat
          t:=term(k); c+:=t; k+:=1;
        until IZ(t,k-1);
        return c;
      end function;
    RealField(101)!
    ExhaustSum(2,
      func
        : IZ:=func
        )>);
    // Jason Kimberley, Jan 20 2017
  • Maple
    interface(quiet=true):
    read("transforms") ;
    Digits := 300 ;
    ZetaM := proc(s,M)
        local v,p;
        v := Zeta(s) ;
        p := 2;
        while p <= M do
            v := v*(1-1/p^s) ;
            p := nextprime(p) ;
        end do:
        v ;
    end proc:
    Hurw := proc(a)
            local T,p,x,L,i,Le,pre,preT,v,t,M ;
        T := 40 ;
        preT := 0.0 ;
        while true do
                1/p/(p+a) ;
                subs(p=1/x,%) ;
                exp(%) ;
                t := taylor(%,x=0,T) ;
                L := [] ;
                for i from 1 to T-1 do
                        L := [op(L),evalf(coeftayl(t,x=0,i))] ;
                end do:
                Le := EULERi(L) ;
            M := -a ;
                v := 1.0 ;
                pre := 0.0 ;
                for i from 2 to nops(Le) do
                        pre := log(v) ;
                        v := v*evalf(ZetaM(i,M))^op(i,Le) ;
                        v := evalf(v) ;
                end do:
            pre := (log(v)+pre)/2. ;
            printf("%.105f\n",%) ;
            if abs(1.0-preT/pre)  < 10^(-Digits/3) then
                break;
            end if;
            preT := pre ;
            T := T+10 ;
        end do:
            pre ;
    end proc:
    A179119 := proc()
        Hurw(1) ;
    end proc:
    A179119() ;
  • Mathematica
    digits = 101; S = NSum[(-1)^n PrimeZetaP[n], {n, 2, Infinity}, Method -> "AlternatingSigns", WorkingPrecision -> digits + 5]; RealDigits[S, 10, digits] // First (* Jean-François Alcover, Sep 11 2015 *)
  • PARI
    eps()=2.>>bitprecision(1.)
    primezeta(s)=my(t=s*log(2)); sum(k=1, lambertw(t/eps())\t, moebius(k)/k*log(abs(zeta(k*s))))
    sumalt(k=2,(-1)^k*primezeta(k)) \\ Charles R Greathouse IV, Aug 03 2016
    
  • PARI
    sumeulerrat(1/(p*(p+1))) \\ Amiram Eldar, Mar 18 2021
    

Formula

P(2) - P(3) + P(4) - P(5) + ..., where P is the prime zeta function. - Charles R Greathouse IV, Aug 03 2016

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

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

A143610 Numbers of the form p^2 * q^3, where p,q are distinct primes.

Original entry on oeis.org

72, 108, 200, 392, 500, 675, 968, 1125, 1323, 1352, 1372, 2312, 2888, 3087, 3267, 4232, 4563, 5324, 6125, 6728, 7688, 7803, 8575, 8788, 9747, 10952, 11979, 13448, 14283, 14792, 15125, 17672, 19652, 19773, 21125, 22472, 22707, 25947, 27436
Offset: 1

Views

Author

M. F. Hasler, Aug 27 2008

Keywords

Comments

Also: numbers with prime signature {3,2}.
This is a subsequence of A114128. [Hasler]
Every a(n) is an Achilles number (A052486). They are minimal, meaning no proper divisor is an Achilles number. - Antonio Roldán, Dec 27 2011

Examples

			The first three terms of this sequence are 3^2 * 2^3 = 72, 2^2 * 3^3 = 108, 5^2 * 2^3 = 200.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Sort[Last/@FactorInteger[n]] == {2, 3}; Select[Range[30000], f] (* Vladimir Joseph Stephan Orlovsky, Oct 09 2009 *)
  • PARI
    for(n=1, 10^5, omega(n)==2 || next; vecsort(factor(n)[,2])==[2,3]~ && print1(n","))
    
  • PARI
    list(lim)=my(v=List(),t);forprime(p=2, (lim\4)^(1/3), t=p^3;forprime(q=2, sqrt(lim\t), if(p==q, next);listput(v,t*q^2))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011
    
  • Python
    from math import isqrt
    from sympy import primepi, primerange, integer_nthroot
    def A143610(n):
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            kmin = kmax >> 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x): return n+x-sum(primepi(isqrt(x//p**3)) for p in primerange(integer_nthroot(x,3)[0]+1))+primepi(integer_nthroot(x,5)[0])
        return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Formula

Sum_{n>=1} 1/a(n) = P(2)*P(3) - P(5) = A085548 * A085541 - A085965 = 0.043280..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020
Showing 1-10 of 71 results. Next