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 11-20 of 26 results. Next

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

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

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

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 5, 4, 6, 5, 12, 14, 19, 22, 27, 36, 38, 51, 77, 86, 128, 141, 163, 163, 207, 233, 259, 260, 514, 657, 813, 983, 1010, 1215, 1255, 1720, 2112, 2256, 3171, 3370, 3499, 3864, 4103, 6292, 7313, 7620, 8374, 10650, 17579, 18462, 23034, 25180
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2020

Keywords

Comments

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

Examples

			The a(n) partitions for n = 1, 5, 7, 9, 10, 11, 13:
  6  21    26       34          35        38           46
     15,6  14,6,6   22,6,6      21,14     26,6,6       34,6,6
           10,10,6  14,14,6     15,14,6   22,10,6      26,14,6
                    14,10,10    15,10,10  14,14,10     21,15,10
                    10,6,6,6,6            14,6,6,6,6   22,14,10
                                          10,10,6,6,6  26,10,10
                                                       15,15,10,6
                                                       22,6,6,6,6
                                                       14,14,6,6,6
                                                       14,10,10,6,6
                                                       10,10,10,10,6
                                                       10,6,6,6,6,6,6
		

Crossrefs

A002100 counts partitions into squarefree semiprimes.
A056768 uses primes instead of squarefree semiprimes.
A101048 counts partitions into semiprimes.
A338902 is the not necessarily squarefree version.
A339113 includes the Heinz numbers of these partitions.
A001358 lists semiprimes, with odd and even terms A046315 and A100484.
A006881 lists squarefree semiprimes, with odd and even terms A046388 and A100484.
A320656 counts factorizations into squarefree semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes, with sum/difference/product A176504/A176506/A087794.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes.

Programs

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

Formula

a(n) = A002100(A006881(n)).

A339839 Number of factorizations of n into distinct primes or semiprimes.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 0, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 0, 2, 2, 2, 2, 1, 2, 2, 2, 1, 4, 1, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 0, 2, 4, 1, 2, 2, 4, 1, 2, 1, 2, 2, 2, 2, 4, 1, 1, 0, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 0, 1, 2, 2, 2, 1, 4, 1, 2, 4
Offset: 1

Views

Author

Gus Wiseman, Dec 20 2020

Keywords

Comments

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

Examples

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

Crossrefs

A008966 allows only primes.
A320732 is the non-strict version.
A339742 does not allow squares of primes.
A339840 lists the positions of zeros.
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.
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A322353 into distinct semiprimes.
- A339839 [this sequence] into distinct primes or semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A000569 counts graphical partitions (A320922).
- A339656 counts loop-graphical partitions (A339658).

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],UnsameQ@@#&&SubsetQ[{1,2},PrimeOmega/@#]&]],{n,100}]
  • PARI
    A339839(n, u=(1+n)) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1) && (dA339839(n/d, d))); (s)); \\ Antti Karttunen, Feb 10 2023

Formula

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

Extensions

Data section extended up to a(105) by Antti Karttunen, Feb 10 2023

A117825 Distance from n-th highly composite number (cf. A002182) to nearest prime.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 11, 13, 1, 11, 1, 17, 1, 1, 13, 13, 1, 1, 17, 1, 17, 1, 1, 17, 17, 17, 1, 1, 19, 37, 37, 1, 17, 23, 1, 29, 1, 1, 19, 1, 19, 23, 1, 19, 31, 1, 19, 1, 1, 1, 1, 23, 1, 29, 23, 23, 1, 23, 71, 37, 1, 1, 31, 1, 23, 53, 1, 31
Offset: 1

Views

Author

Bill McEachen, May 01 2006

Keywords

Comments

a) Conjecture: entries > 1 will always be prime. The entry will be larger than the largest prime factor of the highly composite number.
b) Will 1 always be the most common entry?
c) While a prime may always be located close to each highly composite number, is the converse false?
d) Is there always a prime between successive highly composite numbers?
From Antti Karttunen, Feb 26 2019: (Start)
The second sentence of point (a) follows as both gcd(n, A151799(n)) = 1 and gcd(A151800(n), n) = 1 for all n > 2 and the fact that the highly composite numbers are products of primorials, A002110 (with the least coprime prime > the largest prime factor). See also the conjectures and notes in A129912 and A141345. (End)

Examples

			a(5) = abs(12-11) = 1.
		

Crossrefs

Sequences tied to conjecture a): A228943, A228945.
Cf. also A005235, A060270.

Programs

Formula

a(1) = 1; for n > 1, a(n) = min(A141345(n), A324385(n)). - Antti Karttunen, Feb 26 2019

Extensions

More terms from Don Reble, May 02 2006

A339740 Non-products of distinct primes or squarefree semiprimes.

Original entry on oeis.org

4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336
Offset: 1

Views

Author

Gus Wiseman, Dec 20 2020

Keywords

Comments

Differs from A293243 and A212164 in having 1080, with prime indices {1,1,1,2,2,2,3} and factorization into distinct squarefree numbers 2*3*6*30.

Examples

			The sequence of terms together with their prime indices begins:
      4: {1,1}             80: {1,1,1,1,3}
      8: {1,1,1}           81: {2,2,2,2}
      9: {2,2}             88: {1,1,1,5}
     16: {1,1,1,1}         96: {1,1,1,1,1,2}
     24: {1,1,1,2}        104: {1,1,1,6}
     25: {3,3}            108: {1,1,2,2,2}
     27: {2,2,2}          112: {1,1,1,1,4}
     32: {1,1,1,1,1}      121: {5,5}
     40: {1,1,1,3}        125: {3,3,3}
     48: {1,1,1,1,2}      128: {1,1,1,1,1,1,1}
     49: {4,4}            135: {2,2,2,3}
     54: {1,2,2,2}        136: {1,1,1,7}
     56: {1,1,1,4}        144: {1,1,1,1,2,2}
     64: {1,1,1,1,1,1}    152: {1,1,1,8}
     72: {1,1,1,2,2}      160: {1,1,1,1,1,3}
For example, a complete list of strict factorizations of 72 is: (2*3*12), (2*4*9), (2*36), (3*4*6), (3*24), (4*18), (6*12), (8*9), (72); but since none of these consists of only primes or squarefree semiprimes, 72 is in the sequence.
		

Crossrefs

A013929 allows only primes.
A320894 does not allow primes (but omega is assumed even).
A339741 is the complement.
A339742 has zeros at these positions.
A339840 allows squares of primes.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A320663 counts non-isomorphic multiset partitions into singletons or pairs.
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A050326 into distinct squarefree numbers.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A320732 into primes or semiprimes.
- A322353 into distinct semiprimes.
- A339661 into distinct squarefree semiprimes.
- A339839 into distinct primes or 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).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339617 counts non-graphical partitions of 2n (A339618).
- A321728 is conjectured to count non-half-loop-graphical partitions of n.
The following count partitions/factorizations of even length and give their Heinz numbers:
- A027187/A339846 counts all of even length (A028260).
- A096373/A339737 cannot be partitioned into strict pairs (A320891).
- A338915/A339662 cannot be partitioned into distinct pairs (A320892).
- A339559/A339564 cannot be partitioned into distinct strict pairs (A320894).

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&]}]];
    Select[Range[100],sqps[#]=={}&]

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

A339840 Numbers that cannot be factored into distinct primes or semiprimes.

Original entry on oeis.org

16, 32, 64, 81, 96, 128, 160, 192, 224, 243, 256, 288, 320, 352, 384, 416, 448, 486, 512, 544, 576, 608, 625, 640, 704, 729, 736, 768, 800, 832, 864, 896, 928, 960, 972, 992, 1024, 1088, 1152, 1184, 1215, 1216, 1280, 1312, 1344, 1376, 1408, 1458, 1472, 1504
Offset: 1

Views

Author

Gus Wiseman, Dec 20 2020

Keywords

Comments

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

Examples

			The sequence of terms together with their prime indices begins:
    16: {1,1,1,1}
    32: {1,1,1,1,1}
    64: {1,1,1,1,1,1}
    81: {2,2,2,2}
    96: {1,1,1,1,1,2}
   128: {1,1,1,1,1,1,1}
   160: {1,1,1,1,1,3}
   192: {1,1,1,1,1,1,2}
   224: {1,1,1,1,1,4}
   243: {2,2,2,2,2}
   256: {1,1,1,1,1,1,1,1}
   288: {1,1,1,1,1,2,2}
   320: {1,1,1,1,1,1,3}
   352: {1,1,1,1,1,5}
   384: {1,1,1,1,1,1,1,2}
   416: {1,1,1,1,1,6}
   448: {1,1,1,1,1,1,4}
   486: {1,2,2,2,2,2}
For example, a complete list of all factorizations of 192 into primes or semiprimes is:
  (2*2*2*2*2*2*3)
  (2*2*2*2*2*6)
  (2*2*2*2*3*4)
  (2*2*2*4*6)
  (2*2*3*4*4)
  (2*4*4*6)
  (3*4*4*4)
Since none of these is strict, 192 is in the sequence.
		

Crossrefs

Allowing only primes gives A013929.
Removing all squares of primes gives A339740.
These are the positions of zeros in A339839.
The complement is A339889.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A293511 are a product of distinct squarefree numbers in exactly one way.
A320663 counts non-isomorphic multiset partitions into singletons or pairs.
A338915 cannot be partitioned into distinct pairs (A320892).
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A320732 into primes or semiprimes.
- A322353 into distinct semiprimes.
- A339661 into distinct squarefree semiprimes.
- A339742 into distinct primes or squarefree semiprimes.
- A339839 into distinct primes or semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A321728 is conjectured to count non-half-loop-graphical partitions of n.
- A339617 counts non-graphical partitions of 2n, ranked by A339618.
- A339655 counts non-loop-graphical partitions of 2n (A339657).

Programs

  • Maple
    filter:= proc(n)
      g(map(t -> t[2], ifactors(n)[2]))
    end proc;
    g:= proc(L) option remember; local x,i,j,t,s,Cons,R;
      if nops(L) = 1 then return L[1] > 3
      elif nops(L) = 2 then return max(L) > 4
      fi;
      Cons:= {seq(x[i] + x[i,i] + add(x[j,i], j=1..i-1)
         + add(x[i,j],j=i+1..nops(L)) = L[i], i=1..nops(L))};
      R:= traperror(Optimization:-LPSolve(0,Cons, assume=binary));
      type(R,string)
    end proc:
    select(filter, [$2..2000]); # Robert Israel, Dec 28 2020
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[1000],Select[facs[#],UnsameQ@@#&&SubsetQ[{1,2},PrimeOmega/@#]&]=={}&]

A002121 a(0) = 1, a(1) = 0, a(2) = -1; for n >= 3, a(n) = - a(n-2) + Sum_{ primes p with 3 <= p <= n} a(n-p).

Original entry on oeis.org

1, 0, -1, 1, 1, -1, 0, 2, 0, -2, 2, 4, -3, -2, 8, 1, -8, 8, 12, -11, -4, 25, 4, -24, 21, 40, -31, -16, 82, 14, -81, 71, 131, -99, -48, 258, 46, -249, 223, 422, -303, -162, 825, 169, -791, 714, 1360, -955, -503, 2641, 573, -2479, 2263, 4365, -2941, -1592, 8436, 1978, -7830, 7212, 14083, -9133, -4992, 26970, 6688, -24590
Offset: 0

Views

Author

Keywords

Comments

Arises in studying the Goldbach conjecture.
The last negative term appears to be a(303). - T. D. Noe, Dec 05 2006

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    import Data.List (genericIndex)
    a002121 n = genericIndex a002121_list n
    a002121_list = 1 : 0 : -1 : f 0 (-1) 3 where
       f v w x = y : f w y (x + 1) where
         y = sum (map (a002121 . (x -)) $ takeWhile (<= x) a065091_list) - v
    -- Reinhard Zumkeller, Mar 21 2014
  • Mathematica
    CoefficientList[Series[1/(1+Sum[(-x)^Prime[k],{k,20}]),{x,0,70}],x] (* Harvey P. Dale, Aug 26 2020 *)

Formula

G.f.: 1/(1+Sum_{k>0} (-x)^prime(k)). - Vladeta Jovovic, Mar 29 2003

Extensions

More terms from Vladeta Jovovic, Mar 29 2003
Entry revised by N. J. A. Sloane, Dec 04 2006

A339887 Number of factorizations of n into primes or squarefree semiprimes.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 4, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 4, 1, 2, 2, 4, 1, 3, 1, 2, 2, 2, 2, 4, 1, 2, 1, 2, 1, 5, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Dec 22 2020

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
Conjecture: also the number of semistandard Young tableaux whose entries are the prime indices of n (A323437).
Is this a duplicate of A323437? - R. J. Mathar, Jan 05 2021

Examples

			The a(n) factorizations for n = 36, 60, 180, 360, 420, 840:
  6*6       6*10      5*6*6       6*6*10        2*6*35      6*10*14
  2*3*6     2*5*6     2*6*15      2*5*6*6       5*6*14      2*2*6*35
  2*2*3*3   2*2*15    3*6*10      2*2*6*15      6*7*10      2*5*6*14
            2*3*10    2*3*5*6     2*3*6*10      2*10*21     2*6*7*10
            2*2*3*5   2*2*3*15    2*2*3*5*6     2*14*15     2*2*10*21
                      2*3*3*10    2*2*2*3*15    2*5*6*7     2*2*14*15
                      2*2*3*3*5   2*2*3*3*10    3*10*14     2*2*5*6*7
                                  2*2*2*3*3*5   2*2*3*35    2*3*10*14
                                                2*2*5*21    2*2*2*3*35
                                                2*2*7*15    2*2*2*5*21
                                                2*3*5*14    2*2*2*7*15
                                                2*3*7*10    2*2*3*5*14
                                                2*2*3*5*7   2*2*3*7*10
                                                            2*2*2*3*5*7
		

Crossrefs

See link for additional cross-references.
Only allowing only primes gives A008966.
Not allowing primes gives A320656.
Unlabeled multiset partitions of this type are counted by A320663/A339888.
Allowing squares of primes gives A320732.
The strict version is A339742.
A001055 counts factorizations.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.

Programs

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

Formula

a(A002110(n)) = A000085(n), and in general if n is a product of k distinct primes, a(n) = A000085(k).
a(n) = Sum_{d|n} A320656(n/d), so A320656 is the Moebius transform of this sequence.
Previous Showing 11-20 of 26 results. Next