A320655 -id:A320655 - cp's OEIS frontend

A176506 Difference between the prime indices of the two factors of the n-th semiprime.

0, 1, 0, 2, 3, 1, 2, 4, 0, 5, 3, 6, 1, 7, 4, 8, 0, 5, 2, 6, 9, 10, 3, 7, 11, 1, 12, 4, 13, 8, 2, 9, 14, 5, 15, 10, 6, 16, 3, 0, 17, 11, 12, 4, 18, 13, 19, 1, 7, 20, 8, 21, 14, 5, 22, 0, 15, 23, 16, 9, 2, 24, 17, 25, 6, 10, 26, 3, 18, 27, 11, 7, 28, 19, 1, 29, 12, 20, 2, 21, 4, 30, 8, 31, 13, 22

Offset: 1

Juri-Stepan Gerasimov, Apr 19 2010

Are there no adjacent equal terms? I have verified this up to n = 10^6. - Gus Wiseman, Dec 04 2020

			From _Gus Wiseman_, Dec 04 2020: (Start)
The sequence of semiprimes together with the corresponding differences begins:
   4: 1 - 1 = 0
   6: 2 - 1 = 1
   9: 2 - 2 = 0
  10: 3 - 1 = 2
  14: 4 - 1 = 3
  15: 3 - 2 = 1
  21: 4 - 2 = 2
  22: 5 - 1 = 4
  25: 3 - 3 = 0
  26: 6 - 1 = 5
  33: 5 - 2 = 3
(End)

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A109313.
Cf. A001358, A049084, A084126, A084127.
A087794 is product of the same indices.
A176504 is the sum of the same indices.
A115392 lists positions of first appearances.
A128301 lists positions of 0's.
A172348 lists positions of 1's.
A338898 has this sequence as row differences.
A338900 is the squarefree case.
A338912/A338913 give the two prime indices of semiprimes.
A006881 lists squarefree semiprimes.
A024697 is the sum of semiprimes of weight n.
A056239 gives sum of prime indices (Heinz weight).
A087112 groups semiprimes by greater factor.
A270650/A270652/A338899 give the prime indices of squarefree semiprimes.
A338904 groups semiprimes by weight.
A338907/A338906 list semiprimes of odd/even weight.
A339114/A339115 give the least/greatest semiprime of weight n.
Cf. A000040, A001222, A062198, A112141, A112798, A325699, A320655, A339005.

isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:
A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:
A084126 := proc(n) min(op(numtheory[factorset](A001358(n)))) ; end proc:
A084127 := proc(n) max(op(numtheory[factorset](A001358(n)))) ; end proc:
A176506 := proc(n) numtheory[pi](A084127(n)) - numtheory[pi](A084126(n)) ; end proc: seq(A176506(n),n=1..120) ; # R. J. Mathar, Apr 22 2010
# Alternative:
N:= 500: # to use the first N semiprimes
Primes:= select(isprime, [2,seq(i,i=3..N/2,2)]):
SP:= NULL:
for i from 1 to nops(Primes) do
  for j from 1 to i do
    sp:= Primes[i]*Primes[j];
    if sp > N then break fi;
    SP:= SP, [sp, i-j]
od od:
SP:= sort([SP],(s,t) -> s[1] t[2], SP); # Robert Israel, Jan 17 2019

M = 500; (* to use the first M semiprimes *)
primes = Select[Join[{2}, Range[3, M/2, 2]], PrimeQ];
SP = {};
For[i = 1, i <= Length[primes], i++,
  For[j = 1, j <= i, j++,
    sp = primes[[i]] primes[[j]];
    If[sp > M, Break []];
    AppendTo[SP, {sp, i - j}]
]];
SortBy[SP, First][[All, 2]] (* Jean-François Alcover, Jul 18 2020, after Robert Israel *)
Table[If[!SquareFreeQ[n],0,-Subtract@@PrimePi/@First/@FactorInteger[n]],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

lista(nn) = {my(vsp = select(x->(bigomega(x)==2), [1..nn])); vector(#vsp, k, my(f=factor(vsp[k])[,1]); primepi(vecmax(f)) - primepi(vecmin(f)));} \\ Michel Marcus, Jul 18 2020

a(n) = A049084(A084127(n)) - A049084(A084126(n)). [corrected by R. J. Mathar, Apr 22 2010]

a(n) = A338913(n) - A338912(n). - Gus Wiseman, Dec 04 2020

a(51) and a(69) corrected by R. J. Mathar, Apr 22 2010

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

Gus Wiseman, Oct 23 2018

nonn

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

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

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

Cf. A001055, A001222, A001358, A007717, A028260, A320655, A320892, A320893, A320894, A320911.

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[#]!={}]&]

A087794 Products of prime-indices of factors of semiprimes.

Original entry on oeis.org

1, 2, 4, 3, 4, 6, 8, 5, 9, 6, 10, 7, 12, 8, 12, 9, 16, 14, 15, 16, 10, 11, 18, 18, 12, 20, 13, 21, 14, 20, 24, 22, 15, 24, 16, 24, 27, 17, 28, 25, 18, 26, 28, 32, 19, 30, 20, 30, 30, 21, 33, 22, 32, 36, 23, 36, 34, 24, 36, 36, 35, 25, 38, 26, 40, 39, 27, 40, 40, 28, 42, 44, 29

Offset: 1

Reinhard Zumkeller, Oct 09 2003

nonn

A semiprime (A001358) is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798. - Gus Wiseman, Dec 04 2020

			A001358(20)=57=3*19=A000040(2)*A000040(8), therefore a(20)=2*8=16.
From _Gus Wiseman_, Dec 04 2020: (Start)
The sequence of all semiprimes together with the products of their prime indices begins:
   4: 1 * 1 = 1
   6: 1 * 2 = 2
   9: 2 * 2 = 4
  10: 1 * 3 = 3
  14: 1 * 4 = 4
  15: 2 * 3 = 6
  21: 2 * 4 = 8
  22: 1 * 5 = 5
  25: 3 * 3 = 9
  26: 1 * 6 = 6
(End)

A003963 is the version for not just semiprimes.
A176504 gives the sum of the same two indices.
A176506 gives the difference of the same two indices.
A339361 is the squarefree case.
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A289182/A115392 list the positions of odd/even terms of A001358.
A338898/A338912/A338913 give the prime indices of semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
A338904 groups semiprimes by weight.
Cf. A001222, A046315, A056239, A065516, A084126, A084127, A112798, A128301, A300912, A318990, A320655, A338909, A339362.

Table[If[SquareFreeQ[n],Times@@PrimePi/@First/@FactorInteger[n],PrimePi[Sqrt[n]]^2],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

a(n) = A049084(A084126(n))*A049084(A084127(n)) = A003963(A001358(n)).

a(n) = A003963(A001358(n)) = A338912(n) * A338913(n). - Gus Wiseman, Dec 04 2020

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

Gus Wiseman, Oct 18 2018

nonn

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

Andrew Howroyd, Table of n, a(n) for n = 0..50

Cf. A001055, A001222, A001358, A005117, A006881, A007716, A007717, A037143, A320462, A320655, A320656, A320664, A320665.

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

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

Gus Wiseman, Oct 20 2018

nonn

			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).

Cf. A001055, A001222, A001358, A007716, A007717, A037143, A045778, A320663, A320655.

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

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

Gus Wiseman, Oct 23 2018

nonn

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.

			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.

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

Cf. A001055, A001358, A005117, A006881, A007717, A028260, A318871, A318953, A320655, A320656, A320892, A320893, A320894.

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[#]=={}]&]

A320894 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into distinct squarefree semiprimes.

Original entry on oeis.org

4, 9, 16, 24, 25, 36, 40, 49, 54, 56, 64, 81, 88, 96, 100, 104, 121, 135, 136, 144, 152, 160, 169, 184, 189, 196, 216, 224, 225, 232, 240, 248, 250, 256, 289, 296, 297, 324, 328, 336, 344, 351, 352, 360, 361, 375, 376, 384, 400, 416, 424, 441, 459, 472, 484

Offset: 1

Gus Wiseman, Oct 23 2018

nonn

A squarefree semiprime (A006881) is a product of any two distinct primes.

			A complete list of all strict factorizations of 24 is: (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.

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

Cf. A001055, A001358, A005117, A006881, A007717, A028260, A318871, A318953, A320655, A320656, A320891, A320892, A320893.

strsqfsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsqfsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
Select[Range[100],And[EvenQ[PrimeOmega[#]],strsqfsemfacs[#]=={}]&]

A338907 Semiprimes whose prime indices sum to an odd number.

Original entry on oeis.org

6, 14, 15, 26, 33, 35, 38, 51, 58, 65, 69, 74, 77, 86, 93, 95, 106, 119, 122, 123, 141, 142, 143, 145, 158, 161, 177, 178, 185, 201, 202, 209, 214, 215, 217, 219, 221, 226, 249, 262, 265, 278, 287, 291, 299, 302, 305, 309, 319, 323, 326, 327, 329, 346, 355

Offset: 1

Gus Wiseman, Nov 28 2020

nonn

All terms are squarefree (A005117).
A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.
The semiprimes in A300063; the semiprimes in A332820. - Peter Munn, Dec 25 2020

			The sequence of terms together with their prime indices begins:
      6: {1,2}      95: {3,8}     202: {1,26}
     14: {1,4}     106: {1,16}    209: {5,8}
     15: {2,3}     119: {4,7}     214: {1,28}
     26: {1,6}     122: {1,18}    215: {3,14}
     33: {2,5}     123: {2,13}    217: {4,11}
     35: {3,4}     141: {2,15}    219: {2,21}
     38: {1,8}     142: {1,20}    221: {6,7}
     51: {2,7}     143: {5,6}     226: {1,30}
     58: {1,10}    145: {3,10}    249: {2,23}
     65: {3,6}     158: {1,22}    262: {1,32}
     69: {2,9}     161: {4,9}     265: {3,16}
     74: {1,12}    177: {2,17}    278: {1,34}
     77: {4,5}     178: {1,24}    287: {4,13}
     86: {1,14}    185: {3,12}    291: {2,25}
     93: {2,11}    201: {2,19}    299: {6,9}

A031368 looks at primes instead of semiprimes.
A098350 has this as union of odd-indexed antidiagonals.
A300063 looks at all numbers (not just semiprimes).
A338904 has this as union of odd-indexed rows.
A338906 is the even version.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A056239 gives the sum of prime indices (Heinz weight).
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A289182/A115392 list the positions of odd/even terms in A001358.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338908 lists squarefree semiprimes of even weight.
A339114/A339115 give the least/greatest semiprime of weight n.
Cf. A000040, A001222, A014342, A024697, A062198, A112798, A300061, A319242, A320655, A338910, A339003.
Subsequence of A332820.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],PrimeOmega[#]==2&&OddQ[Total[primeMS[#]]]&]

from math import isqrt
from sympy import primepi, primerange
def A338907(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum((primepi(x//p)-a>>1) for a,p in enumerate(primerange(isqrt(x)+1)))
    return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

Complement of A338906 in A001358.

A174725 a(n) = (A074206(n) + A008683(n))/2.

Original entry on oeis.org

1, 0, 0, 1, 0, 2, 0, 2, 1, 2, 0, 4, 0, 2, 2, 4, 0, 4, 0, 4, 2, 2, 0, 10, 1, 2, 2, 4, 0, 6, 0, 8, 2, 2, 2, 13, 0, 2, 2, 10, 0, 6, 0, 4, 4, 2, 0, 24, 1, 4, 2, 4, 0, 10, 2, 10, 2, 2, 0, 22, 0, 2, 4, 16, 2, 6, 0, 4, 2, 6, 0, 38, 0, 2, 4, 4, 2

Offset: 1

Mats Granvik, Mar 28 2010

nonn

From Mats Granvik, May 25 2017: (Start)
A074206(n) = A002033(n-1) = a(n) + A174726(n).
A008683(n) = a(n) - A174726(n).
Let m = size of matrix a matrix T, and let T be defined as follows:
T(n,k) = if m = 1 then 1 else if mod(n, k) = 0 then if and(n = k, n = m) then 0 else 1 else if and(n = 1, k = m) then 1 else 0
a(n) is then the number of permutation matrices with a positive contribution in the determinant of matrix T. The determinant of T is equal to the Möbius function A008683, see Mathematica program below for how to compute the determinant.
A174726 is the number of permutation matrices with a negative contribution in the determinant of matrix T.
(End)
From Gus Wiseman, Jan 04 2021: (Start)
Also the number of ordered factorizations of n into an even number of factors > 1. The non-ordered case is A339846. For example, the a(n) factorizations for n = 12, 24, 30, 32, 36 are:
  (2*6)  (3*8)      (5*6)   (4*8)      (4*9)
  (3*4)  (4*6)      (6*5)   (8*4)      (6*6)
  (4*3)  (6*4)      (10*3)  (16*2)     (9*4)
  (6*2)  (8*3)      (15*2)  (2*16)     (12*3)
         (12*2)     (2*15)  (2*2*2*4)  (18*2)
         (2*12)     (3*10)  (2*2*4*2)  (2*18)
         (2*2*2*3)          (2*4*2*2)  (3*12)
         (2*2*3*2)          (4*2*2*2)  (2*2*3*3)
         (2*3*2*2)                     (2*3*2*3)
         (3*2*2*2)                     (2*3*3*2)
                                       (3*2*2*3)
                                       (3*2*3*2)
                                       (3*3*2*2)
(End)

Mats Granvik, Table of n, a(n) for n = 1..10000

Cf. A008683, A051731.
The odd version is A174726.
The unordered version is A339846.
A001055 counts factorizations, with strict case A045778.
A058696 counts partitions of even numbers, ranked by A300061.
A074206 counts ordered factorizations, with strict case A254578.
A251683 counts ordered factorizations by product and length.
Other cases of even length:
- A024430 counts set partitions of even length.
- A027187 counts partitions of even length.
- A034008 counts compositions of even length.
- A052841 counts ordered set partitions of even length.
- A067661 counts strict partitions of even length.
- A332305 counts strict compositions of even length
Cf. A002033, A027193, A028260, A050320, A058695, A236913, A316439, A320655, A320656, A339890, A378216 (Dirichlet inverse).

(* From Mats Granvik, May 25 2017: (Start) *)
Clear[t, nn]; nn = 77; t[1, 1] = 1; t[n_, k_] := t[n, k] = If[k == 1, Sum[t[n, k + i], {i, 1, n - 1}], If[Mod[n, k] == 0, t[n/k, 1], 0], 0]; Monitor[Table[Sum[If[Mod[n, k] == 0, MoebiusMu[k]*t[n/k, 1], 0], {k, 1, 77}], {n, 1, nn}], n]
(* The Möbius function as a determinant *) Table[Det[Table[Table[If[m == 1, 1, If[Mod[n, k] == 0, If[And[n == k, n == m], 0, 1], If[And[n == 1, k == m], 1, 0]]], {k, 1, m}], {n, 1, m}]], {m, 1, 42}]
(* (End) *)
ordfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@ordfacs[n/d],{d,Rest[Divisors[n]]}]];
Table[Length[Select[ordfacs[n],EvenQ@*Length]],{n,100}] (* Gus Wiseman, Jan 04 2021 *)

a(n) = (Mobius transform of a(n)) + (Mobius transform of A174726). - Mats Granvik, Apr 04 2010
From Mats Granvik, May 25 2017: (Start)
This sequence is the Moebius transform of A074206.
a(n) = (A074206(n) + A008683(n))/2.
(End)
G.f. A(x) satisfies: A(x) = x + Sum_{i>=2} Sum_{j>=2} A(x^(i*j)). - Ilya Gutkovskiy, May 11 2019

References to A002033(n-1) changed to A074206(n) by Antti Karttunen, Nov 23 2024

A339560 Number of integer partitions of n that can be partitioned into distinct pairs of distinct parts, i.e., into a set of edges.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 2, 4, 5, 8, 8, 13, 17, 22, 28, 39, 48, 62, 81, 101, 127, 167, 202, 253, 318, 395, 486, 608, 736, 906, 1113, 1353, 1637, 2011, 2409, 2922, 3510, 4227, 5060, 6089, 7242, 8661, 10306, 12251, 14503, 17236, 20345, 24045, 28334, 33374, 39223, 46076

Offset: 0

Gus Wiseman, Dec 10 2020

nonn

Naturally, such a partition must have an even number of parts. Its multiplicities form a graphical partition (A000569, A320922), and vice versa.

			The a(3) = 1 through a(11) = 13 partitions (A = 10):
  (21)  (31)  (32)  (42)  (43)    (53)    (54)    (64)    (65)
              (41)  (51)  (52)    (62)    (63)    (73)    (74)
                          (61)    (71)    (72)    (82)    (83)
                          (3211)  (3221)  (81)    (91)    (92)
                                  (4211)  (3321)  (4321)  (A1)
                                          (4221)  (5221)  (4322)
                                          (4311)  (5311)  (4331)
                                          (5211)  (6211)  (4421)
                                                          (5321)
                                                          (5411)
                                                          (6221)
                                                          (6311)
                                                          (7211)
For example, the partition y = (4,3,3,2,1,1) can be partitioned into a set of edges in two ways:
  {{1,2},{1,3},{3,4}}
  {{1,3},{1,4},{2,3}},
so y is counted under a(14).

Eric Weisstein's World of Mathematics, Graphical partition.

A338916 allows equal pairs (x,x).
A339559 counts the complement in even-length partitions.
A339561 gives the Heinz numbers of these partitions.
A339619 counts factorizations of the same type.
A000070 counts non-multigraphical partitions of 2n, ranked by A339620.
A000569 counts graphical partitions, ranked by A320922.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A058696 counts partitions of even numbers, ranked by A300061.
A209816 counts multigraphical partitions, ranked by A320924.
A320655 counts factorizations into semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339617 counts non-graphical partitions of 2n, ranked by A339618.
A339655 counts non-loop-graphical partitions of 2n, ranked by A339657.
A339656 counts loop-graphical partitions, ranked by A339658.
A339659 counts graphical partitions of 2n into k parts.
The following count partitions of even length and give their Heinz numbers:
- A027187 has no additional conditions (A028260).
- 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).
Cf. A001055, A001221, A005117, A007717, A025065, A030229, A320893, A338899, A338903, A339564.

strs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
Table[Length[Select[IntegerPartitions[n],strs[Times@@Prime/@#]!={}&]],{n,0,15}]

A027187(n) = a(n) + A339559(n).

More terms from Jinyuan Wang, Feb 14 2025

cp's OEIS Frontend

A176506 Difference between the prime indices of the two factors of the n-th semiprime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A087794 Products of prime-indices of factors of semiprimes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

Extensions

A320732 Number of factorizations of n into primes or semiprimes.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A320891 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into squarefree semiprimes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A320894 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into distinct squarefree semiprimes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A338907 Semiprimes whose prime indices sum to an odd number.

Views

Author

Keywords

Comments

Examples

Crossrefs