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-18 of 18 results.

A339661 Number of factorizations of n into distinct squarefree semiprimes.

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, 0, 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, 0, 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, Dec 19 2020

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
Also the number of strict multiset partitions of the multiset of prime factors of n, into distinct strict pairs.

Examples

			The a(n) factorizations for n = 210, 1260, 4620, 30030, 69300 are respectively 3, 2, 6, 15, 7:
  (6*35)   (6*10*21)  (6*10*77)   (6*55*91)    (6*10*15*77)
  (10*21)  (6*14*15)  (6*14*55)   (6*65*77)    (6*10*21*55)
  (14*15)             (6*22*35)   (10*33*91)   (6*10*33*35)
                      (10*14*33)  (10*39*77)   (6*14*15*55)
                      (10*21*22)  (14*33*65)   (6*15*22*35)
                      (14*15*22)  (14*39*55)   (10*14*15*33)
                                  (15*22*91)   (10*15*21*22)
                                  (15*26*77)
                                  (21*22*65)
                                  (21*26*55)
                                  (22*35*39)
                                  (26*33*35)
                                  (6*35*143)
                                  (10*21*143)
                                  (14*15*143)
		

Crossrefs

Dirichlet convolution of A008836 (Liouville's lambda) with A339742.
A050326 allows all squarefree numbers, non-strict case A050320.
A320656 is the not necessarily strict version.
A320911 lists all (not just distinct) products of squarefree semiprimes.
A322794 counts uniform factorizations, such as these.
A339561 lists positions of nonzero terms.
A001055 counts factorizations, with strict case A045778.
A001358 lists semiprimes, with squarefree case A006881.
A320655 counts factorizations into semiprimes, with strict case A322353.
The following count vertex-degree partitions and give their Heinz numbers:
- A000070 counts non-multigraphical partitions of 2n (A339620).
- A209816 counts multigraphical partitions (A320924).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).
- A339617 counts non-graphical partitions of 2n (A339618).
- A000569 counts graphical partitions (A320922).
The following count partitions of even length and give their Heinz numbers:
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).
- A339560 can be partitioned into distinct strict pairs (A339561).

Programs

  • Mathematica
    bfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[bfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Table[Length[bfacs[n]],{n,100}]
  • PARI
    A280710(n) = (bigomega(n)==2*issquarefree(n)); \\ From A280710.
    A339661(n, u=(1+n)) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1) && (dA280710(d), s += A339661(n/d, d))); (s)); \\ Antti Karttunen, May 02 2022

Formula

a(n) = Sum_{d|n} (-1)^A001222(d) * A339742(n/d).

Extensions

More terms and secondary offset added by Antti Karttunen, May 02 2022

A336137 Number of set partitions of the binary indices of n with equal block-sums.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 12 2020

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The a(n) set partitions for n = 7, 59, 119, 367, 127:
  {123}    {12456}      {123567}      {1234679}    {1234567}
  {12}{3}  {126}{45}    {1236}{57}    {12346}{79}  {1247}{356}
           {15}{24}{6}  {156}{237}    {1249}{367}  {1256}{347}
                        {17}{26}{35}  {1267}{349}  {1346}{257}
                                      {169}{2347}  {167}{2345}
                                                   {16}{25}{34}{7}
The binary indices of 382 are {2,3,4,5,6,7,9}, with equal block-sum set partitions:
  {{2,7},{3,6},{4,5},{9}}
  {{2,4,6},{3,9},{5,7}}
  {{2,7,9},{3,4,5,6}}
  {{2,3,4,9},{5,6,7}}
  {{2,3,6,7},{4,5,9}}
  {{2,4,5,7},{3,6,9}}
  {{2,3,4,5,6,7,9}}
so a(382) = 7.
		

Crossrefs

These set partitions are counted by A035470.
The version for twice-partitions is A279787.
The version for partitions of partitions is A305551.
The version for factorizations is A321455.
The version for normal multiset partitions is A326518.
The version for distinct block-sums is A336138.
Set partitions of binary indices are A050315.
Normal multiset partitions with equal lengths are A317583.
Normal multiset partitions with equal averages are A326520.
Multiset partitions with equal block-sums are ranked by A326534.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[bpe[n]],SameQ@@Total/@#&]],{n,0,100}]

A066637 Total number of elements in all factorizations of n with all factors > 1.

Original entry on oeis.org

0, 1, 1, 3, 1, 3, 1, 6, 3, 3, 1, 8, 1, 3, 3, 12, 1, 8, 1, 8, 3, 3, 1, 17, 3, 3, 6, 8, 1, 10, 1, 20, 3, 3, 3, 22, 1, 3, 3, 17, 1, 10, 1, 8, 8, 3, 1, 34, 3, 8, 3, 8, 1, 17, 3, 17, 3, 3, 1, 27, 1, 3, 8, 35, 3, 10, 1, 8, 3, 10, 1, 46, 1, 3, 8, 8, 3, 10, 1, 34, 12, 3, 1, 27, 3, 3, 3, 17, 1, 27, 3, 8, 3, 3, 3
Offset: 1

Views

Author

Amarnath Murthy, Dec 28 2001

Keywords

Comments

From Gus Wiseman, Apr 18 2021: (Start)
Number of ways to choose a factor index or position in a factorization of n. The version selecting a factor value is A339564. For example, the factorizations of n = 2, 4, 8, 12, 16, 24, 30 with a selected position (in parentheses) are:
((2)) ((4)) ((8)) ((12)) ((16)) ((24)) ((30))
((2)*2) ((2)*4) ((2)*6) ((2)*8) ((3)*8) ((5)*6)
(2*(2)) (2*(4)) (2*(6)) (2*(8)) (3*(8)) (5*(6))
((2)*2*2) ((3)*4) ((4)*4) ((4)*6) ((2)*15)
(2*(2)*2) (3*(4)) (4*(4)) (4*(6)) (2*(15))
(2*2*(2)) ((2)*2*3) ((2)*2*4) ((2)*12) ((3)*10)
(2*(2)*3) (2*(2)*4) (2*(12)) (3*(10))
(2*2*(3)) (2*2*(4)) ((2)*2*6) ((2)*3*5)
((2)*2*2*2) (2*(2)*6) (2*(3)*5)
(2*(2)*2*2) (2*2*(6)) (2*3*(5))
(2*2*(2)*2) ((2)*3*4)
(2*2*2*(2)) (2*(3)*4)
(2*3*(4))
((2)*2*2*3)
(2*(2)*2*3)
(2*2*(2)*3)
(2*2*2*(3))
(End)

Examples

			a(12) = 8: there are 4 factorizations of 12: (12), (6*2), (4*3), (3*2*2) having 1, 2, 2, 3 elements respectively, a total of 8.
		

References

  • Amarnath Murthy, Generalization of Partition function, Introducing Smarandache Factor partitions, Smarandache Notions Journal, Vol. 11, 1-2-3, Spring 2000.
  • Amarnath Murthy, Length and extent of Smarandache Factor partitions, Smarandache Notions Journal, Vol. 11, 1-2-3, Spring 2000.

Crossrefs

The version for normal multisets is A001787.
The version for compositions is A001792.
The version for partitions is A006128 (strict: A015723).
Choosing a value instead of position gives A339564.
A000070 counts partitions with a selected part.
A001055 counts factorizations.
A002033 and A074206 count ordered factorizations.
A067824 counts strict chains of divisors starting with n.
A336875 counts compositions with a selected part.

Programs

  • Maple
    # Return a list of lists which are factorizations (product representations)
    # of n. Within each sublist, the factors are sorted. A minimum factor in
    # each element of sublists returned can be specified with 'mincomp'.
    # If mincomp=2, the number of sublists contained in the list returned is A001055(n).
    # Example:
    # n=8 and mincomp=2 return [[2,2,2],[4,8],[8]]
    listProdRep := proc(n,mincomp)
        local dvs,resul,f,i,j,rli,tmp ;
        resul := [] ;
        # list returned is empty if n < mincomp
        if n >= mincomp then
            if n = 1 then
                RETURN([1]) ;
            else
                # compute the divisors, and take each divisor
                # as a head element (minimum element) of one of the
                # sublists. Example: for n=8 use {1,2,4,8}, and consider
                # (for mincomp=2) sublists [2,...], [4,...] and [8].
                dvs := numtheory[divisors](n) ;
                for i from 1 to nops(dvs) do
                    # select the head element 'f' from the divisors
                    f := op(i,dvs) ;
                    # if this is already the maximum divisor n
                    # itself, this head element is the last in
                    # the sublist
                    if f =n and f >= mincomp then
                        resul := [op(resul),[f]] ;
                    elif f >= mincomp then
                        # if this is not the maximum element
                        # n itself, produce all factorizations
                        # of the remaining factor recursively.
                        rli := procname(n/f,f) ;
                        # Prepend all the results produced
                        # from the recursion with the head
                        # element for the result.
                        for j from 1 to nops(rli) do
                            tmp := [f,op(op(j,rli))] ;
                            resul := [op(resul),tmp] ;
                        od ;
                    fi ;
                od ;
            fi ;
        fi ;
        resul ;
    end:
    A066637 := proc(n)
        local f,d;
        a := 0 ;
        for d in listProdRep(n,2) do
            a := a+nops(d) ;
        end do:
        a ;
    end proc: # R. J. Mathar, Jul 11 2013
    # second Maple program:
    with(numtheory):
    b:= proc(n, k) option remember; `if`(n>k, 0, [1$2])+
          `if`(isprime(n), 0, (p-> p+[0, p[1]])(add(
          `if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n})))
        end:
    a:= n-> `if`(n<2, 0, b(n$2)[2]):
    seq(a(n), n=1..120); # Alois P. Heinz, Feb 12 2019
  • Mathematica
    g[1, r_] := g[1, r]={1, 0}; g[n_, r_] := g[n, r]=Module[{ds, i, val}, ds=Select[Divisors[n], 1<#<=r&]; val={0, 0}+Sum[g[n/ds[[i]], ds[[i]]], {i, 1, Length[ds]}]; val+{0, val[[1]]}]; a[n_] := g[n, n][[2]]; a/@Range[95] (* g[n, r] = {c, f}, where c is the number of factorizations of n with factors <= r and f is the total number of factors in them. - Dean Hickerson, Oct 28 2002 *)
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];Table[Sum[Length[fac],{fac,facs[n]}],{n,50}] (* Gus Wiseman, Apr 18 2021 *)

A371733 Maximal length of a factorization of n into factors > 1 all having the same sum of prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 13 2024

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. Sum of prime indices is given by A056239.
Factorizations into factors all having the same sum of prime indices are counted by A321455.

Examples

			The factorizations of 588 of this type are (7*7*12), (21*28), (588), so a(588) = 3.
The factorizations of 900 of this type are (5*5*6*6), (9*10*10), (25*36), (30*30), (900), so a(900) = 4.
		

Crossrefs

Positions of 1's are A321453, counted by A321451.
Positions of terms > 1 are A321454, counted by A321452.
Factorizations of this type are counted by A321455, different sums A321469.
For different sums instead of same sums we have A371734.
For set partitions of binary indices we have A371735.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A321142 and A371794 count non-biquanimous strict partitions.
A371789 counts non-quanimous sets, differences A371790.
A371796 counts quanimous sets, differences A371797.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&, Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    Table[Max[Length/@Select[facs[n],SameQ@@hwt/@#&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
    all_have_same_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == 1));
    A371733(n, m=n, facs=List([])) = if(1==n, if(all_have_same_sum_of_pis(facs),#facs,0), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s = max(s, A371733(n/d, d, newfacs)))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(108) by Antti Karttunen, Jan 20 2025

A322788 Irregular triangle read by rows where T(n,k) is the number of uniform multiset partitions of a multiset with d = A027750(n,k) copies of each integer from 1 to n/d.

Original entry on oeis.org

1, 2, 2, 2, 2, 5, 4, 3, 2, 2, 27, 11, 6, 4, 2, 2, 142, 29, 8, 4, 282, 12, 3, 1073, 101, 8, 4, 2, 2, 32034, 1581, 234, 75, 20, 6, 2, 2, 136853, 2660, 10, 4, 1527528, 1985, 91, 4, 4661087, 64596, 648, 20, 5, 2, 2, 227932993, 1280333, 41945, 231, 28, 6
Offset: 1

Views

Author

Gus Wiseman, Dec 26 2018

Keywords

Comments

A multiset partition is uniform if all parts have the same size.

Examples

			Triangle begins:
     1
     2    2
     2    2
     5    4    3
     2    2
    27   11    6    4
     2    2
   142   29    8    4
   282   12    3
  1073  101    8    4
The multiset partitions counted under row 6:
  {123456}          {112233}          {111222}          {111111}
  {123}{456}        {112}{233}        {111}{222}        {111}{111}
  {124}{356}        {113}{223}        {112}{122}        {11}{11}{11}
  {125}{346}        {122}{133}        {11}{12}{22}      {1}{1}{1}{1}{1}{1}
  {126}{345}        {123}{123}        {12}{12}{12}
  {134}{256}        {11}{22}{33}      {1}{1}{1}{2}{2}{2}
  {135}{246}        {11}{23}{23}
  {136}{245}        {12}{12}{33}
  {145}{236}        {12}{13}{23}
  {146}{235}        {13}{13}{22}
  {156}{234}        {1}{1}{2}{2}{3}{3}
  {12}{34}{56}
  {12}{35}{46}
  {12}{36}{45}
  {13}{24}{56}
  {13}{25}{46}
  {13}{26}{45}
  {14}{23}{56}
  {14}{25}{36}
  {14}{26}{35}
  {15}{23}{46}
  {15}{24}{36}
  {15}{26}{34}
  {16}{23}{45}
  {16}{24}{35}
  {16}{25}{34}
  {1}{2}{3}{4}{5}{6}
		

Crossrefs

Programs

  • Mathematica
    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[Select[mps[Join@@Table[Range[n/d],{d}]],SameQ@@Length/@#&]],{n,10},{d,Divisors[n]}]

Formula

T(n,k) = A322794(A002110(n/d)^d), where d = A027750(n,k).

Extensions

More terms from Alois P. Heinz, Jan 30 2019
Terms a(38) and beyond from Andrew Howroyd, Feb 03 2022
Edited by Peter Munn, Mar 05 2025

A371734 Maximal length of a factorization of n into factors > 1 all having different sums of prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 13 2024

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. Sum of prime indices is given by A056239.
Factorizations into factors > 1 all having different sums of prime indices are counted by A321469.

Examples

			The factorizations of 90 of this type are (2*3*15), (2*5*9), (2*45), (3*30), (5*18), (6*15), (90), so a(90) = 3.
		

Crossrefs

For set partitions of binary indices we have A000120, same sums A371735.
Positions of 1's are A000430.
Positions of terms > 1 are A080257.
Factorizations of this type are counted by A321469, same sums A321455.
For same instead of different sums we have A371733.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A321451 counts non-quanimous partitions, ranks A321453.
A321452 counts quanimous partitions, ranks A321454.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&, Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    Table[Max[Length/@Select[facs[n],UnsameQ@@hwt/@#&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
    all_have_different_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == #facs));
    A371734(n, m=n, facs=List([])) = if(1==n, if(all_have_different_sum_of_pis(facs),#facs,0), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s = max(s,A371734(n/d, d, newfacs)))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A371735 Maximal length of a set partition of the binary indices of n into blocks all having the same sum.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 14 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
If a(n) = k then the binary indices of n (row n of A048793) are k-quanimous (counted by A371783).

Examples

			The binary indices of 119 are {1,2,3,5,6,7}, and the set partitions into blocks with the same sum are:
  {{1,7},{2,6},{3,5}}
  {{1,5,6},{2,3,7}}
  {{1,2,3,6},{5,7}}
  {{1,2,3,5,6,7}}
So a(119) = 3.
		

Crossrefs

Set partitions of this type are counted by A035470, A336137.
A version for factorizations is A371733.
Positions of 1's are A371738.
Positions of terms > 1 are A371784.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A070939 gives length of binary expansion.
A321452 counts quanimous partitions, ranks A321454.
A326031 gives weight of the set-system with BII-number n.
A371783 counts k-quanimous partitions.
A371789 counts non-quanimous sets, differences A371790.
A371796 counts quanimous sets, differences A371797.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Table[Max[Length/@Select[sps[bix[n]],SameQ@@Total/@#&]],{n,0,100}]

A327900 Nonprime squarefree numbers whose prime indices all have the same Omega (number of prime factors counted with multiplicity).

Original entry on oeis.org

1, 15, 33, 51, 55, 85, 91, 93, 123, 155, 161, 165, 177, 187, 201, 203, 205, 249, 255, 295, 299, 301, 327, 329, 335, 341, 377, 381, 415, 451, 465, 471, 511, 527, 537, 545, 553, 559, 561, 573, 611, 615, 633, 635, 649, 667, 679, 697, 703, 707, 723, 737, 785, 831
Offset: 1

Views

Author

Gus Wiseman, Sep 30 2019

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.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
   15: {2,3}
   33: {2,5}
   51: {2,7}
   55: {3,5}
   85: {3,7}
   91: {4,6}
   93: {2,11}
  123: {2,13}
  155: {3,11}
  161: {4,9}
  165: {2,3,5}
  177: {2,17}
  187: {5,7}
  201: {2,19}
  203: {4,10}
  205: {3,13}
  249: {2,23}
  255: {2,3,7}
  295: {3,17}
		

Crossrefs

The case including primes and nonsquarefree numbers is A320324.
The version for sum of prime indices is A327901.
The version for mean of prime indices is A327902.

Programs

  • Mathematica
    Select[Range[1000],!PrimeQ[#]&&SquareFreeQ[#]&&SameQ@@PrimeOmega/@PrimePi/@First/@FactorInteger[#]&]
Previous Showing 11-18 of 18 results.