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.

Showing 1-10 of 10 results.

A338900 Difference between the two prime indices of the n-th squarefree semiprime.

Original entry on oeis.org

1, 2, 3, 1, 2, 4, 5, 3, 6, 1, 7, 4, 8, 5, 2, 6, 9, 10, 3, 7, 11, 1, 12, 4, 13, 8, 2, 9, 14, 5, 15, 10, 6, 16, 3, 17, 11, 12, 4, 18, 13, 19, 1, 7, 20, 8, 21, 14, 5, 22, 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
Offset: 1

Views

Author

Gus Wiseman, Nov 16 2020

Keywords

Comments

A squarefree semiprime is a product of any two distinct 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.
Is this sequence an anti-run, i.e., are there no adjacent equal parts? I have verified this conjecture up to n = 10^6. - Gus Wiseman, Nov 18 2020

Crossrefs

A176506 is the not necessarily squarefree version.
A338899 has row-differences equal to this sequence.
A338901 gives positions of first appearances.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes.
A002100 and A338903 count partitions using squarefree semiprimes.
A004526 counts 2-part partitions, with strict case A140106 (shifted left).
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odds A046388 and evens A100484.
A065516 gives first differences of semiprimes.
A166237 gives first differences of squarefree semiprimes.
A270650 and A270652 give the prime indices of squarefree semiprimes.
A338912 and A338913 give the prime indices of semiprimes.

Programs

  • Mathematica
    -Subtract@@PrimePi/@First/@FactorInteger[#]&/@Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&]

Formula

If the n-th squarefree semiprime is prime(x) * prime(y) with x < y, then a(n) = y - x.
a(n) = A270652(n) - A270650(n).

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

Views

Author

Gus Wiseman, Dec 10 2020

Keywords

Comments

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

Examples

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

Crossrefs

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

Programs

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

Formula

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

Extensions

More terms from Jinyuan Wang, Feb 14 2025

A339561 Products of distinct squarefree semiprimes.

Original entry on oeis.org

1, 6, 10, 14, 15, 21, 22, 26, 33, 34, 35, 38, 39, 46, 51, 55, 57, 58, 60, 62, 65, 69, 74, 77, 82, 84, 85, 86, 87, 90, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 126, 129, 132, 133, 134, 140, 141, 142, 143, 145, 146, 150, 155, 156, 158, 159, 161, 166
Offset: 1

Views

Author

Gus Wiseman, Dec 13 2020

Keywords

Comments

First differs from A320911 in lacking 36.
A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
The following are equivalent characteristics for any positive integer n:
(1) the prime factors of n can be partitioned into distinct strict pairs (a set of edges);
(2) n can be factored into distinct squarefree semiprimes;
(3) the prime signature of n is graphical.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}        55: {3,5}         91: {4,6}
      6: {1,2}     57: {2,8}         93: {2,11}
     10: {1,3}     58: {1,10}        94: {1,15}
     14: {1,4}     60: {1,1,2,3}     95: {3,8}
     15: {2,3}     62: {1,11}       106: {1,16}
     21: {2,4}     65: {3,6}        111: {2,12}
     22: {1,5}     69: {2,9}        115: {3,9}
     26: {1,6}     74: {1,12}       118: {1,17}
     33: {2,5}     77: {4,5}        119: {4,7}
     34: {1,7}     82: {1,13}       122: {1,18}
     35: {3,4}     84: {1,1,2,4}    123: {2,13}
     38: {1,8}     85: {3,7}        126: {1,2,2,4}
     39: {2,6}     86: {1,14}       129: {2,14}
     46: {1,9}     87: {2,10}       132: {1,1,2,5}
     51: {2,7}     90: {1,2,2,3}    133: {4,8}
For example, the number 1260 can be factored into distinct squarefree semiprimes in two ways, (6*10*21) or (6*14*15), so 1260 is in the sequence. The number 69300 can be factored into distinct squarefree semiprimes in seven ways:
  (6*10*15*77)
  (6*10*21*55)
  (6*10*33*35)
  (6*14*15*55)
  (6*15*22*35)
  (10*14*15*33)
  (10*15*21*22),
so 69300 is in the sequence. A complete list of all strict factorizations of 24 is: (2*3*4), (2*12), (3*8), (4*6), (24), all of which contain at least one number that is not a squarefree semiprime, so 24 is not in the sequence.
		

Crossrefs

A309356 is a kind of universal embedding.
A320894 is the complement in A028260.
A320911 lists all (not just distinct) products of squarefree semiprimes.
A339560 counts the partitions with these Heinz numbers.
A339661 has nonzero terms at these positions.
A001358 lists semiprimes, with squarefree case A006881.
A005117 lists squarefree numbers.
A320656 counts factorizations into squarefree semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A058696 counts partitions of 2n (A300061).
- A000070 counts non-multigraphical partitions of 2n (A339620).
- A209816 counts multigraphical partitions (A320924).
- A320921 counts connected graphical partitions (A320923).
- 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:
- 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).
- A339560 can be partitioned into distinct strict pairs (A339561 [this sequence]).

Programs

  • Mathematica
    sqs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqs[n/d],Min@@#>d&]],{d,Select[Divisors[n],SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Select[Range[100],sqs[#]!={}&]

Formula

A338914 Number of integer partitions of n of even length whose greatest multiplicity is at most half their length.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 3, 4, 6, 9, 11, 16, 23, 29, 39, 53, 69, 90, 118, 150, 195, 249, 315, 398, 506, 629, 789, 982, 1219, 1504, 1860, 2277, 2798, 3413, 4161, 5051, 6137, 7406, 8948, 10765, 12943, 15503, 18571, 22153, 26432, 31432, 37352, 44268, 52444, 61944, 73141
Offset: 0

Views

Author

Gus Wiseman, Dec 09 2020

Keywords

Comments

These are also integer partitions that can be partitioned into not necessarily distinct edges (pairs of distinct parts). For example, (3,3,2,2) can be partitioned as {{2,3},{2,3}}, so is counted under a(10), but (4,2,2,2) and (4,2,1,1,1,1) cannot be partitioned into edges. The multiplicities of such a partition form a multigraphical partition (A209816, A320924).

Examples

			The a(3) = 1 through a(10) = 11 partitions:
  (21)  (31)  (32)  (42)    (43)    (53)    (54)      (64)
              (41)  (51)    (52)    (62)    (63)      (73)
                    (2211)  (61)    (71)    (72)      (82)
                            (3211)  (3221)  (81)      (91)
                                    (3311)  (3321)    (3322)
                                    (4211)  (4221)    (4321)
                                            (4311)    (4411)
                                            (5211)    (5221)
                                            (222111)  (5311)
                                                      (6211)
                                                      (322111)
		

Crossrefs

A096373 counts the complement in even-length partitions.
A320911 gives the Heinz numbers of these partitions.
A339560 is the strict case.
A339562 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.
A320656 counts factorizations into squarefree semiprimes.
A320921 counts connected graphical partitions, ranked by A320923.
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.
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).
- 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
    Table[Length[Select[IntegerPartitions[n],EvenQ[Length[#]]&&Max@@Length/@Split[#]<=Length[#]/2&]],{n,0,30}]

Formula

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

A338901 Position of the first appearance of prime(n) as a factor in the list of squarefree semiprimes.

Original entry on oeis.org

1, 1, 2, 3, 6, 7, 9, 11, 13, 17, 18, 21, 23, 25, 29, 31, 34, 36, 40, 42, 45, 47, 50, 52, 56, 58, 61, 64, 67, 70, 76, 78, 81, 82, 86, 89, 93, 97, 100, 104, 106, 107, 112, 113, 116, 118, 125, 129, 133, 134, 135, 139, 141, 147, 150, 154, 159, 160, 165, 167, 169
Offset: 1

Views

Author

Gus Wiseman, Nov 16 2020

Keywords

Comments

The a(n)-th squarefree semiprime is the first divisible by prime(n).
After a(1) = 1, these are the positions of even terms in the list of all squarefree semiprimes A006881.

Crossrefs

A001358 lists semiprimes, with odds A046315 and evens A100484.
A004526 counts 2-part partitions, with strict case A140106 (shifted left).
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odds A046388 and evens A100484.
A115392 is the not necessarily squarefree version.
A166237 gives the first differences of squarefree semiprimes.
A270650 and A270652 give the prime indices of squarefree semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A338898 gives prime indices of semiprimes, with differences A176506.
A338899 gives prime indices of squarefree semiprimes, differences A338900.
A338912 and A338913 give the prime indices of semiprimes.

Programs

  • Mathematica
    rs=First/@FactorInteger[#]&/@Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&];
    Table[Position[rs,i][[1,1]],{i,Union@@rs}]

Formula

A006881(a(n)) = A100484(n).

A339559 Number of integer partitions of n that have an even number of parts and cannot be partitioned into distinct pairs of distinct parts, i.e., that are not the multiset union of any set of edges.

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 4, 3, 7, 6, 14, 14, 23, 27, 41, 47, 70, 84, 114, 141, 190, 225, 303, 370, 475, 578, 738, 890, 1131, 1368, 1698, 2058, 2549, 3048, 3759, 4505, 5495, 6574, 7966, 9483, 11450, 13606, 16307, 19351, 23116, 27297, 32470, 38293, 45346, 53342, 62939
Offset: 0

Views

Author

Gus Wiseman, Dec 10 2020

Keywords

Comments

The multiplicities of such a partition form a non-graphical partition.

Examples

			The a(2) = 1 through a(10) = 14 partitions (empty column indicated by dot):
  11   .   22     2111   33       2221     44         3222       55
           1111          2211     4111     2222       6111       3322
                         3111     211111   3311       222111     3331
                         111111            5111       321111     4222
                                           221111     411111     4411
                                           311111     21111111   7111
                                           11111111              222211
                                                                 322111
                                                                 331111
                                                                 421111
                                                                 511111
                                                                 22111111
                                                                 31111111
                                                                 1111111111
For example, the partition y = (4,4,3,3,2,2,1,1,1,1) can be partitioned into a multiset of edges in just three ways:
  {{1,2},{1,2},{1,3},{1,4},{3,4}}
  {{1,2},{1,3},{1,3},{1,4},{2,4}}
  {{1,2},{1,3},{1,4},{1,4},{2,3}}
None of these are strict, so y is counted under a(22).
		

Crossrefs

A320894 ranks these partitions (using Heinz numbers).
A338915 allows equal pairs (x,x).
A339560 counts the complement in even-length partitions.
A339564 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.
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).
- A339560 can be partitioned into distinct strict pairs (A339561).

Programs

  • Mathematica
    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],EvenQ[Length[#]]&&strs[Times@@Prime/@#]=={}&]],{n,0,15}]

Formula

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

Extensions

More terms from Jinyuan Wang, Feb 14 2025

A339003 Numbers of the form prime(x) * prime(y) where x and y are distinct and both odd.

Original entry on oeis.org

10, 22, 34, 46, 55, 62, 82, 85, 94, 115, 118, 134, 146, 155, 166, 187, 194, 205, 206, 218, 235, 253, 254, 274, 295, 298, 314, 334, 335, 341, 358, 365, 382, 391, 394, 415, 422, 451, 454, 466, 482, 485, 514, 515, 517, 527, 538, 545, 554, 566, 614, 626, 635, 649
Offset: 1

Views

Author

Gus Wiseman, Nov 21 2020

Keywords

Comments

The squarefree semiprimes in A332822. - Peter Munn, Dec 25 2020

Examples

			The sequence of terms together with their prime indices begins:
     10: {1,3}     187: {5,7}     358: {1,41}    527: {7,11}
     22: {1,5}     194: {1,25}    365: {3,21}    538: {1,57}
     34: {1,7}     205: {3,13}    382: {1,43}    545: {3,29}
     46: {1,9}     206: {1,27}    391: {7,9}     554: {1,59}
     55: {3,5}     218: {1,29}    394: {1,45}    566: {1,61}
     62: {1,11}    235: {3,15}    415: {3,23}    614: {1,63}
     82: {1,13}    253: {5,9}     422: {1,47}    626: {1,65}
     85: {3,7}     254: {1,31}    451: {5,13}    635: {3,31}
     94: {1,15}    274: {1,33}    454: {1,49}    649: {5,17}
    115: {3,9}     295: {3,17}    466: {1,51}    662: {1,67}
    118: {1,17}    298: {1,35}    482: {1,53}    685: {3,33}
    134: {1,19}    314: {1,37}    485: {3,25}    694: {1,69}
    146: {1,21}    334: {1,39}    514: {1,55}    697: {7,13}
    155: {3,11}    335: {3,19}    515: {3,27}    706: {1,71}
    166: {1,23}    341: {5,11}    517: {5,15}    713: {9,11}
		

Crossrefs

A338910 is the not necessarily squarefree version.
A339004 is the even instead of odd version.
A001358 lists semiprimes, with odd and even terms A046315 and A100484.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odd and even terms A046388 and A100484.
A289182/A115392 list the positions of odd/even terms of A001358.
A300912 lists products of two primes of relatively prime index.
A320656 counts factorizations into squarefree semiprimes.
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.
A338904 groups semiprimes by weight.
A338906/A338907 list semiprimes of even/odd weight.
A339002 lists products of two distinct primes of non-relatively prime index.
A339005 lists products of two distinct primes of divisible index.
Subsequence of A332822.

Programs

  • Mathematica
    Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&&OddQ[Times@@PrimePi/@First/@FactorInteger[#]]&]
  • Python
    from math import isqrt
    from sympy import primepi, primerange
    def A339003(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),1) if a&1)
        return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

Formula

Numbers m such that A001221(m) = A001222(m) = A195017(m) = 2. - Peter Munn, Dec 31 2020

A339004 Numbers of the form prime(x) * prime(y) where x and y are distinct and both even.

Original entry on oeis.org

21, 39, 57, 87, 91, 111, 129, 133, 159, 183, 203, 213, 237, 247, 259, 267, 301, 303, 321, 339, 371, 377, 393, 417, 427, 453, 481, 489, 497, 519, 543, 551, 553, 559, 579, 597, 623, 669, 687, 689, 703, 707, 717, 749, 753, 789, 791, 793, 813, 817, 843, 879, 917
Offset: 1

Views

Author

Gus Wiseman, Nov 22 2020

Keywords

Comments

The squarefree semiprimes in A332821. - Peter Munn, Dec 25 2020

Examples

			The sequence of terms together with their prime indices begins:
     21: {2,4}     267: {2,24}    543: {2,42}
     39: {2,6}     301: {4,14}    551: {8,10}
     57: {2,8}     303: {2,26}    553: {4,22}
     87: {2,10}    321: {2,28}    559: {6,14}
     91: {4,6}     339: {2,30}    579: {2,44}
    111: {2,12}    371: {4,16}    597: {2,46}
    129: {2,14}    377: {6,10}    623: {4,24}
    133: {4,8}     393: {2,32}    669: {2,48}
    159: {2,16}    417: {2,34}    687: {2,50}
    183: {2,18}    427: {4,18}    689: {6,16}
    203: {4,10}    453: {2,36}    703: {8,12}
    213: {2,20}    481: {6,12}    707: {4,26}
    237: {2,22}    489: {2,38}    717: {2,52}
    247: {6,8}     497: {4,20}    749: {4,28}
    259: {4,12}    519: {2,40}    753: {2,54}
		

Crossrefs

A338911 is the not necessarily squarefree version.
A339003 is the odd instead of even version, with not necessarily squarefree version A338910.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A289182/A115392 list the positions of odd/even terms in A001358.
A300912 lists products of pairs of primes with relatively prime indices.
A318990 lists products of pairs of primes with divisible indices.
A320656 counts factorizations into squarefree semiprimes.
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.
A338904 groups semiprimes by weight.
A338906/A338907 list semiprimes of even/odd weight.
Subsequence of A332821.

Programs

  • Mathematica
    Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&&OddQ[Times@@(1+ PrimePi/@First/@FactorInteger[#])]&]
  • Python
    from math import isqrt
    from sympy import primepi, primerange
    def A339004(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),1) if a&1^1)
        return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

Formula

Numbers m such that A001221(m) = A001222(m) = 2 and A195017(m) = -2. - Peter Munn, Dec 31 2020

A338902 Number of integer partitions of the n-th semiprime into semiprimes.

Original entry on oeis.org

1, 1, 1, 2, 3, 2, 4, 7, 7, 10, 17, 25, 21, 34, 34, 73, 87, 103, 149, 176, 206, 281, 344, 479, 725, 881, 1311, 1597, 1742, 1841, 2445, 2808, 3052, 3222, 6784, 9298, 11989, 14533, 15384, 17414, 18581, 19680, 28284, 35862, 38125, 57095, 60582, 64010, 71730, 76016
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2020

Keywords

Comments

A semiprime (A001358) is a product of any two prime numbers.

Examples

			The a(1) = 1 through a(33) = 17 partitions of 4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, where A-Z = 10-35:
  4  6  9  A   E    F   L     M      P      Q       X
           64  A4   96  F6    994    FA     M4      EA9
               644      966   A66    L4     AA6     F99
                        9444  E44    A96    E66     FE4
                              6664   F64    9944    L66
                              A444   9664   A664    P44
                              64444  94444  E444    9996
                                            66644   AA94
                                            A4444   E964
                                            644444  F666
                                                    FA44
                                                    L444
                                                    96666
                                                    A9644
                                                    F6444
                                                    966444
                                                    9444444
		

Crossrefs

A002100 counts partitions into squarefree semiprimes.
A056768 uses primes instead of semiprimes.
A101048 counts partitions into semiprimes.
A338903 is the squarefree version.
A339112 includes the Heinz numbers of these partitions.
A001358 lists semiprimes, with odd and even terms A046315 and A100484.
A037143 lists primes and semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A320655 counts factorizations into semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes, with sum/difference/product A176504/A176506/A087794.
A338899/A270650/A270652 give prime indices of squarefree semiprimes.

Programs

  • Mathematica
    nn=100;Table[Length[IntegerPartitions[n,All,Select[Range[nn],PrimeOmega[#]==2&]]],{n,Select[Range[nn],PrimeOmega[#]==2&]}]

Formula

a(n) = A101048(A001358(n)).

A339191 Partial products of squarefree semiprimes (A006881).

Original entry on oeis.org

6, 60, 840, 12600, 264600, 5821200, 151351200, 4994589600, 169816046400, 5943561624000, 225855341712000, 8808358326768000, 405184483031328000, 20664408634597728000, 1136542474902875040000, 64782921069463877280000, 3757409422028904882240000
Offset: 1

Views

Author

Gus Wiseman, Nov 30 2020

Keywords

Comments

A squarefree semiprime is a product of any two distinct prime numbers.
Do all terms belong to A242031 (weakly decreasing prime signature)?

Examples

			The sequence of terms together with their prime indices begins:
          6: {1,2}
         60: {1,1,2,3}
        840: {1,1,1,2,3,4}
      12600: {1,1,1,2,2,3,3,4}
     264600: {1,1,1,2,2,2,3,3,4,4}
    5821200: {1,1,1,1,2,2,2,3,3,4,4,5}
  151351200: {1,1,1,1,1,2,2,2,3,3,4,4,5,6}
The sequence of terms together with their prime signatures begins:
                   6: (1,1)
                  60: (2,1,1)
                 840: (3,1,1,1)
               12600: (3,2,2,1)
              264600: (3,3,2,2)
             5821200: (4,3,2,2,1)
           151351200: (5,3,2,2,1,1)
          4994589600: (5,4,2,2,2,1)
        169816046400: (6,4,2,2,2,1,1)
       5943561624000: (6,4,3,3,2,1,1)
     225855341712000: (7,4,3,3,2,1,1,1)
    8808358326768000: (7,5,3,3,2,2,1,1)
  405184483031328000: (8,5,3,3,2,2,1,1,1)
		

Crossrefs

A000040 lists the primes, with partial products A002110 (primorials).
A001358 lists semiprimes, with partial products A112141.
A002100 counts partitions into squarefree semiprimes (restricted: A338903)
A000142 lists factorial numbers, with partial products A000178.
A005117 lists squarefree numbers, with partial products A111059.
A006881 lists squarefree semiprimes, with partial sums A168472.
A166237 gives first differences of squarefree semiprimes.
A320655 counts factorizations into semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes.
A338899/A270650/A270652 give prime indices of squarefree semiprimes.
A338901 gives first appearances in the list of squarefree semiprimes.
A339113 gives products of primes of squarefree semiprime index.

Programs

  • Mathematica
    FoldList[Times,Select[Range[20],SquareFreeQ[#]&&PrimeOmega[#]==2&]]
Showing 1-10 of 10 results.