A303386 -id:A303386 - cp's OEIS frontend

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

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

A316439 Irregular triangle where T(n,k) is the number of factorizations of n into k factors > 1, with k ranging from 1 to Omega(n).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 03 2018

			The factorizations of 24 are (2*2*2*3), (2*2*6), (2*3*4), (2*12), (3*8), (4*6), (24) so the 24th row is {1, 3, 2, 1}.
Triangle begins:
  {}
  1
  1
  1  1
  1
  1  1
  1
  1  1  1
  1  1
  1  1
  1
  1  2  1
  1
  1  1
  1  1
  1  2  1  1
  1
  1  2  1
  1
  1  2  1
  1  1
  1  1
  1
  1  3  2  1
  1  1
  1  1
  1  1  1
  1  2  1
  1
  1  3  1

Alois P. Heinz, Rows n = 1..20000, flattened

Cf. A001222 (row lengths), A001055 (row sums), A001970, A007716, A045778, A162247, A259936, A281116, A303386.

g:= proc(n, k) option remember; `if`(n>k, 0, x)+
      `if`(isprime(n), 0, expand(x*add(`if`(d>k, 0,
      g(n/d, d)), d=numtheory[divisors](n) minus {1, n})))
    end:
T:= n-> `if`(n=1, [][], (p-> seq(coeff(p, x, i)
        , i=1..degree(p)))(g(n$2))):
seq(T(n), n=1..50);  # Alois P. Heinz, Aug 11 2019

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[#]==k&]],{n,100},{k,PrimeOmega[n]}]

A301700 Number of aperiodic rooted trees with n nodes.

Original entry on oeis.org

1, 1, 1, 2, 4, 10, 21, 52, 120, 290, 697, 1713, 4200, 10446, 26053, 65473, 165257, 419357, 1068239, 2732509, 7013242, 18059960, 46641983, 120790324, 313593621, 816046050, 2128101601, 5560829666, 14557746453, 38177226541, 100281484375, 263815322761, 695027102020

Offset: 1

Gus Wiseman, Apr 23 2018

nonn

An unlabeled rooted tree is aperiodic if the multiset of branches of the root is an aperiodic multiset, meaning it has relatively prime multiplicities, and each branch is also aperiodic.

			The a(6) = 10 aperiodic trees are (((((o))))), (((o(o)))), ((o((o)))), ((oo(o))), (o(((o)))), (o(o(o))), ((o)((o))), (oo((o))), (o(o)(o)), (ooo(o)).

Andrew Howroyd, Table of n, a(n) for n = 1..500

Cf. A000081, A000740, A000837, A001678, A003238, A004111, A007716, A007916, A100953, A276625, A284639, A290689, A298422, A303386, A303431.

arut[n_]:=arut[n]=If[n===1,{{}},Join@@Function[c,Select[Union[Sort/@Tuples[arut/@c]],GCD@@Length/@Split[#]===1&]]/@IntegerPartitions[n-1]];
Table[Length[arut[n]],{n,20}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
MoebiusT(v)={vector(#v, n, sumdiv(n,d,moebius(n/d)*v[d]))}
seq(n)={my(v=[1]); for(n=2, n, v=concat([1], MoebiusT(EulerT(v)))); v} \\ Andrew Howroyd, Sep 01 2018

Terms a(21) and beyond from Andrew Howroyd, Sep 01 2018

A303431 Aperiodic tree numbers. Matula-Goebel numbers of aperiodic rooted trees.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 12, 13, 15, 18, 20, 22, 24, 26, 29, 30, 31, 33, 37, 39, 40, 41, 44, 45, 47, 48, 50, 52, 54, 55, 58, 60, 61, 62, 65, 66, 71, 72, 74, 75, 78, 79, 80, 82, 87, 88, 89, 90, 93, 94, 96, 99, 101, 104, 108, 109, 110, 111, 113, 116, 117, 120, 122

Offset: 1

Gus Wiseman, Apr 23 2018

nonn

A positive integer is an aperiodic tree number iff either it is equal to 1 or it belongs to A007916 (numbers that are not perfect powers, or numbers whose prime multiplicities are relatively prime) and all of its prime indices are also aperiodic tree numbers, where a prime index of n is a number m such that prime(m) divides n.

			Sequence of aperiodic rooted trees begins:
01 o
02 (o)
03 ((o))
05 (((o)))
06 (o(o))
10 (o((o)))
11 ((((o))))
12 (oo(o))
13 ((o(o)))
15 ((o)((o)))
18 (o(o)(o))
20 (oo((o)))
22 (o(((o))))
24 (ooo(o))
26 (o(o(o)))
29 ((o((o))))
30 (o(o)((o)))
31 (((((o)))))
33 ((o)(((o))))

Cf. A000081, A000740, A000837, A007097, A007916, A052409, A052410, A061775, A111299, A214577, A275024, A276625, A290760, A290822, A291442, A298533, A298536, A301700, A302242, A303386.

zapQ[1]:=True;zapQ[n_]:=And[GCD@@FactorInteger[n][[All,2]]===1,And@@zapQ/@PrimePi/@FactorInteger[n][[All,1]]];
Select[Range[100],zapQ]

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 18 2018

nonn

In a factorization, two primes are equivalent if each factor has in its prime factorization the same multiplicity of both primes.

			The a(36) = 7 factorizations are (2*2*3*3), (2*2*9), (2*3*6), (3*3*4), (2*18), (3*12), (4*9). Missing from this list are (6*6) and (36).

Cf. A001055, A007716, A007717, A020555, A045778, A162247, A281116, A303386.

Cf. A316974, A316979, A316980, A316981.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
dual[eds_]:=Table[First/@Position[eds,x],{x,Union@@eds}];
Table[Length[Select[facs[n],UnsameQ@@dual[primeMS/@#]&]],{n,100}]

a(prime^n) = A000041(n).

a(squarefree) = 1.

A303546 Number of non-isomorphic aperiodic multiset partitions of weight n.

Original entry on oeis.org

1, 3, 9, 29, 90, 285, 909, 2984, 9935, 34113, 119368, 428923, 1574223, 5915235, 22699730, 89000042, 356058539, 1453069854, 6044132793, 25612564200, 110503626702, 485161228675, 2166488899641, 9835209480533, 45370059225227

Offset: 1

Gus Wiseman, Apr 26 2018

nonn

A multiset is aperiodic if its multiplicities are relatively prime. For this sequence neither the parts nor their multiset union are required to be aperiodic, only the multiset of parts.

			Non-isomorphic representatives of the a(3) = 9 aperiodic multiset partitions are:
  {{1,1,1}}, {{1,2,2}}, {{1,2,3}},
  {{1},{1,1}}, {{1},{2,2}}, {{1},{2,3}}, {{2},{1,2}},
  {{1},{2},{2}}, {{1},{2},{3}}.

Andrew Howroyd, Table of n, a(n) for n = 1..50

Cf. A000740, A000837, A001055, A007716, A007916, A049311, A100953, A255906, A269134, A275024, A293993, A301700, A303386, A303431, A303547.

a(n) = Sum_{d|n} mu(d) * A007716(n/d).

A114592 Sum_{n>=1} a(n)/n^s = Product_{k>=2} (1 - 1/k^s).

Original entry on oeis.org

1, -1, -1, -1, -1, 0, -1, 0, -1, 0, -1, 1, -1, 0, 0, 0, -1, 1, -1, 1, 0, 0, -1, 1, -1, 0, 0, 1, -1, 1, -1, 1, 0, 0, 0, 1, -1, 0, 0, 1, -1, 1, -1, 1, 1, 0, -1, 1, -1, 1, 0, 1, -1, 1, 0, 1, 0, 0, -1, 1, -1, 0, 1, 0, 0, 1, -1, 1, 0, 1, -1, 1, -1, 0, 1, 1, 0, 1

Offset: 1

Leroy Quet, Dec 11 2005

sign

For n >= 2, Sum_{k|n} A001055(n/k) * a(k) = 0. A114591(n) = Sum_{k|n} a(k).

First entry greater than 1 in absolute value is a(360) = -2. - Gus Wiseman, Sep 15 2018

			24 can be factored into distinct integers (each >= 2) as 24; as 4*6, 3*8 and 2*12; and as 2*3*4. (A045778(24) = 5).
So a(24) = (-1)^1 + 3*(-1)^2 + (-1)^3 = 1, where the 1 exponent is due to the 1 factor of the 24 = 24 factorization and the 2 exponent is due to the 3 cases of 2 factors each of the 24 = 4*6 = 3*8 = 2*12 factorizations and the 3 exponent is due to the 24 = 2*3*4 factorization.

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

Cf. A001055, A045778, A114591.

Cf. A001222, A162247, A190938, A259936, A281116, A303386, A316441, A319237, A319238.

strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
Table[Sum[(-1)^Length[f],{f,strfacs[n]}],{n,100}] (* Gus Wiseman, Sep 15 2018 *)

A114592aux(n, k) = if(1==n, 1, sumdiv(n, d, if(d > 1 && d <= k && d < n, (-1)*A114592aux(n/d, d-1))) - (n<=k)); \\ After code in A045778.
A114592(n) = A114592aux(n,n); \\ Antti Karttunen, Jul 23 2017

a(1) = 1; for n>= 2, a(n) = sum, over ways to factor n into any number of distinct integers >= 2, of (-1)^(number of integers in a factorization). (See example.)

More terms from Antti Karttunen, Jul 23 2017

A316441 a(n) = Sum (-1)^k where the sum is over all factorizations of n into factors > 1 and k is the number of factors.

Original entry on oeis.org

1, -1, -1, 0, -1, 0, -1, -1, 0, 0, -1, 0, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, -1, 1, 0, 0, -1, 0, -1, 1, -1, -1, 0, 0, 0, 1, -1, 0, 0, 1, -1, 1, -1, 0, 0, 0, -1, 0, 0, 0, 0, 0, -1, 1, 0, 1, 0, 0, -1, 1, -1, 0, 0, 1, 0, 1, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, 0, 1, -1

Offset: 1

Gus Wiseman, Jul 03 2018

sign

First term greater than 1 in absolute value is a(256) = 2.

			The factorizations of 24 are (2*2*2*3), (2*2*6), (2*3*4), (2*12), (3*8), (4*6), (24); so a(24) = 1 - 2 + 3 - 1 = 1.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Cf. A001222, A001055, A045778, A114592, A162247, A190938, A259936, A281116, A303386.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Sum[(-1)^Length[f],{f,facs[n]}],{n,200}]

A316441(n, m=n, k=0) = if(1==n, (-1)^k, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A316441(n/d, d, k+1))); (s)); \\ Antti Karttunen, Sep 08 2018, after Michael B. Porter's code for A001055

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

Secondary offset added by Antti Karttunen, Sep 08 2018

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

A305149 Number of factorizations of n whose distinct factors are pairwise indivisible and greater than 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 3, 1, 3, 1, 3, 2, 2, 1, 4, 2, 2, 2, 3, 1, 5, 1, 2, 2, 2, 2, 6, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 5, 2, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 8, 1, 2, 3, 4, 2, 5, 1, 3, 2, 5, 1, 6, 1, 2, 3, 3, 2, 5, 1, 5, 3, 2, 1, 8, 2, 2, 2, 4, 1, 8, 2, 3, 2, 2, 2, 6, 1, 3, 3, 6, 1, 5, 1, 4, 5

Offset: 1

Gus Wiseman, May 26 2018

nonn

			The a(60) = 8 factorizations are (2*2*3*5), (2*2*15), (3*4*5), (3*20), (4*15), (5*12), (6*10), (60). Missing from this list are (2*3*10), (2*5*6), (2*30).

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

Cf. A001055, A001970, A007716, A034444, A045778, A259936, A281116, A285572, A302242, A303386, A303431, A305001, A305148, A305150.

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],Select[Tuples[Union[#],2],UnsameQ@@#&&Divisible@@#&]=={}&]],{n,100}]

pairwise_indivisible(v) = { for(i=1,#v,for(j=i+1,#v,if(!(v[j]%v[i]),return(0)))); (1); };
A305149(n, m=n, facs=List([])) = if(1==n, pairwise_indivisible(Set(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A305149(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Oct 08 2018

More terms from Antti Karttunen, Oct 08 2018

cp's OEIS Frontend

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

A316439 Irregular triangle where T(n,k) is the number of factorizations of n into k factors > 1, with k ranging from 1 to Omega(n).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A301700 Number of aperiodic rooted trees with n nodes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A303431 Aperiodic tree numbers. Matula-Goebel numbers of aperiodic rooted trees.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A303546 Number of non-isomorphic aperiodic multiset partitions of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A114592 Sum_{n>=1} a(n)/n^s = Product_{k>=2} (1 - 1/k^s).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula