A089723 -id:A089723 - cp's OEIS frontend

A346402 Decimal expansion of Sum_{k>=1} A089723(k)/2^k.

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

Offset: 1

Amiram Eldar, Jul 15 2021

This constant is irrational (Golomb, 1973).

			1.06838993007841232002084765901758238224728960865528...

Solomon W. Golomb, A new arithmetic function of combinatorial significance, Journal of Number Theory, Vol. 5, No. 3 (1973), pp. 218-223.

Cf. A089723.

f[n_] := DivisorSigma[0, GCD @@ FactorInteger[n][[;; , 2]]]; RealDigits[Sum[f[n]/2^n, {n, 1, 1000}], 10, 100][[1]]

1/2 + suminf(k=2, numdiv(gcd(factor(k)[, 2]))/2^k) \\ Michel Marcus, Jul 16 2021

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

A052409 a(n) = largest integer power m for which a representation of the form n = k^m exists (for some k).

Original entry on oeis.org

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

Eric W. Weisstein

nonn

Greatest common divisor of all prime-exponents in canonical factorization of n for n>1: a(n)>1 iff n is a perfect power; a(A001597(k))=A025479(k). - Reinhard Zumkeller, Oct 13 2002
a(1) set to 0 since there is no largest finite integer power m for which a representation of the form 1 = 1^m exists (infinite largest m). - Daniel Forgues, Mar 06 2009
A052410(n)^a(n) = n. - Reinhard Zumkeller, Apr 06 2014
Positions of 1's are A007916. Smallest base is given by A052410. - Gus Wiseman, Jun 09 2020

			n = 72 = 2*2*2*3*3: GCD[exponents] = GCD[3,2] = 1. This is the least n for which a(n) <> A051904(n), the minimum of exponents.
For n = 10800 = 2^4 * 3^3 * 5^2, GCD[4,3,2] = 1, thus a(10800) = 1.

Daniel Forgues, Table of n, a(n) for n = 1..100000
N. J. A. Sloane, Maple programs for A007916, A278028, A278029, A052409, A089723, A277564
Eric Weisstein's World of Mathematics, Power
Eric Weisstein's World of Mathematics, Perfect Power
Index entries for sequences computed from exponents in factorization of n

Cf. A052410, A005361, A051903, A072411-A072414, A124010, A075802, A072776, A270492.
Apart from the initial term essentially the same as A253641.
Differs from A051904 for the first time at n=72, where a(72) = 1, while A051904(72) = 2.
Differs from A158378 for the first time at n=10800, where a(10800) = 1, while A158378(10800) = 2.
Cf. A000005, A000961, A001597, A052410, A303386, A327501.

a052409 1 = 0
a052409 n = foldr1 gcd $ a124010_row n
-- Reinhard Zumkeller, May 26 2012

# See link.
#
a:= n-> igcd(map(i-> i[2], ifactors(n)[2])[]):
seq(a(n), n=1..120);  # Alois P. Heinz, Oct 20 2019

Table[GCD @@ Last /@ FactorInteger[n], {n, 100}] (* Ray Chandler, Jan 24 2006 *)

a(n)=my(k=ispower(n)); if(k, k, n>1) \\ Charles R Greathouse IV, Oct 30 2012

from math import gcd
from sympy import factorint
def A052409(n): return gcd(*factorint(n).values()) # Chai Wah Wu, Aug 31 2022

(define (A052409 n) (if (= 1 n) 0 (gcd (A067029 n) (A052409 (A028234 n))))) ;; Antti Karttunen, Aug 07 2017

a(1) = 0; for n > 1, a(n) = gcd(A067029(n), a(A028234(n))). - Antti Karttunen, Aug 07 2017

More terms from Labos Elemer, Jun 17 2002

A295935 Number of twice-factorizations of n where the latter factorizations are constant, i.e., type (P,P,R).

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 5, 3, 2, 1, 5, 1, 2, 2, 12, 1, 5, 1, 5, 2, 2, 1, 10, 3, 2, 5, 5, 1, 5, 1, 18, 2, 2, 2, 15, 1, 2, 2, 10, 1, 5, 1, 5, 5, 2, 1, 22, 3, 5, 2, 5, 1, 10, 2, 10, 2, 2, 1, 13, 1, 2, 5, 40, 2, 5, 1, 5, 2, 5, 1, 28, 1, 2, 5, 5, 2, 5, 1, 22, 12, 2, 1

Offset: 1

Gus Wiseman, Nov 29 2017

nonn

a(n) is also the number of ways to choose a perfect divisor of each factor in a factorization of n.

			The a(24) = 10 twice-factorizations are:
(2)*(2)*(2)*(3), (2)*(3)*(2*2), (3)*(2*2*2)
(2)*(2)*(6), (2*2)*(6),
(2)*(3)*(4),
(2)*(12),
(3)*(8),
(4)*(6),
(24).

Cf. A000005, A001055, A052409, A052410, A089723, A279784, A281113, A284639, A295923, A295924, A295931.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Sum[Product[Length[Divisors[GCD@@FactorInteger[d][[All,2]]]],{d,f}],{f,facs[n]}],{n,100}]

Dirichlet g.f.: 1/Product_{n > 1}(1 - A089723(n)/n^s).

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5

Offset: 1

Geoffrey Critzer, Jul 09 2015

nonn

Equivalently, a(n) is the number of ways to express the cyclic group Z_n as a direct sum of its Hall subgroups.  A Hall subgroup of a finite group G is a subgroup whose order is coprime to its index.
a(n) is the number of ways to partition the set of distinct prime factors of n.
Also the number of singleton or pairwise coprime factorizations of n. - Gus Wiseman, Sep 24 2019

			a(60) = 5 because we have: 60 = 4*3*5 = 4*15 = 3*20 = 5*12.
For n = 36, its unitary divisors are 1, 4, 9, 36. From these we obtain 36 either as 1*36 or 4*9, thus a(36) = 2. - _Antti Karttunen_, Oct 21 2017

Alois P. Heinz, Table of n, a(n) for n = 1..20000
Wikipedia, Hall subgroup
Index entries for sequences computed from exponents in factorization of n

Cf. A000110, A001055, A001221, A034444, A089233, A258466, A281116, A285572.
Differs from A050320 for the first time at n=36.
Differs from A354870 for the first time at n=210, where a(210) = 15, while A354870(210) = 12.
Cf. A304716, A302569, A304711, A305079.
Related classes of factorizations:
- No conditions: A001055
- Strict: A045778
- Constant: A089723
- Distinct multiplicities: A255231
- Singleton or coprime: A259936
- Relatively prime: A281116
- Aperiodic: A303386
- Stable (indivisible): A305149
- Connected: A305193
- Strict relatively prime: A318721
- Uniform: A319269
- Intersecting: A319786
- Constant or distinct factors coprime: A327399
- Constant or relatively prime: A327400
- Coprime: A327517
- Not relatively prime: A327658
- Distinct factors coprime: A327695

map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015

Table[BellB[PrimeNu[n]], {n, 1, 75}]
(* second program *)
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],Length[#]==1||CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Sep 24 2019 *)

a(n) = my(t=omega(n), x='x, m=contfracpnqn(matrix(2, t\2, y, z, if( y==1, -z*x^2, 1 - (z+1)*x)))); polcoeff(1/(1 - x + m[2, 1]/m[1, 1]) + O(x^(t+1)), t) \\ Charles R Greathouse IV, Jun 30 2017

a(n) = A000110(A001221(n)).

a(n > 1) = A327517(n) + 1. - Gus Wiseman, Sep 24 2019

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A169594 Number of divisors of n, counting divisor multiplicity in n.

Original entry on oeis.org

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

Offset: 1

Joseph L. Pe, Dec 02 2009

The multiplicity of a divisor d > 1 in n is defined as the largest power i for which d^i divides n; and for d = 1 it is defined as 1.
a(n) is also the sum of the multiplicities of the divisors of n.
In other words, a(n) = 1 + sum of the highest exponents e_i for which each number k_i in range 2 .. n divide n, as {k_i}^{e_i} | n. For nondivisors of n this exponent e_i is 0, for n itself it is 1. - Antti Karttunen, May 20 2017
From Gus Wiseman, Mar 25 2021: (Start)
Also the number of strict chains of divisors ending with n and having constant (equal) first quotients. The case starting with 1 is A089723. For example, the a(1) = 1 through a(12) = 7 chains are:
  1  2    3    4      5    6    7    8        9      10    11    12
     1|2  1|3  1|4    1|5  1|6  1|7  1|8      1|9    1|10  1|11  1|12
               2|4         2|6       2|8      3|9    2|10        2|12
               1|2|4       3|6       4|8      1|3|9  5|10        3|12
                                     2|4|8                       4|12
                                     1|2|4|8                     6|12
                                                                 3|6|12
(End)
a(n) depends only on the prime signature of n. - David A. Corneth, Mar 28 2021

			The divisors of 8 are 1, 2, 4, 8 of multiplicity 1, 3, 1, 1, respectively. So a(8) = 1 + 3 + 1 + 1 = 6.

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

Cf. A168512.
Row sums of A286561, A286563 and A286564.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A057567 counts chains of divisors with weakly increasing first quotients.
A067824 counts strict chains of divisors ending with n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.
A342086 counts chains of divisors with strictly increasing first quotients.
A342496 counts partitions with equal first quotients (strict: A342515, ranking: A342522, ordered: A342495).
A342530 counts chains of divisors with distinct first quotients.
Cf. A003238, A003242, A069916, A122651, A309891, A318991, A318992, A325545.
First differences of A078651.

a := n -> ifelse(n < 2, 1, 1 + add(padic:-ordp(n, k), k = 2..n)):
seq(a(n), n = 1..98);  # Peter Luschny, Apr 10 2025

divmult[d_, n_] := Module[{output, i}, If[d == 1, output = 1, If[d == n, output = 1, i = 0; While[Mod[n, d^(i + 1)] == 0, i = i + 1]; output = i]]; output]; dmt0[n_] := Module[{divs, l}, divs = Divisors[n]; l = Length[divs]; Sum[divmult[divs[[i]], n], {i, 1, l}]]; Table[dmt0[i], {i, 1, 40}]
Table[1 + DivisorSum[n, IntegerExponent[n, #] &, # > 1 &], {n, 98}] (* Michael De Vlieger, May 20 2017 *)

A286561(n,k) = { my(i=1); if(1==k, 1, while(!(n%(k^i)), i = i+1); (i-1)); };
A169594(n) = sumdiv(n,d,A286561(n,d)); \\ Antti Karttunen, May 20 2017

a(n) = { if(n == 1, return(1)); my(f = factor(n), u = vecmax(f[, 2]), cf = f, res = numdiv(f) - u + 1); for(i = 2, u, cf[, 2] = f[, 2]\i; res+=numdiv(factorback(cf)) ); res } \\ David A. Corneth, Mar 29 2021

A169594(n) = {my(s=0, k=2); while(k<=n, s+=valuation(n, k); k=k+1); s + 1} \\ Zhuorui He, Aug 28 2025

def a286561(n, k):
    i=1
    if k==1: return 1
    while n%(k**i)==0:
        i+=1
    return i-1
def a(n): return sum([a286561(n, d) for d in divisors(n)]) # Indranil Ghosh, May 20 2017

(define (A169594 n) (add (lambda (k) (A286561bi n k)) 1 n))
;; Implements sum_{i=lowlim..uplim} intfun(i)
(define (add intfun lowlim uplim) (let sumloop ((i lowlim) (res 0)) (cond ((> i uplim) res) (else (sumloop (1+ i) (+ res (intfun i)))))))
;; For A286561bi see A286561. - Antti Karttunen, May 20 2017

From Friedjof Tellkamp, Feb 29 2024: (Start)
a(n) = A309891(n) + 1.
G.f.: x/(1-x) + Sum_{k>=2, j>=1} x^(k^j)/(1-x^(k^j)).
Dirichlet g.f.: zeta(s) * (1 + Sum_{k>=1} (zeta(k*s) - 1)).
Sum_{n>=1} a(n)/n^2 = (7/24) * Pi^2. (End)

Extended by Ray Chandler, Dec 08 2009

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]

A284639 Number of ways to write n > 1 as a power of a product n = (x_1..x_k)^y, where y >= 1, x_i > 1, and x_i <= x_j for i < j.

Original entry on oeis.org

1, 1, 3, 1, 2, 1, 4, 3, 2, 1, 4, 1, 2, 2, 8, 1, 4, 1, 4, 2, 2, 1, 7, 3, 2, 4, 4, 1, 5, 1, 8, 2, 2, 2, 11, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 3, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 17, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 8, 2, 1, 11, 2

Offset: 2

Gus Wiseman, Mar 30 2017

nonn

a(1)=infinity is not an integer so is necessarily omitted from the data.

			The a(64)=17 ways to write 64 as a power of a product:
(2*2*2*2*2*2)^1   (2*2*2)^2  (2*2)^3  (2)^6
(2*2*2*2*4)^1     (2*4)^2    (4)^3
(2*2*2*8)^1       (8)^2
(2*2*4*4)^1
(2*2*16)^1
(2*4*8)^1
(2*32)^1
(4*4*4)^1
(4*16)^1
(8*8)^1
(64)^1.

Vincenzo Librandi, Table of n, a(n) for n = 2..1000

Cf. A001055, A001597, A007916, A089723.

nn=85;
postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[DivisorSum[GCD@@FactorInteger[n][[All,2]],Length[postfacs[n^(1/#)]]&],{n,2,nn}]

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.

cp's OEIS Frontend

A346402 Decimal expansion of Sum_{k>=1} A089723(k)/2^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A007916 Numbers that are not perfect powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A052409 a(n) = largest integer power m for which a representation of the form n = k^m exists (for some k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Scheme

Formula

Extensions

A295935 Number of twice-factorizations of n where the latter factorizations are constant, i.e., type (P,P,R).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A169594 Number of divisors of n, counting divisor multiplicity in n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A284639 Number of ways to write n > 1 as a power of a product n = (x_1..x_k)^y, where y >= 1, x_i > 1, and x_i <= x_j for i < j.