A086971 -id:A086971 - cp's OEIS frontend

A367093 Least positive integer with n more semiprime divisors than semi-sums of prime indices.

1, 90, 630, 2310, 6930, 34650, 30030, 90090, 450450, 570570, 510510, 1531530, 7657650, 14804790, 11741730, 9699690, 29099070, 145495350

Offset: 0

Gus Wiseman, Nov 05 2023

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.
We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.
Are all primorials after 210 included?

			The terms together with their prime indices begin:
       1: {}
      90: {1,2,2,3}
     630: {1,2,2,3,4}
    2310: {1,2,3,4,5}
    6930: {1,2,2,3,4,5}
   34650: {1,2,2,3,3,4,5}
   30030: {1,2,3,4,5,6}
   90090: {1,2,2,3,4,5,6}
  450450: {1,2,2,3,3,4,5,6}
  570570: {1,2,3,4,5,6,8}
  510510: {1,2,3,4,5,6,7}

The first part (semiprime divisors) is A086971, firsts A220264.
The second part (semi-sums of prime indices) is A366739, firsts A367097.
All sums of pairs of prime indices are counted by A367095.
The non-binary version is A367105.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, conjugate A065119.
A056239 adds up prime indices, row sums of A112798.
A299701 counts subset-sums of prime indices, positive A304793.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264
Cf. A000720, A001248, A008967, A117855, A304792, A365541, A365920, A366738, A366740, A366753.

nn=10000;
w=Table[Length[Union[Subsets[prix[n],{2}]]]-Length[Union[Total/@Subsets[prix[n],{2}]]],{n,nn}];
spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
Table[Position[w,k][[1,1]],{k,0,spnm[w]}]

from itertools import count
from sympy import factorint, primepi
from sympy.utilities.iterables import multiset_combinations
def A367093(n):
    for k in count(1):
        c, a = 0, set()
        for s in (sum(p) for p in multiset_combinations({primepi(i):j for i,j in factorint(k).items()},2)):
            if s not in a:
                a.add(s)
            else:
                c += 1
            if c > n:
                break
        if c == n:
            return k # Chai Wah Wu, Nov 13 2023

a(n) is the least positive integer such that A086971(a(n)) - A366739(a(n)) = n.

a(12)-a(16) from Chai Wah Wu, Nov 13 2023

a(17) from Chai Wah Wu, Nov 18 2023

A366753 Number of integer partitions of n without all different sums of two-element submultisets.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 3, 4, 9, 11, 22, 27, 48, 61, 98, 123, 188, 237, 345, 435, 611, 765, 1046, 1305, 1741, 2165, 2840, 3502, 4527, 5562, 7083, 8650, 10908, 13255, 16545, 20016, 24763, 29834, 36587, 43911, 53514, 63964, 77445, 92239, 111015, 131753

Offset: 0

Gus Wiseman, Nov 07 2023

nonn

			The two-element submultisets of y = {1,1,1,2,2,3} are {1,1}, {1,2}, {1,3}, {2,2}, {2,3}, with sums 2, 3, 4, 4, 5, which are not all different, so y is counted under a(10).
The a(8) = 1 through a(13) = 11 partitions:
  (3221)  (32211)  (4321)    (33221)    (4332)      (43321)
                   (32221)   (43211)    (5331)      (53221)
                   (322111)  (322211)   (5421)      (53311)
                             (3221111)  (43221)     (54211)
                                        (322221)    (332221)
                                        (332211)    (432211)
                                        (432111)    (3222211)
                                        (3222111)   (3322111)
                                        (32211111)  (4321111)
                                                    (32221111)
                                                    (322111111)

Semiprime divisors are counted by A086971, distinct sums A366739.
The non-binary complement is A108917, strict A275972, ranks A299702.
These partitions have ranks A366740.
The non-binary version is A366754, strict A316402, ranks A299729.
A276024 counts positive subset-sums of partitions, strict A284640.
A304792 counts subset-sum of partitions, strict A365925.
A365543 counts partitions with a subset-sum k, complement A046663.
A365661 counts strict partitions with a subset-sum k, complement A365663.
A366738 counts semi-sums of partitions, strict A366741.
A367096 lists semiprime divisors, row sums A076290.
Cf. A002033, A001358, A006827, A008967, A122768, A126796, A365923, A365924, A367093, A367095.

Table[Length[Select[IntegerPartitions[n],!UnsameQ@@Total/@Union[Subsets[#,{2}]]&]],{n,0,30}]

A367096 Irregular triangle read by rows where row n lists the semiprime divisors of n. Alternatively, row n lists the semiprime divisors of A002808(n).

Original entry on oeis.org

4, 6, 4, 9, 10, 4, 6, 14, 15, 4, 6, 9, 4, 10, 21, 22, 4, 6, 25, 26, 9, 4, 14, 6, 10, 15, 4, 33, 34, 35, 4, 6, 9, 38, 39, 4, 10, 6, 14, 21, 4, 22, 9, 15, 46, 4, 6, 49, 10, 25, 51, 4, 26, 6, 9, 55, 4, 14, 57, 58, 4, 6, 10, 15, 62, 9, 21, 4, 65, 6, 22, 33, 4, 34

Offset: 1

Gus Wiseman, Nov 08 2023

On the first interpretation, the first three rows are empty. On the second, the first row is (4).

			The semiprime divisors of 30 are {6,10,15}, so row 30 is (6,10,15). Without empty rows, this is row 19.
Triangle begins (empty rows indicated by dots):
   1: .
   2: .
   3: .
   4: 4
   5: .
   6: 6
   7: .
   8: 4
   9: 9
  10: 10
  11: .
  12: 4,6
Without empty rows:
   1: 4
   2: 6
   3: 4
   4: 9
   5: 10
   6: 4,6
   7: 14
   8: 15
   9: 4
  10: 6,9
  11: 4,10
  12: 21

For all divisors we have A027750.
Square terms are counted by A056170.
Row sums are A076290.
Squarefree terms are counted by A079275.
Row lengths are A086971, firsts A220264.
A000005 counts divisors.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, complement A100959.
Cf. A008967, A365829, A366738, A366739, A366740, A367093, A367098.

Table[Select[Divisors[n],PrimeOmega[#]==2&],{n,100}]

row(n) = select(x -> bigomega(x) == 2, divisors(n)); \\ Amiram Eldar, May 02 2025

A346041 Numbers with exactly 1 semiprime divisor.

Original entry on oeis.org

4, 6, 8, 9, 10, 14, 15, 16, 21, 22, 25, 26, 27, 32, 33, 34, 35, 38, 39, 46, 49, 51, 55, 57, 58, 62, 64, 65, 69, 74, 77, 81, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 121, 122, 123, 125, 128, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 169

Offset: 1

Wesley Ivan Hurt, Jul 02 2021

nonn

Numbers of the form p*q or p^k, where p and q are prime and k >= 2.

Numbers k such that A086971(k) = 1. - Wesley Ivan Hurt, Jun 21 2024

			6 is in the sequence since it has exactly 1 semiprime divisor, namely 6.
16 is in the sequence since it has exactly 1 semiprime divisor, namely 4.

Cf. A001358 (semiprimes), A086971.

Select[Range@200,Length@Select[Divisors@#,PrimeOmega@#==2&]==1&] (* Giorgos Kalogeropoulos, Jul 03 2021 *)

isok(k) = sumdiv(k, d, bigomega(d)==2) == 1; \\ Michel Marcus, Jul 03 2021

from sympy import factorint
def ok(n):
    f = factorint(n); w = len(f); W = sum(f.values())
    return (w == 1 and W >= 2) or (w == 2 and W == 2)
print(list(filter(ok, range(170)))) # Michael S. Branicky, Jul 03 2021

A367097 Least positive integer whose multiset of prime indices has exactly n distinct semi-sums.

Original entry on oeis.org

1, 4, 12, 30, 60, 210, 330, 660, 2730, 3570, 6270, 12540, 53130, 79170, 110670, 221340, 514140, 1799490, 2284590, 4196010, 6750870, 13501740, 37532220, 97350330, 131362770, 189620970, 379241940, 735844830, 1471689660

Offset: 0

Gus Wiseman, Nov 09 2023

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.
We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.
From David A. Corneth, Nov 15 2023: (Start)
Terms are cubefree.
bigomega(a(n)) = A001222(a(n)) >= A002024(n) + 1 = floor(sqrt(2n) + 1/2) + 1 for n > 0. (End)

			The prime indices of 60 are {1,1,2,3}, with four semi-sums {2,3,4,5}, and 60 is the first number whose prime indices have four semi-sums, so a(4) = 60.
The terms together with their prime indices begin:
       1: {}
       4: {1,1}
      12: {1,1,2}
      30: {1,2,3}
      60: {1,1,2,3}
     210: {1,2,3,4}
     330: {1,2,3,5}
     660: {1,1,2,3,5}
    2730: {1,2,3,4,6}
    3570: {1,2,3,4,7}
    6270: {1,2,3,5,8}
   12540: {1,1,2,3,5,8}
   53130: {1,2,3,4,5,9}
   79170: {1,2,3,4,6,10}
  110670: {1,2,3,4,7,11}
  221340: {1,1,2,3,4,7,11}
  514140: {1,1,2,3,5,8,13}

The non-binary version is A259941, firsts of A299701.
These are the positions of first appearances in A366739.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, complement A100959.
A056239 adds up prime indices, row sums of A112798.
A299702 ranks knapsack partitions, counted by A108917.
A366738 counts semi-sums of partitions, strict A366741.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264
Cf. A000720, A001248, A002024, A004709, A304793, A365920, A366740, A366753, A367093, A367095.

nn=1000;
w=Table[Length[Union[Total/@Subsets[prix[n],{2}]]],{n,nn}];
spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
v=Table[Position[w,k][[1,1]],{k,0,spnm[w]}]

from itertools import count
from sympy import factorint, primepi
from sympy.utilities.iterables import multiset_combinations
def A367097(n): return next(k for k in count(1) if len({sum(s) for s in multiset_combinations({primepi(i):j for i,j in factorint(k).items()},2)}) == n) # Chai Wah Wu, Nov 13 2023

2 | a(n) for n > 0. - David A. Corneth, Nov 13 2023

a(17)-a(22) from Chai Wah Wu, Nov 13 2023

a(23)-a(28) from David A. Corneth, Nov 13 2023

A102466 Numbers such that the number of divisors is the sum of numbers of prime factors with and without repetitions.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 31, 32, 33, 34, 35, 37, 38, 39, 41, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 67, 69, 71, 73, 74, 77, 79, 81, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109

Offset: 1

Reinhard Zumkeller, Jan 09 2005

nonn

A000005(a(n)) = A001221(a(n)) + A001222(a(n)); prime powers are a subsequence (A000961); complement of A102467; not the same as A085156.
Equals { n | omega(n)=1 or Omega(n)=2 }, that is, these are exactly the prime powers (>1) and semiprimes. - M. F. Hasler, Jan 14 2008
For n > 1: A086971(a(n)) <= 1. - Reinhard Zumkeller, Dec 14 2012

T. D. Noe, Table of n, a(n) for n = 1..1000

Cf. A000005, A001221, A001222, A000961, A102467, A085156, A086971, A135767.

a102466 n = a102466_list !! (n-1)
a102466_list = [x | x <- [1..], a000005 x == a001221 x + a001222 x]
-- Reinhard Zumkeller, Dec 14 2012

with(numtheory):
q:= n-> is(tau(n)=bigomega(n)+nops(factorset(n))):
select(q, [$1..200])[];  # Alois P. Heinz, Jul 14 2023

Select[Range[110],DivisorSigma[0,#]==PrimeOmega[#]+PrimeNu[#]&] (* Harvey P. Dale, Mar 09 2016 *)

is(n)=my(f=factor(n)[,2]); #f==1 || f==[1,1]~ \\ Charles R Greathouse IV, Oct 19 2015

def is_A102466(n) :
    return bool(sloane.A001221(n) == 1 or sloane.A001222(n) == 2)
def A102466_list(n) :
    return [k for k in (1..n) if is_A102466(k)]
A102466_list(109)  # Peter Luschny, Feb 08 2012

A106404 Number of even semiprimes dividing n.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 2, 0, 1, 0, 2, 0, 2, 0, 1, 0, 1, 0, 2, 0, 1, 0, 2, 0, 2, 0, 2, 0, 1, 0, 2, 0, 1, 0, 2, 0, 1, 0, 2, 0, 1, 0, 3, 0, 1, 0, 1, 0, 2, 0, 2, 0, 2, 0, 2, 0, 1, 0, 2, 0, 2, 0, 2, 0, 1, 0, 3, 0, 1, 0, 2, 0, 2, 0, 2, 0, 1, 0, 2, 0, 1, 0, 2, 0, 2, 0, 2, 0

Offset: 1

Reinhard Zumkeller, May 02 2005

nonn

Also the number of prime divisors p|n such that n/p is even. - Gus Wiseman, Jun 06 2018

			a(60) = #{4, 6, 10} = #{2*2, 2*3, 2*5} = 3.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A000005, A000607, A001221, A001227, A008683, A083399, A088705, A100484, A205745, A305614.

a106404 n = length [d | d <- takeWhile (<= n) a100484_list, mod n d == 0]
-- Reinhard Zumkeller, Jan 31 2012

Table[Length[Select[Divisors[n],PrimeQ[#]&&EvenQ[n/#]&]],{n,100}] (* Gus Wiseman, Jun 06 2018 *)
Table[Count[Divisors[n],?(EvenQ[#]&&PrimeOmega[#]==2&)],{n,110}] (* _Harvey P. Dale, May 04 2021 *)
a[n_] := If[EvenQ[n], PrimeNu[n/2], 0]; Array[a, 100] (* Amiram Eldar, Jun 26 2022 *)

a(n)=if(n%2,0,omega(n/2)) \\ Charles R Greathouse IV, Jan 30 2012

def A106404(n):
    return add(1-(n/d)%2 for d in divisors(n) if is_prime(d))
print([A106404(n) for n in (1..105)])  # Peter Luschny, Jan 30 2012

a(n) = A086971(n) - A106405(n).
a(A100484(n)) = 1.
a(A005408(n)) = 0.
a(A005843(n)) > 0 for n>1.
a(2n) = omega(n), a(2n+1) = 0, where omega(n) is the number of distinct prime divisors of n, A001221. - Franklin T. Adams-Watters, Jun 09 2006
a(n) = card { d | d*p = n, d even, p prime }. - Peter Luschny, Jan 30 2012
O.g.f.: Sum_{p prime} x^(2p)/(1 - x^(2p)). - Gus Wiseman, Jun 06 2018

A367402 Number of integer partitions of n whose semi-sums cover an interval of positive integers.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 9, 10, 13, 17, 20, 26, 31, 38, 44, 58, 64, 81, 95, 116, 137, 166, 192, 233, 278, 330, 385, 459, 542, 636, 759, 879, 1038, 1211, 1418, 1656, 1942, 2242, 2618, 3029, 3535, 4060, 4735, 5429, 6299, 7231, 8346, 9556, 11031, 12593, 14482, 16525

Offset: 0

Gus Wiseman, Nov 17 2023

nonn

We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.

			The partition y = (3,2,1,1) has semi-sums {2,3,4,5}, which is an interval, so y is counted under a(7).
The a(1) = 1 through a(8) = 13 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (21)   (22)    (32)     (33)      (43)       (44)
             (111)  (31)    (41)     (42)      (52)       (53)
                    (211)   (221)    (51)      (61)       (62)
                    (1111)  (2111)   (222)     (322)      (71)
                            (11111)  (321)     (2221)     (332)
                                     (2211)    (3211)     (2222)
                                     (21111)   (22111)    (3221)
                                     (111111)  (211111)   (22211)
                                               (1111111)  (32111)
                                                          (221111)
                                                          (2111111)
                                                          (11111111)

For parts instead of sums we have A034296, ranks A073491.
For all subset-sums we have A126796, ranks A325781, strict A188431.
The complement for parts instead of sums is A239955, ranks A073492.
The complement for all sub-sums is A365924, ranks A365830, strict A365831.
The complement is counted by A367403.
The strict case is A367410, complement A367411.
A000009 counts partitions covering an initial interval, ranks A055932.
A086971 counts semi-sums of prime indices.
A261036 counts complete partitions by maximum.
A276024 counts positive subset-sums of partitions, strict A284640.
Cf. A000041, A002033, A046663, A108917, A264401, A304792, A365543, A365661, A365918, A365921.

Table[Length[Select[IntegerPartitions[n], (d=Total/@Subsets[#,{2}];If[d=={}, {}, Range[Min@@d,Max@@d]]==Union[d])&]], {n,0,15}]

A367403 Number of integer partitions of n whose semi-sums do not cover an interval of positive integers.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 2, 5, 9, 13, 22, 30, 46, 63, 91, 118, 167, 216, 290, 374, 490, 626, 810, 1022, 1297, 1628, 2051, 2551, 3176, 3929, 4845, 5963, 7311, 8932, 10892, 13227, 16035, 19395, 23397, 28156, 33803, 40523, 48439, 57832, 68876, 81903, 97212, 115198

Offset: 0

Gus Wiseman, Nov 17 2023

nonn

We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.

			The a(0) = 0 through a(9) = 13 partitions:
  .  .  .  .  .  (311)  (411)   (331)    (422)     (441)
                        (3111)  (421)    (431)     (522)
                                (511)    (521)     (531)
                                (4111)   (611)     (621)
                                (31111)  (3311)    (711)
                                         (4211)    (4311)
                                         (5111)    (5211)
                                         (41111)   (6111)
                                         (311111)  (33111)
                                                   (42111)
                                                   (51111)
                                                   (411111)
                                                   (3111111)

The complement for parts instead of sums is A034296, ranks A073491.
The complement for all sub-sums is A126796, ranks A325781, strict A188431.
For parts instead of sums we have A239955, ranks A073492.
For all subset-sums we have A365924, ranks A365830, strict A365831.
The complement is counted by A367402.
The strict case is A367411, complement A367410.
A000009 counts partitions covering an initial interval, ranks A055932.
A086971 counts semi-sums of prime indices.
A261036 counts complete partitions by maximum.
A276024 counts positive subset-sums of partitions, strict A284640.
Cf. A000041, A002033, A046663, A108917, A264401, A304792, A365543, A365658.

Table[Length[Select[IntegerPartitions[n], (d=Total/@Subsets[#,{2}];If[d=={}, {}, Range[Min@@d,Max@@d]]!=Union[d])&]], {n,0,15}]

A106405 Number of odd semiprimes dividing n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 2, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 2, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 2, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 2, 1, 0, 1, 0, 1, 0, 0, 1, 2, 1, 0, 1, 0, 0, 3

Offset: 1

Reinhard Zumkeller, May 02 2005

nonn

			a(105) = #{15, 21, 35} = #{3*5, 3*7, 5*7} = 3.

Harvey P. Dale, Table of n, a(n) for n = 1..1000

Cf. A000005, A001227, A046315, A086971, A093641, A105441, A106404.

Table[Count[Divisors[n],?(OddQ[#]&&PrimeOmega[#]==2&)],{n,120}] (* _Harvey P. Dale, May 05 2015 *)
a[n_] := Count[e = FactorInteger[n][[;; , 2]], ?(# > 1 &)] + (o = Length[e])*(o - 1)/2 - If[EvenQ[n], If[e[[1]] > 1, o, o - 1], 0]; Array[a, 100] (* _Amiram Eldar, Jun 30 2022 *)

a(n) = A086971(n) - A106404(n);

a(A046315(n)) = 1; a(A093641(n)) = 0; a(A105441(n)) > 0.

cp's OEIS Frontend

A367093 Least positive integer with n more semiprime divisors than semi-sums of prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A366753 Number of integer partitions of n without all different sums of two-element submultisets.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A367096 Irregular triangle read by rows where row n lists the semiprime divisors of n. Alternatively, row n lists the semiprime divisors of A002808(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

A346041 Numbers with exactly 1 semiprime divisor.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

A367097 Least positive integer whose multiset of prime indices has exactly n distinct semi-sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A102466 Numbers such that the number of divisors is the sum of numbers of prime factors with and without repetitions.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Sage

A106404 Number of even semiprimes dividing n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Sage