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 71-80 of 485 results. Next

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

A339617 Number of non-graphical integer partitions of 2n.

Original entry on oeis.org

0, 1, 3, 6, 13, 25, 46, 81, 141, 234, 383, 615, 968, 1503, 2298, 3468, 5176, 7653, 11178, 16212, 23290, 33218, 46996, 66091, 92277, 128122, 176787, 242674, 331338, 450279, 608832, 819748, 1098907, 1467122, 1951020, 2584796, 3411998
Offset: 0

Views

Author

Gus Wiseman, Dec 13 2020

Keywords

Comments

An integer partition is graphical if it comprises the multiset of vertex-degrees of some graph. See A209816 for multigraphical partitions, A000070 for non-multigraphical partitions. Graphical partitions are counted by A000569.
The following are equivalent characteristics for any positive integer n:
(1) the prime indices 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 a(1) = 1 through a(4) = 13 partitions:
  (2)  (4)    (6)      (8)
       (2,2)  (3,3)    (4,4)
       (3,1)  (4,2)    (5,3)
              (5,1)    (6,2)
              (3,2,1)  (7,1)
              (4,1,1)  (3,3,2)
                       (4,2,2)
                       (4,3,1)
                       (5,2,1)
                       (6,1,1)
                       (3,3,1,1)
                       (4,2,1,1)
                       (5,1,1,1)
For example, the partition (2,2,2,2) is not counted under a(4) because there are three possible graphs with the prescribed degrees:
  {{1,2},{1,3},{2,4},{3,4}}
  {{1,2},{1,4},{2,3},{3,4}}
  {{1,3},{1,4},{2,3},{2,4}}
		

Crossrefs

A006881 lists squarefree semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339659 counts graphical partitions of 2n into k parts.
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).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).
- A339617 [this sequence] 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
    prptns[m_]:=Union[Sort/@If[Length[m]==0,{{}},Join@@Table[Prepend[#,m[[ipr]]]&/@prptns[Delete[m,List/@ipr]],{ipr,Select[Prepend[{#},1]&/@Select[Range[2,Length[m]],m[[#]]>m[[#-1]]&],UnsameQ@@m[[#]]&]}]]];
    strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
    Table[Length[Select[strnorm[2*n],Select[prptns[#],UnsameQ@@#&]=={}&]],{n,0,5}]

Formula

a(n) + A000569(n) = A000041(2*n).

A338915 Number of integer partitions of n that have an even number of parts and cannot be partitioned into distinct pairs of not necessarily distinct parts.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 1, 1, 4, 2, 6, 6, 12, 12, 20, 22, 38, 42, 60, 73, 101, 124, 164, 203, 266, 319, 415, 507, 649, 786, 983, 1198, 1499, 1797, 2234, 2673, 3303, 3952, 4826, 5753, 6999, 8330, 10051, 11943, 14357, 16956, 20322, 23997, 28568, 33657, 39897, 46879
Offset: 0

Views

Author

Gus Wiseman, Dec 10 2020

Keywords

Comments

The multiplicities of such a partition form a non-loop-graphical partition (A339655, A339657).

Examples

			The a(7) = 1 through a(12) = 12 partitions:
  211111  2222      411111    222211      222221      3333
          221111    21111111  331111      611111      222222
          311111              511111      22211111    441111
          11111111            22111111    32111111    711111
                              31111111    41111111    22221111
                              1111111111  2111111111  32211111
                                                      33111111
                                                      42111111
                                                      51111111
                                                      2211111111
                                                      3111111111
                                                      111111111111
For example, the partition y = (3,2,2,1,1,1,1,1) can be partitioned into pairs in just three ways:
  {{1,1},{1,1},{1,2},{2,3}}
  {{1,1},{1,1},{1,3},{2,2}}
  {{1,1},{1,2},{1,2},{1,3}}
None of these is strict, so y is counted under a(12).
		

Crossrefs

The Heinz numbers of these partitions are A320892.
The complement in even-length partitions is A338916.
A000070 counts non-multigraphical partitions of 2n, ranked by A339620.
A000569 counts graphical partitions, ranked by A320922.
A001358 lists semiprimes, with squarefree case A006881.
A058696 counts partitions of even numbers, ranked by A300061.
A209816 counts multigraphical partitions, ranked by A320924.
A320655 counts factorizations into semiprimes.
A322353 counts factorizations into distinct 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.
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).
- 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
    smcs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[smcs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],PrimeOmega[#]==2&]}]];
    Table[Length[Select[IntegerPartitions[n],EvenQ[Length[#]]&&smcs[Times@@Prime/@#]=={}&]],{n,0,10}]

Formula

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

Extensions

More terms from Jinyuan Wang, Feb 14 2025

A339618 Heinz numbers of non-graphical integer partitions of even numbers.

Original entry on oeis.org

3, 7, 9, 10, 13, 19, 21, 22, 25, 28, 29, 30, 34, 37, 39, 43, 46, 49, 52, 53, 55, 57, 61, 62, 63, 66, 70, 71, 75, 76, 79, 82, 84, 85, 87, 88, 89, 91, 94, 100, 101, 102, 107, 111, 113, 115, 116, 117, 118, 121, 129, 130, 131, 133, 134, 136, 138, 139, 146, 147
Offset: 1

Views

Author

Gus Wiseman, Dec 18 2020

Keywords

Comments

An integer partition is graphical if it comprises the multiset of vertex-degrees of some graph. Graphical partitions are counted by A000569.
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 distinct strict pairs (a set of edges);
(2) n can be factored into distinct squarefree semiprimes;
(3) the unordered prime signature of n is graphical.

Examples

			The sequence of terms together with their prime indices begins:
      3: {2}         43: {14}        79: {22}
      7: {4}         46: {1,9}       82: {1,13}
      9: {2,2}       49: {4,4}       84: {1,1,2,4}
     10: {1,3}       52: {1,1,6}     85: {3,7}
     13: {6}         53: {16}        87: {2,10}
     19: {8}         55: {3,5}       88: {1,1,1,5}
     21: {2,4}       57: {2,8}       89: {24}
     22: {1,5}       61: {18}        91: {4,6}
     25: {3,3}       62: {1,11}      94: {1,15}
     28: {1,1,4}     63: {2,2,4}    100: {1,1,3,3}
     29: {10}        66: {1,2,5}    101: {26}
     30: {1,2,3}     70: {1,3,4}    102: {1,2,7}
     34: {1,7}       71: {20}       107: {28}
     37: {12}        75: {2,3,3}    111: {2,12}
     39: {2,6}       76: {1,1,8}    113: {30}
For example, there are three possible multigraphs with degrees (1,1,3,3):
  {{1,2},{1,2},{1,2},{3,4}}
  {{1,2},{1,2},{1,3},{2,4}}
  {{1,2},{1,2},{1,4},{2,3}}.
Since none of these is a graph, the Heinz number 100 belongs to the sequence.
		

Crossrefs

A181819 applied to A320894 gives this sequence.
A300061 is a superset.
A339617 counts these partitions.
A320922 ranks the complement, counted by A000569.
A006881 lists squarefree semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339659 counts graphical partitions of 2n into k parts.
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).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).
- A339617 counts non-graphical partitions of 2n (A339618 [this sequence]).
- 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
    strs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strs[n/d],Min@@#>d&]],{d,Select[Divisors[n],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    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[#]]]&&strs[Times@@Prime/@nrmptn[#]]=={}&]

Formula

Equals A300061 \ A320922.
For all n, A181821(a(n)) and A304660(a(n)) belong to A320894.

A056595 Number of nonsquare divisors of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 3, 1, 2, 1, 3, 1, 4, 1, 3, 3, 2, 1, 4, 1, 4, 3, 3, 1, 6, 1, 3, 2, 4, 1, 7, 1, 3, 3, 3, 3, 5, 1, 3, 3, 6, 1, 7, 1, 4, 4, 3, 1, 7, 1, 4, 3, 4, 1, 6, 3, 6, 3, 3, 1, 10, 1, 3, 4, 3, 3, 7, 1, 4, 3, 7, 1, 8, 1, 3, 4, 4, 3, 7, 1, 7, 2, 3, 1, 10, 3, 3, 3, 6, 1, 10, 3, 4, 3, 3, 3, 9, 1, 4, 4, 5, 1, 7, 1
Offset: 1

Views

Author

Labos Elemer, Jul 21 2000

Keywords

Comments

a(A000430(n))=1; a(A030078(n))=2; a(A030514(n))=2; a(A006881(n))=3; a(A050997(n))=3; a(A030516(n))=3; a(A054753(n))=4; a(A000290(n))=A055205(n). - Reinhard Zumkeller, Aug 15 2011

Examples

			a(36)=5 because the set of divisors of 36 has tau(36)=nine elements, {1, 2, 3, 4, 6, 9, 12, 18, 36}, five of which, that is {2, 3, 6, 12, 18}, are not perfect squares.
		

Crossrefs

See A194095 and A194096 for record values and where they occur.

Programs

Formula

a(n) = A000005(n) - A046951(n) = tau(n) - tau(A000188(n)).
Sum_{k=1..n} a(k) ~ n*log(n) + (2*gamma - zeta(2) - 1)*n, where gamma is Euler's constant (A001620). - Amiram Eldar, Dec 01 2023

A065119 Numbers k such that the k-th cyclotomic polynomial is a trinomial.

Original entry on oeis.org

3, 6, 9, 12, 18, 24, 27, 36, 48, 54, 72, 81, 96, 108, 144, 162, 192, 216, 243, 288, 324, 384, 432, 486, 576, 648, 729, 768, 864, 972, 1152, 1296, 1458, 1536, 1728, 1944, 2187, 2304, 2592, 2916, 3072, 3456, 3888, 4374, 4608, 5184, 5832, 6144, 6561, 6912, 7776, 8748, 9216
Offset: 1

Views

Author

Len Smiley, Nov 12 2001

Keywords

Comments

Appears to be numbers of form 2^a * 3^b, a >= 0, b > 0. - Lekraj Beedassy, Sep 10 2004
This is true: see link "Cyclotomic trinomials". - Robert Israel, Jul 14 2015
3-smooth numbers (A003586) which are not powers of 2 (A000079). - Amiram Eldar, Nov 10 2020
These are the conjugates of semiprimes, where conjugation is A122111; or Heinz numbers of conjugates of length-2 partitions. - Gus Wiseman, Nov 09 2023
A multiplicative semigroup: if m and n are in the sequence, then so is m*n. - Antti Karttunen, Jul 13 2024

Examples

			The 54th cyclotomic polynomial is x^18 - x^9 + 1 which is trinomial, so 54 is in the sequence.
From _Gus Wiseman_, Nov 09 2023: (Start)
The terms and conjugate semiprimes, showing their respective Heinz partitions, begin:
    3: (2)              4: (1,1)
    6: (2,1)            6: (2,1)
    9: (2,2)            9: (2,2)
   12: (2,1,1)         10: (3,1)
   18: (2,2,1)         15: (3,2)
   24: (2,1,1,1)       14: (4,1)
   27: (2,2,2)         25: (3,3)
   36: (2,2,1,1)       21: (4,2)
   48: (2,1,1,1,1)     22: (5,1)
   54: (2,2,2,1)       35: (4,3)
   72: (2,2,1,1,1)     33: (5,2)
   81: (2,2,2,2)       49: (4,4)
   96: (2,1,1,1,1,1)   26: (6,1)
(End)
		

References

  • Jean-Marie De Koninck and Armel Mercier, 1001 Problèmes en Théorie Classique Des Nombres, Problem 733, pp. 74 and 310, Ellipses Paris, 2004.

Crossrefs

Differs at the 18th term from A063996.
For primes (A008578) we have conjugates A000079.
For triprimes (A014612) we have conjugates A080193.
A001358 lists semiprimes, squarefree A006881, complement A100959.

Programs

  • Maple
    with(numtheory): a := []; for m from 1 to 3000 do if nops([coeffs(cyclotomic(m,x))])=3 then a := [op(a),m] fi od; print(a);
  • Mathematica
    max = 5000; Sort[Flatten[Table[2^a 3^b, {a, 0, Floor[Log[2, max]]}, {b, Floor[Log[3, max/2^a]]}]]] (* Alonso del Arte, May 19 2016 *)
  • PARI
    isok(n)=my(vp = Vec(polcyclo(n))); sum(k=1, #vp, vp[k] != 0) == 3; \\ Michel Marcus, Jul 11 2015
    
  • PARI
    list(lim)=my(v=List(),N); for(n=1,logint(lim\1,3), N=3^n; while(N<=lim, listput(v,N); N<<=1)); Set(v) \\ Charles R Greathouse IV, Aug 07 2015

Formula

A206787(a(n)) = 4. - Reinhard Zumkeller, Feb 12 2012
a(n) = A033845(n)/2 = 3 * A003586(n). - Robert Israel, Jul 14 2015
Sum_{n>=1} 1/a(n) = 1. - Amiram Eldar, Nov 10 2020

Extensions

Offset set to 1 and more terms from Michel Marcus, Jul 11 2015

A068443 Triangular numbers which are the product of two primes.

Original entry on oeis.org

6, 10, 15, 21, 55, 91, 253, 703, 1081, 1711, 1891, 2701, 3403, 5671, 12403, 13861, 15931, 18721, 25651, 34453, 38503, 49141, 60031, 64261, 73153, 79003, 88831, 104653, 108811, 114481, 126253, 146611, 158203, 171991, 188191, 218791, 226801, 258121, 269011
Offset: 1

Views

Author

Stephan Wagler (stephanwagler(AT)aol.com), Mar 09 2002

Keywords

Comments

These triangular numbers are equal to p * (2p +- 1).
All terms belong to A006987. For n>2 all terms are odd and belong to A095147. - Alexander Adamchuk, Oct 31 2006
A156592 is a subsequence. - Reinhard Zumkeller, Feb 10 2009
Triangular numbers with exactly 4 divisors. - Jon E. Schoenfield, Sep 05 2018

Examples

			Triangular numbers begin 0, 1, 3, 6, 10, ...; 6=2*3, and 2 and 3 are two distinct primes; 10=2*5, and 2 and 5 are two distinct primes, etc. - _Vladimir Joseph Stephan Orlovsky_, Feb 27 2009
a(11) = 1891 and 1891 = 31 * 61.
		

Crossrefs

Programs

  • Maple
    q:= n-> is(numtheory[bigomega](n)=2):
    select(q, [i*(i+1)/2$i=0..1000])[];  # Alois P. Heinz, Mar 27 2024
  • Mathematica
    Select[ Table[ n(n + 1)/2, {n, 1000}], Apply[Plus, Transpose[ FactorInteger[ # ]] [[2]]] == 2 &]
    Select[Accumulate[Range[1000]],PrimeOmega[#]==2&] (* Harvey P. Dale, Apr 03 2016 *)
  • PARI
    list(lim)=my(v=List());forprime(p=2,(sqrtint(lim\1*8+1)+1)\4, if(isprime(2*p-1),listput(v,2*p^2-p)); if(isprime(2*p+1), listput(v,2*p^2+p))); Vec(v) \\ Charles R Greathouse IV, Jun 13 2013

Formula

A010054(a(n))*A064911(a(n)) = 1. - Reinhard Zumkeller, Dec 03 2009
a(n) = A000217(A164977(n)). - Zak Seidov, Feb 16 2015

Extensions

Edited by Robert G. Wilson v, Jul 08 2002
Definition corrected by Zak Seidov, Mar 09 2008

A090076 a(n) = prime(n)*prime(n+2).

Original entry on oeis.org

10, 21, 55, 91, 187, 247, 391, 551, 713, 1073, 1271, 1591, 1927, 2279, 2773, 3233, 3953, 4331, 4891, 5609, 6059, 7031, 8051, 8989, 9991, 10807, 11227, 12091, 13843, 14803, 17399, 18209, 20413, 20989, 23393, 24613, 26219, 28199, 29893, 31313
Offset: 1

Views

Author

Felix Tubiana, Jan 21 2004

Keywords

Comments

Subsequence of A192133. - Reinhard Zumkeller, Jun 26 2011
For n > 1: A078898(a(n)) = 4. - Reinhard Zumkeller, Apr 06 2015

Examples

			a(5) = prime(5)*prime(7) = 11*17 = 187.
		

Crossrefs

Subset of the squarefree semiprimes, A006881. Cf. A006094, A090090.
Cf. A078898.

Programs

  • Haskell
    a090076 n = a090076_list !! (n-1)
    a090076_list = zipWith (*) a000040_list $ drop 2 a000040_list
    -- Reinhard Zumkeller, Dec 17 2014
  • Mathematica
    Table[Prime[n] Prime[n + 2], {n, 1, 40}] (* Robert G. Wilson v, Jan 22 2004 *)
    Last[#]First[#]&/@Partition[Prime[Range[50]],3,1] (* Harvey P. Dale, May 08 2013 *)
  • MuPAD
    ithprime(i)*ithprime(i+2) $ i = 1..40 // Zerinvary Lajos, Feb 26 2007
    
  • Sage
    def prime_gaps(n):
        primegaps = []
        nprimes = primes_first_n(n+1)
        for i in range(2, n+1):
            primegaps.append(nprimes[i]*nprimes[i-2])
        return primegaps
    print(prime_gaps(60)) # Zerinvary Lajos, Jul 08 2008
    

Extensions

Extended by Robert G. Wilson v, Jan 22 2004

A112141 Product of the first n semiprimes.

Original entry on oeis.org

4, 24, 216, 2160, 30240, 453600, 9525600, 209563200, 5239080000, 136216080000, 4495130640000, 152834441760000, 5349205461600000, 203269807540800000, 7927522494091200000, 364666034728195200000, 17868635701681564800000, 911300420785759804800000
Offset: 1

Views

Author

Jonathan Vos Post, Nov 28 2005

Keywords

Comments

Semiprime analog of primorial (A002110). Equivalent for product of what A062198 is for sum.

Examples

			a(10) = 4*6*9*10*14*15*21*22*25*26 = 136216080000, the product of the first 10 semiprimes.
From _Gus Wiseman_, Dec 06 2020: (Start)
The sequence of terms together with their prime signatures begins:
                        4: (2)
                       24: (3,1)
                      216: (3,3)
                     2160: (4,3,1)
                    30240: (5,3,1,1)
                   453600: (5,4,2,1)
                  9525600: (5,5,2,2)
                209563200: (6,5,2,2,1)
               5239080000: (6,5,4,2,1)
             136216080000: (7,5,4,2,1,1)
            4495130640000: (7,6,4,2,2,1)
          152834441760000: (8,6,4,2,2,1,1)
         5349205461600000: (8,6,5,3,2,1,1)
       203269807540800000: (9,6,5,3,2,1,1,1)
      7927522494091200000: (9,7,5,3,2,2,1,1)
    364666034728195200000: (10,7,5,3,2,2,1,1,1)
  17868635701681564800000: (10,7,5,5,2,2,1,1,1)
(End)
		

Crossrefs

Partial sums of semiprimes are A062198.
First differences of semiprimes are A065516.
A000040 lists primes, with partial products A002110 (primorials).
A000142 lists factorials, with partial products A000178 (superfactorials).
A001358 lists semiprimes, with partial products A112141 (this sequence).
A005117 lists squarefree numbers, with partial products A111059.
A006881 lists squarefree semiprimes, with partial products A339191.
A101048 counts partitions into semiprimes (restricted: A338902).
A320655 counts factorizations into semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes, with product/sum/difference A087794/A176504/A176506.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.

Programs

  • Maple
    A112141 := proc(n)
        mul(A001358(i),i=1..n) ;
    end proc:
    seq(A112141(n),n=1..10) ; # R. J. Mathar, Jun 30 2020
  • Mathematica
    NextSemiPrime[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sp = n + sgn; While[c < Abs[k], While[ PrimeOmega[sp] != 2, If[sgn < 0, sp--, sp++]]; If[sgn < 0, sp--, sp++]; c++]; sp + If[sgn < 0, 1, -1]]; f[n_] := Times @@ NestList[ NextSemiPrime@# &, 2^2, n - 1]; Array[f, 18] (* Robert G. Wilson v, Jun 13 2013 *)
    FoldList[Times,Select[Range[30],PrimeOmega[#]==2&]] (* Gus Wiseman, Dec 06 2020 *)
  • PARI
    a(n)=my(v=vector(n),i,k=3);while(iCharles R Greathouse IV, Apr 04 2013
    
  • Python
    from sympy import factorint
    def aupton(terms):
        alst, k, p = [], 1, 1
        while len(alst) < terms:
            if sum(factorint(k).values()) == 2:
                p *= k
                alst.append(p)
            k += 1
        return alst
    print(aupton(18)) # Michael S. Branicky, Aug 31 2021

Formula

a(n) = Product_{i=1..n} A001358(i).
A001222(a(n)) = 2*n.

A179643 Products of exactly 2 distinct squares of primes and a different prime (p^2 * q^2 * r).

Original entry on oeis.org

180, 252, 300, 396, 450, 468, 588, 612, 684, 700, 828, 882, 980, 1044, 1100, 1116, 1300, 1332, 1452, 1476, 1548, 1575, 1692, 1700, 1900, 1908, 2028, 2124, 2156, 2178, 2196, 2205, 2300, 2412, 2420, 2450, 2475, 2548, 2556, 2628, 2844, 2900, 2925, 2988
Offset: 1

Views

Author

Keywords

Comments

A050326(a(n)) = 5, subsequence of A225228. - Reinhard Zumkeller, May 03 2013

Examples

			180 = 2^2 * 3^2 * 5, 252 = 2^2 * 3^2 * 7, 300 = 2^2 * 3 * 5^2, ...
		

Crossrefs

Programs

  • Mathematica
    f[n_]:=Sort[Last/@FactorInteger[n]]=={1,2,2}; Select[Range[3000], f]
  • PARI
    list(lim)=my(v=List(),t);forprime(p=2,sqrt(lim\12),forprime(q=p+1,sqrt(lim\p^2\2),t=(p*q)^2;forprime(r=2,lim\t,if(p==r||q==r,next);listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 19 2011
    
  • Python
    from math import isqrt
    from sympy import primepi, primerange, integer_nthroot
    def A179643(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((t:=primepi(s:=isqrt(y:=isqrt(x//r))))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(x+1))+sum(primepi(isqrt(x//p**3)) for p in primerange(integer_nthroot(x,3)[0]+1))-primepi(integer_nthroot(x,5)[0])
        return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025
Previous Showing 71-80 of 485 results. Next