A052409 -id:A052409 - cp's OEIS frontend

A303710 Number of factorizations of numbers that are not perfect powers using only numbers that are not perfect powers.

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

Offset: 1

Gus Wiseman, Apr 29 2018

nonn

Note that a factorization of a number that is not a perfect power (A007916) is always itself aperiodic, meaning the multiplicities of its factors are relatively prime.

			The a(19) = 4 factorizations of 24 are (2*2*2*3), (2*2*6), (2*12), (24).
The a(23) = 5 factorizations of 30 are (2*3*5), (2*15), (3*10), (5*6), (30).

Cf. A001055, A001597, A007716, A007916, A052409, A052410, A281116, A303365, A303386, A303707, A303708, A303709.

radQ[n_] := And[n > 1, GCD@@FactorInteger[n][[All, 2]] === 1]; facsr[n_] := If[n <= 1, {{}}, Join@@Table[Map[Prepend[#, d] &, Select[facsr[n/d], Min@@# >= d &]], {d, Select[Divisors[n], radQ]}]]; Table[Length[facsr[n]], {n, Select[Range[100], radQ]}]

A304250 Perfect powers whose prime factors span an initial interval of prime numbers.

Original entry on oeis.org

4, 8, 16, 32, 36, 64, 128, 144, 216, 256, 324, 512, 576, 900, 1024, 1296, 1728, 2048, 2304, 2916, 3600, 4096, 5184, 5832, 7776, 8100, 8192, 9216, 11664, 13824, 14400, 16384, 20736, 22500, 26244, 27000, 32400, 32768, 36864, 44100, 46656, 57600, 65536, 72900

Offset: 1

Gus Wiseman, May 13 2018

nonn

The multiset of prime indices of a(n) is the a(n)-th row of A112798. This multiset is normal, meaning it spans an initial interval of positive integers, and periodic, meaning its multiplicities have a common divisor greater than 1.

			Sequence of all normal periodic multisets begins
4:    {1,1}
8:    {1,1,1}
16:   {1,1,1,1}
32:   {1,1,1,1,1}
36:   {1,1,2,2}
64:   {1,1,1,1,1,1}
128:  {1,1,1,1,1,1,1}
144:  {1,1,1,1,2,2}
216:  {1,1,1,2,2,2}
256:  {1,1,1,1,1,1,1,1}
324:  {1,1,2,2,2,2}
512:  {1,1,1,1,1,1,1,1,1}
576:  {1,1,1,1,1,1,2,2}
900:  {1,1,2,2,3,3}
1024: {1,1,1,1,1,1,1,1,1,1}

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

Cf. A000837, A000961 A001597, A007916, A018783, A052409, A052410, A055932, A056239, A112798, A178472, A303431, A303709, A304450, A368682.

Select[Range[1000],FactorInteger[#][[-1,1]]==Prime[Length[FactorInteger[#]]]&&GCD@@FactorInteger[#][[All,2]]>1&]

Intersection of A001597 and A055932.

A322901 Numbers whose prime indices are all powers of the same number.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 34, 36, 37, 38, 40, 41, 42, 43, 44, 46, 47, 48, 49, 50, 52, 53, 54, 56, 57, 58, 59, 61, 62, 63, 64, 67, 68, 71, 72, 73, 74, 76, 79, 80, 81, 82, 83

Offset: 1

Gus Wiseman, Dec 30 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

			The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k). The sequence of all integer partitions whose Heinz numbers belong to the sequence begins: (), (1), (2), (11), (3), (21), (4), (111), (22), (31), (5), (211), (6), (41), (1111), (7), (221), (8), (311), (42), (51), (9), (2111), (33), (61), (222), (411).

Cf. A001597, A018819, A052409, A052410, A056239, A072720, A072721, A302242, A302593, A322900, A322902, A322903.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
radbase[n_]:=n^(1/GCD@@FactorInteger[n][[All,2]]);
Select[Range[100],SameQ@@radbase/@DeleteCases[primeMS[#],1]&]

A333941 Triangle read by rows where T(n,k) is the number of compositions of n with rotational period k.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 0, 2, 2, 0, 0, 3, 2, 3, 0, 0, 2, 4, 6, 4, 0, 0, 4, 6, 9, 8, 5, 0, 0, 2, 6, 15, 20, 15, 6, 0, 0, 4, 8, 24, 32, 35, 18, 7, 0, 0, 3, 10, 27, 56, 70, 54, 28, 8, 0, 0, 4, 12, 42, 84, 125, 120, 84, 32, 9, 0, 0, 2, 10, 45, 120, 210, 252, 210, 120, 45, 10, 0

Offset: 0

Gus Wiseman, Apr 16 2020

A composition of n is a finite sequence of positive integers summing to n.

			Triangle begins:
   1
   0   1
   0   2   0
   0   2   2   0
   0   3   2   3   0
   0   2   4   6   4   0
   0   4   6   9   8   5   0
   0   2   6  15  20  15   6   0
   0   4   8  24  32  35  18   7   0
   0   3  10  27  56  70  54  28   8   0
   0   4  12  42  84 125 120  84  32   9   0
   0   2  10  45 120 210 252 210 120  45  10   0
   0   6  18  66 168 335 450 462 320 162  50  11   0
Row n = 6 counts the following compositions (empty columns indicated by dots):
  .  (6)       (15)    (114)  (1113)  (11112)  .
     (33)      (24)    (123)  (1122)  (11121)
     (222)     (42)    (132)  (1131)  (11211)
     (111111)  (51)    (141)  (1221)  (12111)
               (1212)  (213)  (1311)  (21111)
               (2121)  (231)  (2112)
                       (312)  (2211)
                       (321)  (3111)
                       (411)

Andrew Howroyd, Table of n, a(n) for n = 0..1325 (rows 0..50)

Column k = 1 is A000005.
Row sums are A011782.
Diagonal T(2n,n) is A045630(n).
The strict version is A072574.
A version counting runs is A238279.
Column k = n - 1 is A254667.
Aperiodic compositions are counted by A000740.
Aperiodic binary words are counted by A027375.
The orderless period of prime indices is A052409.
Numbers whose binary expansion is periodic are A121016.
Periodic compositions are counted by A178472.
Period of binary expansion is A302291.
Numbers whose prime signature is aperiodic are A329139.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Rotational symmetries are counted by A138904.
- Constant compositions are A272919.
- Lyndon compositions are A275692.
- Co-Lyndon compositions are A326774.
- Aperiodic compositions are A328594.
- Rotational period is A333632.
- Co-necklaces are A333764.
- Reversed necklaces are A333943.
Cf. A000031, A001037, A008965, A019536, A211100, A291166, A328595, A328596, A329312, A329313, A329326.

Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Function[c,Length[Union[Array[RotateRight[c,#]&,Length[c]]]]==k]]],{n,0,10},{k,0,n}]

T(n,k)=if(n==0, k==0, sumdiv(n, m, sumdiv(gcd(k,m), d, moebius(d)*binomial(m/d-1, k/d-1)))) \\ Andrew Howroyd, Jan 19 2023

T(n,k) = Sum_{m|n} Sum_{d|gcd(k,m)} mu(d)*binomial(m/d-1, k/d-1) for n > 0. - Andrew Howroyd, Jan 19 2023

A367685 Numbers divisible by their multiset multiplicity kernel.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 31, 32, 36, 37, 40, 41, 43, 44, 45, 47, 48, 49, 50, 52, 53, 54, 56, 59, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 88, 89, 92, 96, 97, 98, 99, 100, 101, 103, 104, 107

Offset: 1

Gus Wiseman, Nov 30 2023

nonn

First differs from A344586 in lacking 120.
We define the multiset multiplicity kernel (MMK) of a positive integer n to be the product of (least prime factor with exponent k)^(number of prime factors with exponent k) over all distinct exponents k appearing in the prime factorization of n. For example, 90 has prime factorization 2^1 * 3^2 * 5^1, so for k = 1 we have 2^2, and for k = 2 we have 3^1, so MMK(90) = 12. As an operation on multisets MMK is represented by A367579, and as an operation on their ranks it is represented by A367580.
First differs from A212165 at n=73: A212165(73)=120 is not a term of this. - Amiram Eldar, Dec 04 2023

			The terms together with their prime indices begin:
    1: {}
    2: {1}
    3: {2}
    4: {1,1}
    5: {3}
    7: {4}
    8: {1,1,1}
    9: {2,2}
   11: {5}
   12: {1,1,2}
   13: {6}
   16: {1,1,1,1}
   17: {7}
   18: {1,2,2}
   19: {8}
   20: {1,1,3}
   23: {9}
   24: {1,1,1,2}

Includes all prime-powers A000961.
The only squarefree terms are the primes A008578.
Partitions of this type are counted by A367684.
A007947 gives squarefree kernel.
A027746 lists prime factors, length A001222, indices A112798.
A027748 lists distinct prime factors, length A001221, indices A304038.
A071625 counts distinct prime exponents.
A124010 gives multiset of multiplicities (prime signature), sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A367579 lists MMK, ranks A367580, sum A367581, max A367583.
Cf. A005117, A039956, A051904, A052409, A072774, A130091, A367582, A367584, A367586, A367682.

mmk[n_Integer]:= Product[Min[#]^Length[#]&[First/@Select[FactorInteger[n], Last[#]==k&]], {k,Union[Last/@FactorInteger[n]]}];
Select[Range[100], Divisible[#,mmk[#]]&]

A367859 Multiset multiplicity cokernel (MMC) of n. Product of (greatest prime factor with exponent k)^(number of prime factors with exponent k) over all distinct exponents k appearing in the prime factorization of n.

Original entry on oeis.org

1, 2, 3, 2, 5, 9, 7, 2, 3, 25, 11, 6, 13, 49, 25, 2, 17, 6, 19, 10, 49, 121, 23, 6, 5, 169, 3, 14, 29, 125, 31, 2, 121, 289, 49, 9, 37, 361, 169, 10, 41, 343, 43, 22, 15, 529, 47, 6, 7, 10, 289, 26, 53, 6, 121, 14, 361, 841, 59, 50, 61, 961, 21, 2, 169, 1331

Offset: 1

Gus Wiseman, Dec 03 2023

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.

We define the multiset multiplicity cokernel MMC(m) of a multiset m by the following property, holding for all distinct multiplicities k >= 1. If S is the set of elements of multiplicity k in m, then max(S) has multiplicity |S| in MMC(m). For example, MMC({1,1,2,2,3,4,5}) = {2,2,5,5,5}, and MMC({1,2,3,4,5,5,5,5}) = {4,4,4,4,5}. As an operation on multisets MMC is represented by A367858, and as an operation on their ranks it is represented by A367859.

			90 has prime factorization 2^1*3^2*5^1, so for k = 1 we have 5^2, and for k = 2 we have 3^1, so a(90) = 75.

Positions of 2's are A000079 without 1.
Positions of 3's are A000244 without 1.
Positions of primes (including 1) are A000961.
Depends only on rootless base A052410, see A007916.
Positions of prime powers are A072774.
Positions of squarefree numbers are A130091.
For kernel instead of cokernel we have A367580, ranks of A367579.
Rows of A367858 have this rank, sum A367860, max A061395, min A367587.
A007947 gives squarefree kernel.
A027746 lists prime factors, length A001222, indices A112798.
A027748 lists distinct prime factors, length A001221, indices A304038.
A071625 counts distinct prime exponents.
A124010 gives multiset of multiplicities (prime signature), sorted A118914.
Cf. A005117, A020639, A052409, A175781, A181819, A353742, A367584, A367585.

mmc[q_]:=With[{mts=Length/@Split[q]}, Sort[Table[Max@@Select[q,Count[q,#]==i&], {i,mts}]]];
Table[Times@@mmc[Join@@ConstantArray@@@FactorInteger[n]], {n,30}]

a(n^k) = a(n) for all positive integers n and k.
If n is squarefree, a(n) = A006530(n)^A001222(n).
A055396(a(n)) = A367587(n).
A056239(a(n)) = A367860(n).
A061395(a(n)) = A061395(n).
A001222(a(n)) = A001221(n).
A001221(a(n)) = A071625(n).
A071625(a(n)) = A323022(n).

A072414 Non-Achilles numbers for which LCM of the exponents in the prime factorization of n is not equal to the maximum of the same exponents.

Original entry on oeis.org

360, 504, 540, 600, 756, 792, 936, 1176, 1188, 1224, 1350, 1368, 1400, 1404, 1440, 1500, 1656, 1836, 1960, 2016, 2052, 2088, 2160, 2200, 2232, 2250, 2400, 2484, 2520, 2600, 2646, 2664, 2904, 2952, 3024, 3096, 3132, 3168, 3240, 3348, 3384, 3400, 3500

Offset: 1

Labos Elemer, Jun 17 2002

nonn

Most members of this sequence fail to be Achilles numbers because they have at least one prime factor with multiplicity 1. There are also numbers in the sequence that fail to be Achilles numbers because they are perfect powers: these are precisely the proper powers of members of A072412, so the smallest such is 5184 = 2^6*3^4 = 72^2. - Franklin T. Adams-Watters, Oct 09 2006

			m = 504 = 2*2*2*3*3*7: exponent-set = E = {3,2,1}, max(E) = 3 < lcm(E) = 6, gcd(E) = min(E) = 1.

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

Cf. A005361, A051903, A051904, A052409, A072411, A072412, A052486.

Select[Range@ 3500, And[LCM @@ # != Max@ #, GCD @@ # == Min@ # == 1] &[FactorInteger[#][[All, -1]] ] &] (* Michael De Vlieger, Jul 18 2017 *)

is(n)=my(f=factor(n)[,2]); n>9 && lcm(f)!=vecmax(f) && (#f==1 || vecmin(f)<2) \\ Charles R Greathouse IV, Oct 16 2015

A051903(a(n)) is not equal A072411(a(n)) but the numbers are not in A052486.

A239728 Perfect power but neither square nor cube.

Original entry on oeis.org

32, 128, 243, 2048, 2187, 3125, 7776, 8192, 16807, 78125, 100000, 131072, 161051, 177147, 248832, 279936, 371293, 524288, 537824, 759375, 823543, 1419857, 1594323, 1889568, 2476099, 3200000, 4084101, 5153632, 6436343, 7962624, 8388608, 10000000, 11881376, 17210368

Offset: 1

Jeppe Stig Nielsen, Mar 25 2014

			279936 is included since 279936 = 6^7 is a power and this is not a square or a cube.
59049 = 9^5 not included since this is a square 243^2 = 59049.
32768 = 8^5 not included since this is a cube 32^3 = 32768.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000

Cf. A001597 (perfect powers), A097054 (nonsquare perfect powers), A340585 (noncube perfect powers).

Cf. A008683, A052409.

for(i=1, 2^25, if(gcd(ispower(i), 6) == 1, print(i)))

from sympy import mobius, integer_nthroot
def A239728(n):
    def f(x): return int(n+x-integer_nthroot(x,5)[0]+sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(7,x.bit_length())))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 14 2024

GCD(A052409(a(n)), 6) = 1. - Reinhard Zumkeller, Mar 28 2014

Sum_{n>=1} 1/a(n) = 1 - zeta(2) - zeta(3) + zeta(6) + Sum_{k>=2} mu(k)*(1-zeta(k)) = 0.0448164603... - Amiram Eldar, Dec 21 2020

A294873 a(n) = Product_{d|n, d>1, d = x^(2k-1) for some maximal k >= 1} prime(k).

Original entry on oeis.org

1, 2, 2, 2, 2, 8, 2, 6, 2, 8, 2, 16, 2, 8, 8, 6, 2, 16, 2, 16, 8, 8, 2, 96, 2, 8, 6, 16, 2, 128, 2, 30, 8, 8, 8, 32, 2, 8, 8, 96, 2, 128, 2, 16, 16, 8, 2, 192, 2, 16, 8, 16, 2, 96, 8, 96, 8, 8, 2, 1024, 2, 8, 16, 30, 8, 128, 2, 16, 8, 128, 2, 384, 2, 8, 16, 16, 8, 128, 2, 192, 6, 8, 2, 1024, 8, 8, 8, 96, 2, 1024, 8, 16, 8, 8, 8, 1920, 2, 16, 16, 32, 2, 128, 2

Offset: 1

Antti Karttunen, Nov 11 2017

nonn

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

Cf. A000040, A056595, A052409, A183096.

Cf. A293524, A294874, A294875.

A294873(n) = { my(m=1,e); fordiv(n,d, if(d>1, e = ispower(d); if(!e, m += m, if((e>1)&&(e%2), m *= prime((e+1)/2))))); m; };

a(n) = Product_{d|n, d>1, r = A052409(d) is odd} A000040((r+1)/2).
Other identities. For all n >= 1:
A001222(a(n)) = A056595(n).
A007814(a(n)) = A183096(n).

A295920 Number of twice-factorizations of n of type (P,R,R).

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 8, 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, 17, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Nov 30 2017

nonn

a(n) is also the number of ways to choose a perfect divisor d|n and then a sequence of log_d(n) perfect divisors of d.

			The a(64) = 17 twice-factorizations are:
(2)*(2)*(2)*(2)*(2)*(2)  (2*2)*(2*2)*(2*2)  (2*2*2)*(2*2*2)  (2*2*2*2*2*2)
(2*2)*(2*2)*(4)          (2*2)*(4)*(2*2)    (4)*(2*2)*(2*2)
(2*2)*(4)*(4)            (4)*(2*2)*(4)      (4)*(4)*(2*2)
(2*2*2)*(8)              (8)*(2*2*2)
(4)*(4)*(4)              (4*4*4)
(8)*(8)                  (8*8)
(64)

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

Cf. A000005, A001055, A052409, A052410, A089723, A279789, A281113, A295923, A295924, A295931, A295935.

Table[Sum[Length[Divisors[GCD@@FactorInteger[n^(1/d)][[All,2]]]]^d,{d,Divisors[GCD@@FactorInteger[n][[All,2]]]}],{n,100}]

A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
A295920(n) = if(1==n,n,my(r); sumdiv(A052409(n), d, if(!ispower(n,d,&r),(1/0),numdiv(A052409(r))^d))); \\ Antti Karttunen, Dec 06 2018, after Mathematica-code

a(n) = Sum_{d|A052409(n)} A000005(A052409(n^(1/d)))^d. - Antti Karttunen, Dec 06 2018, after Mathematica-code

More terms from Antti Karttunen, Dec 06 2018

cp's OEIS Frontend

A303710 Number of factorizations of numbers that are not perfect powers using only numbers that are not perfect powers.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A304250 Perfect powers whose prime factors span an initial interval of prime numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A322901 Numbers whose prime indices are all powers of the same number.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A333941 Triangle read by rows where T(n,k) is the number of compositions of n with rotational period k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A367685 Numbers divisible by their multiset multiplicity kernel.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A367859 Multiset multiplicity cokernel (MMC) of n. Product of (greatest prime factor with exponent k)^(number of prime factors with exponent k) over all distinct exponents k appearing in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A072414 Non-Achilles numbers for which LCM of the exponents in the prime factorization of n is not equal to the maximum of the same exponents.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A239728 Perfect power but neither square nor cube.

Views

Author

Keywords

Examples

Links

Crossrefs