A277562 -id:A277562 - cp's OEIS frontend

A007916 Numbers that are not perfect powers.

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

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

A288636 Height of power-tower factorization of n. Row lengths of A278028.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 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, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Jun 12 2017

nonn

After a(1)=0 this sequence has many terms equal to A089723. It first differs at a(64)=2, A089723(64)=4.

First positions of a(n) = {0, 1, 2, 3, 4} are n = {1, 2, 4, 16, 65536}. - Michael De Vlieger, Nov 24 2017

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

Cf. A007916, A089723, A164337, A277562, A277564, A278028.

a[n_] := If[n == 1, 0, 1+a[GCD @@ FactorInteger[n][[All, 2]]]];
Array[a,100]

A278028 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 x_1, ..., x_k.

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane, Nov 09 2016

Row lengths are A288636(n). - Gus Wiseman, Jun 12 2017

			Rows 2 through 32 are:
1,
2,
1, 1,
3,
4,
5,
1, 2,
2, 1,
6,
7,
8,
9,
10,
11,
1, 1, 1,
12,
13,
14,
15,
16,
17,
18,
19,
3, 1,
20,
2, 2,
21,
22,
23,
24,
1, 3,
...

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

See A277564 for another version.

Cf. A007916, A089723, A277562, A288636.

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

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

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 2, 6, 4, 4, 2, 7, 2, 4, 4, 10, 2, 7, 2, 7, 4, 4, 2, 10, 4, 4, 6, 7, 2, 8, 2, 12, 4, 4, 4, 12, 2, 4, 4, 10, 2, 8, 2, 7, 7, 4, 2, 15, 4, 7, 4, 7, 2, 10, 4, 10, 4, 4, 2, 13, 2, 4, 7, 16, 4, 8, 2, 7, 4, 8, 2, 16, 2, 4, 7, 7, 4, 8, 2, 15, 10, 4, 2, 13, 4, 4, 4, 10, 2, 13, 4, 7, 4, 4, 4, 18, 2, 7, 7, 12, 2, 8, 2, 10, 8

Offset: 1

Gus Wiseman, Oct 28 2017

nonn

			The a(12) = 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. A001597, A007916, A052409, A052410, A089723, A164336, A277562, A284639, A288636, A294336, A294338, A294339.

a[n_]:=1+Sum[a[g],{g,Rest[Divisors[GCD@@FactorInteger[n][[All,2]]]]}];
Table[a[2^n],{n,100}]

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)));
A294337(n) = sumdiv(n,d,A294336(d));
\\ Or alternatively, after Mathematica-code as:
A294337(n) = A294336(2^n); \\ Antti Karttunen, Jun 12 2018

a(n) = Sum_{d|n} A294336(d) = A294336(A000079(n)). - Antti Karttunen, Jun 12 2018

More terms from Antti Karttunen, Jun 12 2018

A295931 Number of ways to write n in the form n = (x^y)^z where x, y, and z are positive integers.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Nov 29 2017

nonn

By convention a(1) = 1.

Values can be 1, 3, 6, 9, 10, 15, 18, 21, 27, 28, 30, 36, 45, 54, 60, 63, 84, 90, etc. - Robert G. Wilson v, Dec 10 2017

			The a(256) = 10 ways are:
(2^1)^8    (2^2)^4   (2^4)^2  (2^8)^1
(4^1)^4    (4^2)^2   (4^4)^1
(16^1)^2   (16^2)^1
(256^1)^1

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A000005, A007425, A052409, A052410, A089723, A093771, A175082, A277562, A281113, A284639, A294786, A294336, A294338, A295920, A295923, A295924, A295935.

f:= proc(n) local m,d,t;
  m:= igcd(seq(t[2],t=ifactors(n)[2]));
  add(numtheory:-tau(d),d=numtheory:-divisors(m))
end proc:
f(1):= 1:
map(f, [$1..100]); # Robert Israel, Dec 19 2017

Table[Sum[DivisorSigma[0,d],{d,Divisors[GCD@@FactorInteger[n][[All,2]]]}],{n,100}]

a(A175082(k)) = 1, a(A093771(k)) = 3.

a(n) = Sum_{d|A052409(n)} A000005(d).

A294338 Number of ways to write n as a finite power-tower of positive integers greater than one, allowing both left and right nesting of parentheses.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 5, 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, 6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Oct 28 2017

nonn

			The a(16) = 5 ways are: 16, 4^2, (2^2)^2, 2^4, 2^(2^2).

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

Cf. A001597, A007916, A052409, A052410, A089723, A164336, A277562, A284639, A288636, A294336, A294337, A294339.

A294338 := proc(n)
    local expo,g,a,d ;
    if n =1 then
        return 1;
    end if;
    # compute gcd of the set of prime power exponents (A052409)
    ifactors(n)[2] ;
    [ seq(op(2,ep),ep=%)] ;
    igcd(op(%)) ;
    # set of divisors of A052409 (without the 1)
    g := numtheory[divisors](%) minus {1} ;
    a := 0 ;
    for d in g do
        # recursive (sort of convolution) call
        a := a+ procname(d)*procname(root[d](n)) ;
    end do:
    1+a ;
end proc:
seq(A294338(n),n=1..120) ; # R. J. Mathar, Nov 27 2017

a[n_]:=1+Sum[a[n^(1/g)]*a[g],{g,Rest[Divisors[GCD@@FactorInteger[n][[All,2]]]]}];
Array[a,100]

A323520 Numbers of the form p^(k^2) where p is prime and k >= 0.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jan 17 2019

nonn

Cf. A000290, A000961, A001248, A001597, A026478, A050376, A277562, A323519, A323526, A323528.

Select[Range[100],#==1||PrimePowerQ[#]&&IntegerQ[Sqrt[FactorInteger[#][[1,2]]]]&]

A294339 Number of ways to write 2^n as a finite power-tower of positive integers greater than one, allowing both left and right nesting of parentheses.

Original entry on oeis.org

1, 2, 2, 5, 2, 6, 2, 12, 5, 6, 2, 19, 2, 6, 6, 32, 2, 19, 2, 19, 6, 6, 2, 56, 5, 6, 12, 19, 2, 26, 2, 79, 6, 6, 6, 71, 2, 6, 6, 56, 2, 26, 2, 19, 19, 6, 2, 169, 5, 19, 6, 19, 2, 56, 6, 56, 6, 6, 2, 101, 2, 6, 19, 203, 6, 26, 2, 19, 6, 26, 2, 237, 2, 6, 19, 19

Offset: 1

Gus Wiseman, Oct 28 2017

nonn

			The a(6) = 6 ways are 64, 8^2, (2^3)^2, 4^3, (2^2)^3, 2^6.

Hans Havermann, Table of n, a(n) for n = 1..10000

Cf. A001055, A001597, A007916, A052409, A052410, A089723, A164336, A277562, A284639, A288636, A294336, A294337, A294338.

f:= proc(n) option remember; local F,t,s,g,a;
  F:= ifactors(n)[2];
  g:= igcd(op(map(t -> t[2],F)));
  t:= 1;
  for s in numtheory:-divisors(g) minus {1} do
    t:= t + procname(mul(a[1]^(a[2]/s),a=F))*procname(s)
  od;
  t
end proc:
seq(f(2^n),n=1..100); # Robert Israel, Dec 01 2017

a[n_]:=1+Sum[a[n^(1/g)]*a[g],{g,Rest[Divisors[GCD@@FactorInteger[n][[All,2]]]]}];
Table[a[2^n],{n,100}]

a(n) = A294338(2^n). - R. J. Mathar, Nov 27 2017

cp's OEIS Frontend

A007916 Numbers that are not perfect powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

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

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A288636 Height of power-tower factorization of n. Row lengths of A278028.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

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

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A295931 Number of ways to write n in the form n = (x^y)^z where x, y, and z are positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple