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.

Showing 1-6 of 6 results.

A304647 Smallest term of A304636 that requires exactly n iterations to reach a fixed point under the x -> A181819(x) map.

Original entry on oeis.org

5, 8, 30, 360, 1801800, 2746644314348614680000, 13268350773236509446586539974366689358164301703214270074935844483572035447570761114173070859428708074413696096366645684575600000000
Offset: 0

Views

Author

Gus Wiseman, May 15 2018

Keywords

Comments

The first entry 5 is optional so has offset 0.

Examples

			The list of multisets with Heinz numbers in the sequence is the following. The number of k's in row n + 1 is equal to the k-th largest term of row n.
                     5: {3}
                     8: {1,1,1}
                    30: {1,2,3}
                   360: {1,1,1,2,2,3}
               1801800: {1,1,1,2,2,3,3,4,5,6}
2746644314348614680000: {1,1,1,1,1,1,2,2,2,2,2,3,3,3,3,4,4,4,5,5,5,6,6,7,7,8,9,10}
		

Crossrefs

Programs

  • Mathematica
    Function[m,Times@@Prime/@m]/@NestList[Join@@Table[Table[i,{Reverse[#][[i]]}],{i,Length[#]}]&,{3},6]

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

A307534 Heinz numbers of strict integer partitions with 3 parts, all of which are odd.

Original entry on oeis.org

110, 170, 230, 310, 374, 410, 470, 506, 590, 670, 682, 730, 782, 830, 902, 935, 970, 1030, 1034, 1054, 1090, 1265, 1270, 1298, 1370, 1394, 1426, 1474, 1490, 1570, 1598, 1606, 1670, 1705, 1790, 1826, 1886, 1910, 1955, 1970, 2006, 2110, 2134, 2162, 2255, 2266
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2019

Keywords

Comments

The enumeration of these partitions by sum is given by A001399.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   110: {1,3,5}
   170: {1,3,7}
   230: {1,3,9}
   310: {1,3,11}
   374: {1,5,7}
   410: {1,3,13}
   470: {1,3,15}
   506: {1,5,9}
   590: {1,3,17}
   670: {1,3,19}
   682: {1,5,11}
   730: {1,3,21}
   782: {1,7,9}
   830: {1,3,23}
   902: {1,5,13}
   935: {3,5,7}
   970: {1,3,25}
  1030: {1,3,27}
  1034: {1,5,15}
  1054: {1,7,11}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1000],SquareFreeQ[#]&&PrimeNu[#]==3&&OddQ[Times@@PrimePi/@First/@FactorInteger[#]]&]
  • Python
    from math import isqrt
    from sympy import primepi, primerange, integer_nthroot, nextprime
    def A307534(n):
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x): return int(n+x-sum((primepi(x//(k*m))+1>>1)-(b+1>>1) for a,k in filter(lambda x:x[0]&1,enumerate(primerange(2,integer_nthroot(x,3)[0]+1),1)) for b,m in filter(lambda x:x[0]&1,enumerate(primerange(nextprime(k)+1,isqrt(x//k)+1),a+2))))
        return bisection(f,n,n) # Chai Wah Wu, Oct 20 2024

A304687 Start with the multiset of prime multiplicities of n. Given a multiset, take the multiset of its multiplicities. Repeat until a constant multiset {k,k,...,k} is reached, and set a(n) to the sum of this multiset (k times the length).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 16 2018

Keywords

Examples

			The following are examples showing the reduction of a multiset starting with the multiset of prime multiplicities of n.
         a(60) = 2: {1,1,2} -> {1,2} -> {1,1}.
        a(360) = 3: {1,2,3} -> {1,1,1}.
       a(1260) = 4: {1,1,2,2} -> {2,2}.
a(21492921450) = 6: {1,1,2,2,3,3} -> {2,2,2}.
		

Crossrefs

Programs

  • Maple
    a:= proc(n) map(i-> i[2], ifactors(n)[2]);
          while nops({%[]})>1 do [coeffs(add(x^i, i=%))] od;
          add(i, i=%)
        end:
    seq(a(n), n=1..100);  # Alois P. Heinz, May 17 2018
  • Mathematica
    Table[If[n==1,0,NestWhile[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],!SameQ@@#&]//Total],{n,360}]

A304634 Numbers n with prime omicron 2, meaning A304465(n) = 2.

Original entry on oeis.org

4, 6, 9, 10, 12, 14, 15, 18, 20, 21, 22, 24, 25, 26, 28, 33, 34, 35, 36, 38, 39, 40, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 65, 68, 69, 72, 74, 75, 76, 77, 80, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 98, 99, 100, 104
Offset: 1

Views

Author

Gus Wiseman, May 15 2018

Keywords

Comments

If n > 1 is not a prime number, we have A056239(n) >= Omega(n) >= omega(n) >= A071625(n) >= ... >= omicron(n) > 1 where Omega = A001222, omega = A001221, and omicron = A304465.

Examples

			This is a list of normalized factorizations (see A112798) of selected entries:
    4: {1,1}
    6: {1,2}
   12: {1,1,2}
   24: {1,1,1,2}
   36: {1,1,2,2}
   48: {1,1,1,1,2}
   60: {1,1,2,3}
   72: {1,1,1,2,2}
   96: {1,1,1,1,1,2}
  120: {1,1,1,2,3}
  144: {1,1,1,1,2,2}
  180: {1,1,2,2,3}
  192: {1,1,1,1,1,1,2}
  216: {1,1,1,2,2,2}
  240: {1,1,1,1,2,3}
  288: {1,1,1,1,1,2,2}
  384: {1,1,1,1,1,1,1,2}
  420: {1,1,2,3,4}
  432: {1,1,1,1,2,2,2}
  480: {1,1,1,1,1,2,3}
  576: {1,1,1,1,1,1,2,2}
  768: {1,1,1,1,1,1,1,1,2}
  840: {1,1,1,2,3,4}
  864: {1,1,1,1,1,2,2,2}
  960: {1,1,1,1,1,1,2,3}
		

Crossrefs

Programs

  • Mathematica
    Join@@Position[Table[Switch[n,1,0,?PrimeQ,1,,NestWhile[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Length[#]>1&]//First],{n,100}],2]

A325267 Number of integer partitions of n with omicron 2.

Original entry on oeis.org

0, 0, 1, 1, 3, 5, 7, 12, 17, 24, 33, 44, 57, 76, 100, 129, 168, 214, 282, 355, 462, 586, 755, 937, 1202, 1493, 1900, 2349, 2944, 3621, 4520, 5514, 6813, 8298, 10150, 12240, 14918, 17931, 21654, 25917, 31081, 37029, 44256, 52474, 62405, 73724, 87378, 102887
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

The Heinz numbers of these partitions are given by A304634.
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. We define the omicron of an integer partition to be 0 if the partition is empty, 1 if it is a singleton, and otherwise the second-to-last part of its omega-sequence. 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

			The a(1) = 1 through a(8) = 17 partitions:
  (11)  (21)  (22)   (32)    (33)     (43)      (44)
              (31)   (41)    (42)     (52)      (53)
              (211)  (221)   (51)     (61)      (62)
                     (311)   (411)    (322)     (71)
                     (2111)  (2211)   (331)     (332)
                             (3111)   (511)     (422)
                             (21111)  (2221)    (611)
                                      (3211)    (3221)
                                      (4111)    (3311)
                                      (22111)   (4211)
                                      (31111)   (5111)
                                      (211111)  (22211)
                                                (32111)
                                                (41111)
                                                (221111)
                                                (311111)
                                                (2111111)
		

Crossrefs

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]==2&]],{n,0,30}]
Showing 1-6 of 6 results.