A007916 -id:A007916 - cp's OEIS frontend

A367580 Multiset multiplicity kernel (MMK) of n. 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.

1, 2, 3, 2, 5, 4, 7, 2, 3, 4, 11, 6, 13, 4, 9, 2, 17, 6, 19, 10, 9, 4, 23, 6, 5, 4, 3, 14, 29, 8, 31, 2, 9, 4, 25, 4, 37, 4, 9, 10, 41, 8, 43, 22, 15, 4, 47, 6, 7, 10, 9, 26, 53, 6, 25, 14, 9, 4, 59, 18, 61, 4, 21, 2, 25, 8, 67, 34, 9, 8, 71, 6, 73, 4, 15, 38

Offset: 1

Gus Wiseman, Nov 26 2023

nonn

As an operation on multisets, this is represented by A367579.

			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 a(90) = 12.

Positions of 2's are A000079 without 1.
Positions of 3's are A000244 without 1.
Positions of primes (including 1) are A000961.
Positions of prime(k) are prime powers prime(k)^i, rows of A051128.
Depends only on rootless base A052410, see A007916.
Positions of prime powers are A072774.
Positions of squarefree numbers are A130091.
Agrees with A181819 at positions A367683, counted by A367682.
Rows of A367579 have this rank, sum A367581, max A367583, min A055396.
Positions of first appearances are A367584, sorted A367585.
Positions of powers of 2 are A367586.
Divides n at positions A367685, counted by A367684.
The opposite version (cokernel) is A367859.
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. A001597, A005117, A020639, A051904, A052409, A062770, A175781, A238747, A288636, A353742, A367582.

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

a(n^k) = a(n) for all positive integers n and k.
A001221(a(n)) = A071625(n).
A001222(a(n)) = A001221(n).
If n is squarefree, a(n) = A020639(n)^A001222(n).
A056239(a(n)) = A367581(n).

A375702 Length of the n-th maximal run of adjacent (increasing by one at a time) non-perfect-powers.

Original entry on oeis.org

2, 3, 6, 8, 1, 4, 3, 12, 14, 16, 18, 20, 3, 2, 15, 24, 26, 19, 8, 17, 12, 32, 34, 18, 17, 38, 40, 42, 27, 16, 46, 48, 50, 52, 54, 56, 58, 60, 38, 23, 64, 66, 68, 70, 34, 37, 74, 76, 78, 80, 46, 35, 84, 86, 88, 22, 67, 70, 9, 11, 94, 96, 98, 100, 102, 39, 64

Offset: 1

Gus Wiseman, Aug 27 2024

nonn

Non-perfect-powers A007916 are numbers with no proper integer roots.

			The list of all non-perfect-powers, split into runs, begins:
   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
Row n has length a(n), first A375703, last A375704, sum A375705.

For nonsquarefree numbers we have A053797, anti-runs A373409.
For squarefree numbers we have A120992, anti-runs A373127.
For nonprime numbers we have A176246, anti-runs A373403.
For prime-powers we have A373675, anti-runs A373576.
For non-prime-powers we have A373678, anti-runs A373679.
The anti-run version is A375736, sum A375737.
For runs of non-perfect-powers (A007916):
- length: A375702 (this).
- first: A375703
- last: A375704
- sum: A375705
A001597 lists perfect-powers, differences A053289.
A007916 lists non-perfect-powers, differences A375706.
A046933 counts composite numbers between primes.
Cf. A007674, A045542, A061399, A216765, A251092, A375708, A375714.

radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
Length/@Split[Select[Range[100],radQ],#1+1==#2&]//Most

For n > 2 we have a(n) = A053289(n+1) - 1.

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}]

A304464 Start with the normalized multiset of prime factors of n > 1. Given a multiset, take the multiset of its multiplicities. Repeat this until a multiset of size 1 is obtained. a(n) is the unique element of this multiset.

Original entry on oeis.org

0, 1, 2, 2, 3, 2, 4, 3, 2, 2, 5, 2, 6, 2, 2, 4, 7, 2, 8, 2, 2, 2, 9, 2, 2, 2, 3, 2, 10, 3, 11, 5, 2, 2, 2, 2, 12, 2, 2, 2, 13, 3, 14, 2, 2, 2, 15, 2, 2, 2, 2, 2, 16, 2, 2, 2, 2, 2, 17, 2, 18, 2, 2, 6, 2, 3, 19, 2, 2, 3, 20, 2, 21, 2, 2, 2, 2, 3, 22, 2, 4, 2, 23

Offset: 1

Gus Wiseman, May 13 2018

nonn

a(1) = 0 by convention.

			Starting with the normalized multiset of prime factors of 360, we obtain {1,1,1,2,2,3} -> {1,2,3} -> {1,1,1} -> {3}, so a(360) = 3.

Cf. A000005, A001222, A001597, A005117, A007916, A055932, A056239, A112798, A181819, A182850, A182857, A275870, A296150, A303945, A304465.

Table[If[n===1,0,NestWhile[Sort[Length/@Split[#]]&,If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]],Length[#]>1&]//First],{n,100}]

a(prime(n)) = n.
a(p^n) = n where p is any prime number and n > 1.
a(product of n > 1 distinct primes) = n.

A317658 Number of positions in the n-th free pure symmetric multifunction (with empty expressions allowed) with one atom.

Original entry on oeis.org

1, 2, 3, 3, 4, 4, 5, 4, 4, 5, 6, 5, 5, 6, 7, 4, 6, 6, 7, 8, 5, 7, 7, 8, 5, 9, 5, 6, 8, 8, 9, 5, 6, 10, 6, 5, 7, 9, 9, 10, 6, 7, 11, 7, 6, 8, 10, 10, 6, 11, 7, 8, 12, 8, 7, 9, 11, 11, 7, 12, 8, 9, 13, 5, 9, 8, 10, 12, 12, 8, 13, 9, 10, 14, 6, 10, 9, 11, 13, 13

Offset: 1

Gus Wiseman, Aug 03 2018

nonn

Given a positive integer n > 1 we construct a unique free pure symmetric multifunction e(n) by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)].

Also the number of positions in the orderless Mathematica expression with e-number n.

			The first twenty Mathematica expressions:
   1: o
   2: o[]
   3: o[][]
   4: o[o]
   5: o[][][]
   6: o[o][]
   7: o[][][][]
   8: o[o[]]
   9: o[][o]
  10: o[o][][]
  11: o[][][][][]
  12: o[o[]][]
  13: o[][o][]
  14: o[o][][][]
  15: o[][][][][][]
  16: o[o,o]
  17: o[o[]][][]
  18: o[][o][][]
  19: o[o][][][][]
  20: o[][][][][][][]

Mathematica Reference, Orderless

First differs from A277615 at a(128) = 5, A277615(128) = 6.
Cf. A001003, A052893, A053492, A255906, A277996, A279944, A280000.
Cf. A317652, A317653, A317654, A317655, A317656.

nn=100;
radQ[n_]:=If[n===1,False,GCD@@FactorInteger[n][[All,2]]===1];
rad[n_]:=rad[n]=If[n===0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
Clear[radPi];Set@@@Array[radPi[rad[#]]==#&,nn];
exp[n_]:=If[n===1,x,With[{g=GCD@@FactorInteger[n][[All,2]]},Apply[exp[radPi[Power[n,1/g]]],exp/@Flatten[Cases[FactorInteger[g],{p_?PrimeQ,k_}:>ConstantArray[PrimePi[p],k]]]]]];
Table[exp[n],{n,1,nn}]

a(rad(x)^(prime(y_1) * ... * prime(y_k))) = a(x) + a(y_1) + ... + a(y_k).
e(2^(2^n)) = o[o,...,o].
e(2^prime(2^prime(2^...))) = o[o[...o[o]]].
e(rad(rad(rad(...)^2)^2)^2) = o[o][o]...[o].

A336415 Number of divisors of n! with equal prime multiplicities.

Original entry on oeis.org

1, 1, 2, 4, 6, 10, 13, 21, 24, 28, 33, 49, 53, 85, 94, 100, 104, 168, 173, 301, 307, 317, 334, 590, 595, 603, 636, 642, 652, 1164, 1171, 2195, 2200, 2218, 2283, 2295, 2301, 4349, 4478, 4512, 4519, 8615, 8626, 16818, 16836, 16844, 17101, 33485, 33491, 33507, 33516, 33582

Offset: 0

Gus Wiseman, Jul 22 2020

nonn

A number k has "equal prime multiplicities" (or is "uniform") iff its prime signature is constant, meaning that k is a power of a squarefree number.

			The a(n) uniform divisors of n for n = 1, 2, 6, 8, 30, 36 are the columns:
  1  2  6  8  30  36
     1  3  6  15  30
        2  4  10  16
        1  3   8  15
           2   6  10
           1   5   9
               4   8
               3   6
               2   5
               1   4
                   3
                   2
                   1
In 20!, the multiplicity of the third prime (5) is 4 but the multiplicity of the fourth prime (7) is 2. Hence there are 2^3 - 1 = 3 divisors with all exponents 3 (we subtract |{1}| = 1 from that count as 1 has no exponent 3). - _David A. Corneth_, Jul 27 2020

David A. Corneth, Table of n, a(n) for n = 0..9999
Jon Maiga, Computer-generated formulas for A336415, Sequence Machine.

The version for distinct prime multiplicities is A336414.
The version for nonprime perfect powers is A336416.
Uniform partitions are counted by A047966.
Uniform numbers are A072774, with nonprime terms A182853.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
Maximum divisor with distinct prime multiplicities is A327498.
Uniform divisors are counted by A327527.
Maximum uniform divisor is A336618.
1st differences are given by A048675.
Cf. A000005, A000961, A001222, A001597, A007916, A112798, A124010, A327499.
Factorial numbers: A000142, A007489, A022559, A027423, A048656, A071626, A108731, A325272, A325273, A325617.

Table[Length[Select[Divisors[n!],SameQ@@Last/@FactorInteger[#]&]],{n,0,15}]

a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); (#ex==0) || (vecmin(ex) == vecmax(ex))); \\ Michel Marcus, Jul 24 2020

a(n) = {if(n<2, return(1)); my(f = primes(primepi(n)), res = 1, t = #f); f = vector(#f, i, val(n, f[i])); for(i = 1, f[1], while(f[t] < i, t--; ); res+=(1<David A. Corneth, Jul 27 2020

a(n) = A327527(n!).

Terms a(31) and onwards from David A. Corneth, Jul 27 2020

A375714 Positions of non-successions of consecutive non-perfect-powers. Numbers k such that the k-th non-perfect-power is at least two fewer than the next.

Original entry on oeis.org

2, 5, 11, 19, 20, 24, 27, 39, 53, 69, 87, 107, 110, 112, 127, 151, 177, 196, 204, 221, 233, 265, 299, 317, 334, 372, 412, 454, 481, 497, 543, 591, 641, 693, 747, 803, 861, 921, 959, 982, 1046, 1112, 1180, 1250, 1284, 1321, 1395, 1471, 1549, 1629, 1675, 1710

Offset: 1

Gus Wiseman, Sep 10 2024

nonn

Non-perfect-powers (A007916) are numbers with no proper integer roots.

			The initial non-perfect-powers are 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, which increase by more than one after term 2, term 5, term 11, etc.

First differences are A375702.
Positions of terms > 1 in A375706 (differences of A007916).
The complement for non-prime-powers is A375713, differences A373672.
The complement is A375740.
The version for non-prime-powers is A375928, differences A110969.
Prime-powers inclusive:
- terms: A000961
- differences: A057820
- runs: A373675, A373673, A373674, A174965
- antiruns: A373576, A120430, A006549, A373671
Non-prime-powers inclusive:
- terms: A361102
- differences: A375708
- runs: A373678, A373676, A373677, A110969
- antiruns: A373679, A373575, A255346, A373672
Cf. A000040, A001223, A046933, A053289, A073783, A246655, A251092, A375734.

radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
ce=Select[Range[100],radQ];
Select[Range[Length[ce]-1],!ce[[#+1]]==ce[[#]]+1&]

from itertools import count, islice
from sympy import perfect_power
def A375714_gen(): # generator of terms
    a, b = -1, 0
    for n in count(1):
        c = not perfect_power(n)
        if c:
            a += 1
        if b&(c^1):
            yield a
        b = c
A375714_list = list(islice(A375714_gen(),52)) # Chai Wah Wu, Sep 11 2024

A007916(a(n)+1) - A007916(a(n)) > 1.

A376340 Sorted positions of first appearances in A057820, the sequence of first differences of prime-powers.

Original entry on oeis.org

1, 4, 9, 12, 18, 24, 34, 47, 60, 79, 117, 178, 198, 206, 215, 244, 311, 402, 465, 614, 782, 1078, 1109, 1234, 1890, 1939, 1961, 2256, 2290, 3149, 3377, 3460, 3502, 3722, 3871, 4604, 4694, 6634, 8073, 8131, 8793, 12370, 12661, 14482, 14990, 15912, 17140, 19166

Offset: 1

Gus Wiseman, Sep 22 2024

nonn

			The terms together with their prime indices begin:
     1: {}
     4: {1,1}
     9: {2,2}
    12: {1,1,2}
    18: {1,2,2}
    24: {1,1,1,2}
    34: {1,7}
    47: {15}
    60: {1,1,2,3}
    79: {22}
   117: {2,2,6}
   178: {1,24}
   198: {1,2,2,5}
   206: {1,27}
   215: {3,14}
   244: {1,1,18}

For compression instead of sorted firsts we have A376308.
For run-lengths instead of sorted firsts we have A376309.
For run-sums instead of sorted firsts we have A376310.
The version for squarefree numbers is the unsorted version of A376311.
The unsorted version is A376341.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, first differences A057820.
A003242 counts compressed compositions, ranks A333489.
A005117 lists squarefree numbers, differences A076259.
A024619 and A361102 list non-prime-powers, first differences A375708.
A116861 counts partitions by compressed sum, by compressed length A116608.
Cf. A001597, A006549, A007916, A025475, A037201, A053289, A078147, A110969, A120430, A174965, A373948, A375706.

q=Differences[Select[Range[100],PrimePowerQ]];
Select[Range[Length[q]],!MemberQ[Take[q,#-1],q[[#]]]&]

A377054 First term of the n-th differences of the powers of primes. Inverse zero-based binomial transform of A000961.

Original entry on oeis.org

1, 1, 0, 0, 0, 1, -5, 15, -34, 63, -97, 115, -54, -251, 1184, -3536, 8736, -18993, 37009, -64545, 98442, -121393, 82008, 147432, -860818, 2710023, -7110594, 17077281, -38873146, 85085287, -179965647, 367885014, -725051280, 1372311999, -2481473550, 4257624252

Offset: 0

Gus Wiseman, Oct 22 2024

sign

			The sixth differences of A000961 begin: -5, 10, -9, 1, 6, -10, 16, -18, ..., so a(6) = -5.

The version for primes is A007442, noncomposites A030016, composites A377036.
For squarefree numbers we have A377041, nonsquarefree A377049.
This is the first column of the array A377051.
For antidiagonal-sums we have A377052, absolute A377053.
For positions of first zeros we have A377055.
A000040 lists the primes, differences A001223, seconds A036263.
A000961 lists the powers of primes, differences A057820.
A001597 lists perfect-powers, complement A007916.
A008578 lists the noncomposites, differences A075526.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.
Cf. A025475, A053707, A093555, A174965, A175804, A246655, A361102, A376340, A376596, A377056.

q=Select[Range[100],#==1||PrimePowerQ[#]&];
Table[Sum[(-1)^(j-k)*Binomial[j,k]*q[[1+k]],{k,0,j}],{j,0,Length[q]/2}]

The inverse zero-based binomial transform of a sequence (q(0), q(1), q(2), ...) is the sequence p given by:

p(j) = sum_{k=0..j} (-1)^(j-k)*binomial(j,k)*q(k)

A377703 First differences of the sequence A345531(k) = least prime-power greater than the k-th prime.

Original entry on oeis.org

1, 3, 1, 5, 3, 3, 4, 2, 6, 1, 9, 2, 4, 2, 10, 2, 3, 7, 2, 6, 2, 8, 8, 4, 2, 4, 2, 4, 8, 7, 9, 2, 10, 2, 6, 6, 4, 2, 10, 2, 10, 2, 4, 2, 12, 12, 4, 2, 4, 6, 2, 2, 13, 7, 6, 2, 6, 4, 2, 6, 18, 4, 2, 4, 14, 6, 6, 6, 4, 6, 2, 12, 6, 4, 6, 8, 4, 8, 10, 2, 10, 2, 6

Offset: 1

Gus Wiseman, Nov 07 2024

nonn

What is the union of this sequence? In particular, does it contain 17?

First differences of A345531.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect-powers, differences A053289, seconds A376559.
A007916 lists the non-perfect-powers, differences A375706, seconds A376562.
A024619 lists the non-prime-powers, differences A375735, seconds A376599.
A080101 counts prime-powers between primes (exclusive).
A246655 lists the prime-powers, differences A057820 without first term.
A361102 lists the non-powers of primes, differences A375708.
A366833 counts prime-powers between primes, see A053607, A304521, A377057 (positive), A377286 (zero), A377287 (one), A377288 (two).
A377432 counts perfect-powers between primes, see A377434 (one), A377436 (zero), A377466 (multiple).
Cf. A000040, A008864, A031218, A377281, A377282, A377286, A377289, A377468.

Differences[Table[NestWhile[#+1&, Prime[n]+1,!PrimePowerQ[#]&],{n,100}]]

from sympy import factorint, prime, nextprime
def A377703(n): return -next(filter(lambda m:len(factorint(m))<=1, count((p:=prime(n))+1)))+next(filter(lambda m:len(factorint(m))<=1, count(nextprime(p)+1))) # Chai Wah Wu, Nov 14 2024

cp's OEIS Frontend

A367580 Multiset multiplicity kernel (MMK) of n. 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.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A375702 Length of the n-th maximal run of adjacent (increasing by one at a time) non-perfect-powers.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A304464 Start with the normalized multiset of prime factors of n > 1. Given a multiset, take the multiset of its multiplicities. Repeat this until a multiset of size 1 is obtained. a(n) is the unique element of this multiset.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A317658 Number of positions in the n-th free pure symmetric multifunction (with empty expressions allowed) with one atom.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A336415 Number of divisors of n! with equal prime multiplicities.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A375714 Positions of non-successions of consecutive non-perfect-powers. Numbers k such that the k-th non-perfect-power is at least two fewer than the next.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

Formula

A376340 Sorted positions of first appearances in A057820, the sequence of first differences of prime-powers.

Views

Author

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.