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-10 of 19 results. Next

A320655 Number of factorizations of n into semiprimes. Number of multiset partitions of the multiset of prime factors of n, into pairs.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Comments

The characteristic function of nonzero terms is A065043. - R. J. Mathar, Jan 18 2021

Examples

			The a(900) = 5 factorizations into semiprimes:
  900 = (4*9*25)
  900 = (4*15*15)
  900 = (6*6*25)
  900 = (6*10*15)
  900 = (9*10*10)
The a(900) = 5 multiset partitions into pairs:
  {{1,1},{2,2},{3,3}}
  {{1,1},{2,3},{2,3}}
  {{1,2},{1,2},{3,3}}
  {{1,2},{1,3},{2,3}}
  {{2,2},{1,3},{1,3}}
		

Crossrefs

Programs

  • Mathematica
    semfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[semfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Table[Length[semfacs[n]],{n,100}]
  • PARI
    A320655(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((2==bigomega(d)&&(d<=m)), s += A320655(n/d, d))); (s)); \\ Antti Karttunen, Dec 06 2020

Extensions

Data section extended up to 105 terms by Antti Karttunen, Dec 06 2020

A320656 Number of factorizations of n into squarefree semiprimes. Number of multiset partitions of the multiset of prime factors of n, into strict pairs.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Examples

			The a(4620) = 6 factorizations into squarefree semiprimes:
  4620 = (6*10*77)
  4620 = (6*14*55)
  4620 = (6*22*35)
  4620 = (10*14*33)
  4620 = (10*21*22)
  4620 = (14*15*22)
The a(4620) = 6 multiset partitions into strict pairs:
  {{1,2},{1,3},{4,5}}
  {{1,2},{1,4},{3,5}}
  {{1,2},{1,5},{3,4}}
  {{1,3},{1,4},{2,5}}
  {{1,3},{2,4},{1,5}}
  {{1,4},{2,3},{1,5}}
The a(69300) = 10 factorizations into squarefree semiprimes:
  69300 = (6*6*35*55)
  69300 = (6*10*15*77)
  69300 = (6*10*21*55)
  69300 = (6*10*33*35)
  69300 = (6*14*15*55)
  69300 = (6*15*22*35)
  69300 = (10*10*21*33)
  69300 = (10*14*15*33)
  69300 = (10*15*21*22)
  69300 = (14*15*15*22)
The a(69300) = 10 multiset partitions into strict pairs:
  {{1,2},{1,2},{3,4},{3,5}}
  {{1,2},{1,3},{2,3},{4,5}}
  {{1,2},{1,3},{2,4},{3,5}}
  {{1,2},{1,3},{2,5},{3,4}}
  {{1,2},{1,4},{2,3},{3,5}}
  {{1,2},{2,3},{1,5},{3,4}}
  {{1,3},{1,3},{2,4},{2,5}}
  {{1,3},{1,4},{2,3},{2,5}}
  {{1,3},{2,3},{2,4},{1,5}}
  {{1,4},{2,3},{2,3},{1,5}}.
The a(210) = 3 factorizations into squarefree semiprimes: 210 = (6*35) = (10*21) = (14*15). - _Antti Karttunen_, Nov 02 2022
		

Crossrefs

Programs

  • Mathematica
    bepfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[bepfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Table[Length[bepfacs[n]],{n,100}]
  • PARI
    A320656(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&issquarefree(d)&&2==bigomega(d), s += A320656(n/d, d))); (s)); \\ Antti Karttunen, Nov 02 2022

Formula

a(A002110(n)) = A123023(n). - Antti Karttunen, Nov 02 2022

Extensions

Data section extended up to a(120) and the secondary offset added by Antti Karttunen, Nov 02 2022

A320892 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into distinct semiprimes.

Original entry on oeis.org

16, 64, 81, 96, 144, 160, 224, 256, 324, 352, 384, 400, 416, 486, 544, 576, 608, 625, 640, 729, 736, 784, 864, 896, 928, 960, 992, 1024, 1184, 1215, 1296, 1312, 1344, 1376, 1408, 1440, 1504, 1536, 1600, 1664, 1696, 1701, 1888, 1936, 1944, 1952, 2016, 2025
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

A semiprime (A001358) is a product of any two not necessarily distinct primes.
If A025487(k) is in the sequence then so is every number with the same prime signature. - David A. Corneth, Oct 23 2018
Numbers for which A001222(n) is even and A322353(n) is zero. - Antti Karttunen, Dec 06 2018

Examples

			A complete list of all factorizations of 1296 into semiprimes is:
  1296 = (4*4*9*9)
  1296 = (4*6*6*9)
  1296 = (6*6*6*6)
None of these is strict, so 1296 belongs to the sequence.
		

Crossrefs

Programs

  • Mathematica
    strsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Select[Range[1000],And[EvenQ[PrimeOmega[#]],strsemfacs[#]=={}]&]
  • PARI
    A322353(n, m=n, facs=List([])) = if(1==n, my(u=apply(bigomega,Vec(facs))); (0==length(u)||(2==vecmin(u)&&2==vecmax(u))), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A322353(n/d, d-1, newfacs))); (s));
    isA300892(n) = if(bigomega(n)%2,0,(0==A322353(n))); \\ Antti Karttunen, Dec 06 2018

A319057 Minimum sum of a strict factorization of n into factors > 1.

Original entry on oeis.org

0, 2, 3, 4, 5, 5, 7, 6, 9, 7, 11, 7, 13, 9, 8, 10, 17, 9, 19, 9, 10, 13, 23, 9, 25, 15, 12, 11, 29, 10, 31, 12, 14, 19, 12, 11, 37, 21, 16, 11, 41, 12, 43, 15, 14, 25, 47, 12, 49, 15, 20, 17, 53, 14, 16, 13, 22, 31, 59, 12, 61, 33, 16, 14, 18, 16, 67, 21, 26
Offset: 1

Views

Author

Gus Wiseman, Sep 09 2018

Keywords

Comments

a(n) >= A001414(n), with equality iff n is squarefree or four times a squarefree number (i.e., A000188(n) <= 2). - Charlie Neder, Sep 10 2018

Examples

			The strict factorizations of 48 are (48), (2*24), (3*16), (4*12), (6*8), (2*3*8), (2*4*6), with sums 48, 26, 19, 16, 14, 13, 12 respectively, so a(48) = 12.
		

Crossrefs

Programs

  • Mathematica
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[strfacs[n/d],Min@@#>d&],{d,Rest[Divisors[n]]}]];
    Table[Min[Total/@strfacs[n]],{n,100}]

A320891 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into squarefree semiprimes.

Original entry on oeis.org

4, 9, 16, 24, 25, 40, 49, 54, 56, 64, 81, 88, 96, 104, 121, 135, 136, 144, 152, 160, 169, 184, 189, 224, 232, 240, 248, 250, 256, 289, 296, 297, 324, 328, 336, 344, 351, 352, 361, 375, 376, 384, 400, 416, 424, 459, 472, 486, 488, 513, 528, 529, 536, 544, 560
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct primes.
Also numbers with an even number x of prime factors, whose greatest prime multiplicity exceeds x/2.

Examples

			A complete list of all factorizations of 24 is:
  (2*2*2*3),
  (2*2*6), (2*3*4),
  (2*12), (3*8), (4*6),
  (24).
All of these contain at least one number that is not a squarefree semiprime, so 24 belongs to the sequence.
		

Crossrefs

Programs

  • Mathematica
    semfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[semfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    Select[Range[100],And[EvenQ[PrimeOmega[#]],semfacs[#]=={}]&]

A320894 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into distinct squarefree semiprimes.

Original entry on oeis.org

4, 9, 16, 24, 25, 36, 40, 49, 54, 56, 64, 81, 88, 96, 100, 104, 121, 135, 136, 144, 152, 160, 169, 184, 189, 196, 216, 224, 225, 232, 240, 248, 250, 256, 289, 296, 297, 324, 328, 336, 344, 351, 352, 360, 361, 375, 376, 384, 400, 416, 424, 441, 459, 472, 484
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct primes.

Examples

			A complete list of all strict factorizations of 24 is: (2*3*4), (2*12), (3*8), (4*6), (24). All of these contain at least one number that is not a squarefree semiprime, so 24 belongs to the sequence.
		

Crossrefs

Programs

  • Mathematica
    strsqfsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsqfsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    Select[Range[100],And[EvenQ[PrimeOmega[#]],strsqfsemfacs[#]=={}]&]

A320893 Numbers with an even number of prime factors (counted with multiplicity) that can be factored into squarefree semiprimes (A320911) but cannot be factored into distinct semiprimes (A320892).

Original entry on oeis.org

1296, 7776, 10000, 12960, 18144, 19440, 21600, 27216, 28512, 33696, 36000, 38416, 42336, 42768, 44064, 46656, 48600, 49248, 50544, 50625, 59616, 60000, 66096, 73872, 75168, 77760, 80352, 89424, 95256, 95904, 98784, 100000
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

A semiprime (A001358) is a product of any two not necessarily distinct primes.

Crossrefs

Programs

  • Mathematica
    sqfsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfsemfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    strsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Select[Range[10000],And[EvenQ[PrimeOmega[#]],strsemfacs[#]=={},sqfsemfacs[#]!={}]&]

A320665 Number of non-isomorphic multiset partitions of weight n with no singletons or vertices that appear only once.

Original entry on oeis.org

1, 0, 1, 1, 5, 6, 27, 47, 169, 406, 1327, 3790, 12560, 39919, 136821, 470589, 1687981, 6162696, 23173374, 88981796, 349969596, 1405386733, 5764142220, 24111709328, 102825231702, 446665313598, 1975339030948, 8888051121242, 40667889052853, 189126710033882, 893526261542899
Offset: 0

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Comments

The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}. This sequence counts non-isomorphic multiset partitions with no singletons whose dual also has no singletons.

Examples

			Non-isomorphic representatives of the a(2) = 1 through a(6) = 27 multiset partitions:
  {{1,1}}  {{1,1,1}}  {{1,1,1,1}}    {{1,1,1,1,1}}    {{1,1,1,1,1,1}}
                      {{1,1,2,2}}    {{1,1,2,2,2}}    {{1,1,1,2,2,2}}
                      {{1,1},{1,1}}  {{1,1},{1,1,1}}  {{1,1,2,2,2,2}}
                      {{1,1},{2,2}}  {{1,1},{1,2,2}}  {{1,1,2,2,3,3}}
                      {{1,2},{1,2}}  {{1,1},{2,2,2}}  {{1,1},{1,1,1,1}}
                                     {{1,2},{1,2,2}}  {{1,1,1},{1,1,1}}
                                                      {{1,1},{1,2,2,2}}
                                                      {{1,1,1},{2,2,2}}
                                                      {{1,1,2},{1,2,2}}
                                                      {{1,1},{2,2,2,2}}
                                                      {{1,1,2},{2,2,2}}
                                                      {{1,1},{2,2,3,3}}
                                                      {{1,1,2},{2,3,3}}
                                                      {{1,2},{1,1,2,2}}
                                                      {{1,2},{1,2,2,2}}
                                                      {{1,2},{1,2,3,3}}
                                                      {{1,2,2},{1,2,2}}
                                                      {{1,2,3},{1,2,3}}
                                                      {{2,2},{1,1,2,2}}
                                                      {{1,1},{1,1},{1,1}}
                                                      {{1,1},{1,2},{2,2}}
                                                      {{1,1},{2,2},{2,2}}
                                                      {{1,1},{2,2},{3,3}}
                                                      {{1,1},{2,3},{2,3}}
                                                      {{1,2},{1,2},{1,2}}
                                                      {{1,2},{1,2},{2,2}}
                                                      {{1,2},{1,3},{2,3}}
		

Crossrefs

Programs

  • PARI
    \\ See links in A339645 for combinatorial species functions.
    seq(n)={my(A=symGroupSeries(n)); NumUnlabeledObjsSeq(sCartProd(sExp(A-x*sv(1)), sExp(A-x*sv(1))))} \\ Andrew Howroyd, Jan 17 2023
    
  • PARI
    Vec(G(20,1)) \\ G defined in A369287. - Andrew Howroyd, Jan 28 2024

Extensions

Terms a(11) and beyond from Andrew Howroyd, Jan 17 2023

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

A080688 Resort the index of A064553 using A080444 and maintaining ascending order within each grouping: seen as a triangle read by rows, the n-th row contains the A001055(n) numbers m with A064553(m)=n.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 6, 11, 13, 8, 10, 17, 9, 19, 14, 23, 29, 12, 15, 22, 31, 37, 26, 41, 21, 43, 16, 20, 25, 34, 47, 53, 18, 33, 38, 59, 61, 28, 35, 46, 67, 39, 71, 58, 73, 79, 24, 30, 44, 51, 55, 62, 83, 49, 89, 74, 97, 27, 57, 101, 52, 65, 82
Offset: 1

Views

Author

Alford Arnold, Mar 23 2003

Keywords

Comments

The number 12 can be written as 3*2*2, 4*3, 6*2 and 12 corresponding to each of the four values (12,15,22,31) in the example. Note that A001055(12) = 4. Since A001055(n) depends only on the least prime signature, the values 1,2,4,6,8,12,16,24,30,32,36,... A025487 are of special interest when counting multisets. (see for example, A035310 and A035341).
A064553(T(n,k)) = A080444(n,k) = n for k=1..A001055(n); T(n,1) = A064554(n); T(n,A001055(n)) = A064554(n). - Reinhard Zumkeller, Oct 01 2012
Row n is the sorted list of shifted Heinz numbers of factorizations of n into factors > 1, where the shifted Heinz number of a factorization (y_1, ..., y_k) is prime(y_1 - 1) * ... * prime(y_k - 1). - Gus Wiseman, Sep 05 2018

Examples

			a(18),a(19),a(20) and a(21) are 12,15,22 and 31 because A064553(12,15,22,31) = (12,12,12,12) similarly, A064553(36,45,66,76,93,95,118,121,149) = (36,36,36,36,36,36,36,36,36)
From _Gus Wiseman_, Sep 05 2018: (Start)
Triangle begins:
   1
   2
   3
   4  5
   7
   6 11
  13
   8 10 17
   9 19
  14 23
  29
  12 15 22 31
  37
  26 41
  21 43
  16 20 25 34 47
Corresponding triangle of factorizations begins:
  (),
  (2),
  (3),
  (2*2), (4),
  (5),
  (2*3), (6),
  (7),
  (2*2*2), (2*4), (8),
  (3*3), (9),
  (2*5), (10),
  (11),
  (2*2*3), (3*4), (2*6), (12).
(End)
		

Crossrefs

Programs

  • Haskell
    a080688 n k = a080688_row n !! (k-1)
    a080688_row n = map (+ 1) $ take (a001055 n) $
                    elemIndices n $ map fromInteger a064553_list
    a080688_tabl = map a080688_row [1..]
    a080688_list = concat a080688_tabl
    -- Reinhard Zumkeller, Oct 01 2012
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
    Table[Sort[Table[Times@@Prime/@(f-1),{f,facs[n]}]],{n,20}] (* Gus Wiseman, Sep 05 2018 *)

Extensions

More terms from Sean A. Irvine, Oct 05 2011
Keyword tabf added and definition complemented accordingly by Reinhard Zumkeller, Oct 01 2012
Showing 1-10 of 19 results. Next