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

A320924 Heinz numbers of multigraphical partitions.

Original entry on oeis.org

1, 4, 9, 12, 16, 25, 27, 30, 36, 40, 48, 49, 63, 64, 70, 75, 81, 84, 90, 100, 108, 112, 120, 121, 144, 147, 154, 160, 165, 169, 175, 189, 192, 196, 198, 210, 220, 225, 243, 250, 252, 256, 264, 270, 273, 280, 286, 289, 300, 324, 325, 336, 343, 351, 352, 360
Offset: 1

Views

Author

Gus Wiseman, Oct 24 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
An integer partition is multigraphical if it comprises the multiset of vertex-degrees of some multigraph.
Also Heinz numbers of integer partitions of even numbers whose greatest part is less than or equal to half the sum of parts, i.e., numbers n whose sum of prime indices A056239(n) is even and at least twice the greatest prime index A061395(n). - Gus Wiseman, May 23 2021

Examples

			The sequence of all multigraphical partitions begins: (), (11), (22), (211), (1111), (33), (222), (321), (2211), (3111), (21111), (44), (422), (111111), (431), (332), (2222), (4211), (3221), (3311), (22211), (41111), (32111), (55), (221111).
From _Gus Wiseman_, May 23 2021: (Start)
The sequence of terms together with their prime indices and a multigraph realizing each begins:
    1:      () | {}
    4:    (11) | {{1,2}}
    9:    (22) | {{1,2},{1,2}}
   12:   (112) | {{1,3},{2,3}}
   16:  (1111) | {{1,2},{3,4}}
   25:    (33) | {{1,2},{1,2},{1,2}}
   27:   (222) | {{1,2},{1,3},{2,3}}
   30:   (123) | {{1,3},{2,3},{2,3}}
   36:  (1122) | {{1,2},{3,4},{3,4}}
   40:  (1113) | {{1,4},{2,4},{3,4}}
   48: (11112) | {{1,2},{3,5},{4,5}}
   49:    (44) | {{1,2},{1,2},{1,2},{1,2}}
   63:   (224) | {{1,3},{1,3},{2,3},{2,3}}
(End)
		

Crossrefs

These partitions are counted by A209816.
The case with odd weights is A322109.
The conjugate case of equality is A340387.
The conjugate version with odd weights allowed is A344291.
The conjugate opposite version is A344292.
The opposite version with odd weights allowed is A344296.
The conjugate version is A344413.
The conjugate opposite version with odd weights allowed is A344414.
The case of equality is A344415.
The opposite version is A344416.
A000070 counts non-multigraphical partitions.
A025065 counts palindromic partitions.
A035363 counts partitions into even parts.
A056239 adds up prime indices, row sums of A112798.
A110618 counts partitions that are the vertex-degrees of some set multipartition with no singletons.
A334201 adds up all prime indices except the greatest.

Programs

  • Mathematica
    prptns[m_]:=Union[Sort/@If[Length[m]==0,{{}},Join@@Table[Prepend[#,m[[ipr]]]&/@prptns[Delete[m,List/@ipr]],{ipr,Select[Prepend[{#},1]&/@Select[Range[2,Length[m]],m[[#]]>m[[#-1]]&],UnsameQ@@m[[#]]&]}]]];
    Select[Range[1000],prptns[Flatten[MapIndexed[Table[#2,{#1}]&,If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]]]]!={}&]

Formula

Members m of A300061 such that A061395(m) <= A056239(m)/2. - Gus Wiseman, May 23 2021

A320922 Heinz numbers of graphical partitions.

Original entry on oeis.org

1, 4, 12, 16, 27, 36, 40, 48, 64, 81, 90, 108, 112, 120, 144, 160, 192, 225, 243, 252, 256, 270, 300, 324, 336, 352, 360, 400, 432, 448, 480, 567, 576, 625, 630, 640, 675, 729, 750, 756, 768, 792, 810, 832, 840, 900, 972, 1000, 1008, 1024, 1056, 1080, 1120
Offset: 1

Views

Author

Gus Wiseman, Oct 24 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
An integer partition is graphical if it comprises the vertex-degrees of some simple graph.

Examples

			The sequence of all graphical partitions begins: (), (11), (211), (1111), (222), (2211), (3111), (21111), (111111), (2222), (3221), (22211), (41111), (32111), (221111), (311111), (2111111), (3322), (22222), (42211).
		

Crossrefs

Programs

  • Mathematica
    prptns[m_]:=Union[Sort/@If[Length[m]==0,{{}},Join@@Table[Prepend[#,m[[ipr]]]&/@prptns[Delete[m,List/@ipr]],{ipr,Select[Prepend[{#},1]&/@Select[Range[2,Length[m]],m[[#]]>m[[#-1]]&],UnsameQ@@m[[#]]&]}]]];
    Select[Range[1000],Select[prptns[Flatten[MapIndexed[Table[#2,{#1}]&,If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]]]],UnsameQ@@#&]!={}&]

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

A368095 Number of non-isomorphic set-systems of weight n satisfying a strict version of the axiom of choice.

Original entry on oeis.org

1, 1, 2, 4, 8, 17, 39, 86, 208, 508, 1304
Offset: 0

Views

Author

Gus Wiseman, Dec 24 2023

Keywords

Comments

A set-system is a finite set of finite nonempty sets. The weight of a set-system is the sum of cardinalities of its elements.
The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 17 set-systems:
  {1}  {12}    {123}      {1234}        {12345}
       {1}{2}  {1}{23}    {1}{234}      {1}{2345}
               {2}{12}    {12}{34}      {12}{345}
               {1}{2}{3}  {13}{23}      {14}{234}
                          {3}{123}      {23}{123}
                          {1}{2}{34}    {4}{1234}
                          {1}{3}{23}    {1}{2}{345}
                          {1}{2}{3}{4}  {1}{23}{45}
                                        {1}{24}{34}
                                        {1}{4}{234}
                                        {2}{13}{23}
                                        {2}{3}{123}
                                        {3}{13}{23}
                                        {4}{12}{34}
                                        {1}{2}{3}{45}
                                        {1}{2}{4}{34}
                                        {1}{2}{3}{4}{5}
		

Crossrefs

For labeled graphs we have A133686, complement A367867.
For unlabeled graphs we have A134964, complement A140637.
For set-systems we have A367902, complement A367903.
These set-systems have BII-numbers A367906, complement A367907.
The complement is A368094, connected A368409.
Repeats allowed: A368098, ranks A368100, complement A368097, ranks A355529.
Minimal multiset partitions not of this type are counted by A368187.
The connected case is A368410.
Factorizations of this type are counted by A368414, complement A368413.
Allowing repeated edges gives A368422, complement A368421.
A000110 counts set-partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

Programs

  • Mathematica
    Table[Length[Select[bmp[n], UnsameQ@@#&&And@@UnsameQ@@@#&&Select[Tuples[#], UnsameQ@@#&]!={}&]], {n,0,10}]

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

Original entry on oeis.org

1, 4, 6, 9, 10, 14, 15, 21, 22, 24, 25, 26, 33, 34, 35, 36, 38, 39, 40, 46, 49, 51, 54, 55, 56, 57, 58, 60, 62, 65, 69, 74, 77, 82, 84, 85, 86, 87, 88, 90, 91, 93, 94, 95, 100, 104, 106, 111, 115, 118, 119, 121, 122, 123, 126, 129, 132, 133, 134, 135, 136, 140
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Comments

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

Examples

			9000 is in the sequence and can be factored in either of two ways: (4*6*15*25) or (4*9*10*25).
		

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[100],And[EvenQ[PrimeOmega[#]],strsemfacs[#]!={}]&]

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 4, 1, 1, 1, 5, 1, 4, 1, 4, 1, 1, 1, 7, 2, 1, 3, 4, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 4, 4, 1, 1, 12, 2, 4, 1, 4, 1, 7, 1, 7, 1, 1, 1, 7, 1, 1, 4, 11, 1, 1, 1, 4, 1, 1, 1, 16, 1, 1, 4, 4, 1, 1, 1, 12, 5, 1, 1, 7, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Jul 18 2018

Keywords

Comments

In a factorization, two primes are equivalent if each factor has in its prime factorization the same multiplicity of both primes.

Examples

			The a(36) = 7 factorizations are (2*2*3*3), (2*2*9), (2*3*6), (3*3*4), (2*18), (3*12), (4*9). Missing from this list are (6*6) and (36).
		

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]]}]];
    dual[eds_]:=Table[First/@Position[eds,x],{x,Union@@eds}];
    Table[Length[Select[facs[n],UnsameQ@@dual[primeMS/@#]&]],{n,100}]

Formula

a(prime^n) = A000041(n).
a(squarefree) = 1.

A320663 Number of non-isomorphic multiset partitions of weight n using singletons or pairs.

Original entry on oeis.org

1, 1, 4, 7, 21, 40, 106, 216, 534, 1139, 2715, 5962, 14012, 31420, 73484, 167617, 392714, 908600, 2140429, 5015655, 11905145, 28228533, 67590229, 162067916, 391695348, 949359190, 2316618809, 5673557284, 13979155798, 34583650498, 86034613145, 214948212879
Offset: 0

Views

Author

Gus Wiseman, Oct 18 2018

Keywords

Examples

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

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    gs(v) = {sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(2*v[i]*v[j]/g))) + sum(i=1, #v, my(r=v[i]); (1 + (1+r)%2)*x^r + ((1+r)\2)*x^(2*r))}
    a(n)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(gs(p) + O(x*x^n), -n))[n]); s/n!} \\ Andrew Howroyd, Oct 26 2018

Extensions

Terms a(11) and beyond from Andrew Howroyd, Oct 26 2018

A320732 Number of factorizations of n into primes or semiprimes.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 20 2018

Keywords

Examples

			The a(60) = 5 factorizations are (2*2*3*5), (2*2*15), (2*3*10), (2*5*6), (3*4*5), (4*15), (6*10).
		

Crossrefs

Programs

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

A020555 Number of multigraphs on n labeled edges (with loops). Also number of genetically distinct states amongst n individuals.

Original entry on oeis.org

1, 2, 9, 66, 712, 10457, 198091, 4659138, 132315780, 4441561814, 173290498279, 7751828612725, 393110572846777, 22385579339430539, 1419799938299929267, 99593312799819072788, 7678949893962472351181, 647265784993486603555551, 59357523410046023899154274
Offset: 0

Views

Author

Gilbert Labelle (gilbert(AT)lacim.uqam.ca), Simon Plouffe, N. J. A. Sloane

Keywords

Comments

Also the number of factorizations of (p_n#)^2. - David W. Wilson, Apr 30 2001
Also the number of multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}. - Gus Wiseman, Jul 18 2018
a(n) gives the number of genetically distinct states for n diploid individuals in the case that maternal and paternal alleles transmitted to the individuals are not distinguished (if maternal and paternal alleles are distinguished, then the number of states is A000110(2n)). - Noah A Rosenberg, Aug 23 2022

Examples

			From _Gus Wiseman_, Jul 18 2018: (Start)
The a(2) = 9 multiset partitions of {1, 1, 2, 2}:
  (1122),
  (1)(122), (2)(112), (11)(22), (12)(12),
  (1)(1)(22), (1)(2)(12), (2)(2)(11),
  (1)(1)(2)(2).
(End)
		

References

  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018
  • E. Keith Lloyd, Math. Proc. Camb. Phil. Soc., vol. 103 (1988), 277-284.
  • A. Murthy, Generalization of partition function, introducing Smarandache factor partitions. Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000.
  • G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004.

Crossrefs

Row n=2 of A219727. - Alois P. Heinz, Nov 26 2012
See also A322764. Row 0 of the array in A322765.
Main diagonal of A346500.

Programs

  • Maple
    B := n -> combinat[bell](n):
    P := proc(m,n) local k; global B; option remember;
    if n = 0 then B(m)  else
    (1/2)*( P(m+2,n-1) + P(m+1,n-1) + add( binomial(n-1,k)*P(m,k), k=0..n-1) ); fi; end;
    r:=m->[seq(P(m,n),n=0..20)]; r(0); # N. J. A. Sloane, Dec 30 2018
  • Mathematica
    max = 16; s = Series[Exp[-3/2 + Exp[x]/2]*Sum[Exp[Binomial[n+1, 2]*x]/n!, {n, 0, 3*max }], {x, 0, max}] // Normal; a[n_] := SeriesCoefficient[s, {x, 0, n}]*n!; Table[a[n] // Round, {n, 0, max} ] (* Jean-François Alcover, Apr 23 2014, after Vladeta Jovovic *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Length[mps[Ceiling[Range[1/2,n,1/2]]]],{n,5}] (* Gus Wiseman, Jul 18 2018 *)

Formula

Lloyd's article gives a complicated explicit formula.
E.g.f.: exp(-3/2 + exp(x)/2)*Sum_{n>=0} exp(binomial(n+1, 2)*x)/n! [probably in the Labelle paper]. - Vladeta Jovovic, Apr 27 2004
a(n) = A001055(A002110(n)^2). - Alois P. Heinz, Aug 23 2022

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[#]=={}]&]
Previous Showing 11-20 of 82 results. Next