A164336 -id:A164336 - cp's OEIS frontend

A164337 Those powers of primes that are not in sequence A164336.

64, 729, 1024, 4096, 15625, 16384, 32768, 59049, 117649, 262144, 531441, 1048576, 1771561, 2097152, 4194304, 4782969, 4826809, 9765625, 14348907, 16777216, 24137569, 47045881, 67108864, 148035889, 244140625, 268435456, 282475249

Offset: 1

Leroy Quet, Aug 13 2009

nonn

Michael De Vlieger, Table of n, a(n) for n = 1..10000

Cf. A164336.

Block[{a = {1}, k = 50, b, nn = 10^9}, Do[If[Length@ # == 1 && MemberQ[a, First@ #], AppendTo[a, i]] &[FactorInteger[i][[All, -1]]], {i, 2, Prime@ k}]; b = Complement[Range@ Prime@ k, a]; Union@ Flatten@ MapIndexed[Prime[First@ #2]^#1 &, DeleteCases[#, {}]] &@ Map[Function[p, TakeWhile[b, # <= Floor@ Log[p, nn] &]], Prime@ Range@ k]] (* Michael De Vlieger, Aug 31 2017 *)

More terms from Hagen von Eitzen, Oct 03 2009

A007916 Numbers that are not perfect powers.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 82, 83

Offset: 1

R. Muller

From Gus Wiseman, Oct 23 2016: (Start)
There is a 1-to-1 correspondence between integers N >= 2 and sequences a(x_1),a(x_2),...,a(x_k) of terms from this sequence. Every N >= 2 can be written uniquely as a "power tower"
N = a(x_1)^a(x_2)^a(x_3)^...^a(x_k),
where the exponents are to be nested from the right.
Proof: If N is not a perfect power then N = a(x) for some x, and we are done. Otherwise, write N = a(x_1)^M for some M >=2, and repeat the process. QED
Of course, prime numbers also have distinct power towers (see A164336). (End)
These numbers can be computed with a modified Sieve of Eratosthenes: (1) start at n=2; (2) if n is not crossed out, then append n to the sequence and cross out all powers of n; (3) set n = n+1 and go to step 2. - Sam Alexander, Dec 15 2003
These are all numbers such that the multiplicities of the prime factors have no common divisor. The first number in the sequence whose prime multiplicities are not coprime is 180 = 2 * 2 * 3 * 3 * 5. Mathematica: CoprimeQ[2,2,1]->False. - Gus Wiseman, Jan 14 2017

			Example of the power tower factorizations for the first nine positive integers: 1=1, 2=a(1), 3=a(2), 4=a(1)^a(1), 5=a(3), 6=a(4), 7=a(5), 8=a(1)^a(2), 9=a(2)^a(1). - _Gus Wiseman_, Oct 20 2016

N. J. A. Sloane, Table of n, a(n) for n = 1..9875
Joakim Munkhammar, The Riemann zeta function as a sum of geometric series, The Mathematical Gazette (2020) Vol. 104, Issue 561, 527-530.
N. J. A. Sloane, Maple programs for A007916, A278028, A278029, A052409, A089723, A277564
F. Smarandache, Only Problems, Not Solutions!, Xiquan Publ., Phoenix-Chicago, 1993
Index entries for sequences generated by sieves

Complement of A001597. Union of A052485 and A052486.
Cf. A144338, A277562, A277564, A075802.
Cf. A153158 (squares of these numbers).
See A277562, A277564, A277576, A277615 for more about the power towers.
A278029 is a left inverse.
Cf. A052409.

a007916 n = a007916_list !! (n-1)
a007916_list = filter ((== 1) . foldl1 gcd . a124010_row) [2..]
-- Reinhard Zumkeller, Apr 13 2012

[n : n in [2..1000] | not IsPower(n) ];

Maple
```
See link.
```

a = {}; Do[If[Apply[GCD, Transpose[FactorInteger[n]][[2]]] == 1, a = Append[a, n]], {n, 2, 200}];
Select[Range[2,200],GCD@@FactorInteger[#][[All,-1]]===1&] (* Michael De Vlieger, Oct 21 2016. Corrected by Gus Wiseman, Jan 14 2017 *)

is(n)=!ispower(n)&&n>1 \\ Charles R Greathouse IV, Jul 01 2013

from sympy import mobius, integer_nthroot
def A007916(n):
    def f(x): return int(n+1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Aug 13 2024

A075802(a(n)) = 0. - Reinhard Zumkeller, Mar 19 2009
Gcd(exponents in prime factorization of a(n)) = 1, cf. A124010. - Reinhard Zumkeller, Apr 13 2012
a(n) ~ n. - Charles R Greathouse IV, Jul 01 2013
A052409(a(n)) = 1. - Ridouane Oudra, Nov 23 2024

More terms from Henry Bottomley, Sep 12 2000
Edited by Charles R Greathouse IV, Mar 18 2010
Further edited by N. J. A. Sloane, Nov 09 2016

A050361 Number of factorizations into distinct prime powers greater than 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1

Offset: 1

Christian G. Bower, Oct 15 1999

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).

The number of unordered factorizations of n into 1 and exponentially odd prime powers, i.e., p^e where p is a prime and e is odd (A246551). - Amiram Eldar, Jun 12 2025

			From _Gus Wiseman_, Jul 30 2022: (Start)
The A000688(216) = 9 factorizations of 216 into prime powers are:
  (2*2*2*3*3*3)
  (2*2*2*3*9)
  (2*2*2*27)
  (2*3*3*3*4)
  (2*3*4*9)
  (2*4*27)
  (3*3*3*8)
  (3*8*9)
  (8*27)
Of these, the a(216) = 4 strict cases are:
  (2*3*4*9)
  (2*4*27)
  (3*8*9)
  (8*27)
(End)

Antti Karttunen, Table of n, a(n) for n = 1..100000 (first 10000 terms from Reinhard Zumkeller)
Index entries for sequences computed from exponents in factorization of n.

Cf. A000009, A050360, A050362, A050363, A050364, A101296, A246551.
Cf. A124010.
This is the strict case of A000688.
Positions of 1's are A004709, complement A046099.
The case of primes (instead of prime-powers) is A008966, non-strict A000012.
The non-strict additive version allowing 1's A023893, ranked by A302492.
The non-strict additive version is A023894, ranked by A355743.
The additive version (partitions) is A054685, ranked by A356065.
The additive version allowing 1's is A106244, ranked by A302496.
A001222 counts prime-power divisors.
A005117 lists all squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists all prime-powers (A000961 includes 1), towers A164336.
A296131 counts twice-factorizations of type PQR, non-strict A295935.
Cf. A001970, A002110, A025487, A055887, A063834, A076610, A085970, A279786, A302590, A302601, A354911, A355742.

a050361 = product . map a000009 . a124010_row
-- Reinhard Zumkeller, Aug 28 2014

A050361 := proc(n)
    local a,f;
    if n = 1 then
        1;
    else
        a := 1 ;
        for f in ifactors(n)[2] do
            a := a*A000009(op(2,f)) ;
        end do:
    end if;
end proc: # R. J. Mathar, May 25 2017

Table[Times @@ PartitionsQ[Last /@ FactorInteger[n]], {n, 99}] (* Arkadiusz Wesolowski, Feb 27 2017 *)

A000009(n,k=(n-!(n%2))) = if(!n,1,my(s=0); while(k >= 1, if(k<=n, s += A000009(n-k,k)); k -= 2); (s));
A050361(n) = factorback(apply(A000009,factor(n)[,2])); \\ Antti Karttunen, Nov 17 2019

Dirichlet g.f.: Product_{n is a prime power >1}(1 + 1/n^s).
Multiplicative with a(p^e) = A000009(e).
a(A002110(k))=1.
a(n) = A050362(A101296(n)). - R. J. Mathar, May 26 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.26020571070524171076..., where f(x) = (1-x) * Product_{k>=1} (1 + x^k). - Amiram Eldar, Oct 03 2023

A023893 Number of partitions of n into prime power parts (1 included); number of nonisomorphic Abelian subgroups of symmetric group S_n.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 14, 20, 27, 36, 48, 63, 82, 105, 134, 171, 215, 269, 335, 415, 511, 626, 764, 929, 1125, 1356, 1631, 1953, 2333, 2776, 3296, 3903, 4608, 5427, 6377, 7476, 8744, 10205, 11886, 13818, 16032, 18565, 21463, 24768, 28536

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(6) = 10 partitions:
  ()  (1)  (2)   (3)    (4)     (5)      (33)
           (11)  (21)   (22)    (32)     (42)
                 (111)  (31)    (41)     (51)
                        (211)   (221)    (222)
                        (1111)  (311)    (321)
                                (2111)   (411)
                                (11111)  (2211)
                                         (3111)
                                         (21111)
                                         (111111)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Index entries for sequences related to groups

Cf. A009490, A023894 (first differences), A062297 (number of Abelian subgroups).
Cf. A018819, A062051, A131995.
The multiplicative version (factorizations) is A000688.
Not allowing 1's gives A023894, strict A054685, ranked by A355743.
The version for just primes (not prime-powers) is A034891, strict A036497.
The strict version is A106244.
These partitions are ranked by A302492.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970, A279784, A295935, A356065.

Table[Length[Select[IntegerPartitions[n],Count[Map[Length,FactorInteger[#]], 1] == Length[#] &]], {n, 0, 35}] (* Geoffrey Critzer, Oct 25 2015 *)
nmax = 50; Clear[P]; P[m_] := P[m] = Product[Product[1/(1-x^(p^k)), {k, 1, m}], {p, Prime[Range[PrimePi[nmax]]]}]/(1-x)+O[x]^nmax // CoefficientList[ #, x]&; P[1]; P[m=2]; While[P[m] != P[m-1], m++]; P[m] (* Jean-François Alcover, Aug 31 2016 *)

lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (isprimepower(k), 1/(1-x^k), 1))/(1-x); for (n=0, m, print1(polcoeff(gf, n, t), ", "));} \\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023893(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())+1
    return (c(n)+sum(c(k)*A023893(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: (Product_{p prime} Product_{k>=1} 1/(1-x^(p^k))) / (1-x).

A023894 Number of partitions of n into prime power parts (1 excluded).

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 3, 4, 6, 7, 9, 12, 15, 19, 23, 29, 37, 44, 54, 66, 80, 96, 115, 138, 165, 196, 231, 275, 322, 380, 443, 520, 607, 705, 819, 950, 1099, 1268, 1461, 1681, 1932, 2214, 2533, 2898, 3305, 3768, 4285, 4872, 5530, 6267, 7094, 8022, 9060

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(9) = 7 partitions:
  ()  .  (2)  (3)  (4)   (5)   (33)   (7)    (8)     (9)
                   (22)  (32)  (42)   (43)   (44)    (54)
                               (222)  (52)   (53)    (72)
                                      (322)  (332)   (333)
                                             (422)   (432)
                                             (2222)  (522)
                                                     (3222)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
E. Grosswald, Partitions into prime powers

The multiplicative version (factorizations) is A000688, coprime A354911.
Allowing 1's gives A023893, strict A106244, ranked by A302492.
The strict version is A054685.
The version for just primes is ranked by A076610, squarefree A356065.
Twice-partitions of this type are counted by A279784, factorizations A295935.
These partitions are ranked by A355743.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970.

Table[Length[Select[IntegerPartitions[n],And@@PrimePowerQ/@#&]],{n,0,30}] (* Gus Wiseman, Jul 28 2022 *)

is_primepower(n)= {ispower(n, , &n); isprime(n)}
lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (is_primepower(k), 1/(1-x^k), 1)); for (n=0, m, print1(polcoeff(gf, n, t), ", "));}
\\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023894(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())
    return (c(n)+sum(c(k)*A023894(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: Prod(p prime, Prod(k >= 1, 1/(1-x^(p^k))))

A355743 Numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 15, 17, 19, 21, 23, 25, 27, 31, 33, 35, 41, 45, 49, 51, 53, 55, 57, 59, 63, 67, 69, 75, 77, 81, 83, 85, 93, 95, 97, 99, 103, 105, 109, 115, 119, 121, 123, 125, 127, 131, 133, 135, 147, 153, 155, 157, 159, 161, 165, 171, 175, 177, 179, 187

Offset: 1

Gus Wiseman, Jul 24 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Also MM-numbers of multiset partitions into constant multisets, where the multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

			The terms together with their prime indices begin:
   1: {}
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  11: {5}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  31: {11}
  33: {2,5}
  35: {3,4}
  41: {13}
  45: {2,2,3}

Robert Price, Table of n, a(n) for n = 1..1410

The multiplicative version is A000688, strict A050361, coprime A354911.
The case of only primes (not all prime-powers) is A076610, strict A302590.
Allowing prime index 1 gives A302492.
These are the products of elements of A302493.
Requiring n to be a prime-power gives A302601.
These are the positions of 1's in A355741.
The squarefree case is A356065.
The complement is A356066.
A001222 counts prime-power divisors.
A023894 counts ptns into prime-powers, strict A054685, with 1's A023893.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.
Cf. A085970, A106244, A279784, A295935, A355731, A356064.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],And@@PrimePowerQ/@primeMS[#]&]

A294336 Number of ways to write n as a finite power-tower a^(b^(c^...)) of positive integers greater than one.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Oct 28 2017

nonn

Möbius-transform of A294337. - Antti Karttunen, Jun 12 2018

			The a(4096) = 7 ways are: 2^12, 4^6, 8^4, 8^(2^2), 16^3, 64^2, 4096.

Antti Karttunen, Table of n, a(n) for n = 1..65537

Cf. A000041, A001055, A001597, A007916, A052409, A052410, A089723, A164336, A277562, A284639, A288636, A294337, A294338, A294339, A375598, A375599.

Array[1+Sum[#0[g],{g,Rest[Divisors[GCD@@FactorInteger[#1][[All,2]]]]}]&,200]

A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
A294336(n) = if(1==n,n,sumdiv(A052409(n),d,A294336(d))); \\ Antti Karttunen, Jun 12 2018, after Mathematica-code.

a(1) = 1; for n > 1, a(n) = Sum_{d|A052409(n)} a(d). - Antti Karttunen, Jun 12 2018, after Mathematica-code.

a(n) = A294337(A052409(n)) for n >= 2. - Pontus von Brömssen, Aug 20 2024

More terms from Antti Karttunen, Jun 12 2018

A085970 Number of integers ranging from 2 to n that are not prime-powers.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 5, 5, 6, 6, 7, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 13, 14, 15, 16, 17, 17, 18, 19, 20, 20, 21, 21, 22, 23, 24, 24, 25, 25, 26, 27, 28, 28, 29, 30, 31, 32, 33, 33, 34, 34, 35, 36, 36, 37, 38, 38, 39, 40, 41, 41, 42, 42, 43

Offset: 1

Reinhard Zumkeller, Jul 06 2003

nonn

For n > 2, a(n) gives the number of duplicate eliminations performed by the Sieve of Eratosthenes when sieving the interval [2, n]. - Felix Fröhlich, Dec 10 2016
Number of terms of A024619 <= n. - Felix Fröhlich, Dec 10 2016
First differs from A082997 at n = 30. - Gus Wiseman, Jul 28 2022

			The a(30) = 13 numbers: 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30. - _Gus Wiseman_, Jul 28 2022

Felix Fröhlich, Table of n, a(n) for n = 1..10000

Cf. A000720, A024619, A085972.
The complement is counted by A065515, without 1's A025528.
For primes instead of prime-powers we have A065855, with 1's A062298.
Partial sums of A143731.
The version not treating 1 as a prime-power is A356068.
A000688 counts factorizations into prime-powers.
A001222 counts prime-power divisors.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A023893, A023894, A036234, A069513, A207375, A354911, A355743, A356064, A356065.

With[{nn = 75}, Table[n - Count[#, k_ /; k < n] - 1, {n, nn}] &@ Join[{1}, Select[Range@ nn, PrimePowerQ]]] (* Michael De Vlieger, Dec 11 2016 *)

a(n) = my(i=0); forcomposite(c=4, n, if(!isprimepower(c), i++)); i \\ Felix Fröhlich, Dec 10 2016

from sympy import primepi, integer_nthroot
def A085970(n): return n-1-sum(primepi(integer_nthroot(n,k)[0]) for k in range(1,n.bit_length())) # Chai Wah Wu, Aug 20 2024

a(n) = Max{A024619(k)<=n} k;

a(n) = n - A065515(n) = A085972(n) - A000720(n).

Name modified by Gus Wiseman, Jul 28 2022. Normally 1 is not considered a prime-power, cf. A000961, A246655.

A277562 Numbers of the form c(x_1)^c(x_2)^...^c(x_k) where each c(i) = A007916(i) is a non-perfect-power, k >= 2, and the exponents are nested from the right.

Original entry on oeis.org

16, 81, 256, 512, 625, 1296, 2401, 6561, 10000, 14641, 19683, 20736, 28561, 38416, 50625, 65536, 83521, 104976, 130321, 160000, 194481, 234256, 279841, 331776, 390625, 456976, 614656, 707281, 810000, 923521, 1185921, 1336336, 1500625, 1679616, 1874161, 1953125, 2085136, 2313441, 2560000, 2825761, 3111696, 3418801

Offset: 1

Gus Wiseman, Oct 19 2016

nonn

Non-perfect-powers, or NPPs (A007916), are numbers whose prime multiplicities are relatively prime. As discussed in A007916, the expansion of a positive integer into a tower of NPPs is unique and always possible. 65536=2^2^2^2 is the smallest number that requires a tower of height more than 3.

			       16 = 2^2^2        81 = 3^2^2       256 = 2^2^3       512 = 2^3^2
      625 = 5^2^2      1296 = 6^2^2      2401 = 7^2^2      6561 = 3^2^3
    10000 = 10^2^2    14641 = 11^2^2    19683 = 3^3^2     20736 = 12^2^2
    28561 = 13^2^2    38416 = 14^2^2    50625 = 15^2^2
    65536 = 2^2^2^2   83521 = 17^2^2   104976 = 18^2^2   130321 = 19^2^2
   160000 = 20^2^2   194481 = 21^2^2   234256 = 22^2^2   279841 = 23^2^2
   331776 = 24^2^2   390625 = 5^2^3    456976 = 26^2^2   614656 = 28^2^2
   707281 = 29^2^2   810000 = 30^2^2   923521 = 31^2^2  1185921 = 33^2^2
  1336336 = 34^2^2  1500625 = 35^2^2  1679616 = 6^2^3   1874161 = 37^2^2
  1953125 = 5^3^2   2085136 = 38^2^2  2313441 = 39^2^2  2560000 = 40^2^2
  2825761 = 41^2^2  3111696 = 42^2^2  3418801 = 43^2^2  3748096 = 44^2^2
  4100625 = 45^2^2  4477456 = 46^2^2  4879681 = 47^2^2  5308416 = 48^2^2
  5764801 = 7^2^3   6250000 = 50^2^2  6765201 = 51^2^2  7311616 = 52^2^2
  7890481 = 53^2^2  8503056 = 54^2^2  9150625 = 55^2^2  9834496 = 56^2^2

David A. Corneth, Table of n, a(n) for n = 1..10025 (terms <= 10^16)
R. K. Guy and J. L. Selfridge, The nesting and roosting habits of the laddered parenthesis, Amer. Math. Monthly 80 (8) (1973), 868-876.
R. K. Guy and J. L. Selfridge, The nesting and roosting habits of the laddered parenthesis (annotated cached copy)

Cf. A007916, A001597, A164336, A164337, A106490 (Quetian Superfactorization).

radicalQ[1]:=False;
radicalQ[n_]:=SameQ[GCD@@FactorInteger[n][[All,2]],1];
hyperfactor[1]:={};
hyperfactor[n_?radicalQ]:={n};
hyperfactor[n_]:=With[{g=GCD@@FactorInteger[n][[All,2]]},Prepend[hyperfactor[g],Product[Apply[Power[#1,#2/g]&,r],{r,FactorInteger[n]}]]];
Select[Range[10^6],Length[hyperfactor[#]]>2&]

Edited by N. J. A. Sloane, Nov 09 2016

Offset changed to 1 by David A. Corneth, Apr 30 2024

A277564 Let {c(i)} = A007916 denote the sequence of numbers > 1 which are not perfect powers. Every positive integer n has a unique representation as a tower n = c(x_1)^c(x_2)^c(x_3)^...^c(x_k), where the exponents are nested from the right. The sequence is an irregular triangle read by rows, where the n-th row lists n followed by x_1, ..., x_k.

Original entry on oeis.org

1, 2, 1, 3, 2, 4, 1, 1, 5, 3, 6, 4, 7, 5, 8, 1, 2, 9, 2, 1, 10, 6, 11, 7, 12, 8, 13, 9, 14, 10, 15, 11, 16, 1, 1, 1, 17, 12, 18, 13, 19, 14, 20, 15, 21, 16, 22, 17, 23, 18, 24, 19, 25, 3, 1, 26, 20, 27, 2, 2, 28, 21, 29, 22, 30, 23, 31, 24, 32, 1, 3, 33, 25, 34, 26, 35, 27, 36, 4, 1, 37, 28, 38, 29, 39, 30, 40, 31

Offset: 1

Gus Wiseman, Oct 20 2016

The row lengths are A288636(n) + 1. - Gus Wiseman, Jun 12 2017

See A278028 for a version in which row n simply lists x_1, x_2, ..., x_k (omitting the initial n).

			1 is represented by the empty sequence (), by convention.
Successive rows of the triangle are as follows (c(k) denotes the k-th non-prime-power, A007916(k)):
2, 1,
3, 2,
4, 1, 1,
5, 3,
6, 4, because 6 = c(4)
7, 5,
8, 1, 2, because 8 = 2^3 = c(1)^c(2)
9, 2, 1,
10, 6,
11, 7,
...
16, 1, 1, 1, because 16 = 2^4 = c(1)^4 = c(1)^(c(1)^2) = c[1]^(c[1]^c[1])
17, 12,
...
This sequence represents a bijection N -> Q where Q is the set of all finite sequences of positive integers: 1->(), 2->(1), 3->(2), 4->(1 1), 5->(3), 6->(4), 7->(5), 8->(1 2), 9->(2 1), ...

Gus Wiseman, Table of n, a(n) for n = 1..20131
N. J. A. Sloane, Maple programs for A007916, A278028, A278029, A052409, A089723, A277564

Cf. A007916, A277562, A164336, A000961, A164337, A278028, A089723.

Maple
```
See link.
```

nn=10000;radicalQ[1]:=False;radicalQ[n_]:=SameQ[GCD@@FactorInteger[n][[All,2]],1];
hyperfactor[1]:={};hyperfactor[n_?radicalQ]:={n};hyperfactor[n_]:=With[{g=GCD@@FactorInteger[n][[All,2]]},Prepend[hyperfactor[g],Product[Apply[Power[#1,#2/g]&,r],{r,FactorInteger[n]}]]];
rad[0]:=1;rad[n_?Positive]:=rad[n]=NestWhile[#+1&,rad[n-1]+1,Not[radicalQ[#]]&];Set@@@Array[radPi[rad[#]]==#&,nn];
Flatten[Join[{#},radPi/@hyperfactor[#]]&/@Range[nn]]

Edited by N. J. A. Sloane, Nov 09 2016

cp's OEIS Frontend

A164337 Those powers of primes that are not in sequence A164336.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Extensions

A007916 Numbers that are not perfect powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A050361 Number of factorizations into distinct prime powers greater than 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A023893 Number of partitions of n into prime power parts (1 included); number of nonisomorphic Abelian subgroups of symmetric group S_n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A023894 Number of partitions of n into prime power parts (1 excluded).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A355743 Numbers whose prime indices are all prime-powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A294336 Number of ways to write n as a finite power-tower a^(b^(c^...)) of positive integers greater than one.

Views

Author

Keywords

Comments