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.

Previous Showing 31-40 of 61 results. Next

A339620 Heinz numbers of non-multigraphical partitions of even numbers.

Original entry on oeis.org

3, 7, 10, 13, 19, 21, 22, 28, 29, 34, 37, 39, 43, 46, 52, 53, 55, 57, 61, 62, 66, 71, 76, 79, 82, 85, 87, 88, 89, 91, 94, 101, 102, 107, 111, 113, 115, 116, 117, 118, 129, 130, 131, 133, 134, 136, 138, 139, 146, 148, 151, 155, 156, 159, 163, 166, 171, 172, 173
Offset: 1

Views

Author

Gus Wiseman, Dec 18 2020

Keywords

Comments

An integer partition is non-multigraphical if it does not comprise the multiset of vertex-degrees of any multigraph (multiset of non-loop edges). Multigraphical partitions are counted by A209816, non-multigraphical partitions by A000070.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.
The following are equivalent characteristics for any positive integer n:
(1) the multiset of prime indices of n can be partitioned into strict pairs (a multiset of edges);
(2) n can be factored into squarefree semiprimes;
(3) the unordered prime signature of n is multigraphical.

Examples

			The sequence of terms together with their prime indices begins:
      3: {2}         53: {16}          94: {1,15}
      7: {4}         55: {3,5}        101: {26}
     10: {1,3}       57: {2,8}        102: {1,2,7}
     13: {6}         61: {18}         107: {28}
     19: {8}         62: {1,11}       111: {2,12}
     21: {2,4}       66: {1,2,5}      113: {30}
     22: {1,5}       71: {20}         115: {3,9}
     28: {1,1,4}     76: {1,1,8}      116: {1,1,10}
     29: {10}        79: {22}         117: {2,2,6}
     34: {1,7}       82: {1,13}       118: {1,17}
     37: {12}        85: {3,7}        129: {2,14}
     39: {2,6}       87: {2,10}       130: {1,3,6}
     43: {14}        88: {1,1,1,5}    131: {32}
     46: {1,9}       89: {24}         133: {4,8}
     52: {1,1,6}     91: {4,6}        134: {1,19}
For example, a complete lists of all loop-multigraphs with degrees (5,2,1) is:
  {{1,1},{1,1},{1,2},{2,3}}
  {{1,1},{1,1},{1,3},{2,2}}
  {{1,1},{1,2},{1,2},{1,3}},
but since none of these is a multigraph (they have loops), the Heinz number 66 belongs to the sequence.
		

Crossrefs

A000070 counts these partitions.
A300061 is a superset.
A320891 has image under A181819 equal to this set of terms.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
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 [this sequence]).
- A209816 counts multigraphical partitions (A320924).
- A147878 counts connected multigraphical partitions (A320925).
- 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).

Programs

  • Mathematica
    prpts[m_]:=If[Length[m]==0,{{}},Join@@Table[Prepend[#,ipr]&/@prpts[Fold[DeleteCases[#1,#2,{1},1]&,m,ipr]],{ipr,Select[Subsets[Union[m],{2}],MemberQ[#,m[[1]]]&]}]];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],EvenQ[Length[nrmptn[#]]]&&prpts[nrmptn[#]]=={}&]

Formula

Equals A300061 \ A320924.
For all n, both A181821(a(n)) and A304660(a(n)) belong to A320891.

A339659 Irregular triangle read by rows where T(n,k) is the number of graphical partitions of 2n into k parts, 0 <= k <= 2n.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 2, 1, 1, 0, 0, 0, 0, 2, 3, 2, 1, 1, 0, 0, 0, 0, 1, 4, 5, 3, 2, 1, 1, 0, 0, 0, 0, 1, 4, 7, 7, 5, 3, 2, 1, 1, 0, 0, 0, 0, 0, 4, 9, 11, 11, 7, 5, 3, 2, 1, 1, 0, 0, 0, 0, 0, 2, 11, 15, 17, 15, 11, 7, 5, 3, 2, 1, 1
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2020

Keywords

Comments

Conjecture: The column sums 1, 0, 1, 2, 7, 20, 67, ... are given by A304787.
An integer partition is graphical if it comprises the multiset of vertex-degrees of some graph. Graphical partitions are counted by A000569.

Examples

			Triangle begins:
  1
  0 0 1
  0 0 0 1 1
  0 0 0 1 2 1 1
  0 0 0 0 2 3 2 1 1
  0 0 0 0 1 4 5 3 2 1 1
  0 0 0 0 1 4 7 7 5 3 2 1 1
For example, row n = 5 counts the following partitions:
  3322  22222  222211  2221111  22111111  211111111  1111111111
        32221  322111  3211111  31111111
        33211  331111  4111111
        42211  421111
               511111
		

Crossrefs

A000569 gives the row sums.
A004250 is the central column.
A005408 gives the row lengths.
A008284/A072233 is the version counting all partitions.
A259873 is the left half of the triangle.
A309356 is a universal embedding.
A027187 counts partitions of even length.
A339559 = partitions that cannot be partitioned into distinct strict pairs.
A339560 = partitions that can be partitioned into distinct strict pairs.
The following count vertex-degree partitions and give their Heinz numbers:
- A000070 counts non-multigraphical partitions of 2n (A339620).
- A000569 counts graphical partitions (A320922).
- A058696 counts partitions of 2n (A300061).
- A147878 counts connected multigraphical partitions (A320925).
- A209816 counts multigraphical partitions (A320924).
- A320921 counts connected graphical partitions (A320923).
- A321728 is conjectured to count non-half-loop-graphical partitions of n.
- A339617 counts non-graphical partitions of 2n (A339618).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).

Programs

  • Mathematica
    prpts[m_]:=If[Length[m]==0,{{}},Join@@Table[Prepend[#,ipr]&/@prpts[Fold[DeleteCases[#1,#2,{1},1]&,m,ipr]],{ipr,Subsets[Union[m],{2}]}]];
    strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
    Table[Length[Select[strnorm[2*n],Length[Union[#]]==k&&Select[prpts[#],UnsameQ@@#&]!={}&]],{n,0,5},{k,0,2*n}]

A322353 Number of factorizations of n into distinct semiprimes; a(1) = 1 by convention.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 2, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 2, 1, 1, 1, 1, 0, 2, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2018

Keywords

Comments

A semiprime (A001358) is a product of any two prime numbers. In the even case, these factorizations have A001222(n)/2 factors. - Gus Wiseman, Dec 31 2020
Records 1, 2, 3, 4, 5, 9, 13, 15, 17, ... occur at 1, 60, 210, 840, 1260, 4620, 27720, 30030, 69300, ...

Examples

			a(4) = 1, as there is just one way to factor 4 into distinct semiprimes, namely as {4}.
From _Gus Wiseman_, Dec 31 2020: (Start)
The a(n) factorizations for n = 60, 210, 840, 1260, 4620, 12600, 18480:
  4*15   6*35    4*6*35    4*9*35    4*15*77    4*6*15*35    4*6*10*77
  6*10   10*21   4*10*21   4*15*21   4*21*55    4*6*21*25    4*6*14*55
         14*15   4*14*15   6*10*21   4*33*35    4*9*10*35    4*6*22*35
                 6*10*14   6*14*15   6*10*77    4*9*14*25    4*10*14*33
                           9*10*14   6*14*55    4*10*15*21   4*10*21*22
                                     6*22*35    6*10*14*15   4*14*15*22
                                     10*14*33                6*10*14*22
                                     10*21*22
                                     14*15*22
(End)
		

Crossrefs

Unlabeled multiset partitions of this type are counted by A007717.
The version for partitions is A112020, or A101048 without distinctness.
The non-strict version is A320655.
Positions of zeros include A320892.
Positions of nonzero terms are A320912.
The case of squarefree factors is A339661, or A320656 without distinctness.
Allowing prime factors gives A339839, or A320732 without distinctness.
A322661 counts loop-graphs, ranked by A320461.
A001055 counts factorizations, with strict case A045778.
A001358 lists semiprimes, with squarefree case A006881.
A027187 counts partitions of even length, ranked by A028260.
A037143 lists primes and semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes.
A339846 counts even-length factorizations, with ordered version A174725.

Programs

Formula

a(n) = Sum_{d|n} (-1)^A001222(d) * A339839(n/d). - Gus Wiseman, Dec 31 2020

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

A339888 Number of non-isomorphic multiset partitions of weight n into singletons or strict pairs.

Original entry on oeis.org

1, 1, 3, 5, 13, 23, 55, 104, 236, 470, 1039, 2140, 4712, 9962, 21961, 47484, 105464, 232324, 521338, 1167825, 2651453, 6031136, 13863054, 31987058, 74448415, 174109134, 410265423, 971839195, 2317827540, 5558092098, 13412360692, 32542049038, 79424450486
Offset: 0

Views

Author

Gus Wiseman, Jan 09 2021

Keywords

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 13 multiset partitions:
  {{1}}  {{1,2}}    {{1},{2,3}}    {{1,2},{1,2}}
         {{1},{1}}  {{2},{1,2}}    {{1,2},{3,4}}
         {{1},{2}}  {{1},{1},{1}}  {{1,3},{2,3}}
                    {{1},{2},{2}}  {{1},{1},{2,3}}
                    {{1},{2},{3}}  {{1},{2},{1,2}}
                                   {{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}}
		

Crossrefs

The version for set partitions is A000085, with ordered version A080599.
The case of integer partitions is 1 + A004526(n), ranked by A003586.
Non-isomorphic multiset partitions are counted by A007716.
The case without singletons is A007717.
The version allowing non-strict pairs (x,x) is A320663.
A001190 counts rooted trees with out-degrees <= 2, ranked by A292050.
A339742 counts factorizations into distinct primes or squarefree semiprimes.
A339887 counts factorizations into primes or squarefree semiprimes.

Programs

  • PARI
    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 + ((r-1)\2)*x^(2*r))}
    a(n)={if(n==0, 1, my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(gs(p) + O(x*x^n), -n))[n]); s/n!)} \\ Andrew Howroyd, Apr 16 2021

Extensions

Terms a(11) and beyond from Andrew Howroyd, Apr 16 2021

A340785 Number of factorizations of 2n into even factors > 1.

Original entry on oeis.org

1, 2, 1, 3, 1, 2, 1, 5, 1, 2, 1, 4, 1, 2, 1, 7, 1, 3, 1, 4, 1, 2, 1, 7, 1, 2, 1, 4, 1, 3, 1, 11, 1, 2, 1, 6, 1, 2, 1, 7, 1, 3, 1, 4, 1, 2, 1, 12, 1, 3, 1, 4, 1, 3, 1, 7, 1, 2, 1, 7, 1, 2, 1, 15, 1, 3, 1, 4, 1, 3, 1, 12, 1, 2, 1, 4, 1, 3, 1, 12, 1, 2, 1, 7, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 30 2021

Keywords

Examples

			The a(n) factorizations for n = 2*2, 2*4, 2*8, 2*12, 2*16, 2*32, 2*36, 2*48 are:
  4    8      16       24     32         64           72      96
  2*2  2*4    2*8      4*6    4*8        8*8          2*36    2*48
       2*2*2  4*4      2*12   2*16       2*32         4*18    4*24
              2*2*4    2*2*6  2*2*8      4*16         6*12    6*16
              2*2*2*2         2*4*4      2*4*8        2*6*6   8*12
                              2*2*2*4    4*4*4        2*2*18  2*6*8
                              2*2*2*2*2  2*2*16               4*4*6
                                         2*2*2*8              2*2*24
                                         2*2*4*4              2*4*12
                                         2*2*2*2*4            2*2*4*6
                                         2*2*2*2*2*2          2*2*2*12
                                                              2*2*2*2*6
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
The version for partitions is A035363 (A066207).
The odd version is A340101.
The even length case is A340786.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
A316439 counts factorizations by product and length
A340102 counts odd-length factorizations of odd numbers into odd factors.
- Even -
A027187 counts partitions of even length/maximum (A028260/A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum.
A340601 counts partitions of even rank (A340602).
Even bisection of A349906.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Select[#,OddQ]=={}&]],{n,2,100,2}]
  • PARI
    A349906(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&!(d%2), s += A349906(n/d, d))); (s));
    A340785(n) = A349906(2*n); \\ Antti Karttunen, Dec 13 2021

Formula

a(n) = A349906(2*n). - Antti Karttunen, Dec 13 2021

A338905 Irregular triangle read by rows where row n lists all squarefree semiprimes with prime indices summing to n.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 33, 35, 34, 39, 55, 38, 51, 65, 77, 46, 57, 85, 91, 58, 69, 95, 119, 143, 62, 87, 115, 133, 187, 74, 93, 145, 161, 209, 221, 82, 111, 155, 203, 247, 253, 86, 123, 185, 217, 299, 319, 323, 94, 129, 205, 259, 341, 377, 391, 106, 141
Offset: 3

Views

Author

Gus Wiseman, Nov 28 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.

Examples

			Triangle begins:
   6
  10
  14  15
  21  22
  26  33  35
  34  39  55
  38  51  65  77
  46  57  85  91
  58  69  95 119 143
  62  87 115 133 187
  74  93 145 161 209 221
  82 111 155 203 247 253
  86 123 185 217 299 319 323
		

Crossrefs

A004526 (shifted right) gives row lengths.
A025129 (shifted right) gives row sums.
A056239 gives sum of prime indices (Heinz weight).
A339116 is a different triangle whose diagonals are these rows.
A338904 is the not necessarily squarefree version, with row sums A024697.
A338907/A338908 are the union of odd/even rows.
A339114/A332765 are the row minima/maxima.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A087112 groups semiprimes by greater factor.
A168472 gives partial sums of 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.

Programs

  • Mathematica
    Table[Sort[Table[Prime[k]*Prime[n-k],{k,(n-1)/2}]],{n,3,10}]

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

A339195 Triangle of squarefree numbers grouped by greatest prime factor, read by rows.

Original entry on oeis.org

1, 2, 3, 6, 5, 10, 15, 30, 7, 14, 21, 35, 42, 70, 105, 210, 11, 22, 33, 55, 66, 77, 110, 154, 165, 231, 330, 385, 462, 770, 1155, 2310, 13, 26, 39, 65, 78, 91, 130, 143, 182, 195, 273, 286, 390, 429, 455, 546, 715, 858, 910, 1001, 1365, 1430, 2002, 2145, 2730, 3003, 4290, 5005, 6006, 10010, 15015, 30030
Offset: 0

Views

Author

Gus Wiseman, Dec 02 2020

Keywords

Comments

Also Heinz numbers of subsets of {1..n} that contain n if n>0, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
A019565 in its triangle form, with each row's terms in increasing order. - Peter Munn, Feb 26 2021
From David James Sycamore, Jan 09 2025: (Start)
Alternative definition, with offset = 1: a(1) = 1. For n>1 if a(n-1) = A002110(k), a(n) = prime(k+1). Otherwise a(n) is the smallest novel squarefree number whose prime factors have already occurred as previous terms.
Permutation of A005117, Squarefree version A379746. (End)

Examples

			Triangle begins:
   1
   2
   3   6
   5  10  15  30
   7  14  21  35  42  70  105  210
		

Crossrefs

A011782 gives row lengths.
A339360 gives row sums.
A008578 (shifted) is column k = 1.
A100484 is column k = 2.
A001748 is column k = 3.
A002110 is column k = 2^(n-1).
A070826 is column k = 2^(n-1) - 1.
A209862 takes prime indices to binary indices in these terms.
A246867 groups squarefree numbers by Heinz weight, with row sums A147655.
A261144 divides the n-th row by prime(n), with row sums A054640.
A339116 is the restriction to semiprimes, with row sums A339194.
A005117 lists squarefree numbers, ordered lexicographically by prime factors: A019565.
A006881 lists squarefree semiprimes.
A072047 counts prime factors of squarefree numbers.
A319246 is the sum of prime indices of the n-th squarefree number.
A329631 lists prime indices of squarefree numbers, reversed: A319247.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
Cf. A379746.

Programs

  • Maple
    T:= proc(n) option remember; `if`(n=0, 1, (p-> map(
          x-> x*p, {seq(T(i), i=0..n-1)})[])(ithprime(n)))
        end:
    seq(T(n), n=0..6);  # Alois P. Heinz, Jan 08 2025
  • Mathematica
    Table[Prime[n]*Sort[Times@@Prime/@#&/@Subsets[Range[n-1]]],{n,5}]

Formula

For n > 1, T(n,k) = prime(n) * A261144(n-1,k).
a(n) = A019565(A379770(n)). - Michael De Vlieger, Jan 08 2025

Extensions

Row n=0 (term 1) prepended by Alois P. Heinz, Jan 08 2025

A339742 Number of factorizations of n into distinct primes or squarefree semiprimes.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 0, 2, 0, 1, 1, 4, 1, 0, 2, 2, 2, 1, 1, 2, 2, 0, 1, 4, 1, 1, 1, 2, 1, 0, 0, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 3, 1, 2, 1, 0, 2, 4, 1, 1, 2, 4, 1, 0, 1, 2, 1, 1, 2, 4, 1, 0, 0, 2, 1, 3, 2, 2, 2, 0, 1, 3, 2, 1, 2, 2, 2, 0, 1, 1, 1, 1, 1, 4, 1, 0, 4
Offset: 1

Views

Author

Gus Wiseman, Dec 20 2020

Keywords

Comments

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 singletons or strict pairs, i.e., into a set of half-loops and edges;
(2) n can be factored into distinct primes or squarefree semiprimes.

Examples

			The a(n) factorizations for n = 6, 30, 60, 210, 420 are respectively 2, 4, 3, 10, 9:
  (6)    (5*6)    (6*10)    (6*35)     (2*6*35)
  (2*3)  (2*15)   (2*5*6)   (10*21)    (5*6*14)
         (3*10)   (2*3*10)  (14*15)    (6*7*10)
         (2*3*5)            (5*6*7)    (2*10*21)
                            (2*3*35)   (2*14*15)
                            (2*5*21)   (2*5*6*7)
                            (2*7*15)   (3*10*14)
                            (3*5*14)   (2*3*5*14)
                            (3*7*10)   (2*3*7*10)
                            (2*3*5*7)
		

Crossrefs

Dirichlet convolution of A008966 with A339661.
A008966 allows only primes.
A339661 does not allow primes, only squarefree semiprimes.
A339740 lists the positions of zeros.
A339741 lists the positions of positive terms.
A339839 allows nonsquarefree semiprimes.
A339887 is the non-strict version.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A013929 cannot be factored into distinct primes.
A293511 are a product of distinct squarefree numbers in exactly one way.
A320663 counts non-isomorphic multiset partitions into singletons or pairs.
A339840 cannot be factored into distinct primes or semiprimes.
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A050320 into squarefree numbers.
- A050326 into distinct squarefree numbers.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A320732 into primes or semiprimes.
- A322353 into distinct semiprimes.
- A339742 [this sequence] into distinct primes or squarefree semiprimes.
- A339839 into distinct primes or semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A000569 counts graphical partitions (A320922).
- A058696 counts all partitions of 2n (A300061).
- A209816 counts multigraphical partitions (A320924).
- A339656 counts loop-graphical partitions (A339658).
-
The following count partitions/factorizations of even length and give their Heinz numbers:
- A027187/A339846 has no additional conditions (A028260).
- A338914/A339562 can be partitioned into edges (A320911).
- A338916/A339563 can be partitioned into distinct pairs (A320912).
- A339559/A339564 cannot be partitioned into distinct edges (A320894).
- A339560/A339619 can be partitioned into distinct edges (A339561).

Programs

  • Mathematica
    sqps[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqps[n/d],Min@@#>d&]],{d,Select[Divisors[n],PrimeQ[#]||SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Table[Length[sqps[n]],{n,100}]
  • PARI
    A353471(n) = (numdiv(n)==2*omega(n));
    A339742(n, u=(1+n)) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1) && (dA353471(d), s += A339742(n/d, d))); (s)); \\ Antti Karttunen, May 02 2022

Formula

a(n) = Sum_{d|n squarefree} A339661(n/d).

Extensions

More terms from Antti Karttunen, May 02 2022
Previous Showing 31-40 of 61 results. Next