A296119 -id:A296119 - cp's OEIS frontend

A371453 Numbers whose binary indices are all squarefree semiprimes.

32, 512, 544, 8192, 8224, 8704, 8736, 16384, 16416, 16896, 16928, 24576, 24608, 25088, 25120, 1048576, 1048608, 1049088, 1049120, 1056768, 1056800, 1057280, 1057312, 1064960, 1064992, 1065472, 1065504, 1073152, 1073184, 1073664, 1073696, 2097152, 2097184

Offset: 1

Gus Wiseman, Apr 02 2024

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.

			The terms together with their binary expansions and binary indices begin:
       32:                 100000 ~ {6}
      512:             1000000000 ~ {10}
      544:             1000100000 ~ {6,10}
     8192:         10000000000000 ~ {14}
     8224:         10000000100000 ~ {6,14}
     8704:         10001000000000 ~ {10,14}
     8736:         10001000100000 ~ {6,10,14}
    16384:        100000000000000 ~ {15}
    16416:        100000000100000 ~ {6,15}
    16896:        100001000000000 ~ {10,15}
    16928:        100001000100000 ~ {6,10,15}
    24576:        110000000000000 ~ {14,15}
    24608:        110000000100000 ~ {6,14,15}
    25088:        110001000000000 ~ {10,14,15}
    25120:        110001000100000 ~ {6,10,14,15}
  1048576:  100000000000000000000 ~ {21}

Partitions of this type are counted by A002100, squarefree case of A101048.
For primes instead of squarefree semiprimes we get A326782.
For prime indices instead of binary indices we have A339113, A339112.
Allowing any squarefree numbers gives A368533.
This is the squarefree case of A371454.
A001358 lists squarefree semiprimes, squarefree A006881.
A005117 lists squarefree numbers.
A048793 lists binary indices, reverse A272020, length A000120, sum A029931.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.
Cf. A087086, A112798, A296119, A302478, A326031, A367905, A368109, A371450.

M:= 26: # for terms < 2^M
P:= select(isprime, [$2..(M+1)/2]): nP:= nops(P):
S:= select(`<`,{seq(seq(P[i]*P[j],i=1..j-1),j=1..nP)},M+1):
R:= map(proc(s) local i; add(2^(i-1),i=s) end proc, combinat:-powerset(S) minus {{}}):
sort(convert(R,list)); # Robert Israel, Apr 04 2024

bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
sqfsemi[n_]:=SquareFreeQ[n]&&PrimeOmega[n]==2;
Select[Range[10000],And@@sqfsemi/@bix[#]&]

def A371453(n): return sum(1<<A006881(i)-1 for i, j in enumerate(bin(n)[:1:-1],1) if j=='1')

from math import isqrt
from sympy import primepi, primerange
def A371453(n):
    def f(x,n): return int(n+x+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    def A006881(n):
        m, k = n, f(n,n)
        while m != k:
            m, k = k, f(k,n)
        return m
    return sum(1<<A006881(i)-1 for i, j in enumerate(bin(n)[:1:-1],1) if j=='1') # Chai Wah Wu, Aug 16 2024

A371454 Numbers whose binary indices are all semiprimes.

Original entry on oeis.org

8, 32, 40, 256, 264, 288, 296, 512, 520, 544, 552, 768, 776, 800, 808, 8192, 8200, 8224, 8232, 8448, 8456, 8480, 8488, 8704, 8712, 8736, 8744, 8960, 8968, 8992, 9000, 16384, 16392, 16416, 16424, 16640, 16648, 16672, 16680, 16896, 16904, 16928, 16936, 17152

Offset: 1

Gus Wiseman, Apr 02 2024

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.

			The terms together with their binary expansions and binary indices begin:
     8:           1000 ~ {4}
    32:         100000 ~ {6}
    40:         101000 ~ {4,6}
   256:      100000000 ~ {9}
   264:      100001000 ~ {4,9}
   288:      100100000 ~ {6,9}
   296:      100101000 ~ {4,6,9}
   512:     1000000000 ~ {10}
   520:     1000001000 ~ {4,10}
   544:     1000100000 ~ {6,10}
   552:     1000101000 ~ {4,6,10}
   768:     1100000000 ~ {9,10}
   776:     1100001000 ~ {4,9,10}
   800:     1100100000 ~ {6,9,10}
   808:     1100101000 ~ {4,6,9,10}

Partitions of this type are counted by A101048, squarefree case A002100.
For primes instead of semiprimes we get A326782.
For prime indices instead of binary indices we have A339112, A339113.
The squarefree case is A371453.
A001358 lists semiprimes, squarefree A006881.
A005117 lists squarefree numbers.
A048793 lists binary indices, reverse A272020, length A000120, sum A029931.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.
Cf. A087086, A296119, A302478, A326031, A367905, A368109, A368533, A371450.

bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
semi[n_]:=PrimeOmega[n]==2;
Select[Range[10000],And@@semi/@bix[#]&]

from math import isqrt
from sympy import primepi, primerange
def A371454(n):
    def f(x,n): return int(n+x+((t:=primepi(s:=isqrt(x)))*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    def A001358(n):
        m, k = n, f(n,n)
        while m != k:
            m, k = k, f(k,n)
        return m
    return sum(1<<A001358(i)-1 for i, j in enumerate(bin(n)[:1:-1],1) if j=='1') # Chai Wah Wu, Aug 16 2024

A301753 Number of ways to choose a strict rooted partition of each part in a rooted partition of n.

Original entry on oeis.org

1, 1, 2, 3, 6, 9, 16, 25, 43, 66, 108, 166, 269, 408, 643, 975, 1517, 2277, 3497, 5223, 7936, 11803, 17736, 26219, 39174, 57594, 85299, 124957, 183987, 268158, 392685, 569987, 830282, 1200843, 1740422, 2507823, 3620550, 5197885, 7472229, 10694865, 15319700

Offset: 1

Gus Wiseman, Mar 26 2018

nonn

A rooted partition of n is an integer partition of n - 1.

			The a(7) = 16 rooted twice-partitions:
(5), (32), (41),
(2)(2), (3)(1), (4)(), (21)(1), (31)(),
(1)(1)(1), (2)(1)(), (3)()(), (21)()(),
(1)(1)()(), (2)()()(),
(1)()()()(),
()()()()()().

Andrew Howroyd, Table of n, a(n) for n = 1..500

Cf. A002865, A032305, A063834, A093637, A270995, A281113, A296119, A301422, A301462, A301467, A301480, A301706, A301751.

nn=50;
ser=x*Product[1/(1-PartitionsQ[n-1]x^n),{n,nn}];
Table[SeriesCoefficient[ser,{x,0,n}],{n,nn}]

seq(n)={my(u=Vec(prod(k=1, n-1, 1 + x^k + O(x^n)))); Vec(1/prod(k=1, n-1, 1 - u[k]*x^k + O(x^n)))} \\ Andrew Howroyd, Aug 29 2018

O.g.f.: x * Product_{n > 0} 1/(1 - A000009(n-1) x^n).

A383310 Number of ways to choose a strict multiset partition of a factorization of n into factors > 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 5, 2, 3, 1, 8, 1, 3, 3, 9, 1, 8, 1, 8, 3, 3, 1, 20, 2, 3, 5, 8, 1, 12, 1, 19, 3, 3, 3, 24, 1, 3, 3, 20, 1, 12, 1, 8, 8, 3, 1, 46, 2, 8, 3, 8, 1, 20, 3, 20, 3, 3, 1, 38, 1, 3, 8, 37, 3, 12, 1, 8, 3, 12, 1, 67, 1, 3, 8, 8, 3, 12, 1, 46, 9, 3

Offset: 1

Gus Wiseman, Apr 26 2025

nonn

			The a(36) = 24 choices:
  {{2,2,3,3}}  {{2},{2,3,3}}  {{2},{3},{2,3}}
  {{2,2,9}}    {{3},{2,2,3}}  {{2},{3},{6}}
  {{2,3,6}}    {{2,2},{3,3}}
  {{2,18}}     {{2},{2,9}}
  {{3,3,4}}    {{9},{2,2}}
  {{3,12}}     {{2},{3,6}}
  {{4,9}}      {{3},{2,6}}
  {{6,6}}      {{6},{2,3}}
  {{36}}       {{2},{18}}
               {{3},{3,4}}
               {{4},{3,3}}
               {{3},{12}}
               {{4},{9}}

The case of a unique choice (positions of 1) is A008578.
This is the strict case of A050336.
For distinct strict blocks we have A050345.
For integer partitions we have A261049, strict case of A001970.
For strict blocks that are not necessarily distinct we have A296119.
Twice-partitions of this type are counted by A296122.
For normal multisets we have A317776, strict case of A255906.
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers, distinct A050326.
A281113 counts twice-factorizations, strict A296121, see A296118, A296120.
Cf. A000009, A005117, A045782, A050342, A279785, A293243, A293511, A302494, A316439, A358914, A381992, A382201.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
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[Sum[Length[Select[mps[y],UnsameQ@@#&]],{y,facs[n]}],{n,30}]

A383308 Number of integer partitions of n that can be partitioned into sets with a common sum.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 8, 6, 10, 13, 15, 13, 31

Offset: 0

Gus Wiseman, Apr 25 2025

Any strict partition can be partitioned into a single set, so we have a lower bound a(n) >= A000009(n).

			The multiset (3,2,2,1,1) has partition {{3},{1,2},{1,2}}, so is counted under a(9).
The a(1) = 1 through a(9) = 13 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)         (9)
       (11)  (21)   (22)    (32)     (33)      (43)       (44)        (54)
             (111)  (31)    (41)     (42)      (52)       (53)        (63)
                    (1111)  (11111)  (51)      (61)       (62)        (72)
                                     (222)     (421)      (71)        (81)
                                     (321)     (1111111)  (431)       (333)
                                     (2211)               (521)       (432)
                                     (111111)             (2222)      (531)
                                                          (3311)      (621)
                                                          (11111111)  (3321)
                                                                      (32211)
                                                                      (222111)
                                                                      (111111111)

Twice-partitions of this type (into sets with a common sum) are counted by A279788.
Multiset partitions of this type are ranked by A326534 /\ A302478.
For distinct instead of equal sums we have A381992, see also A382077.
The complement is counted by A381994, ranks A381719.
Partitions of prime indices of this type are counted by A382080.
Normal multiset partitions of this type are counted by A382429, see A326518.
For constant instead of strict blocks we have A383093, ranks A383014.
A000041 counts integer partitions, strict A000009.
A001055 counts factorizations, strict A045778.
Set multipartitions: A089259, A116540, A270995, A296119, A318360.
Cf. A050320, A293511, A321455, A358914, A381633, A381717, A381993.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
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[IntegerPartitions[n],Length[Select[mps[#],And@@UnsameQ@@@#&&SameQ@@Total/@#&]]>0&]],{n,0,10}]

A383311 Number of ways to choose a set multipartition (multiset of sets) of a factorization of n into factors > 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 7, 1, 3, 3, 7, 1, 7, 1, 7, 3, 3, 1, 16, 2, 3, 4, 7, 1, 12, 1, 12, 3, 3, 3, 20, 1, 3, 3, 16, 1, 12, 1, 7, 7, 3, 1, 33, 2, 7, 3, 7, 1, 16, 3, 16, 3, 3, 1, 34, 1, 3, 7, 22, 3, 12, 1, 7, 3, 12, 1, 49, 1, 3, 7, 7, 3, 12, 1, 33, 7, 3

Offset: 1

Gus Wiseman, Apr 28 2025

nonn

First differs from A296119 at a(36) = 20, A296119(36) = 21.

			The a(36) = 20 choices are:
  {{2,3,6}}  {{2,3},{2,3}}  {{2},{3},{2,3}}  {{2},{2},{3},{3}}
  {{2,18}}   {{2},{2,9}}    {{2},{2},{9}}
  {{3,12}}   {{2},{3,6}}    {{2},{3},{6}}
  {{4,9}}    {{3},{2,6}}    {{3},{3},{4}}
  {{36}}     {{6},{2,3}}
             {{2},{18}}
             {{3},{3,4}}
             {{3},{12}}
             {{4},{9}}
             {{6},{6}}

The case of a unique choice (positions of 1) is A008578.
For multisets of multisets we have A050336.
For sets of sets we have A050345.
For normal multisets we have A116540, strong A330783.
For integer partitions instead of factorizations we have A089259.
Twice-partitions of this type are counted by A270995.
For sets of multisets we have A383310 (distinct products A296118).
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers, distinct A050326.
A281113 counts twice-factorizations, see A294788, A296120, A296121.
A302478 gives MM-numbers of set multipartitions.
A302494 gives MM-numbers of sets of sets.
Cf. A000009, A001970, A005117, A050342, A116539, A279785, A293243, A293511, A296119, A316439, A381992, A382077.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
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[Sum[Length[Select[mps[y], And@@UnsameQ@@@#&]], {y,facs[n]}],{n,100}]

A368603 Products of odd primes of squarefree index. MM-numbers of set multipartitions.

Original entry on oeis.org

1, 3, 5, 9, 11, 13, 15, 17, 25, 27, 29, 31, 33, 39, 41, 43, 45, 47, 51, 55, 59, 65, 67, 73, 75, 79, 81, 83, 85, 87, 93, 99, 101, 109, 113, 117, 121, 123, 125, 127, 129, 135, 137, 139, 141, 143, 145, 149, 153, 155, 157, 163, 165, 167, 169, 177, 179, 181, 187

Offset: 1

Gus Wiseman, Jan 08 2024

nonn

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. The multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

A set multipartition is a finite multiset of finite nonempty sets.

			The terms together with the corresponding set multipartitions begin:
   1: {}
   3: {{1}}
   5: {{2}}
   9: {{1},{1}}
  11: {{3}}
  13: {{1,2}}
  15: {{1},{2}}
  17: {{4}}
  25: {{2},{2}}
  27: {{1},{1},{1}}
  29: {{1,3}}
  31: {{5}}
  33: {{1},{3}}
  39: {{1},{1,2}}
  41: {{6}}
  43: {{1,4}}
  45: {{1},{1},{2}}

Odd case of A302478.
Products of odd terms of A302491.
A049311 counts non-isomorphic set multipartitions, strict A283877.
A050320 counts set multipartitions of prime indices.
A056239 adds up prime indices, row sums of A112798.
A089259 counts set multipartitions of integer partitions.
A116540 counts set multipartitions covering an initial interval by weight.
A368533 lists numbers with squarefree binary indices.
Cf. A000040, A000720, A001222, A005117, A006450, A076610, A270995, A296119, A302242, A302590, A339113.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],OddQ[#]&&And@@SquareFreeQ/@prix[#]&]

cp's OEIS Frontend

A371453 Numbers whose binary indices are all squarefree semiprimes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Python

Python

A371454 Numbers whose binary indices are all semiprimes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

A301753 Number of ways to choose a strict rooted partition of each part in a rooted partition of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A383310 Number of ways to choose a strict multiset partition of a factorization of n into factors > 1.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A383308 Number of integer partitions of n that can be partitioned into sets with a common sum.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A383311 Number of ways to choose a set multipartition (multiset of sets) of a factorization of n into factors > 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A368603 Products of odd primes of squarefree index. MM-numbers of set multipartitions.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica