A030629 -id:A030629 - 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

A030513 Numbers with 4 divisors.

Original entry on oeis.org

6, 8, 10, 14, 15, 21, 22, 26, 27, 33, 34, 35, 38, 39, 46, 51, 55, 57, 58, 62, 65, 69, 74, 77, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 125, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 177, 178, 183, 185, 187

Offset: 1

Jeff Burch

Essentially the same as A007422.
Numbers which are either the product of two distinct primes (A006881) or the cube of a prime (A030078).
4*a(n) are the solutions to A048272(x) = Sum_{d|x} (-1)^d = 4. - Benoit Cloitre, Apr 14 2002
Since A119479(4)=3, there are never more than 3 consecutive integers in the sequence. Triples of consecutive integers start at 33, 85, 93, 141, 201, ... (A039833). No such triple contains a term of the form p^3. - Ivan Neretin, Feb 08 2016
Numbers that are equal to the product of their proper divisors (A007956) (proof in Sierpiński). - Bernard Schott, Apr 04 2022

Wacław Sierpiński, Elementary Theory of Numbers, Ex. 2 p. 174, Warsaw, 1964.

Cf. A000005, A007422, A007956, A030515, A035533, A048272, A119479.

Equals the disjoint union of A006881 and A030078.

[n: n in [1..200] | DivisorSigma(0, n) eq 4]; // Vincenzo Librandi, Jul 16 2015

Select[Range[200], DivisorSigma[0,#]==4&] (* Harvey P. Dale, Apr 06 2011 *)

is(n)=numdiv(n)==4 \\ Charles R Greathouse IV, May 18 2015

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A030513(n):
    def f(x): return int(n+x-primepi(integer_nthroot(x,3)[0])+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Aug 16 2024

{n : A000005(n) = 4}. - Juri-Stepan Gerasimov, Oct 10 2009

Incorrect comments removed by Charles R Greathouse IV, Mar 18 2010

A030631 Numbers with 13 divisors.

Original entry on oeis.org

4096, 531441, 244140625, 13841287201, 3138428376721, 23298085122481, 582622237229761, 2213314919066161, 21914624432020321, 353814783205469041, 787662783788549761, 6582952005840035281, 22563490300366186081, 39959630797262576401, 116191483108948578241

Offset: 1

Jeff Burch

12th powers of primes. The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. - Omar E. Pol, May 06 2008

R. J. Mathar, Table of n, a(n) for n = 1..1000

Cf. A000040, A013670, A030629.

Subsequence of A008456.

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

(Prime@Range@30)^12  (* Vladimir Joseph Stephan Orlovsky, Apr 11 2011 *)

is(n)=isprimepower(n)==12 \\ Charles R Greathouse IV, Jun 19 2016

a(n) = A000040(n)^(13-1) = A000040(n)^(12). - Omar E. Pol, May 06 2008
From Amiram Eldar, Jan 24 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(12)/zeta(24) = 218517792968475/(236364091*Pi^12).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(12) = 638512875/(691*Pi^12) = 1/A013670. (End)

A139572 Numbers with 37 divisors.

Original entry on oeis.org

68719476736, 150094635296999121, 14551915228366851806640625, 2651730845859653471779023381601, 30912680532870672635673352936887453361, 12646218552730347184269489080961456410641

Offset: 1

Omar E. Pol, May 07 2008

36th powers of primes. The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime.

T. D. Noe, Table of n, a(n) for n = 1..1000
OEIS Wiki, Index entries for number of divisors

Cf. A000005, A000040, A001248, A030514, A030516, A030629, A030631, A030635, A030637, A137486, A137492, A139571, A139573.

A139572 := proc(n) ithprime(n)^36 ; end proc: seq(A139572(n),n=1..10) ; # R. J. Mathar, Feb 05 2010

Prime[Range[10]]^36 (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)

a(n)=prime(n)^36 \\ Charles R Greathouse IV, Jun 19 2016

a(n) = A000040(n)^(37-1) = A000040(n)^36.

More terms from R. J. Mathar, Feb 05 2010

A139573 Numbers with 41 divisors.

Original entry on oeis.org

1099511627776, 12157665459056928801, 9094947017729282379150390625, 6366805760909027985741435139224001, 452592555681759518058893560348969204658401

Offset: 1

Omar E. Pol, May 07 2008

40th powers of primes. The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime.

T. D. Noe, Table of n, a(n) for n = 1..1000
OEIS Wiki, Index entries for number of divisors

Cf. A000005, A000040, A001248, A030514, A030516, A030629, A030631, A030635, A030637, A137486, A137492, A139571, A139572.

Prime[Range[9]]^40 (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)

a(n)=prime(n)^40 \\ Charles R Greathouse IV, Jun 19 2016

a(n)=A000040(n)^(41-1)=A000040(n)^40.

More terms from Jon E. Schoenfield, May 18 2010

A030628 1 together with numbers of the form p*q^4 and p^9, where p and q are distinct primes.

Original entry on oeis.org

1, 48, 80, 112, 162, 176, 208, 272, 304, 368, 405, 464, 496, 512, 567, 592, 656, 688, 752, 848, 891, 944, 976, 1053, 1072, 1136, 1168, 1250, 1264, 1328, 1377, 1424, 1539, 1552, 1616, 1648, 1712, 1744, 1808, 1863, 1875, 2032, 2096, 2192, 2224, 2349, 2384

Offset: 1

Jeff Burch

Also 1 together with numbers with 10 divisors. Also numbers n such that product of all proper divisors of n equals n^4.

If M(n) denotes the product of all divisors of n, then n is said to be k-multiplicatively perfect if M(n)=n^k. All such numbers are of the form p*q^(k-1) or p^(2k-1). The sequence A030628 is therefore 5-multiplicatively perfect. See the Links for A007422. - Walter Kehowski, Sep 13 2005

D. M. Burton, Elementary Number Theory, Allyn and Bacon Inc., Boston, MA, 1976. p. 119.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, entry for 48, page 106, 1997.

R. J. Mathar, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Divisor Product
Index to sequences related to prime signature

Cf. A030515, A030627, A030629.

with(numtheory): k:=5: MPL:=[]: for z from 1 to 1 do for n from 1 to 5000 do if convert(divisors(n),`*`) = n^k then MPL:=[op(MPL),n] fi od; od; MPL; # Walter Kehowski, Sep 13 2005

Join[{1},Select[Range[6000],DivisorSigma[0,#]==10&]] (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)
Select[Range[2500],Times@@Most[Divisors[#]]==#^4&] (* Harvey P. Dale, Nov 04 2024 *)

{v=[]; for(n=1,500,v=concat(v, if(numdiv(n)==10,n,",")); ); v} \\ Jason Earls, Jun 18 2001

list(lim)=my(v=List([1]), t); forprime(p=2, (lim\2+.5)^(1/4), t=p^4; forprime(q=2, lim\t, if(p==q, next); listput(v, t*q))); forprime(p=2,(lim+.5)^(1/9),listput(v,p^9)); vecsort(Vec(v)) \\ Charles R Greathouse IV, Apr 26 2012

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

Union A178739 U A179665 {1}. - R. J. Mathar, Apr 03 2011

Better description from Sharon Sela (sharonsela(AT)hotmail.com), Dec 23 2001

More terms from Walter Kehowski, Sep 13 2005

A067004 Number of numbers <= n with same number of divisors as n.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 4, 2, 2, 3, 5, 1, 6, 4, 5, 1, 7, 2, 8, 3, 6, 7, 9, 1, 3, 8, 9, 4, 10, 2, 11, 5, 10, 11, 12, 1, 12, 13, 14, 3, 13, 4, 14, 6, 7, 15, 15, 1, 4, 8, 16, 9, 16, 5, 17, 6, 18, 19, 17, 1, 18, 20, 10, 1, 21, 7, 19, 11, 22, 8, 20, 2, 21, 23, 12, 13, 24, 9, 22, 2, 2, 25, 23, 3, 26, 27

Offset: 1

Henry Bottomley, Dec 21 2001

nonn

			a(10)=3 since 6,8,10 each have four divisors. a(11)=5 since 2,3,5,7,11 each have two divisors.

Paul Tek, Table of n, a(n) for n = 1..10000

Cf. A000005, A008479, A058933, A067003. Inverse of A000040, A001248, A030513, A030514, A030515, A030516, A030626, A030627, A030628, A030629, A030630, A030631, A030632, A030633, A030634, A030635, A030636, A030637, A030638 etc.
Cf. also A079788, A138009.
A047983(n) = a(n)-1.

N:= 1000: # to get a(1) to a(N)
R:= Vector(N):
for n from 1 to N do
  v:= numtheory:-tau(n);
  R[v]:= R[v]+1;
  A[n]:= R[v];
od:
seq(A[n],n=1..N); # Robert Israel, May 04 2015

b[_] = 0;
a[n_] := a[n] = With[{t = DivisorSigma[0, n]}, b[t] = b[t]+1];
Array[a, 105] (* Jean-François Alcover, Dec 20 2021 *)

a(n)=my(d=numdiv(n)); sum(k=1,n,numdiv(k)==d) \\ Charles R Greathouse IV, Sep 02 2015

Ordinal transform of A000005. - Franklin T. Adams-Watters, Aug 28 2006

a(A000040(n)^(p-1)) = n if p is prime. - Robert Israel, May 04 2015

A173533 Numbers with 53 divisors.

Original entry on oeis.org

4503599627370496, 6461081889226673298932241, 2220446049250313080847263336181640625, 88124787089723195184393736687912818113311201, 1420429319844313329730664601483335671261683881745483121, 8415003868347247618489696679505181495471801448798649088081

Offset: 1

Omar E. Pol, Oct 16 2010

52nd powers of primes.

The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime.

Omar E. Pol, Determinacion geometrica de los numeros primos y perfectos

Cf. A000005, A000040, A001248, A030514, A030516, A030629, A030631, A030635, A030637, A137486, A137492, A139571, A139572, A139573, A139574, A139575.

Prime[Range[9]]^52 (* Vladimir Joseph Stephan Orlovsky, May 06 2011 *)

a(n)=prime(n)^52 \\ Charles R Greathouse IV, Jun 19 2016

a(n) = A000040(n)^(53-1) = A000040(n)^52.

A182944 Square array A(i,j), i >= 1, j >= 1, of prime powers prime(i)^j, by descending antidiagonals.

Original entry on oeis.org

2, 4, 3, 8, 9, 5, 16, 27, 25, 7, 32, 81, 125, 49, 11, 64, 243, 625, 343, 121, 13, 128, 729, 3125, 2401, 1331, 169, 17, 256, 2187, 15625, 16807, 14641, 2197, 289, 19, 512, 6561, 78125, 117649, 161051, 28561, 4913, 361, 23

Offset: 1

Clark Kimberling, Dec 14 2010

We alternatively refer to this sequence as a triangle T(.,.), with T(n,k) = A(k,n-k+1) = prime(k)^(n-k+1).
The monotonic ordering of this sequence, prefixed by 1, is A000961.
The joint-rank array of this sequence is A182869.
Main diagonal gives A062457. - Omar E. Pol, Sep 11 2018

			Square array A(i,j) begins:
  i \ j: 1      2      3      4      5  ...
  ---\-------------------------------------
  1:     2,     4,     8,    16,    32, ...
  2:     3,     9,    27,    81,   243, ...
  3:     5,    25,   125,   625,  3125, ...
  4:     7,    49,   343,  2401, 16807, ...
  ...
The triangle T(n,k) begins:
  n\k:  1     2     3     4     5     6  ...
  1:    2
  2:    4     3
  3:    8     9     5
  4:   16    27    25     7
  5:   32    81   125    49    11
  6:   64   243   625   343   121    13
  ...

Michael De Vlieger, Table of n, a(n) for n = 1..11325 (rows 1..150, flattened)
Michael De Vlieger, Diagram showing row n of triangle in a semicircle as noted, with a color function associated with the magnitude of T(n,k) compared to 2^n in light blue, where prime(n) is the smallest and the prime power indicated in red the largest in the row.

Cf. A000961, A006939 (row products of triangle), A062457, A182945, A332979 (row maxima of triangle).
Columns: A000040 (1), A001248 (2), A030078 (3), A030514 (4), A050997 (5), A030516 (6), A092759 (7), A179645 (8), A179665 (9), A030629 (10).
A319075 extends the array with 0th powers.
Subtable of A242378, A284457, A329332.

TableForm[Table[Prime[n]^j,{n,1,14},{j,1,8}]]

From Peter Munn, Dec 29 2019: (Start)
A(i,j) = A182945(j,i) = A319075(j,i).
A(i,j) = A242378(i-1,2^j) = A329332(2^(i-1),j).
A(i,i) = A062457(i).
(End)

Clarified in respect of alternate reading as a triangle by Peter Munn, Aug 28 2022

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

A030513 Numbers with 4 divisors.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Formula

Extensions

A030631 Numbers with 13 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A139572 Numbers with 37 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A139573 Numbers with 41 divisors.

Views

Author

Keywords

Comments

Links

Crossrefs