A297155 -id:A297155 - cp's OEIS frontend

A156552 Unary-encoded compressed factorization of natural numbers.

0, 1, 2, 3, 4, 5, 8, 7, 6, 9, 16, 11, 32, 17, 10, 15, 64, 13, 128, 19, 18, 33, 256, 23, 12, 65, 14, 35, 512, 21, 1024, 31, 34, 129, 20, 27, 2048, 257, 66, 39, 4096, 37, 8192, 67, 22, 513, 16384, 47, 24, 25, 130, 131, 32768, 29, 36, 71, 258, 1025, 65536, 43, 131072, 2049, 38, 63, 68, 69, 262144

Offset: 1

Leonid Broukhis, Feb 09 2009

The primes become the powers of 2 (2 -> 1, 3 -> 2, 5 -> 4, 7 -> 8); the composite numbers are formed by taking the values for the factors in the increasing order, multiplying them by the consecutive powers of 2, and summing. See the Example section.
From Antti Karttunen, Jun 27 2014: (Start)
The odd bisection (containing even terms) halved gives A244153.
The even bisection (containing odd terms), when one is subtracted from each and halved, gives this sequence back.
(End)
Question: Are there any other solutions that would satisfy the recurrence r(1) = 0; and for n > 1, r(n) = Sum_{d|n, d>1} 2^A033265(r(d)), apart from simple variants 2^k * A156552(n)? See also A297112, A297113. - Antti Karttunen, Dec 30 2017

			For 84 = 2*2*3*7 -> 1*1 + 1*2 + 2*4 + 8*8 =  75.
For 105 = 3*5*7 -> 2*1 + 4*2 + 8*4 = 42.
For 137 = p_33 -> 2^32 = 4294967296.
For 420 = 2*2*3*5*7 -> 1*1 + 1*2 + 2*4 + 4*8 + 8*16 = 171.
For 147 = 3*7*7 = p_2 * p_4 * p_4 -> 2*1 + 8*2 + 8*4 = 50.

David A. Corneth, Table of n, a(n) for n = 1..10000 (first 1024 terms from Antti Karttunen)
Hans Havermann, Factorization of the first 10000 terms, in format [[primes], [exponents]]
A puzzle by Sergey Orlov (in Russian)
Index entries for sequences related to binary expansion of n
Index entries for sequences that are permutations of the natural numbers
Index entries for sequences computed from indices in prime factorization

One less than A005941.
Inverse permutation: A005940 with starting offset 0 instead of 1.
Cf. A000079, A000120, A001222, A052126, A054429, A061395, A064216, A064989, A003188, A243071, A243065-A243066, A244153, A243354, A112798, A125106, A056239, A161511.
Cf. also A297106, A297112 (Möbius transform), A297113, A153013, A290308, A300827, A323243, A323244, A323247, A324201, A324812 (n for which a(n) is a square), A324813, A324822, A324823, A324398, A324713, A324815, A324819, A324865, A324866, A324867.
Other related permutations: A253551, A253792, A253564, A253791, A277195, A297163, A297164, A297165, A297166, A302023, A305418, A322863, A322864.

Table[Floor@ Total@ Flatten@ MapIndexed[#1 2^(#2 - 1) &, Flatten[ Table[2^(PrimePi@ #1 - 1), {#2}] & @@@ FactorInteger@ n]], {n, 67}] (* Michael De Vlieger, Sep 08 2016 *)

a(n) = {my(f = factor(n), p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ David A. Corneth, Mar 08 2019

A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n)))); \\ (based on the given recurrence) - Antti Karttunen, Mar 08 2019

# Program corrected per instructions from Leonid Broukhis. - Antti Karttunen, Jun 26 2014
# However, it gives correct answers only up to n=136, before corruption by a wrap-around effect.
# Note that the correct answer for n=137 is A156552(137) = 4294967296.
$max = $ARGV[0];
$pow = 0;
foreach $i (2..$max) {
@a = split(/ /, `factor $i`);
shift @a;
$shift = 0;
$cur = 0;
while ($n = int shift @a) {
$prime{$n} = 1 << $pow++ if !defined($prime{$n});
$cur |= $prime{$n} << $shift++;
}
print "$cur, ";
}
print "\n";
(Scheme, with memoization-macro definec from Antti Karttunen's IntSeq-library, two different implementations)
(definec (A156552 n) (cond ((= n 1) 0) (else (+ (A000079 (+ -2 (A001222 n) (A061395 n))) (A156552 (A052126 n))))))
(definec (A156552 n) (cond ((= 1 n) (- n 1)) ((even? n) (+ 1 (* 2 (A156552 (/ n 2))))) (else (* 2 (A156552 (A064989 n))))))
;; Antti Karttunen, Jun 26 2014

from sympy import primepi, factorint
def A156552(n): return sum((1<Chai Wah Wu, Mar 10 2023

From Antti Karttunen, Jun 26 2014: (Start)
a(1) = 0, a(n) = A000079(A001222(n)+A061395(n)-2) + a(A052126(n)).
a(1) = 0, a(2n) = 1+2*a(n), a(2n+1) = 2*a(A064989(2n+1)). [Compare to the entanglement recurrence A243071].
For n >= 0, a(2n+1) = 2*A244153(n+1). [Follows from the latter clause of the above formula.]
a(n) = A005941(n) - 1.
As a composition of related permutations:
a(n) = A003188(A243354(n)).
a(n) = A054429(A243071(n)).
For all n >= 1, A005940(1+a(n)) = n and for all n >= 0, a(A005940(n+1)) = n. [The offset-0 version of A005940 works as an inverse for this permutation.]
This permutations also maps between the partition-lists A112798 and A125106:
A056239(n) = A161511(a(n)). [The sums of parts of each partition (the total sizes).]
A003963(n) = A243499(a(n)). [And also the products of those parts.]
(End)
From Antti Karttunen, Oct 09 2016: (Start)
A161511(a(n)) = A056239(n).
A029837(1+a(n)) = A252464(n). [Binary width of terms.]
A080791(a(n)) = A252735(n). [Number of nonleading 0-bits.]
A000120(a(n)) = A001222(n). [Binary weight.]
For all n >= 2, A001511(a(n)) = A055396(n).
For all n >= 2, A000120(a(n))-1 = A252736(n). [Binary weight minus one.]
A252750(a(n)) = A252748(n).
a(A250246(n)) = A252754(n).
a(A005117(n)) = A277010(n). [Maps squarefree numbers to a permutation of A003714, fibbinary numbers.]
A085357(a(n)) = A008966(n). [Ditto for their characteristic functions.]
For all n >= 0:
a(A276076(n)) = A277012(n).
a(A276086(n)) = A277022(n).
a(A260443(n)) = A277020(n).
(End)
From Antti Karttunen, Dec 30 2017: (Start)
For n > 1, a(n) = Sum_{d|n, d>1} 2^A033265(a(d)). [See comments.]
More linking formulas:
A106737(a(n)) = A000005(n).
A290077(a(n)) = A000010(n).
A069010(a(n)) = A001221(n).
A136277(a(n)) = A181591(n).
A132971(a(n)) = A008683(n).
A106400(a(n)) = A008836(n).
A268411(a(n)) = A092248(n).
A037011(a(n)) = A010052(n) [conjectured, depends on the exact definition of A037011].
A278161(a(n)) = A046951(n).
A001316(a(n)) = A061142(n).
A277561(a(n)) = A034444(n).
A286575(a(n)) = A037445(n).
A246029(a(n)) = A181819(n).
A278159(a(n)) = A124859(n).
A246660(a(n)) = A112624(n).
A246596(a(n)) = A069739(n).
A295896(a(n)) = A053866(n).
A295875(a(n)) = A295297(n).
A284569(a(n)) = A072411(n).
A286574(a(n)) = A064547(n).
A048735(a(n)) = A292380(n).
A292272(a(n)) = A292382(n).
A244154(a(n)) = A048673(n), a(A064216(n)) = A244153(n).
A279344(a(n)) = A279339(n), a(A279338(n)) = A279343(n).
a(A277324(n)) = A277189(n).
A037800(a(n)) = A297155(n).
For n > 1, A033265(a(n)) = 1+A297113(n).
(End)
From Antti Karttunen, Mar 08 2019: (Start)
a(n) = A048675(n) + A323905(n).
a(A324201(n)) = A000396(n), provided there are no odd perfect numbers.
The following sequences are derived from or related to the base-2 expansion of a(n):
A000265(a(n)) = A322993(n).
A002487(a(n)) = A323902(n).
A005187(a(n)) = A323247(n).
A324288(a(n)) = A324116(n).
A323505(a(n)) = A323508(n).
A079559(a(n)) = A323512(n).
A085405(a(n)) = A323239(n).
The following sequences are obtained by applying to a(n) a function that depends on the prime factorization of its argument, which goes "against the grain" because a(n) is the binary code of the factorization of n, which in these cases is then factored again:
A000203(a(n)) = A323243(n).
A033879(a(n)) = A323244(n) = 2*a(n) - A323243(n),
A294898(a(n)) = A323248(n).
A000005(a(n)) = A324105(n).
A000010(a(n)) = A324104(n).
A083254(a(n)) = A324103(n).
A001227(a(n)) = A324117(n).
A000593(a(n)) = A324118(n).
A001221(a(n)) = A324119(n).
A009194(a(n)) = A324396(n).
A318458(a(n)) = A324398(n).
A192895(a(n)) = A324100(n).
A106315(a(n)) = A324051(n).
A010052(a(n)) = A324822(n).
A053866(a(n)) = A324823(n).
A001065(a(n)) = A324865(n) = A323243(n) - a(n),
A318456(a(n)) = A324866(n) = A324865(n) OR a(n),
A318457(a(n)) = A324867(n) = A324865(n) XOR a(n),
A318458(a(n)) = A324398(n) = A324865(n) AND a(n),
A318466(a(n)) = A324819(n) = A323243(n) OR 2*a(n),
A318467(a(n)) = A324713(n) = A323243(n) XOR 2*a(n),
A318468(a(n)) = A324815(n) = A323243(n) AND 2*a(n).
(End)

More terms from Antti Karttunen, Jun 28 2014

A252464 a(1) = 0, a(2n) = 1 + a(n), a(2n+1) = 1 + a(A064989(2n+1)); also binary width of terms of A156552 and A243071.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Dec 20 2014

nonn

a(n) tells how many iterations of A252463 are needed before 1 is reached, i.e., the distance of n from 1 in binary trees like A005940 and A163511.

Similarly for A253553 in trees A253563 and A253565. - Antti Karttunen, Apr 14 2019

			From _Gus Wiseman_, Apr 02 2019: (Start)
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so a(n) is the size of the inner lining of the integer partition with Heinz number n, which is also the size of the largest hook of the same partition. For example, the partition with Heinz number 715 is (6,5,3), with diagram
  o o o o o o
  o o o o o
  o o o
which has inner lining
          o o
      o o o
  o o o
and largest hook
  o o o o o o
  o
  o
both of which have size 8, so a(715) = 8.
(End)

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

Cf. A000040, A000079, A001221, A001222, A005940, A029837, A061395, A156552 (A005941), A163511, A243071, A252461, A252463, A252735, A252736, A252759, A253553, A253563, A253565, A297113, A297155, A297167, A324870, A324872.
Right edge of irregular triangle A265146.
Cf. also A246348.
Cf. A052126, A056239, A093641, A112798, A257990, A325133, A325134, A325135.

Table[If[n==1,1,PrimeOmega[n]+PrimePi[FactorInteger[n][[-1,1]]]]-1,{n,100}] (* Gus Wiseman, Apr 02 2019 *)

A061395(n) = if(n>1, primepi(vecmax(factor(n)[, 1])), 0);
A252464(n) = (bigomega(n) + A061395(n) - 1); \\ Antti Karttunen, Apr 14 2019

from sympy import primepi, primeomega, primefactors
def A252464(n): return primeomega(n)+primepi(max(primefactors(n)))-1 if n>1 else 0 # Chai Wah Wu, Jul 17 2023

a(1) = 0; for n > 1: a(n) = 1 + a(A252463(n)).
a(n) = A029837(1+A243071(n)). [a(n) = binary width of terms of A243071.]
a(n) = A029837(A005941(n)) = A029837(1+A156552(n)). [Also binary width of terms of A156552.]
Other identities. For all n >= 1:
a(A000040(n)) = n.
a(A001248(n)) = n+1.
a(A030078(n)) = n+2.
And in general, a(prime(n)^k) = n+k-1.
a(A000079(n)) = n. [I.e., a(2^n) = n.]
For all n >= 2:
a(n) = A001222(n) + A061395(n) - 1 = A001222(n) + A252735(n) = A061395(n) + A252736(n) = 1 + A252735(n) + A252736(n).
a(n) = A325134(n) - 1. - Gus Wiseman, Apr 02 2019
From Antti Karttunen, Apr 14 2019: (Start)
a(1) = 0; for n > 1: a(n) = 1 + a(A253553(n)).
a(n) = A001221(n) + A297167(n) = A297113(n) + A297155(n).
(End).

A297113 a(1) = 0, a(2) = 1, after which, a(n) = a(n/2) if n is of the form 4k+2, and otherwise a(n) = 1+a(A252463(n)).

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Dec 26 2017

nonn

From Gus Wiseman, Apr 06 2019: (Start)
Also the number of squares in the Young diagram of the integer partition with Heinz number n that are graph-distance 1 from the lower-right boundary, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). For example, the partition (6,5,5,3) with Heinz number 7865 has diagram
  o o o o o o
  o o o o o
  o o o o o
  o o o
with inner rim
            o
          o
        o o
  o o o
of size 7, so a(7867) = 7.
(End)

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

One more than A297167 (after the initial term).
Cf. A001511, A008683, A033265, A064989, A156552, A252463, A252464, A297112, A297155.
Cf. also A297157, A297161, A297162.
Cf. A052126, A065770, A112798, A115994, A174090, A325166, A325167, A325169.

Table[If[n==1,0,PrimePi[FactorInteger[n][[-1,1]]]+PrimeOmega[n]-PrimeNu[n]],{n,100}] (* Gus Wiseman, Apr 06 2019 *)

A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A297113(n) = if(n<=2,n-1,if(n%2,1+A297113(A064989(n)), !(n%4)+A297113(n/2)));

\\ More complex way, after Moebius transform:
A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));
A297112(n) = sumdiv(n,d,moebius(n/d)*A156552(d));
A297113(n) = if(1==n,0,1+valuation(A297112(n),2));

;; With memoization-macro definec.
(definec (A297113 n) (cond ((<= n 2) (- n 1)) ((= 2 (modulo n 4)) (A297113 (/ n 2))) (else (+ 1 (A297113 (A252463 n))))))

a(1) = 0, a(2) = 1, after which, a(n) = a(n/2) if n is of the form 4k+2, and otherwise a(n) = 1+a(A252463(n)) .
For n > 1, a(n) = A001511(A297112(n)), where A297112(n) = Sum_{d|n} moebius(n/d)*A156552(d).
a(n) = A252464(n) - A297155(n).
For n > 1, a(n) = 1+A033265(A156552(n)) = 1+A297167(n) = A046660(n) + A061395(n). - Last two sums added by Antti Karttunen, Sep 02 2018
Other identities. For all n >= 1:
a(A000040(n)) = n. [Each n occurs for the first time at the n-th prime.]

A386577 Irregular triangle read by rows where T(n,k) is the number of permutations of the multiset of prime factors of n with k adjacent equal terms.

Original entry on oeis.org

1, 1, 0, 1, 1, 2, 0, 1, 0, 0, 1, 0, 1, 2, 0, 1, 1, 2, 0, 1, 2, 0, 2, 0, 0, 0, 0, 1, 1, 1, 2, 0, 1, 1, 2, 0, 2, 0, 2, 0, 1, 0, 2, 2, 0, 0, 1, 2, 0, 0, 0, 1, 1, 2, 0, 1, 6, 0, 0, 1, 0, 0, 0, 0, 1, 2, 0, 2, 0, 2, 0, 2, 2, 2, 0, 1, 2, 0, 2, 0, 0, 2, 2, 0, 1

Offset: 1

Gus Wiseman, Aug 01 2025

Are the rows all unimodal?

Counts permutations of prime factors by "inseparability". For "separability" we have A374252.

			The prime indices of 12 are {1,1,2}, and we have:
- 1 permutation (1,2,1) with 0 adjacent equal parts
- 2 permutations (1,1,2), (2,1,1) with 1 adjacent equal part
- 0 permutations with 2 adjacent equal parts
so row 12 is (1,2,0).
Row 48 counts the following permutations:
  .  .  (1,1,1,2,1)  (1,1,1,1,2)  .
        (1,1,2,1,1)  (2,1,1,1,1)
        (1,2,1,1,1)
Row 144 counts the following permutations:
  .  (1,1,2,1,2,1)  (1,1,1,2,1,2)  (1,1,1,2,2,1)  (1,1,1,1,2,2)  .
     (1,2,1,1,2,1)  (1,1,2,1,1,2)  (1,1,2,2,1,1)  (2,2,1,1,1,1)
     (1,2,1,2,1,1)  (1,2,1,1,1,2)  (1,2,2,1,1,1)
                    (2,1,1,1,2,1)  (2,1,1,1,1,2)
                    (2,1,1,2,1,1)
                    (2,1,2,1,1,1)
Triangle begins:
   1:
   2: 1
   3: 1
   4: 0  1
   6: 1
   6: 2  0
   7: 1
   8: 0  0  1
   9: 0  1
  10: 2  0
  11: 1
  12: 1  2  0
  13: 1
  14: 2  0
  15: 2  0
  16: 0  0  0  1
  17: 1
  18: 1  2  0
  19: 1
  20: 1  2  0
  21: 2  0
  22: 2  0
  23: 1
  24: 0  2  2  0

Row lengths are A001222.
The minima of each row are A010051.
Sorted positions of first appearances appear to be A025487.
Column k = last is A069513.
Row sums are A168324 or A008480.
The number of trailing zeros in each row is A297155 = A001221-1.
Column k = 1 is A335452.
The number of leading zeros in each row is A374246.
For separability instead of inseparability we have A374252.
For a multiset with prescribed multiplicities we have A386578, separability A386579.
A003242 and A335452 count anti-runs, ranks A333489, patterns A005649.
A025065(n - 2) counts partitions of inseparable type, ranks A335126, sums of A386586.
A124762 gives inseparability of standard compositions, separability A333382.
A325534 counts separable multisets, ranks A335433, sums of A386583.
A325535 counts inseparable multisets, ranks A335448, sums of A386584.
A336106 counts partitions of separable type, ranks A335127, sums of A386585.
Cf. A001221, A051903, A051904, A106351, A238130, A336102, A373957, A382525, A386581, A386632.

Table[Length[Select[Permutations[Flatten[ConstantArray@@@FactorInteger[n]]],Function[q,Length[Select[Range[Length[q]-1],q[[#]]==q[[#+1]]&]]==k]]],{n,30},{k,0,PrimeOmega[n]-1}]

A355538 Partial sum of A001221 (number of distinct prime factors) minus 1, ranging from 2 to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 5, 5, 6, 6, 7, 8, 9, 9, 10, 10, 11, 11, 12, 12, 14, 14, 14, 15, 16, 17, 18, 18, 19, 20, 21, 21, 23, 23, 24, 25, 26, 26, 27, 27, 28, 29, 30, 30, 31, 32, 33, 34, 35, 35, 37, 37, 38, 39, 39, 40, 42, 42, 43, 44, 46, 46

Offset: 1

Gus Wiseman, Jul 23 2022

nonn

For initial terms up to 30 we have a(n) = Log_2 A355537(n).

The sum of the same range is A000096.
The product of the same range is A000142, Heinz number A070826.
For divisors (not just prime factors) we get A002541, also A006218, A077597.
A shifted variation is A013939.
The unshifted version is A022559, product A327486, w/o multiplicity A355537.
The ranges themselves are the rows of A131818 (shifted).
Partial sums of A297155 (shifted).
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A066843 gives partial sums of A000005.
Cf. A000720, A002110, A076610, A305054, A355538, A355731, A355733, A355741, A355744, A355745, A355746, A355747.

Table[Total[(PrimeNu[#]-1)&/@Range[2,n]],{n,1,100}]

a(n) = A013939(n) - n + 1.

A385574 Number of integer partitions of n with the same number of adjacent equal parts as adjacent unequal parts.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 2, 4, 5, 6, 10, 11, 13, 17, 20, 30, 36, 44, 55, 70, 86, 98, 128, 156, 190, 235, 288, 351, 409, 499, 603, 722, 863, 1025, 1227, 1461, 1757, 2061, 2444, 2892, 3406, 3996, 4708, 5497, 6430, 7595, 8835, 10294, 12027, 13971, 16252, 18887, 21878

Offset: 0

Gus Wiseman, Jul 04 2025

nonn

These are also integer partitions of n with the same number of distinct parts as maximal anti-runs of parts.

			The partition (5,3,2,1,1,1,1) has 4 runs ((5),(3),(2),(1,1,1,1)) and 4 anti-runs ((5,3,2,1),(1),(1),(1)) so is counted under a(14).
The a(1) = 1 through a(10) = 10 reversed partitions (A = 10):
  (1)  (2)  (3)  (4)    (5)    (6)    (7)    (8)      (9)      (A)
                 (112)  (113)  (114)  (115)  (116)    (117)    (118)
                        (122)         (133)  (224)    (144)    (226)
                                      (223)  (233)    (225)    (244)
                                             (11123)  (11124)  (334)
                                                      (11223)  (11125)
                                                               (11134)
                                                               (11224)
                                                               (11233)
                                                               (12223)

Christian Sievers, Table of n, a(n) for n = 0..1000

The RHS is counted by A116608, rank statistic A297155.
The LHS is counted by A133121, rank statistic A046660.
For related inequalities see A212165, A212168, A361204.
For subsets instead of partitions see A217615, A385572, A385575.
These partitions are ranked by A385576.
A000041 counts integer partitions, strict A000009.
A007690 counts partitions with no singletons, complement A183558.
A034296 counts flat or gapless partitions, ranks A066311 or A073491.
A034839 counts subsets by number maximal runs, for partitions A384881, strict A116674.
A047993 counts partitions with max part = length (A106529).
A098859 counts Wilf partitions (complement A336866), compositions A242882.
A268193 counts partitions by maximal anti-runs, strict A384905, subsets A384893.
A355394 counts partitions with neighbors, complement A356236.
Cf. A000071, A003114, A008284, A010027, A047966, A210034, A325324, A325325, A356606, A384882, A384885.

Table[Length[Select[Reverse/@IntegerPartitions[n],Length[Union[#]]==Length[Split[#,#2!=#1&]]&]],{n,0,30}]

lista(n)=Vec(polcoef((prod(i=1,n,1+x^i/(t*(1-t*x^i))+O(x*x^n))-1)*t+1,0,t)) \\ Christian Sievers, Jul 18 2025

For a partition p, let s(p) be its sum, e(p) the number of equal adjacent pairs, and d(p) the number of distinct adjacent pairs. Then Sum_{p partition} x^s(p) * t^(e(p)-d(p)) = (Product_{i>=1} (1 + x^i/(t*(1-t*x^i))) - 1) * t + 1, so a(n) is the coefficient of x^n*t^0 of this expression. - Christian Sievers, Jul 18 2025

A385576 Numbers whose prime indices have the same number of distinct elements as maximal anti-runs.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 12, 13, 17, 18, 19, 20, 23, 28, 29, 31, 37, 41, 43, 44, 45, 47, 50, 52, 53, 59, 61, 63, 67, 68, 71, 73, 75, 76, 79, 83, 89, 92, 97, 98, 99, 101, 103, 107, 109, 113, 116, 117, 120, 124, 127, 131, 137, 139, 147, 148, 149, 151, 153, 157, 163

Offset: 1

Gus Wiseman, Jul 04 2025

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.

These are also numbers with the same number of adjacent equal prime indices as adjacent unequal prime indices.

			The prime indices of 2640 are {1,1,1,1,2,3,5}, with 4 distinct parts {1,2,3,5} and 4 maximal anti-runs ((1),(1),(1),(2,3,5)), so 2640 is in the sequence.
The terms together with their prime indices begin:
   1: {}
   2: {1}
   3: {2}
   5: {3}
   7: {4}
  11: {5}
  12: {1,1,2}
  13: {6}
  17: {7}
  18: {1,2,2}
  19: {8}
  20: {1,1,3}
  23: {9}
  28: {1,1,4}
  29: {10}
  31: {11}
  37: {12}
  41: {13}
  43: {14}
  44: {1,1,5}
  45: {2,2,3}
  47: {15}

The LHS is the rank statistic A001221, triangle counted by A116608.
The RHS is the rank statistic A375136, triangle counted by A133121.
These partitions are counted by A385574.
A034296 counts flat or gapless partitions, ranks A066311 or A073491.
A047993 counts partitions with max part = length, ranks A106529.
A356235 counts partitions with a neighborless singleton, ranks A356237.
A384877 gives lengths of maximal anti-runs of binary indices, firsts A384878.
A384893 counts subsets by maximal anti-runs, for partitions A268193, strict A384905.
A385572 counts subsets with the same number of runs as anti-runs, ranks A385575.
Cf. A044813, A046660, A210034, A297155, A356226, A361205, A384889, A385213.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],#==1||PrimeNu[#]==Length[Split[prix[#],UnsameQ]]&]

A001221(a(n)) = A375136(a(n)).

cp's OEIS Frontend

A156552 Unary-encoded compressed factorization of natural numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Perl

Python

Formula

Extensions

A252464 a(1) = 0, a(2n) = 1 + a(n), a(2n+1) = 1 + a(A064989(2n+1)); also binary width of terms of A156552 and A243071.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A297113 a(1) = 0, a(2) = 1, after which, a(n) = a(n/2) if n is of the form 4k+2, and otherwise a(n) = 1+a(A252463(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Scheme

Formula

A386577 Irregular triangle read by rows where T(n,k) is the number of permutations of the multiset of prime factors of n with k adjacent equal terms.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A355538 Partial sum of A001221 (number of distinct prime factors) minus 1, ranging from 2 to n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula

A385574 Number of integer partitions of n with the same number of adjacent equal parts as adjacent unequal parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A385576 Numbers whose prime indices have the same number of distinct elements as maximal anti-runs.

Views

Author

Keywords

Comments

Examples

Crossrefs