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 129 results. Next

A353426 Number of integer partitions of n that are empty or a singleton or whose multiplicities are a sub-multiset that is already counted.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 2, 1, 2, 2, 3, 3, 5, 4, 6, 5, 6, 6, 7, 8, 10, 12, 12, 14, 13, 13, 18, 15, 16, 19, 20, 20, 32, 37, 53, 74, 105
Offset: 0

Views

Author

Gus Wiseman, May 16 2022

Keywords

Comments

a(n) is number of integer partitions of n whose Heinz number belongs to A353393, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The a(n) partitions for selected n (A..M = 10..22):
  n=1: n=4:  n=14:     n=16:     n=17:     n=18:        n=22:
------------------------------------------------------------------
  (1)  (4)   (E)       (G)       (H)       (I)          (M)
       (22)  (5522)    (4444)    (652211)  (7722)       (9922)
             (532211)  (6622)    (742211)  (752211)     (972211)
                       (642211)  (832211)  (842211)     (A62211)
                       (732211)            (932211)     (B52211)
                                           (333222111)  (C42211)
                                                        (D32211)
		

Crossrefs

The non-recursive version is A325702, ranked by A325755.
The version for compositions is A353391, non-recursive A353390.
These partitions are ranked by A353393, nonprime A353389.
A047966 counts uniform partitions, compositions A329738.
A239455 counts Look-and-Say partitions, ranked by A351294.

Programs

  • Mathematica
    oosQ[y_]:=Length[y]<=1||MemberQ[Subsets[Sort[y],{Length[Union[y]]}],Sort[Length/@Split[y]]]&&oosQ[Sort[Length/@Split[y]]];
    Table[Length[Select[IntegerPartitions[n],oosQ]],{n,0,30}]

A353503 Numbers whose product of prime indices equals their product of prime exponents (prime signature).

Original entry on oeis.org

1, 2, 12, 36, 40, 112, 352, 832, 960, 1296, 2176, 2880, 4864, 5376, 11776, 12544, 16128, 29696, 33792, 34560, 38400, 63488, 64000, 101376, 115200, 143360, 151552, 159744, 335872, 479232, 704512, 835584, 1540096, 1658880, 1802240
Offset: 1

Views

Author

Gus Wiseman, May 17 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. A number's prime signature (row n A124010) is the sequence of positive exponents in its prime factorization.

Examples

			The terms together with their prime indices begin:
     1: {}
     2: {1}
    12: {1,1,2}
    36: {1,1,2,2}
    40: {1,1,1,3}
   112: {1,1,1,1,4}
   352: {1,1,1,1,1,5}
   832: {1,1,1,1,1,1,6}
   960: {1,1,1,1,1,1,2,3}
  1296: {1,1,1,1,2,2,2,2}
  2176: {1,1,1,1,1,1,1,7}
  2880: {1,1,1,1,1,1,2,2,3}
  4864: {1,1,1,1,1,1,1,1,8}
  5376: {1,1,1,1,1,1,1,1,2,4}
		

Crossrefs

For shadows instead of exponents we get A003586, counted by A008619.
The LHS (product of prime indices) is A003963, counted by A339095.
The RHS (product of prime exponents) is A005361, counted by A266477.
The version for shadows instead of indices is A353399, counted by A353398.
These partitions are counted by A353506.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A124010 gives prime signature, sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A353394 gives product of shadows of prime indices, firsts A353397.

Programs

  • Mathematica
    Select[Range[1000],Times@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]^k]==Times@@Last/@FactorInteger[#]&]
  • Python
    from itertools import count, islice
    from math import prod
    from sympy import primepi, factorint
    def A353503_gen(startvalue=1): # generator of terms >= startvalue
        return filter(lambda n: n == 1 or prod((f:=factorint(n)).values()) == prod(primepi(p)**e for p,e in f.items()), count(max(startvalue,1)))
    A353503_list = list(islice(A353503_gen(),20)) # Chai Wah Wu, May 20 2022

Formula

A003963(a(n)) = A005361(a(n)).

A367583 Greatest element in row n of A367579 (multiset multiplicity kernel).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 28 2023

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.
We define the multiset multiplicity kernel MMK(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 min(S) has multiplicity |S| in MMK(m). For example, MMK({1,1,2,2,3,4,5}) = {1,1,3,3,3}, and MMK({1,2,3,4,5,5,5,5}) = {1,1,1,1,5}.

Examples

			For 450 = 2^1 * 3^2 * 5^2, we have MMK({1,2,2,3,3}) = {1,2,2} so a(450) = 2.
		

Crossrefs

Positions of first appearances are A008578.
Depends only on rootless base A052410, see A007916, A052409.
For minimum instead of maximum element we have A055396.
Row maxima of A367579.
Greatest prime index of A367580.
Positions of 1's are A367586 (powers of even squarefree numbers).
The opposite version is A367587.
A007947 gives squarefree kernel.
A072774 lists powers of squarefree numbers.
A112798 lists prime indices, length A001222, sum A056239, reverse A296150.
A124010 gives prime signature, sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A238747 gives prime metasignature, reverse A353742.
A304038 lists distinct prime indices, length A001221, sum A066328.
A363486 gives least prime index of greatest exponent.
A363487 gives greatest prime index of greatest exponent.
A364191 gives least prime index of least exponent.
A364192 gives greatest prime index of least exponent.

Programs

  • Mathematica
    mmk[q_]:=With[{mts=Length/@Split[q]},Sort[Table[Min@@Select[q,Count[q,#]==i&],{i,mts}]]];
    Table[If[n==1,0,Max@@mmk[PrimePi/@Join@@ConstantArray@@@If[n==1,{},FactorInteger[n]]]],{n,1,100}]

Formula

a(n) = A061395(A367580(n)).
a(n^k) = a(n) for all positive integers n and k.
If n is a power of a squarefree number, a(n) = A055396(n).

A181820 a(1) = 1; for n > 1, if A025487(n) = Product p(i)^e(i), then a(n) = Product p(e(i)).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 10, 8, 11, 9, 14, 12, 13, 15, 22, 20, 17, 21, 18, 26, 16, 25, 28, 19, 33, 30, 34, 24, 35, 44, 23, 39, 42, 38, 40, 55, 27, 52, 29, 50, 51, 36, 49, 66, 46, 56, 65, 45, 68, 31, 70, 57, 32, 60, 77, 78, 58, 88, 85, 63, 76, 37, 110, 69, 48, 84, 91, 75, 102, 62, 54, 98, 104, 95
Offset: 1

Views

Author

Matthew Vandermast, Dec 07 2010

Keywords

Comments

A permutation of the positive integers.
The partition given by the prime signature of A025487(n) has Heinz number a(n). - Pontus von Brömssen, Mar 25 2023

Examples

			A025487(8) = 24 = 2^3*3 has the exponents (3,1) in its canonical prime factorization. Accordingly, a(8) = prime(3)*prime(1) (i.e., A000040(3)*A000040(1)), which equals 5*2=10.
		

Crossrefs

A181815 is another mapping from the members of A025487 to the positive integers. Also see A181819, A181821.
Cf. A000040, A122111, A361808 (inverse), A361809 (fixed points).

Formula

a(n) = A181819(A025487(n)).
a(n) = A122111(A181815(n)).

A323307 Number of ways to fill a matrix with the parts of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 4, 2, 6, 3, 12, 18, 12, 2, 36, 4, 10, 20, 72, 2, 60, 4, 40, 60, 24, 3, 120, 80, 14, 360, 120, 4, 240, 2, 240, 42, 32, 70, 720, 6, 27, 112, 480, 2, 210, 4, 84, 420, 40, 4, 1440, 280, 280, 108, 224, 5, 1260, 224, 420, 180, 22, 2, 840, 6, 72, 1680, 2880
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(22) = 24 matrices:
  [111112] [111121] [111211] [112111] [121111] [211111]
.
  [111] [111] [111] [112] [121] [211]
  [112] [121] [211] [111] [111] [111]
.
  [11] [11] [11] [11] [12] [21]
  [11] [11] [12] [21] [11] [11]
  [12] [21] [11] [11] [11] [11]
.
  [1] [1] [1] [1] [1] [2]
  [1] [1] [1] [1] [2] [1]
  [1] [1] [1] [2] [1] [1]
  [1] [1] [2] [1] [1] [1]
  [1] [2] [1] [1] [1] [1]
  [2] [1] [1] [1] [1] [1]
		

Crossrefs

Programs

  • Mathematica
    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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Array[Length[ptnmats[Times@@Prime/@nrmptn[#]]]&,30]

Formula

a(n) = A318762(n) * A000005(A056239(n)).

A325258 a(1) = 1; otherwise, first differences of Levine's sequence A011784.

Original entry on oeis.org

1, 1, 1, 1, 3, 7, 28, 171, 2624, 172613, 139584150, 6837485347187, 266437138079023501057, 508009471379222384299345337895696, 37745517525533091954228691786161750063795478326636142, 5347426383812697233786139576220412396732847744407175515852823296919414647252347610750
Offset: 0

Views

Author

Gus Wiseman, Apr 16 2019

Keywords

Comments

a(n) is the number of nonnegative integers k such that the maximum adjusted frequency depth among integer partitions of k is n. For example, the a(5) = 7 numbers are 7, 8, 9, 10, 11, 12, and 13.
The adjusted frequency depth of an integer partition is 0 if the partition is empty, and otherwise it is 1 plus the number of times one must take the multiset of multiplicities to reach a singleton. For example, the partition (32211) has adjusted frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2). The enumeration of integer partitions by adjusted frequency depth is given by A325280. The adjusted frequency depth of the integer partition with Heinz number n is A323014(n). The maximum adjusted frequency depth for partitions of n is A325282(n).

Crossrefs

Programs

  • Mathematica
    grw[q_]:=Join@@Table[ConstantArray[i,q[[Length[q]-i+1]]],{i,Length[q]}];
    ReplacePart[Differences[Last/@NestList[grw,{1,1},9]],2->1]

A332671 Number of non-unimodal permutations of the multiset of prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 22 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.
A sequence of integers is unimodal if it is the concatenation of a weakly increasing and a weakly decreasing sequence.

Examples

			The a(n) permutations for n = 18, 30, 36, 42, 50, 54, 60, 66, 70, 72:
  212  213  1212  214  313  2122  1213  215  314  11212
       312  2112  412       2212  1312  512  413  12112
            2121                  2113            12121
                                  2131            21112
                                  3112            21121
                                  3121            21211
		

Crossrefs

Dominated by A008480.
The complement is counted by A332288.
A more interesting version is A332672.
Unimodal compositions are A001523.
Non-unimodal permutations are A059204.
Non-unimodal compositions are A115981.
Non-unimodal normal sequences are A328509.
Heinz numbers of partitions with non-unimodal run-lengths are A332282.
Compositions whose negation is not unimodal are A332669.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    unimodQ[q_]:=Or[Length[q]<=1,If[q[[1]]<=q[[2]],unimodQ[Rest[q]],OrderedQ[Reverse[q]]]];
    Table[Length[Select[Permutations[primeMS[n]],!unimodQ[#]&]],{n,100}]

Formula

a(n) + A332288(n) = A008480(n).
a(A181821(n)) = A332672(n).

A353397 Replace prime(k) with prime(2^k) in the prime factorization of n.

Original entry on oeis.org

1, 3, 7, 9, 19, 21, 53, 27, 49, 57, 131, 63, 311, 159, 133, 81, 719, 147, 1619, 171, 371, 393, 3671, 189, 361, 933, 343, 477, 8161, 399, 17863, 243, 917, 2157, 1007, 441, 38873, 4857, 2177, 513, 84017, 1113, 180503, 1179, 931, 11013, 386093, 567, 2809, 1083
Offset: 1

Views

Author

Gus Wiseman, May 17 2022

Keywords

Examples

			The terms together with their prime indices begin:
      1: {}
      3: {2}
      7: {4}
      9: {2,2}
     19: {8}
     21: {2,4}
     53: {16}
     27: {2,2,2}
     49: {4,4}
     57: {2,8}
    131: {32}
     63: {2,2,4}
		

Crossrefs

These are the positions of first appearances in A353394.
A001222 counts prime factors with multiplicity, distinct A001221.
A003963 gives product of prime indices, counted by A339095.
A033844 lists primes indexed by powers of 2.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914, product A005361.
A181819 gives prime shadow, firsts A181821, relatively prime A325131.
Equivalent sequence with prime(2*k) instead of prime(2^k): A297002.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Times@@Prime/@(2^primeMS[n]),{n,100}]
  • PARI
    a(n) = my(f=factor(n)); for(k=1, #f~, f[k,1] = prime(2^primepi(f[k,1]))); factorback(f); \\ Michel Marcus, May 20 2022
    
  • Python
    from math import prod
    from sympy import prime, primepi, factorint
    def A353397(n): return prod(prime(2**primepi(p))**e for p, e in factorint(n).items()) # Chai Wah Wu, May 20 2022

Formula

If n = prime(e_1)...prime(e_k), then a(n) = prime(2^(e_1))...prime(2^(e_k)).
Sum_{n>=1} 1/a(n) = 1/Product_{k>=1} (1 - 1/prime(2^k)) = 1.90812936178871496289... . - Amiram Eldar, Dec 09 2022

A382773 Number of ways to permute a multiset whose multiplicities are the prime indices of n so that the run-lengths are all different.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 0, 1, 2, 2, 0, 1, 0, 1, 0, 4, 4, 1, 0, 4, 4, 0, 0, 1, 6, 1, 0, 4, 6, 4, 0, 1, 6, 4, 0, 1, 6, 1, 0, 0, 8, 1, 0, 4, 0, 6, 0, 1, 0, 6, 0, 6, 8, 1, 0, 1, 10, 0, 0, 8, 6, 1, 0, 8, 6, 1, 0, 1, 10, 0, 0, 6, 6, 1, 0, 0, 12, 1, 0, 16
Offset: 1

Views

Author

Gus Wiseman, Apr 09 2025

Keywords

Comments

This described multiset (row n of A305936, Heinz number A181821) is generally not the same as the multiset of prime indices of n (A112798). For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(n) partitions for n = 6, 21, 30, 46:
  (1,1,2)  (1,1,1,1,2,2)  (1,1,1,2,2,3)  (1,1,1,1,1,1,1,1,1,2)
  (2,1,1)  (1,1,1,2,2,1)  (1,1,1,3,2,2)  (1,1,1,1,1,1,1,2,1,1)
           (1,2,2,1,1,1)  (2,2,1,1,1,3)  (1,1,1,1,1,1,2,1,1,1)
           (2,2,1,1,1,1)  (2,2,3,1,1,1)  (1,1,1,1,1,2,1,1,1,1)
                          (3,1,1,1,2,2)  (1,1,1,1,2,1,1,1,1,1)
                          (3,2,2,1,1,1)  (1,1,1,2,1,1,1,1,1,1)
                                         (1,1,2,1,1,1,1,1,1,1)
                                         (2,1,1,1,1,1,1,1,1,1)
		

Crossrefs

Positions of 1 are A008578.
For anti-run permutations we have A335125.
For just prime indices we have A382771, firsts A382772, equal A382857.
These permutations for factorials are counted by A382774, equal A335407.
For equal instead of distinct run-lengths we have A382858.
Positions of 0 are A382912, complement A382913.
A044813 lists numbers whose binary expansion has distinct run-lengths, equal A140690.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A098859 counts partitions with distinct multiplicities, ordered A242882.
A239455 counts Look-and-Say partitions, ranks A351294, conjugate A381432.
A329738 counts compositions with equal run-lengths, ranks A353744.
A329739 counts compositions with distinct run-lengths, ranks A351596.
A351293 counts non-Look-and-Say partitions, ranks A351295, conjugate A381433.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Table[Length[Select[Permutations[nrmptn[n]],UnsameQ@@Length/@Split[#]&]],{n,100}]

Formula

a(n) = A382771(A181821(n)) = A382771(A304660(n)).

A182855 Numbers that require exactly five iterations to reach a fixed point under the x -> A181819(x) map.

Original entry on oeis.org

60, 84, 90, 120, 126, 132, 140, 150, 156, 168, 180, 198, 204, 220, 228, 234, 240, 252, 260, 264, 270, 276, 280, 294, 300, 306, 308, 312, 315, 336, 340, 342, 348, 350, 364, 372, 378, 380, 396, 408, 414, 420, 440, 444, 450, 456, 460, 468, 476, 480, 490, 492, 495
Offset: 1

Views

Author

Matthew Vandermast, Jan 04 2011

Keywords

Comments

In each case, 2 is the fixed point that is reached (1 is the other fixed point of the x -> A181819(x) map).
Includes all integers whose prime signature a) contains two or more distinct numbers, and b) contains no number that occurs the same number of times as any other number. The first member of this sequence that does not fit that description is 75675600, whose prime signature is (4,3,2,2,1,1).
A full characterization is: Numbers whose prime signature (1) has not all equal multiplicities but (2) the numbers of distinct parts appearing with each distinct multiplicity are all equal. For example, the prime signature of 2520 is {1,1,2,3}, which satisfies (1) but fails (2), as the numbers of distinct parts appearing with each distinct multiplicity are 1 (with multiplicity 2, the part being 1) and 2 (with multiplicity 1, the parts being 2 and 3). Hence the sequence does not contain 2520. - Gus Wiseman, Jan 02 2019

Examples

			1. 180 requires exactly five iterations under the x -> A181819(x) map to reach a fixed point (namely, 2).  A181819(180) = 18;  A181819(18) = 6; A181819(6) = 4; A181819(4) = 3;  A181819(3) = 2 (and A181819(2) = 2).
2. The prime signature of 180 (2^2*3^2*5) is (2,2,1).
a. Two distinct numbers appear in (2,2,1) (namely, 1 and 2).
b. Neither 1 nor 2 appears in (2,2,1) the same number of times as any other number that appears there.
		

Crossrefs

Numbers n such that A182850(n) = 5. See also A182853, A182854.
Subsequence of A059404 and A182851. Includes A085987 and A179642 as subsequences.

Programs

  • Mathematica
    Select[Range[1000],With[{sig=Sort[Last/@FactorInteger[#]]},And[!SameQ@@Length/@Split[sig],SameQ@@Length/@Union/@GatherBy[sig,Length[Position[sig,#]]&]]]&] (* Gus Wiseman, Jan 02 2019 *)
Previous Showing 61-70 of 129 results. Next