A367096 -id:A367096 - cp's OEIS frontend

A366740 Positive integers whose semiprime divisors do not all have different Heinz weights (sum of prime indices, A056239).

90, 180, 210, 270, 360, 420, 450, 462, 525, 540, 550, 630, 720, 810, 840, 858, 900, 910, 924, 990, 1050, 1080, 1100, 1155, 1170, 1260, 1326, 1350, 1386, 1440, 1470, 1530, 1575, 1620, 1650, 1666, 1680, 1710, 1716, 1800, 1820, 1848, 1870, 1890, 1911, 1938, 1980

Offset: 1

Gus Wiseman, Nov 05 2023

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.
From Robert Israel, Nov 06 2023: (Start)
Positive integers divisible by the product of four primes, prime(i)*prime(j)*prime(k)*prime(l), i < j <= k < l, with i + l = j + k.
All positive multiples of terms are terms. (End)

			The semiprime divisors of 90 are (6,9,10,15), with prime indices ({1,2},{2,2},{1,3},{2,3}) with sums (3,4,4,5), which are not all different, so 90 is in the sequence.
The terms together with their prime indices begin:
    90: {1,2,2,3}
   180: {1,1,2,2,3}
   210: {1,2,3,4}
   270: {1,2,2,2,3}
   360: {1,1,1,2,2,3}
   420: {1,1,2,3,4}
   450: {1,2,2,3,3}
   462: {1,2,4,5}
   525: {2,3,3,4}
   540: {1,1,2,2,2,3}
   550: {1,3,3,5}
   630: {1,2,2,3,4}
   720: {1,1,1,1,2,2,3}

The complement is too dense.
For all divisors instead of just semiprimes we have A299729, strict A316402.
Distinct semi-sums of prime indices are counted by A366739.
Partitions of this type are counted by A366753, non-binary A366754.
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 distinct subset-sums of prime indices, positive A304793.
A299702 ranks knapsack partitions, counted by A108917, strict A275972.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264
Cf. A000720, A001248, A008967, A365541, A365920, A366737, A366738, A366741, A367093, A367095, A367097.

N:= 10^4: # for terms <= N
P:= select(isprime, [$1..N]): nP:= nops(P):
R:= {}:
for i from 1 while P[i]*P[i+1]^2*P[i+2] < N do
  for j from i+1 while P[i]*P[j]^2 * P[j+1] < N do
    for k from j do
      l:= j+k-i;
      if l <= k or l > nP then break fi;
      v:= P[i]*P[j]*P[k]*P[l];
      if v <= N then
        R:= R union {seq(t,t=v..N,v)};
      fi
od od od:
sort(convert(R,list)); # Robert Israel, Nov 06 2023

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[1000],!UnsameQ@@Total/@Union[Subsets[prix[#],{2}]]&]

These are numbers k such that A086971(k) > A366739(k).

A366739 Number of distinct semi-sums of the multiset of prime indices of n. Number of distinct sums of prime indices of semiprime divisors of n (counted by A086971).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 04 2023

nonn

First differs from A086971 at a(90) = 3, A086971(90) = 4.
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.

			The prime indices of 90 are {1,2,2,3}, with semi-sums
  3 = 1+2
  4 = 1+3 (or 2+2)
  5 = 2+3
so a(90) = 3.
Alternatively, the semiprime divisors of 90 are (6,9,10,15), with prime indices ({1,2},{2,2},{1,3},{2,3}) with sums (3,4,4,5) so a(90) = 3.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences related to prime indices in the factorization of n.

The non-binary version is A299701.
Summing over partitions gives A366738, strict A366741.
For all sums of pairs of elements we have A367095.
Positions of first appearances are A367097.
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.
A299702 ranks knapsack partitions, counted by A108917.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264
Cf. A000005, A000720, A001248, A008967, A117855, A304792, A304793, A365541, A365920, A366740, A366753, A367093.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Union[Total/@Subsets[prix[n],{2}]]],{n,100}]

A366739(n) = #Set(apply(d->((f)->sum(i=1,#f~,f[i,2]*primepi(f[i,1])))(factor(d)), select(d->2==bigomega(d), divisors(n)))); \\ Antti Karttunen, Jan 20 2025

a(n) <= A086971(n). - Antti Karttunen, Jan 20 2025

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

A367095 Number of distinct sums of pairs (repeats allowed) of prime indices of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 3, 1, 3, 3, 1, 1, 3, 1, 3, 3, 3, 1, 3, 1, 3, 1, 3, 1, 5, 1, 1, 3, 3, 3, 3, 1, 3, 3, 3, 1, 6, 1, 3, 3, 3, 1, 3, 1, 3, 3, 3, 1, 3, 3, 3, 3, 3, 1, 5, 1, 3, 3, 1, 3, 6, 1, 3, 3, 6, 1, 3, 1, 3, 3, 3, 3, 6, 1, 3, 1, 3, 1, 6, 3, 3, 3, 3, 1, 5, 3, 3, 3, 3, 3, 3, 1, 3, 3, 3, 1, 6, 1, 3, 5

Offset: 1

Gus Wiseman, Nov 06 2023

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.

Is the image missing only 2 and 4?

			The prime indices of 15 are {2,3}, with sums of pairs:
  2+2 = 4
  2+3 = 5
  3+3 = 6
so a(15) = 3.
The prime indices of 180 are {1,1,2,2,3}, with sums of pairs:
  1+1 = 2
  1+2 = 3
  1+3 = 4
  2+2 = 4
  2+3 = 5
  3+3 = 6
so a(180) = 5.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences related to prime indices in the factorization of n.

Depends only on squarefree kernel A007947. (Even more exactly, on A322591 - Antti Karttunen, Jan 20 2025)
Positions of first appearances appear to be a subset of A325986.
For 2-element submultisets we have A366739, for all submultisets A299701.
A001222 counts prime factors (also indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, conjugate A065119.
A056239 adds up prime indices, row sums of A112798.
A304793 counts positive subset-sums of prime indices.
A367096 lists semiprime divisors, count A086971.
Cf. A001248, A004526, A008967, A366737, A366738, A366740, A367093, A367097.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Union[Total/@Tuples[prix[n],2]]],{n,100}]

A367095(n) = if(1==n, 0, my(pis=apply(primepi,factor(n)[,1]), pairsums = vector(binomial(1+#pis,2)), k=0); for(i=1,#pis,for(j=i,#pis,k++; pairsums[k] = pis[i]+pis[j])); #Set(pairsums)); \\ Antti Karttunen, Jan 20 2025

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

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

Original entry on oeis.org

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}]

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

A367098 Number of divisors of n with exactly two distinct prime factors.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 09 2023

			The a(n) divisors for n = 1, 6, 12, 24, 36, 60, 72, 120, 144, 216, 288, 360:
  .  6  6   6   6   6   6   6   6    6    6    6
        12  12  12  10  12  10  12   12   12   10
            24  18  12  18  12  18   18   18   12
                36  15  24  15  24   24   24   15
                    20  36  20  36   36   36   18
                        72  24  48   54   48   20
                            40  72   72   72   24
                                144  108  96   36
                                     216  144  40
                                          288  45
                                               72

Amiram Eldar, Table of n, a(n) for n = 1..10000

For just one distinct prime factor we have A001222 (prime-power divisors).
This sequence counts divisors belonging to A007774.
Counting all prime factors gives A086971, firsts A220264.
Column k = 2 of A146289.
- Positions of zeros are A000961 (powers of primes), complement A024619.
- Positions of ones are A006881 (squarefree semiprimes).
- Positions of twos are A054753.
- Positions of first appearances are A367099.
A001221 counts distinct prime factors.
A001358 lists semiprimes, complement A100959.
A367096 lists semiprime divisors, sum A076290.
Cf. A000005, A001248, A056170, A079275, A090885, A366740.

Table[Length[Select[Divisors[n], PrimeNu[#]==2&]],{n,100}]
a[1] = 0; a[n_] := (Total[(e = FactorInteger[n][[;; , 2]])]^2 - Total[e^2])/2; Array[a, 100] (* Amiram Eldar, Jan 08 2024 *)

a(n) = {my(e = factor(n)[, 2]); (vecsum(e)^2 - e~*e)/2;} \\ Amiram Eldar, Jan 08 2024

a(n) = (A001222(n)^2 - A090885(n))/2. - Amiram Eldar, Jan 08 2024

A367099 Least positive integer such that the number of divisors having two distinct prime factors is n.

Original entry on oeis.org

1, 6, 12, 24, 36, 60, 72, 120, 144, 216, 288, 360, 432, 960, 720, 864, 1296, 1440, 1728, 2160, 2592, 3456, 7560, 4320, 5184, 7776, 10800, 8640, 10368, 12960, 15552, 17280, 20736, 40320, 25920, 31104, 41472, 60480, 64800, 51840, 62208, 77760, 93312

Offset: 0

Gus Wiseman, Nov 09 2023

nonn

Does this contain every power of six, namely 1, 6, 36, 216, 1296, 7776, ...?

Yes, every power of six is a term, since 6^k = 2^k * 3^k is the least positive integer having n = tau(6^k) - (2k+1) divisors with two distinct prime factors. - Ivan N. Ianakiev, Nov 11 2023

			The divisors of 60 having two distinct prime factors are: 6, 10, 12, 15, 20. Since 60 is the first number having five such divisors, we have a(5) = 60.
The terms together with their prime indices begin:
     1: {}
     6: {1,2}
    12: {1,1,2}
    24: {1,1,1,2}
    36: {1,1,2,2}
    60: {1,1,2,3}
    72: {1,1,1,2,2}
   120: {1,1,1,2,3}
   144: {1,1,1,1,2,2}
   216: {1,1,1,2,2,2}
   288: {1,1,1,1,1,2,2}
   360: {1,1,1,2,2,3}
   432: {1,1,1,1,2,2,2}
   960: {1,1,1,1,1,1,2,3}
   720: {1,1,1,1,2,2,3}
   864: {1,1,1,1,1,2,2,2}

Amiram Eldar, Table of n, a(n) for n = 0..4469

The version for all divisors is A005179 (firsts of A000005).
For all prime factors (A001222) we have A220264, firsts of A086971.
Positions of first appearances in A367098 (counts divisors in A007774).
A000961 lists prime powers, complement A024619.
A001221 counts distinct prime factors.
A001358 lists semiprimes, squarefree A006881, complement A100959.
A367096 lists semiprime divisors, sum A076290.
Cf. A001248, A054753, A056170, A079275, A146289, A366740, A367093.

nn=1000;
w=Table[Length[Select[Divisors[n],PrimeNu[#]==2&]],{n,nn}];
spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
Table[Position[w,k][[1,1]],{k,0,spnm[w]}]

a(n) = my(k=1); while (sumdiv(k, d, omega(d)==2) != n, k++); k; \\ Michel Marcus, Nov 11 2023

cp's OEIS Frontend

A366740 Positive integers whose semiprime divisors do not all have different Heinz weights (sum of prime indices, A056239).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

A366739 Number of distinct semi-sums of the multiset of prime indices of n. Number of distinct sums of prime indices of semiprime divisors of n (counted by A086971).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A367095 Number of distinct sums of pairs (repeats allowed) of prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

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

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

A367098 Number of divisors of n with exactly two distinct prime factors.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula