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

A325239 Irregular triangle read by rows where row 1 is {1} and row n > 1 is the sequence starting with n and repeatedly applying A181819 until 2 is reached.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 15 2019

Keywords

Comments

The function A181819 maps n = p^i*...*q^j to prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n.

Examples

			Triangle begins:
   1              26 4 3 2        51 4 3 2          76 6 4 3 2
   2              27 5 2          52 6 4 3 2        77 4 3 2
   3 2            28 6 4 3 2      53 2              78 8 5 2
   4 3 2          29 2            54 10 4 3 2       79 2
   5 2            30 8 5 2        55 4 3 2          80 14 4 3 2
   6 4 3 2        31 2            56 10 4 3 2       81 7 2
   7 2            32 11 2         57 4 3 2          82 4 3 2
   8 5 2          33 4 3 2        58 4 3 2          83 2
   9 3 2          34 4 3 2        59 2              84 12 6 4 3 2
  10 4 3 2        35 4 3 2        60 12 6 4 3 2     85 4 3 2
  11 2            36 9 3 2        61 2              86 4 3 2
  12 6 4 3 2      37 2            62 4 3 2          87 4 3 2
  13 2            38 4 3 2        63 6 4 3 2        88 10 4 3 2
  14 4 3 2        39 4 3 2        64 13 2           89 2
  15 4 3 2        40 10 4 3 2     65 4 3 2          90 12 6 4 3 2
  16 7 2          41 2            66 8 5 2          91 4 3 2
  17 2            42 8 5 2        67 2              92 6 4 3 2
  18 6 4 3 2      43 2            68 6 4 3 2        93 4 3 2
  19 2            44 6 4 3 2      69 4 3 2          94 4 3 2
  20 6 4 3 2      45 6 4 3 2      70 8 5 2          95 4 3 2
  21 4 3 2        46 4 3 2        71 2              96 22 4 3 2
  22 4 3 2        47 2            72 15 4 3 2       97 2
  23 2            48 14 4 3 2     73 2              98 6 4 3 2
  24 10 4 3 2     49 3 2          74 4 3 2          99 6 4 3 2
  25 3 2          50 6 4 3 2      75 6 4 3 2       100 9 3 2
		

Crossrefs

Row lengths are A182850(n) + 1.
See A353510 for a full square array version of this table.

Programs

  • Mathematica
    red[n_]:=Times@@Prime/@Last/@If[n==1,{},FactorInteger[n]];
    Table[NestWhileList[red,n,#>2&],{n,30}]

Formula

A001222(T(n,k)) = A323023(n,k), n > 2, k <= A182850(n).

A325276 Irregular triangle read by rows where row n is the omega-sequence of n!.

Original entry on oeis.org

1, 2, 2, 1, 4, 2, 2, 1, 5, 3, 2, 2, 1, 7, 3, 3, 1, 8, 4, 3, 2, 2, 1, 11, 4, 3, 2, 2, 1, 13, 4, 3, 2, 2, 1, 15, 4, 4, 1, 16, 5, 4, 2, 2, 1, 19, 5, 4, 2, 2, 1, 20, 6, 4, 2, 2, 1, 22, 6, 4, 2, 1, 24, 6, 5, 2, 2, 1, 28, 6, 5, 2, 2, 1, 29, 7, 5, 2, 2, 1
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).

Examples

			Triangle begins:
  {}
  {}
   1
   2  2  1
   4  2  2  1
   5  3  2  2  1
   7  3  3  1
   8  4  3  2  2  1
  11  4  3  2  2  1
  13  4  3  2  2  1
  15  4  4  1
  16  5  4  2  2  1
  19  5  4  2  2  1
  20  6  4  2  2  1
  22  6  4  2  1
  24  6  5  2  2  1
  28  6  5  2  2  1
  29  7  5  2  2  1
  32  7  5  2  2  1
  33  8  5  2  2  1
  36  8  5  2  2  1
  38  8  5  2  2  1
  40  8  6  2  2  1
  41  9  6  2  2  1
  45  9  6  2  2  1
  47  9  6  2  2  1
  49  9  6  3  2  2  1
  52  9  6  3  2  2  1
  55  9  6  3  2  2  1
  56 10  6  3  2  2  1
  59 10  6  3  2  2  1
		

Crossrefs

Row lengths are A325272. Row sums are A325274. Row n is row A325275(n) of A112798. Second-to-last column is A325273. Column k = 1 is A022559. Column k = 2 is A000720. Column k = 3 is A071626.
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&]];
    Table[omseq[n!],{n,0,30}]

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

A325278 Smallest number with adjusted frequency depth n.

Original entry on oeis.org

1, 2, 4, 6, 12, 60, 2520, 1286485200, 35933692027611398678865941374040400000
Offset: 0

Views

Author

Gus Wiseman, Apr 17 2019

Keywords

Comments

The adjusted frequency depth of a positive integer n is 0 if n = 1, and otherwise it is 1 plus the number of times one must apply A181819 to reach a prime number, where A181819(k = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of k. For example, 180 has adjusted frequency depth 5 because we have: 180 -> 18 -> 6 -> 4 -> 3.
Differs from A182857 in having 2 instead of 3.

Crossrefs

A subsequence of 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
    nn=10000;
    fd[n_]:=Switch[n,1,0,?PrimeQ,1,,1+fd[Times@@Prime/@Last/@FactorInteger[n]]];
    fds=fd/@Range[nn];
    Sort[Table[Position[fds,x][[1,1]],{x,Union[fds]}]]

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 *)

A323024 Numbers with exactly three distinct exponents in their prime factorization, or three distinct parts in their prime signature.

Original entry on oeis.org

360, 504, 540, 600, 720, 756, 792, 936, 1008, 1176, 1188, 1200, 1224, 1350, 1368, 1400, 1404, 1440, 1500, 1584, 1620, 1656, 1836, 1872, 1960, 2016, 2052, 2088, 2160, 2200, 2232, 2250, 2268, 2352, 2400, 2448, 2484, 2520, 2600, 2646, 2664, 2736, 2800, 2880, 2904
Offset: 1

Views

Author

Gus Wiseman, Jan 02 2019

Keywords

Comments

Positions of 3's in A071625.
Numbers k such that A001221(A181819(k)) = 3.
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} r(n)/((n-1)*psi(n)) = 0.030575..., where psi is the Dedekind psi function (A001615), and r(n) = Sum_{d|n, 1Amiram Eldar, Oct 18 2020

Examples

			1500 = 2^2 * 3^1 * 5^3 has three distinct exponents {1, 2, 3}, so belongs to the sequence.
52500 = 2^2 * 3^1 * 5^4 * 7^1 has three distinct exponents {1, 2, 4}, so belongs to the sequence.
		

Crossrefs

Programs

  • Mathematica
    tom[n_]:=Length[Union[Last/@If[n==1,{},FactorInteger[n]]]];
    Select[Range[1000],tom[#]==3&]
  • PARI
    is(n) = #Set(factor(n)[, 2]) == 3 \\ David A. Corneth, Jan 02 2019

A323025 Numbers with exactly four distinct exponents in their prime factorization, or four distinct parts in their prime signature.

Original entry on oeis.org

75600, 105840, 113400, 118800, 126000, 140400, 151200, 158760, 178200, 183600, 198000, 205200, 210600, 211680, 232848, 234000, 237600, 246960, 248400, 252000, 261360, 275184, 275400, 280800, 283500, 294000, 302400, 306000, 307800, 313200, 315000, 334800
Offset: 1

Views

Author

Gus Wiseman, Jan 02 2019

Keywords

Comments

Positions of 4's in A071625.
Numbers k such that A001221(A181819(k)) = 4.
Is a(n) ~ c * n for some c? - David A. Corneth, Jan 09 2019
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} r(n)/((n-1)*psi(n)) = 0.00035750... (corresponding to c = 2797.1... in the question above, whose answer is affirmative), where psi is the Dedekind psi function (A001615), and r(n) = Sum_{d_1|n, 1Amiram Eldar, Oct 18 2020

Examples

			126000 = 2^4 * 3^2 * 5^3 * 7^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
831600 = 2^4 * 3^3 * 5^2 * 7^1 * 11^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
		

Crossrefs

Programs

  • Mathematica
    tom[n_]:=Length[Union[Last/@If[n==1,{},FactorInteger[n]]]];
    Select[Range[100000],tom[#]==4&]
  • PARI
    is(n) = #Set(factor(n)[, 2]) == 4 \\ David A. Corneth, Jan 09 2019

A323055 Numbers with exactly two distinct exponents in their prime factorization, or two distinct parts in their prime signature.

Original entry on oeis.org

12, 18, 20, 24, 28, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 68, 72, 75, 76, 80, 84, 88, 90, 92, 96, 98, 99, 104, 108, 112, 116, 117, 120, 124, 126, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160, 162, 164, 168, 171, 172, 175, 176, 180, 184, 188, 189, 192, 198, 200
Offset: 1

Views

Author

Gus Wiseman, Jan 03 2019

Keywords

Comments

The first term is A006939(2) = 12.
First differs from A059404 in lacking 360, whose prime signature has three distinct parts.
Positions of 2's in A071625.
Numbers k such that A001221(A181819(k)) = 2.
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} 1/((n-1)*psi(n)) = 0.3611398..., where psi is the Dedekind psi function (A001615) (Sanna, 2020). - Amiram Eldar, Oct 18 2020

Examples

			3000 = 2^3 * 3^1 * 5^3 has two distinct exponents {1, 3}, so belongs to the sequence.
		

Crossrefs

One distinct exponent: A062770 or A072774.
Two distinct exponents: this sequence.
Three distinct exponents: A323024.
Four distinct exponents: A323025.
Five distinct exponents: A323056.

Programs

  • Maple
    isA323055 := proc(n)
        local eset;
        eset := {};
        for pf in ifactors(n)[2] do
            eset := eset union {pf[2]} ;
        end do:
        simplify(nops(eset) = 2 ) ;
    end proc:
    for n from 12 to 1000 do
        if isA323055(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Jan 09 2019
  • Mathematica
    Select[Range[100],Length[Union[Last/@FactorInteger[#]]]==2&]

A325250 Number of integer partitions of n whose omega-sequence is strict (no repeated parts).

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 5, 2, 5, 4, 6, 2, 11, 3, 10, 12, 17, 12, 31, 22, 42, 47, 57, 60, 98, 94, 119, 143, 174, 182, 256, 253, 321, 365, 425, 480, 615, 645, 803, 946, 1180, 1341, 1766, 2021, 2607, 3145, 3951, 4727, 6123, 7236, 9136
Offset: 0

Views

Author

Gus Wiseman, Apr 16 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. 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).
The Heinz numbers of these partitions are given by A325247.

Examples

			The a(1) = 1 through a(10) = 6 partitions (A = 10):
  1  2   3    4     5      6       7        8         9          A
     11  111  22    11111  33      1111111  44        333        55
              1111         222              2222      222111     3322
                           2211             3311      111111111  4411
                           111111           11111111             22222
                                                                 1111111111
		

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).
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
    omseq[ptn_List]:=If[ptn=={},{},Length/@NestWhileList[Sort[Length/@Split[#]]&,ptn,Length[#]>1&]];
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@omseq[#]&]],{n,0,30}]

Formula

a(n) + A325262(n) = A000041(n).

A325251 Numbers whose omega-sequence covers an initial interval of positive integers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 82, 83, 84, 85, 86
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 enumeration of these partitions by sum is given by A325260.

Examples

			The sequence of terms together with their omega sequences begins:
   1:              31: 1             63: 3 2 2 1
   2: 1            33: 2 2 1         65: 2 2 1
   3: 1            34: 2 2 1         67: 1
   4: 2 1          35: 2 2 1         68: 3 2 2 1
   5: 1            37: 1             69: 2 2 1
   6: 2 2 1        38: 2 2 1         71: 1
   7: 1            39: 2 2 1         73: 1
   9: 2 1          41: 1             74: 2 2 1
  10: 2 2 1        43: 1             75: 3 2 2 1
  11: 1            44: 3 2 2 1       76: 3 2 2 1
  12: 3 2 2 1      45: 3 2 2 1       77: 2 2 1
  13: 1            46: 2 2 1         79: 1
  14: 2 2 1        47: 1             82: 2 2 1
  15: 2 2 1        49: 2 1           83: 1
  17: 1            50: 3 2 2 1       84: 4 3 2 2 1
  18: 3 2 2 1      51: 2 2 1         85: 2 2 1
  19: 1            52: 3 2 2 1       86: 2 2 1
  20: 3 2 2 1      53: 1             87: 2 2 1
  21: 2 2 1        55: 2 2 1         89: 1
  22: 2 2 1        57: 2 2 1         90: 4 3 2 2 1
  23: 1            58: 2 2 1         91: 2 2 1
  25: 2 1          59: 1             92: 3 2 2 1
  26: 2 2 1        60: 4 3 2 2 1     93: 2 2 1
  28: 3 2 2 1      61: 1             94: 2 2 1
  29: 1            62: 2 2 1         95: 2 2 1
		

Crossrefs

Positions of normal numbers (A055932) in A325248.
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
    normQ[m_]:=Or[m=={},Union[m]==Range[Max[m]]];
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    Select[Range[100],normQ[omseq[#]]&]
Previous Showing 11-20 of 42 results. Next