A092759 -id:A092759 - cp's OEIS frontend

A030078 Cubes of primes.

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

Patrick De Geest

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

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

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.

T. D. Noe, Table of n, a(n) for n = 1..1000
Xavier Gourdon and Pascal Sebah, Some Constants from Number theory.
Eric Weisstein's World of Mathematics, Prime Power.
Wikipedia, p-group, Classification.
Index to sequences related to prime signature

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.
Cf. A036966, A085541, A088453, A157289, A258600.

a030078 = a000578 . a000040
a030078_list = map a000578 a000040_list -- Reinhard Zumkeller, May 26 2012

[p^3: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014

Array[Prime[ # ]^3&, 5! ] (* Vladimir Joseph Stephan Orlovsky, Sep 01 2008 *)

a(n)=prime(n)^3 \\ Charles R Greathouse IV, Mar 20 2013

from sympy import prime, primerange
def aupton(terms): return [p**3 for p in primerange(1, prime(terms)+1)]
print(aupton(35)) # Michael S. Branicky, Aug 27 2021

[p**3 for p in prime_range(100)] # Zerinvary Lajos, May 15 2007

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)

A101296 n has the a(n)-th distinct prime signature.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 6, 2, 4, 4, 7, 2, 6, 2, 6, 4, 4, 2, 8, 3, 4, 5, 6, 2, 9, 2, 10, 4, 4, 4, 11, 2, 4, 4, 8, 2, 9, 2, 6, 6, 4, 2, 12, 3, 6, 4, 6, 2, 8, 4, 8, 4, 4, 2, 13, 2, 4, 6, 14, 4, 9, 2, 6, 4, 9, 2, 15, 2, 4, 6, 6, 4, 9, 2, 12, 7, 4, 2, 13, 4, 4, 4, 8, 2, 13, 4, 6, 4, 4, 4, 16, 2, 6, 6, 11, 2, 9, 2, 8, 9, 4, 2, 15, 2, 9, 4, 12, 2, 9, 4, 6, 6, 4, 4, 17

Offset: 1

David Wasserman, Dec 21 2004

From Antti Karttunen, May 12 2017: (Start)
Restricted growth sequence transform of A046523, the least representative of each prime signature. Thus this partitions the natural numbers to the same equivalence classes as A046523, i.e., for all i, j: a(i) = a(j) <=> A046523(i) = A046523(j), and for that reason satisfies in that respect all the same conditions as A046523. For example, we have, for all i, j: if a(i) = a(j), then:
    A000005(i) = A000005(j), A008683(i) = A008683(j), A286605(i) = A286605(j).
So, this sequence (instead of A046523) can be used for finding sequences where a(n)'s value is dependent only on the prime signature of n, that is, only on the multiset of prime exponents in the factorization of n. (End)
This is also the restricted growth sequence transform of many other sequences, for example, that of A181819. See further comments there. - Antti Karttunen, Apr 30 2022

			From _David A. Corneth_, May 12 2017: (Start)
1 has prime signature (), the first distinct prime signature. Therefore, a(1) = 1.
2 has prime signature (1), the second distinct prime signature after (1). Therefore, a(2) = 2.
3 has prime signature (1), as does 2. Therefore, a(3) = a(2) = 2.
4 has prime signature (2), the third distinct prime signature after () and (1). Therefore, a(4) = 3. (End)
From _Antti Karttunen_, May 12 2017: (Start)
Construction of restricted growth sequences: In this case we start with a(1) = 1 for A046523(1) = 1, and thereafter, for all n > 1, we use the least so far unused natural number k for a(n) if A046523(n) has not been encountered before, otherwise [whenever A046523(n) = A046523(m), for some m < n], we set a(n) = a(m).
For n = 2, A046523(2) = 2, which has not been encountered before (first prime), thus we allot for a(2) the least so far unused number, which is 2, thus a(2) = 2.
For n = 3, A046523(2) = 2, which was already encountered as A046523(1), thus we set a(3) = a(2) = 2.
For n = 4, A046523(4) = 4, not encountered before (first square of prime), thus we allot for a(4) the least so far unused number, which is 3, thus a(4) = 3.
For n = 5, A046523(5) = 2, as for the first time encountered at n = 2, thus we set a(5) = a(2) = 2.
For n = 6, A046523(6) = 6, not encountered before (first semiprime pq with distinct p and q), thus we allot for a(6) the least so far unused number, which is 4, thus a(6) = 4.
For n = 8, A046523(8) = 8, not encountered before (first cube of a prime), thus we allot for a(8) the least so far unused number, which is 5, thus a(8) = 5.
For n = 9, A046523(9) = 4, as for the first time encountered at n = 4, thus a(9) = 3.
(End)
From _David A. Corneth_, May 12 2017: (Start)
(Rough) description of an algorithm of computing the sequence:
Suppose we want to compute a(n) for n in [1..20].
We set up a vector of 20 elements, values 0, and a number m = 1, the minimum number we haven't checked and c = 0, the number of distinct prime signatures we've found so far.
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
We check the prime signature of m and see that it's (). We increase c with 1 and set all elements up to 20 with prime signature () to 1. In the process, we adjust m. This gives:
[1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]. The least number we haven't checked is m = 2. 2 has prime signature (1). We increase c with 1 and set all elements up to 20 with prime signature (1) to 2. In the process, we adjust m. This gives:
[1, 2, 2, 0, 2, 0, 2, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 0, 2, 0]
We check the prime signature of m = 4 and see that its prime signature is (2). We increase c with 1 and set all numbers up to 20 with prime signature (2) to 3. This gives:
[1, 2, 2, 3, 2, 0, 2, 0, 3, 0, 2, 0, 2, 0, 0, 0, 2, 0, 2, 0]
Similarily, after m = 6, we get
[1, 2, 2, 3, 2, 4, 2, 0, 3, 4, 2, 0, 2, 4, 4, 0, 2, 0, 2, 0], after m = 8 we get:
[1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 0, 2, 4, 4, 0, 2, 0, 2, 0], after m = 12 we get:
[1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 6, 2, 4, 4, 0, 2, 6, 2, 0], after m = 16 we get:
[1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 6, 2, 4, 4, 7, 2, 6, 2, 0], after m = 20 we get:
[1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 6, 2, 4, 4, 7, 2, 6, 2, 8]. Now, m > 20 so we stop. (End)
The above method is inefficient, because the step "set all elements a(n) up to n = Nmax with prime signature s(n) = S[c] to c" requires factoring all integers up to Nmax (or at least comparing their signature, once computed, with S[c]) again and again. It is much more efficient to run only once over each m = 1..Nmax, compute its prime signature s(m), add it to an ordered list in case it did not occur earlier, together with its "rank" (= new size of the list), and assign that rank to a(m). The list of prime signatures is much shorter than [1..Nmax]. One can also use m'(m) := the smallest n with the prime signature of m (which is faster to compute than to search for the signature) as representative for s(m), and set a(m) := a(m'(m)). Then it is sufficient to have just one counter (number of prime signatures seen so far) as auxiliary variable, in addition to the sequence to be computed. - _M. F. Hasler_, Jul 18 2019

Michel Marcus (terms 1..10000) & Antti Karttunen, Table of n, a(n) for n = 1..100000

Cf. A025487, A046523, A064839 (ordinal transform of this sequence), A181819, and arrays A095904, A179216.
Cf. A000005, A008683.
Sequences that are unions of finite number (>= 2) of equivalence classes determined by the values that this sequence obtains (i.e., sequences mentioned in David A. Corneth's May 12 2017 formula): A001358 (A001248 U A006881, values 3 & 4), A007422 (values 1, 4, 5), A007964 (2, 3, 4, 5), A014612 (5, 6, 9), A030513 (4, 5), A037143 (1, 2, 3, 4), A037144 (1, 2, 3, 4, 5, 6, 9), A080258 (6, 7), A084116 (2, 4, 5), A167171 (2, 4), A217856 (6, 9).
Cf. also A077462, A305897 (stricter variants, with finer partitioning) and A254524, A286603, A286605, A286610, A286619, A286621, A286622, A286626, A286378 for other similarly constructed sequences.

A101296 := proc(n)
    local a046523, a;
    a046523 := A046523(n) ;
    for a from 1 do
        if A025487(a) = a046523 then
            return a;
        elif A025487(a) > a046523 then
            return -1 ;
        end if;
    end do:
end proc: # R. J. Mathar, May 26 2017

With[{nn = 120}, Function[s, Table[Position[Keys@s, k_ /; MemberQ[k, n]][[1, 1]], {n, nn}]]@ Map[#1 -> #2 & @@ # &, Transpose@ {Values@ #, Keys@ #}] &@ PositionIndex@ Table[Times @@ MapIndexed[Prime[First@ #2]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]] - Boole[n == 1], {n, nn}] ] (* Michael De Vlieger, May 12 2017, Version 10 *)

find(ps, vps) = {for (k=1, #vps, if (vps[k] == ps, return(k)););}
lisps(nn) = {vps = []; for (n=1, nn, ps = vecsort(factor(n)[,2]); ips = find(ps, vps); if (! ips, vps = concat(vps, ps); ips = #vps); print1(ips, ", "););} \\ Michel Marcus, Nov 15 2015; edited by M. F. Hasler, Jul 16 2019

rgs_transform(invec) = { my(occurrences = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(occurrences,invec[i]), my(pp = mapget(occurrences, invec[i])); outvec[i] = outvec[pp] , mapput(occurrences,invec[i],i); outvec[i] = u; u++ )); outvec; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
write_to_bfile(1,rgs_transform(vector(100000,n,A046523(n))),"b101296.txt");
\\ Antti Karttunen, May 12 2017

A025487(a(n)) = A046523(n).
Indices of records give A025487. - Michel Marcus, Nov 16 2015
From David A. Corneth, May 12 2017: (Start) [Corresponding characteristic function in brackets]
a(A000012(n)) = 1 (sig.: ()).      [A063524]
a(A000040(n)) = 2 (sig.: (1)).     [A010051]
a(A001248(n)) = 3 (sig.: (2)).     [A302048]
a(A006881(n)) = 4 (sig.: (1,1)).   [A280710]
a(A030078(n)) = 5 (sig.: (3)).
a(A054753(n)) = 6 (sig.: (1,2)).   [A353472]
a(A030514(n)) = 7 (sig.: (4)).
a(A065036(n)) = 8 (sig.: (1,3)).
a(A007304(n)) = 9 (sig.: (1,1,1)). [A354926]
a(A050997(n)) = 10 (sig.: (5)).
a(A085986(n)) = 11 (sig.: (2,2)).
a(A178739(n)) = 12 (sig.: (1,4)).
a(A085987(n)) = 13 (sig.: (1,1,2)).
a(A030516(n)) = 14 (sig.: (6)).
a(A143610(n)) = 15 (sig.: (2,3)).
a(A178740(n)) = 16 (sig.: (1,5)).
a(A189975(n)) = 17 (sig.: (1,1,3)).
a(A092759(n)) = 18 (sig.: (7)).
a(A189988(n)) = 19 (sig.: (2,4)).
a(A179643(n)) = 20 (sig.: (1,2,2)).
a(A189987(n)) = 21 (sig.: (1,6)).
a(A046386(n)) = 22 (sig.: (1,1,1,1)).
a(A162142(n)) = 23 (sig.: (2,2,2)).
a(A179644(n)) = 24 (sig.: (1,1,4)).
a(A179645(n)) = 25 (sig.: (8)).
a(A179646(n)) = 26 (sig.: (2,5)).
a(A163569(n)) = 27 (sig.: (1,2,3)).
a(A179664(n)) = 28 (sig.: (1,7)).
a(A189982(n)) = 29 (sig.: (1,1,1,2)).
a(A179666(n)) = 30 (sig.: (3,4)).
a(A179667(n)) = 31 (sig.: (1,1,5)).
a(A179665(n)) = 32 (sig.: (9)).
a(A189990(n)) = 33 (sig.: (2,6)).
a(A179669(n)) = 34 (sig.: (1,2,4)).
a(A179668(n)) = 35 (sig.: (1,8)).
a(A179670(n)) = 36 (sig.: (1,1,1,3)).
a(A179671(n)) = 37 (sig.: (3,5)).
a(A162143(n)) = 38 (sig.: (2,2,2)).
a(A179672(n)) = 39 (sig.: (1,1,6)).
a(A030629(n)) = 40 (sig.: (10)).
a(A179688(n)) = 41 (sig.: (1,3,3)).
a(A179689(n)) = 42 (sig.: (2,7)).
a(A179690(n)) = 43 (sig.: (1,1,2,2)).
a(A189991(n)) = 44 (sig.: (4,4)).
a(A179691(n)) = 45 (sig.: (1,2,5)).
a(A179692(n)) = 46 (sig.: (1,9)).
a(A179693(n)) = 47 (sig.: (1,1,1,4)).
a(A179694(n)) = 48 (sig.: (3,6)).
a(A179695(n)) = 49 (sig.: (2,2,3)).
a(A179696(n)) = 50 (sig.: (1,1,7)).
(End)

Data section extended to 120 terms by Antti Karttunen, May 12 2017

Minor edits/corrections by M. F. Hasler, Jul 18 2019

A030626 Numbers with exactly 8 divisors.

Original entry on oeis.org

24, 30, 40, 42, 54, 56, 66, 70, 78, 88, 102, 104, 105, 110, 114, 128, 130, 135, 136, 138, 152, 154, 165, 170, 174, 182, 184, 186, 189, 190, 195, 222, 230, 231, 232, 238, 246, 248, 250, 255, 258, 266, 273, 282, 285, 286, 290, 296, 297, 310, 318, 322, 328, 344, 345, 351, 354, 357, 366, 370, 374, 375, 376, 385, 399, 402

Offset: 1

Jeff Burch

nonn

Since A119479(8)=7, there are never more than 7 consecutive terms. Runs of 7 consecutive terms start at 171897, 180969, 647385, ... (subsequence of A049053). - Ivan Neretin, Feb 08 2016

Jon E. Schoenfield, Table of n, a(n) for n = 1..10000 (first 1000 terms from R. J. Mathar)
Jérôme Germoni, Nombres à huit diviseurs, Images des Mathématiques, CNRS, 2017 (in French).
Eric Weisstein's World of Mathematics, Divisor Product.

Essentially the same as A111398.

Cf. A000005, A007304, A049053, A065036, A092759, A119479.

[n: n in [1..400] | DivisorSigma(0, n) eq 8]; // Vincenzo Librandi, Oct 05 2017

select(numtheory:-tau=8, [$1..1000]); # Robert Israel, Dec 17 2014

Select[Range[400], DivisorSigma[0, #]== 8 &] (* Vincenzo Librandi, Oct 05 2017 *)

Vec(select(x->x==8,vector(500, i, numdiv(i)),1)) \\ Michel Marcus, Dec 17 2014

from sympy import divisor_count
isok = lambda n: divisor_count(n) == 8
print([n for n in range(1, 400) if isok(n)]) # Darío Clavijo, Oct 17 2023

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A030626(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 int(n+x-sum(primepi(x//(k*m))-b for a,k in enumerate(primerange(integer_nthroot(x,3)[0]+1),1) for b,m in enumerate(primerange(k+1,isqrt(x//k)+1),a+1))-sum(primepi(x//p**3) for p in primerange(integer_nthroot(x,3)[0]+1))+primepi(integer_nthroot(x,4)[0])-primepi(integer_nthroot(x,7)[0]))
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

A000005(a(n))=8. - Juri-Stepan Gerasimov, Oct 10 2009

Equals A065036 (p*q^3) U A007304 (p*q*r) U A092759 (p^7). - Amarnath Murthy, Apr 21 2001

A179665 a(n) = prime(n)^9.

Original entry on oeis.org

512, 19683, 1953125, 40353607, 2357947691, 10604499373, 118587876497, 322687697779, 1801152661463, 14507145975869, 26439622160671, 129961739795077, 327381934393961, 502592611936843, 1119130473102767

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 23 2010

Product_{n >= 2, m_n = (a(n) mod 4) - 2} ((a(n) + 1) / (a(n) - 1))^m_n = 209865342976 / 209844223875. - Dimitris Valianatos, May 13 2020

			a(1) = 512 since the ninth power of the first prime is 2^9 = 512. - _Wesley Ivan Hurt_, Mar 27 2014

T. D. Noe, Table of n, a(n) for n = 1..1000
Xavier Gourdon and Pascal Sebah, Some Constants from Number theory.
Will Nicholes, Prime Signatures.

Cf. A001248, A030078, A030514, A050997, A030516, A092759, A179645.

Cf. A013667, A013676, A085969.

[p^9: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014

A179665:=n->ithprime(n)^9; seq(A179665(n), n=1..30); # Wesley Ivan Hurt, Mar 27 2014

Array[Prime[ # ]^9&, 30]
Prime[Range[30]]^9 (* Harvey P. Dale, Jul 20 2024 *)

a(n)=prime(n)^9 \\ Charles R Greathouse IV, Jul 20 2011

a(n) = A000040(n)^9 = A001017(A000040(n)). - Wesley Ivan Hurt, Mar 27 2014
Sum_{n>=1} 1/a(n) = P(9) = 0.0020044675... (A085969). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 24 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(9)/zeta(18) = A013667/A013676.
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(9) = 1/A013667. (End)

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Original entry on oeis.org

1, 8, 27, 32, 125, 128, 216, 243, 343, 512, 864, 1000, 1331, 1944, 2048, 2187, 2197, 2744, 3125, 3375, 3456, 4000, 4913, 6859, 7776, 8192, 9261, 10648, 10976, 12167, 13824, 16000, 16807, 17496, 17576, 19683, 24389, 25000, 27000, 29791, 30375, 31104, 32768, 35937

Offset: 1

Amiram Eldar, Jul 03 2020

nonn

This sequence is a permutation of A355038.
This sequence is also a permutation of the exponentially odd numbers (A268335) multiplied by the square of their squarefree kernel (A007947).
a(n)/rad(a(n)) is a permutation of the squares.
a(n)/rad(a(n))^2 is a permutation of the exponentially odd numbers.

			8 = 2^3 is a term since the exponent of its prime factor 2 is 3 which is odd and larger than 1.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Intersection of A001694 and A268335.
Intersection of A036966 and A268335.
A355038 in ascending order.
A030078, A050997, A092759, A179665, A079395 and A138031 are subsequences.
Cf. A007947, A065487.

Join[{1}, Select[Range[10^5], AllTrue[Last /@ FactorInteger[#], #1 > 1 && OddQ[#1] &] &]]

from math import isqrt, prod
from sympy import factorint
def afind(N): # all terms up to limit N
    cands = (n**2*prod(factorint(n**2)) for n in range(1, isqrt(N//2)+2))
    return sorted(c for c in cands if c <= N)
print(afind(4*10**4)) # Michael S. Branicky, Jun 16 2022

Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/(p*(p^2-1))) = 1.2312911... (A065487).

A085967 Decimal expansion of the prime zeta function at 7.

Original entry on oeis.org

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

Offset: 0

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

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

			0.0082838328561335925351...

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.

Jason Kimberley, Table of n, a(n) for n = 0..1902
Henri Cohen, High Precision Computation of Hardy-Littlewood Constants, Preprint, 1998.
Henri Cohen, High-precision computation of Hardy-Littlewood constants. [pdf copy, with permission]
X. Gourdon and P. Sebah, Some Constants from Number theory
R. J. Mathar, Series of reciprocal powers of k-almost primes, arXiv:0803.0900 [math.NT], 2008-2009. Table 1.
Eric Weisstein's World of Mathematics, Prime Zeta Function

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

Cf. A013665, A092759.

R := RealField(106);
PrimeZeta := func;
[0] cat Reverse(IntegerToSequence(Floor(PrimeZeta(7,47)*10^105)));
// Jason Kimberley, Dec 30 2016

s[n_] := s[n] = Join[{0, 0}, Sum[ MoebiusMu[k]*Log[Zeta[7*k]]/k, {k, 1, n}] // RealDigits[#, 10, 104] & // First]; 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[ 7], 10, 111][[1]] (* Robert G. Wilson v, Sep 03 2014 *)

sumeulerrat(1/p,7) \\ Hugo Pfoertner, Feb 03 2020

P(7) = Sum_{p prime} 1/p^7 = Sum_{n>=1} mobius(n)*log(zeta(7*n))/n.

Equals Sum_{k>=1} 1/A092759(k). - Amiram Eldar, Jul 27 2020

A189987 Numbers with prime factorization p*q^6.

Original entry on oeis.org

192, 320, 448, 704, 832, 1088, 1216, 1458, 1472, 1856, 1984, 2368, 2624, 2752, 3008, 3392, 3645, 3776, 3904, 4288, 4544, 4672, 5056, 5103, 5312, 5696, 6208, 6464, 6592, 6848, 6976, 7232, 8019, 8128, 8384, 8768, 8896, 9477, 9536, 9664, 10048, 10432, 10688

Offset: 1

Vladimir Joseph Stephan Orlovsky, May 03 2011

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures
Index to sequences related to prime signature

Cf. A030632, A092759.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,6}; Select[Range[30000],f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\2)^(1/6), t=p^6;forprime(q=2, lim\t, if(p==q, next);listput(v,t*q))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from sympy import primepi, primerange, integer_nthroot
def A189987(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(x//p**6) for p in primerange(integer_nthroot(x,6)[0]+1))+primepi(integer_nthroot(x,7)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 22 2025

A229125 Numbers of the form p * m^2, where p is prime and m > 0: union of A228056 and A000040.

Original entry on oeis.org

2, 3, 5, 7, 8, 11, 12, 13, 17, 18, 19, 20, 23, 27, 28, 29, 31, 32, 37, 41, 43, 44, 45, 47, 48, 50, 52, 53, 59, 61, 63, 67, 68, 71, 72, 73, 75, 76, 79, 80, 83, 89, 92, 97, 98, 99, 101, 103, 107, 108, 109, 112, 113, 116, 117, 124, 125, 127, 128, 131, 137, 139, 147, 148, 149

Offset: 1

Chris Boyd, Sep 14 2013

nonn

No term is the product of two other terms.
Squares of terms and pairwise products of distinct terms form a subsequence of A028260.
Numbers n such that A162642(n) = 1. - Jason Kimberley, Oct 10 2016
Numbers k such that A007913(k) is a prime number. - Amiram Eldar, Jul 27 2020

Chris Boyd, Table of n, a(n) for n = 1..10000
Eckford Cohen, Arithmetical notes, IX. On the set of integers representable as a product of a prime and square, Acta Arithmetica, Vol. 7 (1962), pp. 417-420.

Subsequence of A026424.
Cf. A028260, A162642, A229153.
Cf. A007913, A013661.
Subsequences: A000040, A030078, A050997, A054753, A092759, A179643, A179665, A246551.

With[{nn=70},Take[Union[Flatten[Table[p*m^2,{p,Prime[Range[nn]]},{m,nn}]]], nn]] (* Harvey P. Dale, Dec 02 2014 *)

test(n)=isprime(core(n))
for(n=1,200,if(test(n), print1(n",")))

from math import isqrt
from sympy import primepi
def A229125(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(x//y**2) for y in range(1,isqrt(x)+1))
    return bisection(f,n,n) # Chai Wah Wu, Jan 30 2025

The number of terms not exceeding x is (Pi^2/6) * x/log(x) + O(x/(log(x))^2) (Cohen, 1962). - Amiram Eldar, Jul 27 2020

cp's OEIS Frontend

A030078 Cubes of primes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Python

Sage

Formula

A101296 n has the a(n)-th distinct prime signature.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A030626 Numbers with exactly 8 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Python

Formula

A179665 a(n) = prime(n)^9.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A085967 Decimal expansion of the prime zeta function at 7.

Views

Author

Keywords

Comments

Examples