cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 61-70 of 255 results. Next

A335451 Number of permutations of the prime indices of n with all equal parts contiguous and none appearing more than twice.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 0, 1, 2, 1, 2, 1, 2, 2, 0, 1, 2, 1, 2, 2, 2, 1, 0, 1, 2, 0, 2, 1, 6, 1, 0, 2, 2, 2, 2, 1, 2, 2, 0, 1, 6, 1, 2, 2, 2, 1, 0, 1, 2, 2, 2, 1, 0, 2, 0, 2, 2, 1, 6, 1, 2, 2, 0, 2, 6, 1, 2, 2, 6, 1, 0, 1, 2, 2, 2, 2, 6, 1, 0, 0, 2, 1, 6, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jun 21 2020

Keywords

Comments

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.

Examples

			The a(90) = 6 permutations are (1,2,2,3), (1,3,2,2), (2,2,1,3), (2,2,3,1), (3,1,2,2), (3,2,2,1).
		

Crossrefs

Separations are counted by A003242 and A335452 and ranked by A333489.
Permutations of prime indices are counted by A008480.
Unsorted prime signature is A124010. Sorted prime signature is A118914.
Permutations of prime indices with equal parts contiguous are A333175.
STC-numbers of permutations of prime indices are A333221.
(1,2,1) and (2,1,2)-avoiding permutations of prime indices are A333175.
Numbers whose prime indices are inseparable are A335448.
(1,2,1) or (2,1,2)-matching permutations of prime indices are A335460.
(1,2,1) and (2,1,2)-matching permutations of prime indices are A335462.
Strict permutations of prime indices are counted by A335489.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Permutations[primeMS[n]],!MatchQ[#,{_,x_,,x_,_}]&]],{n,100}]

Formula

a(n) = A001221(n)! if n is cubefree, otherwise 0.

A352490 Nonexcedance set of A122111. Numbers k > A122111(k), where A122111 represents partition conjugation using Heinz numbers.

Original entry on oeis.org

4, 8, 12, 16, 18, 24, 27, 32, 36, 40, 48, 50, 54, 60, 64, 72, 80, 81, 90, 96, 100, 108, 112, 120, 128, 135, 140, 144, 150, 160, 162, 168, 180, 192, 196, 200, 216, 224, 225, 240, 243, 250, 252, 256, 270, 280, 288, 300, 315, 320, 324, 336, 352, 360, 375, 378
Offset: 1

Views

Author

Gus Wiseman, Mar 20 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The sequence lists all Heinz numbers of partitions whose Heinz number is greater than that of their conjugate.

Examples

			The terms together with their prime indices begin:
    4: (1,1)
    8: (1,1,1)
   12: (2,1,1)
   16: (1,1,1,1)
   18: (2,2,1)
   24: (2,1,1,1)
   27: (2,2,2)
   32: (1,1,1,1,1)
   36: (2,2,1,1)
   40: (3,1,1,1)
   48: (2,1,1,1,1)
   50: (3,3,1)
   54: (2,2,2,1)
   60: (3,2,1,1)
   64: (1,1,1,1,1,1)
For example, the partition (4,4,1,1) has Heinz number 196 and its conjugate (4,2,2,2) has Heinz number 189, and 196 > 189, so 196 is in the sequence, and 189 is not.
		

Crossrefs

These partitions are counted by A000701.
The opposite version is A352487, weak A352489.
The weak version is A352488, counted by A046682.
These are the positions of positive terms in A352491.
A000041 counts integer partitions, strict A000009.
A000700 counts self-conjugate partitions, ranked by A088902 (cf. A258116).
A003963 = product of prime indices, conjugate A329382.
A056239 adds up prime indices, row sums of A112798 and A296150.
A122111 = partition conjugation using Heinz numbers, parts A321649/A321650.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A173018 counts permutations by excedances, weak A123125.
A330644 counts non-self-conjugate partitions, ranked by A352486.
A352521 counts compositions by subdiagonals, rank statistic A352514.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],#>Times@@Prime/@conj[primeMS[#]]&]

Formula

a(n) > A122111(a(n)).

A360460 Two times the median of the unordered prime signature of n; a(1) = 1.

Original entry on oeis.org

1, 2, 2, 4, 2, 2, 2, 6, 4, 2, 2, 3, 2, 2, 2, 8, 2, 3, 2, 3, 2, 2, 2, 4, 4, 2, 6, 3, 2, 2, 2, 10, 2, 2, 2, 4, 2, 2, 2, 4, 2, 2, 2, 3, 3, 2, 2, 5, 4, 3, 2, 3, 2, 4, 2, 4, 2, 2, 2, 2, 2, 2, 3, 12, 2, 2, 2, 3, 2, 2, 2, 5, 2, 2, 3, 3, 2, 2, 2, 5, 8, 2, 2, 2, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Feb 14 2023

Keywords

Comments

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). Since the denominator is always 1 or 2, the median can be represented as an integer by multiplying by 2.
A number's unordered prime signature (row n of A118914) is the multiset of positive exponents in its prime factorization.

Examples

			The unordered prime signature of 2520 is {1,1,2,3}, with median 3/2, so a(2520) = 3.
		

Crossrefs

The version for divisors is A063655.
For mean instead of two times median we have A088529/A088530.
Prime signature is A124010, unordered A118914.
The version for prime indices is A360005.
The version for distinct prime indices is A360457.
The version for distinct prime factors is A360458.
The version for prime factors is A360459.
Positions of even terms are A360553.
Positions of odd terms are A360554.
The version for 0-prepended differences is A360555.
A112798 lists prime indices, length A001222, sum A056239.
A304038 lists distinct prime indices.
A325347 counts partitions w/ integer median, complement A307683.
A329976 counts partitions with median multiplicity 1.
A359893 and A359901 count partitions by median, odd-length A359902.

Programs

  • Mathematica
    Table[If[n==1,1,2*Median[Last/@FactorInteger[n]]],{n,100}]

A329382 Product of exponents of prime factors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 4, 2, 1, 3, 1, 2, 4, 4, 1, 6, 1, 3, 4, 2, 1, 4, 8, 2, 9, 3, 1, 6, 1, 5, 4, 2, 8, 8, 1, 2, 4, 4, 1, 6, 1, 3, 9, 2, 1, 5, 16, 12, 4, 3, 1, 12, 8, 4, 4, 2, 1, 8, 1, 2, 9, 6, 8, 6, 1, 3, 4, 12, 1, 10, 1, 2, 18, 3, 16, 6, 1, 5, 16, 2, 1, 8, 8, 2, 4, 4, 1, 12, 16, 3, 4, 2, 8, 6, 1, 24, 9, 16, 1, 6, 1, 4, 18
Offset: 1

Views

Author

Antti Karttunen, Nov 17 2019

Keywords

Comments

Also the product of parts of the conjugate of the integer partition with Heinz number n, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). For example, the partition (3,2) with Heinz number 15 has conjugate (2,2,1) with product a(15) = 4. - Gus Wiseman, Mar 27 2022

Crossrefs

This is the conjugate version of A003963 (product of prime indices).
The solutions to a(n) = A003963(n) are A325040, counted by A325039.
The Heinz number of the conjugate partition is given by A122111.
These are the row products of A321649 and of A321650.
A000700 counts self-conj partitions, ranked by A088902, complement A330644.
A008480 counts permutations of prime indices, conjugate A321648.
A056239 adds up prime indices, row sums of A112798 and of A296150.
A124010 gives prime signature, sorted A118914, sum A001222.
A238744 gives the conjugate of prime signature, rank A238745.

Programs

  • Mathematica
    Table[Times @@ FactorInteger[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]][[All, -1]], {n, 105}] (* Michael De Vlieger, Jan 21 2020 *)
  • PARI
    A005361(n) = factorback(factor(n)[, 2]); \\ from A005361
    A034386(n) = prod(i=1, primepi(n), prime(i));
    A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
    A329382(n) = A005361(A108951(n));
    
  • PARI
    A329382(n) = if(1==n,1,my(f=factor(n),e=0,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m)); \\ Antti Karttunen, Jan 14 2020

Formula

a(n) = A005361(A108951(n)).
A329605(n) >= a(n) >= A329617(n) >= A329378(n).
a(A019565(n)) = A284001(n).
From Antti Karttunen, Jan 14 2020: (Start)
If n = p(k1)^e(k1) * p(k2)^e(k2) * p(k3)^e(k3) * ... * p(kx)^e(kx), with p(n) = A000040(n) and k1 > k2 > k3 > ... > kx, then a(n) = e(k1)^(k1-k2) * (e(k1)+e(k2))^(k2-k3) * (e(k1)+e(k2)+e(k3))^(k3-k4) * ... * (e(k1)+e(k2)+...+e(kx))^kx.
a(n) = A000005(A331188(n)) = A329605(A052126(n)).
(End)
a(n) = A003963(A122111(n)). - Gus Wiseman, Mar 27 2022

A335446 Number of (1,2,1)-matching permutations of the prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 2, 0, 0, 0, 1, 1, 0, 0, 3, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 6, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 7, 0, 0, 0, 1, 0, 0, 0, 3, 0, 0, 0, 6, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jun 13 2020

Keywords

Comments

Depends only on unsorted prime signature (A124010), but not only on sorted prime signature (A118914).
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 a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

Examples

			The a(n) permutations for n = 12, 24, 36, 60, 72, 90, 120, 144:
  (121)  (1121)  (1212)  (1213)  (11212)  (1232)  (11213)  (111212)
         (1211)  (1221)  (1231)  (11221)  (2132)  (11231)  (111221)
                 (2121)  (1312)  (12112)  (2312)  (11312)  (112112)
                         (1321)  (12121)  (2321)  (11321)  (112121)
                         (2131)  (12211)          (12113)  (112211)
                         (3121)  (21121)          (12131)  (121112)
                                 (21211)          (12311)  (121121)
                                                  (13112)  (121211)
                                                  (13121)  (122111)
                                                  (13211)  (211121)
                                                  (21131)  (211211)
                                                  (21311)  (212111)
                                                  (31121)
                                                  (31211)
		

Crossrefs

Positions of zeros are A065200.
The avoiding version is A335449.
Patterns are counted by A000670.
Permutations of prime indices are counted by A008480.
Unimodal permutations of prime indices are counted by A332288.
(1,2,1) and (2,1,2)-avoiding permutations of prime indices are A333175.
STC-numbers of permutations of prime indices are A333221.
Patterns matched by standard compositions are counted by A335454.
(1,2,1) or (2,1,2)-matching permutations of prime indices are A335460.
(1,2,1) and (2,1,2)-matching permutations of prime indices are A335462.
Dimensions of downsets of standard compositions are A335465.
(1,2,1)-matching compositions are ranked by A335466.
(1,2,1)-matching compositions are counted by A335470.
(1,2,1)-matching patterns are counted by A335509.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Permutations[primeMS[n]],MatchQ[#,{_,x_,_,y_,_,x_,_}/;x
    				

A336417 Number of perfect-power divisors of superprimorials A006939.

Original entry on oeis.org

1, 1, 2, 5, 15, 44, 169, 652, 3106, 15286, 89933, 532476, 3698650, 25749335, 204947216, 1636097441, 14693641859, 132055603656, 1319433514898, 13186485900967, 144978145009105, 1594375302986404, 19128405558986057, 229508085926717076, 2983342885319348522
Offset: 0

Views

Author

Gus Wiseman, Jul 24 2020

Keywords

Comments

A number is a perfect power iff it is 1 or its prime exponents (signature) are not relatively prime.
The n-th superprimorial number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(0) = 1 through a(4) = 15 divisors:
  1  2  12  360  75600
-------------------------
  1  1   1    1      1
         4    4      4
              8      8
              9      9
             36     16
                    25
                    27
                    36
                   100
                   144
                   216
                   225
                   400
                   900
                  3600
		

Crossrefs

A000325 is the uniform version.
A076954 can be used instead of A006939.
A336416 gives the same for factorials instead of superprimorials.
A000217 counts prime power divisors of superprimorials.
A000961 gives prime powers.
A001597 gives perfect powers, with complement A007916.
A006939 gives superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A091050 counts perfect power divisors.
A181818 gives products of superprimorials.
A294068 counts factorizations using perfect powers.
A317829 counts factorizations of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    perpouQ[n_]:=Or[n==1,GCD@@FactorInteger[n][[All,2]]>1];
    Table[Length[Select[Divisors[chern[n]],perpouQ]],{n,0,5}]
  • PARI
    a(n) = {1 + sum(k=2, n, moebius(k)*(1 - prod(i=1, n, 1 + i\k)))} \\ Andrew Howroyd, Aug 30 2020

Formula

a(n) = A091050(A006939(n)).
a(n) = 1 + Sum_{k=2..n} mu(k)*(1 - Product_{i=1..n} 1 + floor(i/k)). - Andrew Howroyd, Aug 30 2020

Extensions

Terms a(10) and beyond from Andrew Howroyd, Aug 30 2020

A336419 Number of divisors d of the n-th superprimorial A006939(n) with distinct prime exponents such that the quotient A006939(n)/d also has distinct prime exponents.

Original entry on oeis.org

1, 2, 4, 10, 24, 64, 184, 536, 1608, 5104, 16448, 55136, 187136, 658624, 2339648, 8618208, 31884640, 121733120, 468209408, 1849540416, 7342849216
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

A number has distinct prime exponents iff its prime signature is strict.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(0) = 1 through a(3) = 10 divisors:
  1  2  12  360
-----------------
  1  1   1    1
     2   3    5
         4    8
        12    9
             18
             20
             40
             45
             72
            360
		

Crossrefs

A000110 shifted once to the left dominates this sequence.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    Table[Length[Select[Divisors[chern[n]],UnsameQ@@Last/@FactorInteger[#]&&UnsameQ@@Last/@FactorInteger[chern[n]/#]&]],{n,0,6}]
  • PARI
    recurse(n,k,b,d)={if(k>n, 1, sum(i=0, k, if((i==0||!bittest(b,i)) && (i==k||!bittest(d,k-i)), self()(n, k+1, bitor(b, 1<Andrew Howroyd, Aug 30 2020

Extensions

a(10)-a(20) from Andrew Howroyd, Aug 31 2020

A352487 Excedance set of A122111. Numbers k < A122111(k), where A122111 represents partition conjugation using Heinz numbers.

Original entry on oeis.org

3, 5, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93, 94
Offset: 1

Views

Author

Gus Wiseman, Mar 19 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The sequence lists all Heinz numbers of partitions whose Heinz number is less than that of their conjugate.

Examples

			The terms together with their prime indices begin:
   3: (2)
   5: (3)
   7: (4)
  10: (3,1)
  11: (5)
  13: (6)
  14: (4,1)
  15: (3,2)
  17: (7)
  19: (8)
  21: (4,2)
  22: (5,1)
  23: (9)
  25: (3,3)
  26: (6,1)
  28: (4,1,1)
For example, the partition (4,1,1) has Heinz number 28 and its conjugate (3,1,1,1) has Heinz number 40, and 28 < 40, so 28 is in the sequence, and 40 is not.
		

Crossrefs

These partitions are counted by A000701.
The weak version is A352489, counted by A046682.
The opposite version is A352490, weak A352488.
These are the positions of negative terms in A352491.
A000041 counts integer partitions, strict A000009.
A000700 counts self-conjugate partitions, ranked by A088902 (cf. A258116).
A003963 = product of prime indices, conjugate A329382.
A008292 is the triangle of Eulerian numbers (version without zeros).
A008480 counts permutations of prime indices, conjugate A321648.
A056239 adds up prime indices, row sums of A112798 and A296150.
A122111 = partition conjugation using Heinz numbers, parts A321649/A321650.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A173018 counts permutations by excedances, weak A123125.
A238744 = partition conjugate of prime signature, ranked by A238745.
A330644 counts non-self-conjugate partitions, ranked by A352486.
A352521 counts compositions by subdiagonals, rank statistic A352514.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],#
    				

Formula

a(n) < A122111(a(n)).

A353507 Product of multiplicities of the prime exponents (signature) of n; a(1) = 0.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 1, 2, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, May 19 2022

Keywords

Comments

Warning: If the prime multiplicities of n are a multiset y, this sequence gives the product of multiplicities in y, not the product of y.
Differs from A351946 at A351946(1260) = 4, a(1260) = 2.
Differs from A327500 at A327500(450) = 3, a(450) = 2.
We set a(1) = 0 so that the positions of first appearances are the primorials A002110.
Also the product of the prime metasignature of n (row n of A238747).

Examples

			The prime signature of 13860 is (2,2,1,1,1), with multiplicities (2,3), so a(13860) = 6.
		

Crossrefs

Positions of first appearances are A002110.
The prime indices themselves have product A003963, counted by A339095.
The prime signature itself has product A005361, counted by A266477.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A071625 counts distinct prime exponents (third omega).
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A238747 gives prime metasignature, sorted A353742.
A323022 gives fourth omega.

Programs

  • Maple
    f:= proc(n) local M,s;
      M:= ifactors(n)[2][..,2];
      mul(numboccur(s,M),s=convert(M,set));
    end proc:
    f(1):= 0:
    map(f, [$1..100]); # Robert Israel, May 19 2023
  • Mathematica
    Table[If[n==1,0,Times@@Length/@Split[Sort[Last/@FactorInteger[n]]]],{n,100}]
    Join[{0},Table[Times@@(Length/@Split[FactorInteger[n][[;;,2]]]),{n,2,100}]] (* Harvey P. Dale, Oct 20 2024 *)
  • Python
    from math import prod
    from itertools import groupby
    from sympy import factorint
    def A353507(n): return 0 if n == 1 else prod(len(list(g)) for k, g in groupby(factorint(n).values())) # Chai Wah Wu, May 20 2022

Formula

A357864 Numbers whose prime indices have strictly decreasing run-sums. Heinz numbers of the partitions counted by A304430.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 24, 25, 27, 29, 31, 32, 37, 41, 43, 45, 47, 48, 49, 53, 59, 61, 64, 67, 71, 73, 79, 80, 81, 83, 89, 96, 97, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 135, 137, 139, 149, 151, 157, 160, 163, 167, 169, 173
Offset: 1

Views

Author

Gus Wiseman, Oct 19 2022

Keywords

Comments

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.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
The sequence of runs of a sequence consists of its maximal consecutive constant subsequences when read left-to-right. For example, the runs of (2,2,1,1,1,3,2,2) are (2,2), (1,1,1), (3), (2,2), with sums (4,3,3,4).

Examples

			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}
   13: {6}
   16: {1,1,1,1}
   17: {7}
   19: {8}
   23: {9}
   24: {1,1,1,2}
   25: {3,3}
   27: {2,2,2}
   29: {10}
For example, the prime indices of 24 are {1,1,1,2}, with run-sums (3,2), which are strictly decreasing, so 24 is in the sequence.
		

Crossrefs

Subsequence of A304686.
These partitions are counted by A304430.
These are the indices of rows in A354584 that are strictly decreasing.
The weakly decreasing version is A357861, counted by A304406.
The opposite version is A357862, counted by A304428, complement A357863.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[300],Greater@@Total/@Split[primeMS[#]]&]
Previous Showing 61-70 of 255 results. Next