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 11-20 of 51 results. Next

A304678 Numbers with weakly increasing prime multiplicities.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 81, 82, 83
Offset: 1

Views

Author

Gus Wiseman, May 16 2018

Keywords

Comments

Complement of A112769.

Examples

			12 = 2*2*3 has prime multiplicities (2,1) so is not in the sequence.
36 = 2*2*3*3 has prime multiplicities (2,2) so is in the sequence.
150 = 2*3*5*5 has prime multiplicities (1,1,2) so is in the sequence.
		

Crossrefs

Programs

  • Maple
    q:= n-> (l-> (t-> andmap(i-> l[i, 2]<=l[i+1, 2],
            [$1..t-1]))(nops(l)))(sort(ifactors(n)[2])):
    select(q, [$1..120])[];  # Alois P. Heinz, Nov 11 2019
  • Mathematica
    Select[Range[200],OrderedQ[FactorInteger[#][[All,2]]]&]
    Select[Range[90],Min[Differences[FactorInteger[#][[;;,2]]]]>=0&] (* Harvey P. Dale, Jan 28 2024 *)
  • PARI
    isok(n) = my(vm = factor(n)[,2]); vm == vecsort(vm); \\ Michel Marcus, May 17 2018

A325238 First positive integer with each omega-sequence.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 96, 120, 128, 192, 210, 216, 240, 256, 360, 384, 420, 480, 512, 720, 768, 840, 900, 960, 1024, 1260, 1296, 1440, 1536, 1680, 1920, 2048, 2310, 2520, 2880, 3072, 3360, 3840, 4096, 4620, 5040, 5760, 6144, 6720
Offset: 1

Views

Author

Gus Wiseman, Apr 14 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = frequency depth of n, and the k-th part is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, given by red(n = p^i*...*q^j) = prime(i)*...*prime(j), i.e., the product of primes indexed by the prime exponents of n.

Examples

			The sequence of terms together with their omega-sequences begins:
    1:
    2: 1
    4: 2 1
    6: 2 2 1
    8: 3 1
   12: 3 2 2 1
   16: 4 1
   24: 4 2 2 1
   30: 3 3 1
   32: 5 1
   36: 4 2 1
   48: 5 2 2 1
   60: 4 3 2 2 1
   64: 6 1
   96: 6 2 2 1
  120: 5 3 2 2 1
  128: 7 1
  192: 7 2 2 1
  210: 4 4 1
  216: 6 2 1
  240: 6 3 2 2 1
  256: 8 1
  360: 6 3 3 1
  384: 8 2 2 1
  420: 5 4 2 2 1
		

Crossrefs

Programs

  • Mathematica
    tomseq[n_]:=If[n<=1,{},Most[FixedPointList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]]]]];
    omseqs=Table[Total/@tomseq[n],{n,1000}];
    Sort[Table[Position[omseqs,x][[1,1]],{x,Union[omseqs]}]]

A325273 Prime omicron of n!.

Original entry on oeis.org

0, 0, 1, 2, 2, 2, 3, 2, 2, 2, 4, 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, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 2, 2, 2, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).
The prime omicron of n (A304465) is 0 if n is 1, 1 if n is prime, and otherwise the second-to-last part of the omega-sequence of n. For example, the prime omicron of 180 is 2.
Conjecture: all terms after a(10) = 4 are less than 4.
From James Rayman, Apr 17 2021: (Start)
The conjecture is false. a(3804) = 4. In fact, there are 91 values of n < 10000 such that a(n) = 4.
The first value of n such that a(n) = 5 is 37934. For any other n < 5*10^5, a(n) < 5. (End)

Crossrefs

a(n) = A055396(A325275(n)/2).
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).

Programs

  • Mathematica
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    omicron[n_]:=Switch[n,1,0,?PrimeQ,1,,omseq[n][[-2]]];
    Table[omicron[n!],{n,0,100}]
  • Python
    from sympy.ntheory import *
    def red(v):
        r = {}
        for i in v: r[i] = r.get(i, 0) + 1
        return r
    def omicron(v):
        if len(v) == 0: return 0
        if len(v) == 1: return v[0]
        else: return omicron(list(red(v).values()))
    f, a_list = {}, []
    for i in range(101):
        a_list.append(omicron(list(f.values())))
        g = factorint(i+1)
        for k in g: f[k] = f.get(k, 0) + g[k]
    print(a_list) # James Rayman, Apr 17 2021

Extensions

More terms from James Rayman, Apr 17 2021

A353859 Triangle read by rows where T(n,k) is the number of integer compositions of n with composition run-sum trajectory of length k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 3, 1, 0, 0, 4, 2, 2, 0, 0, 7, 7, 2, 0, 0, 0, 14, 14, 4, 0, 0, 0, 0, 23, 29, 12, 0, 0, 0, 0, 0, 39, 56, 25, 8, 0, 0, 0, 0, 0, 71, 122, 53, 10, 0, 0, 0, 0, 0, 0, 124, 246, 126, 16, 0, 0, 0, 0, 0, 0, 0, 214, 498, 264, 48, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Jun 02 2022

Keywords

Comments

Every sequence can be uniquely split into a sequence of non-overlapping runs. 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). The run-sum trajectory is obtained by repeatedly taking the run-sums transformation (or condensation, represented by A353847) until an anti-run is reached. For example, the trajectory (2,4,2,1,1) -> (2,4,2,2) -> (2,4,4) -> (2,8) is counted under T(10,4).

Examples

			Triangle begins:
   1
   0   1
   0   1   1
   0   3   1   0
   0   4   2   2   0
   0   7   7   2   0   0
   0  14  14   4   0   0   0
   0  23  29  12   0   0   0   0
   0  39  56  25   8   0   0   0   0
   0  71 122  53  10   0   0   0   0   0
   0 124 246 126  16   0   0   0   0   0   0
   0 214 498 264  48   0   0   0   0   0   0   0
For example, row n = 5 counts the following compositions:
  (5)    (113)    (1121)
  (14)   (122)    (1211)
  (23)   (221)
  (32)   (311)
  (41)   (1112)
  (131)  (2111)
  (212)  (11111)
		

Crossrefs

Column k = 1 is A003242, ranked by A333489, complement A261983.
Row sums are A011782.
Positive row-lengths are A070939.
The version for partitions is A353846, ranked by A353841.
This statistic (trajectory length) is ranked by A353854, firsts A072639.
Counting by length of last part instead of number of parts gives A353856.
A333627 ranks the run-lengths of standard compositions.
A353847 represents the run-sums of a composition, partitions A353832.
A353853-A353859 pertain to composition run-sum trajectory.
A353932 lists run-sums of standard compositions.

Programs

  • Mathematica
    rsc[y_]:=If[y=={},{},NestWhileList[Total/@Split[#]&,y,MatchQ[#,{_,x_,x_,_}]&]];
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Length[rsc[#]]==k&]],{n,0,10},{k,0,n}]

A325268 Triangle read by rows where T(n,k) is the number of integer partitions of n with omicron k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 3, 0, 1, 0, 1, 5, 0, 0, 1, 0, 1, 7, 2, 0, 0, 1, 0, 1, 12, 1, 0, 0, 0, 1, 0, 1, 17, 2, 1, 0, 0, 0, 1, 0, 1, 24, 4, 0, 0, 0, 0, 0, 1, 0, 1, 33, 5, 1, 1, 0, 0, 0, 0, 1, 0, 1, 44, 9, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 57, 14, 3, 0, 1
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

The omega-sequence of an integer partition is the sequence of lengths of the multisets obtained by repeatedly taking the multiset of multiplicities until a singleton is reached. The omicron of the partition is 0 if the omega-sequence is empty, 1 if it is a singleton, and otherwise the second-to-last part. For example, the partition (32211) has chain of multisets of multiplicities {1,1,2,2,3} -> {1,2,2} -> {1,2} -> {1,1} -> {2}, so its omega-sequence is (5,3,2,2,1), and its omicron is 2.

Examples

			Triangle begins:
  1
  0  1
  0  1  1
  0  1  1  1
  0  1  3  0  1
  0  1  5  0  0  1
  0  1  7  2  0  0  1
  0  1 12  1  0  0  0  1
  0  1 17  2  1  0  0  0  1
  0  1 24  4  0  0  0  0  0  1
  0  1 33  5  1  1  0  0  0  0  1
  0  1 44  9  1  0  0  0  0  0  0  1
  0  1 57 14  3  0  1  0  0  0  0  0  1
  0  1 76 20  3  0  0  0  0  0  0  0  0  1
Row n = 8 counts the following partitions.
  (8)  (44)       (431)  (2222)  (11111111)
       (53)       (521)
       (62)
       (71)
       (332)
       (422)
       (611)
       (3221)
       (3311)
       (4211)
       (5111)
       (22211)
       (32111)
       (41111)
       (221111)
       (311111)
       (2111111)
		

Crossrefs

Row sums are A000041. Column k = 2 is A325267.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or A325280 (length/frequency depth).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Switch[#,{},0,{},1,,NestWhile[Sort[Length/@Split[#]]&,#,Length[#]>1&]//First]==k&]],{n,0,10},{k,0,n}]
  • PARI
    omicron(p)={if(!#p, 0, my(r=1); while(#p > 1, my(L=List(), k=0); r=#p; for(i=1, #p, if(i==#p||p[i]<>p[i+1], listput(L,i-k); k=i)); listsort(L); p=L); r)}
    row(n)={my(v=vector(1+n)); forpart(p=n, v[1 + omicron(Vec(p))]++); v}
    { for(n=0, 10, print(row(n))) } \\ Andrew Howroyd, Jan 18 2023

A305563 Number of reducible integer partitions of n.

Original entry on oeis.org

1, 2, 3, 4, 7, 7, 15, 16, 27, 30, 56, 56, 100, 105, 157, 188, 287, 303, 470, 524, 724, 850, 1197, 1339, 1856, 2135, 2814, 3305, 4360, 4951, 6532, 7561, 9563, 11195, 14165, 16328, 20631, 23866, 29471, 34320, 42336, 48672, 59872, 69139, 83625, 96911, 117153
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Comments

A multiset m whose distinct elements are m_1, m_2, ..., m_k with multiplicities y_1, y_2, ..., y_k is reducible if either m is of size 1 or gcd(m_1, ..., m_k) = 1 and the multiset {y_1, ..., y_k} is also reducible.

Examples

			The a(6) = 7 reducible integer partitions are (6), (51), (411), (321), (3111), (21111), (111111). Missing from this list are (42), (33), (222), (2211).
		

Crossrefs

Programs

  • Mathematica
    ptnredQ[y_]:=Or[Length[y]==1,And[GCD@@y==1,ptnredQ[Sort[Length/@Split[y],Greater]]]];
    Table[Length[Select[IntegerPartitions[n],ptnredQ]],{n,20}]

A325249 Sum of the omega-sequence of n.

Original entry on oeis.org

0, 1, 1, 3, 1, 5, 1, 4, 3, 5, 1, 8, 1, 5, 5, 5, 1, 8, 1, 8, 5, 5, 1, 9, 3, 5, 4, 8, 1, 7, 1, 6, 5, 5, 5, 7, 1, 5, 5, 9, 1, 7, 1, 8, 8, 5, 1, 10, 3, 8, 5, 8, 1, 9, 5, 9, 5, 5, 1, 12, 1, 5, 8, 7, 5, 7, 1, 8, 5, 7, 1, 10, 1, 5, 8, 8, 5, 7, 1, 10, 5, 5, 1, 12, 5
Offset: 1

Views

Author

Gus Wiseman, Apr 16 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).

Examples

			The omega-sequence of 180 is (5,3,2,2,1) with sum 13, so a(180) = 13.
		

Crossrefs

Positions of m's are A000040 (m = 1), A001248 (m = 3), A030078 (m = 4), A068993 (m = 5), A050997 (m = 6), A325264 (m = 7).
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number).

Programs

  • Mathematica
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    Table[Total[omseq[n]],{n,100}]

Formula

a(n) = A056239(A325248(n)).
a(n!) = A325274(n).

A325248 Heinz number of the omega-sequence of n.

Original entry on oeis.org

1, 2, 2, 6, 2, 18, 2, 10, 6, 18, 2, 90, 2, 18, 18, 14, 2, 90, 2, 90, 18, 18, 2, 126, 6, 18, 10, 90, 2, 50, 2, 22, 18, 18, 18, 42, 2, 18, 18, 126, 2, 50, 2, 90, 90, 18, 2, 198, 6, 90, 18, 90, 2, 126, 18, 126, 18, 18, 2, 630, 2, 18, 90, 26, 18, 50, 2, 90, 18, 50
Offset: 1

Views

Author

Gus Wiseman, Apr 16 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The omega-sequence of 180 is (5,3,2,2,1) with Heinz number 990, so a(180) = 990.
		

Crossrefs

Positions of squarefree terms are A325247.
Positions of normal numbers (A055932) are A325251.
First positions of each distinct term are A325238.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number).

Programs

  • Mathematica
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    Table[Times@@Prime/@omseq[n],{n,100}]

Formula

A001222(a(n)) = A323014(n).
A061395(a(n)) = A001222(n).
A304465(n) = A055396(a(n)/2).
A325249(n) = A056239(a(n)).
a(n!) = A325275(n).

A317246 Heinz numbers of supernormal integer partitions.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 18, 30, 32, 60, 64, 90, 128, 150, 180, 210, 256, 300, 360, 450, 512, 540, 600, 1024, 1350, 1500, 2048, 2250, 2310, 2520, 3780, 4096, 4200, 5880, 8192, 9450, 10500, 12600, 13230, 15750, 16384, 17640, 18900, 20580, 26460, 29400, 30030
Offset: 1

Views

Author

Gus Wiseman, Jul 24 2018

Keywords

Comments

An integer partition is supernormal if either (1) it is of the form 1^n for some n >= 0, or (2a) it spans an initial interval of positive integers, and (2b) its multiplicities, sorted in weakly decreasing order, are themselves a supernormal integer partition.

Examples

			Sequence of supernormal integer partitions begins: (), (1), (11), (21), (111), (211), (1111), (221), (321), (11111), (3211), (111111), (3221), (1111111), (3321), (32211), (4321).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    supnrm[q_]:=Or[q=={}||Union[q]=={1},And[Union[q]==Range[Max[q]],supnrm[Sort[Length/@Split[q],Greater]]]];
    Select[Range[10000],supnrm[primeMS[#]]&]

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

Views

Author

Gus Wiseman, May 13 2018

Keywords

Comments

a(1) = 0 by convention.

Examples

			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.
		

Crossrefs

Programs

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

Formula

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.
Previous Showing 11-20 of 51 results. Next