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 41-50 of 50 results.

A380692 Numbers k such that the least prime dividing k is larger than the maximum exponent in the prime factorization of k; a(1) = 1 by convention.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 99
Offset: 1

Views

Author

Amiram Eldar, Jan 30 2025

Keywords

Comments

First differs from A368110 at n = 68: A368110(68) = 98 is not a term of this sequence.
First differs from its subsequence A380695 at n = 428: a(428) = 625 is not a term of A380695.
Numbers k such that A020639(k) > A051903(k).
All the squarefree numbers (A005117) are terms, and all the odd terms of A067259 are terms of this sequence.
Disjoint union of the sequences S_k, k >= 1, where S_k is the sequence of p-rough numbers (numbers whose prime factors are all greater than or equal to p), with p = nextprime(k) = A151800(k), whose maximum exponent in their prime factorization is k (i.e., numbers that are (k+1)-free but not k-free, where k-free numbers are numbers whose prime factorization exponents do not exceed k).
The asymptotic density of this sequence is Sum_{i>=1} d(i) = 0.68213349032332767778..., where d(i), the density of S_i, equals f(i+1) * Product_{primes p <= i} ((1-1/p)/(1-1/p^(i+1))) - f(i) * Product_{primes p <= i} ((1-1/p)/(1-1/p^i)), f(i) = 1/zeta(i) if i >= 2, and f(1) = 0.

Examples

			6 = 2^1 * 3^1 is a term since 2 > 1.
8 = 2^3 is not a term since 2 < 3.
		

Crossrefs

Subsequence of A380693.
Subsequences: A005117, A136327, the intersection of A005408 and A067259, A380694, A380695.

Programs

  • Mathematica
    q[k_] := k == 1 || Module[{f = FactorInteger[k]}, f[[1, 1]] > Max[f[[;; , 2]]]]; Select[Range[100], q]
  • PARI
    isok(k) = if(k == 1, 0, my(f = factor(k), e = f[, 2]); f[1, 1] > vecmax(e));

A380695 Numbers k such that the prime index of the least prime dividing k is larger than or equal to the maximum exponent in the prime factorization of k; a(1) = 1 by convention.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 99
Offset: 1

Views

Author

Amiram Eldar, Jan 30 2025

Keywords

Comments

First differs from A368110 at n = 68: A368110(68) = 98 is not a term of this sequence.
Subsequence of A380692 and first differs from it at n = 428: A380692(428) = 625 is not a term of this sequence.
Numbers k such that A055396(k) >= A051903(k).
Disjoint union of the sequences S_k, k >= 1, where S_k is the sequence of p-rough numbers (numbers whose prime factors are all greater than or equal to p), with p = prime(k), whose maximum exponent in their prime factorization is k (i.e., numbers that are (k+1)-free but not k-free, where k-free numbers are numbers whose prime factorization exponents do not exceed k).
The asymptotic density of this sequence is Sum_{i>=1} d(i) = 0.68165919742420048618..., where d(i), the density of S_i, equals f(i+1) * Product_{primes p < prime(i)} ((1-1/p)/(1-1/p^(i+1))) - f(i) * Product_{primes p < prime(i)} ((1-1/p)/(1-1/p^i)), f(i) = 1/zeta(i) if i >= 2, and f(1) = 0.

Examples

			2 = 2^1 is a term since PrimePi(2) = 1 >= 1.
4 = 2^2 is not a term since PrimePi(2) = 1 is smaller than the exponent 2.
25 = 5^2 is a term since PrimePi(5) = 3 >= 2.
		

Crossrefs

Subsequence of A380692 and A380693.
A380694 is a subsequence.

Programs

  • Mathematica
    q[k_] := k == 1 || Module[{f = FactorInteger[k]}, f[[1, 1]] >= Prime[Max[f[[;; , 2]]]]]; Select[Range[100], q]
  • PARI
    isok(k) = if(k == 1, 1, my(f = factor(k), e = f[,2]); f[1,1] >= prime(vecmax(e)));

A387112 Numbers with (strictly) choosable initial intervals of prime indices.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2025

Keywords

Comments

First differs from A371088 in having a(86) = 121.
The initial interval of a nonnegative integer x is the set {1,...,x}.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We say that a set or sequence of nonempty sets is choosable iff it is possible to choose a different element from each set. For example, ({1,2},{1},{1,3}) is choosable because we have the choice (2,1,3), but ({1,2,3},{1},{1,3},{2}) is not.
This sequence lists all numbers k such that if the prime indices of k are (x1,x2,...,xz), then the sequence of sets (initial intervals) ({1,...,x1},{1,...,x2},...,{1,...,xz}) is choosable.

Examples

			The prime indices of 85 are {3,7}, with initial intervals {{1,2,3},{1,2,3,4,5,6,7}}, which are choosable, so 85 is in the sequence
The prime indices of 90 are {1,2,2,3}, with initial intervals {{1},{1,2},{1,2},{1,2,3}}, which are not choosable, so 90 is not in the sequence.
		

Crossrefs

Partitions of this type are counted by A238873, complement A387118.
For partitions instead of initial intervals we have A276078, complement A276079.
For prime factors instead of initial intervals we have A368100, complement A355529.
For divisors instead of initial intervals we have A368110, complement A355740.
These are all the positions of nonzero terms in A387111, complement A387134.
The complement is A387113.
For strict partitions instead of initial intervals we have A387176, complement A387137.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A120383 lists numbers divisible by all of their prime indices.
A367902 counts choosable set-systems, complement A367903.
A370582 counts sets with choosable prime factors, complement A370583.
A370585 counts maximal subsets with choosable prime factors.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Select[Tuples[Range/@prix[#]],UnsameQ@@#&]!={}&]

A387176 Numbers whose prime indices do not have choosable sets of strict integer partitions. Zeros of A387115.

Original entry on oeis.org

4, 8, 9, 12, 16, 18, 20, 24, 27, 28, 32, 36, 40, 44, 45, 48, 52, 54, 56, 60, 63, 64, 68, 72, 76, 80, 81, 84, 88, 90, 92, 96, 99, 100, 104, 108, 112, 116, 117, 120, 124, 125, 126, 128, 132, 135, 136, 140, 144, 148, 152, 153, 156, 160, 162, 164, 168, 171, 172
Offset: 1

Views

Author

Gus Wiseman, Aug 27 2025

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We say that a sequence of nonempty sets is choosable iff it is possible to choose a different element from each set. For example, ({1,2},{1},{1,3}) is choosable because we have the choice (2,1,3), but ({1},{2},{1,3},{2,3}) is not.

Crossrefs

The complement for all partitions appears to be A276078, counted by A052335.
For all partitions we appear to have A276079, counted by A387134.
For divisors instead of strict partitions we have A355740, counted by A370320.
Twice-partitions of this type (into distinct strict partitions) are counted by A358914.
The complement for divisors is A368110, counted by A239312.
The complement for initial intervals is A387112, counted by A238873, see A387111.
For initial intervals instead of strict partitions we have A387113, counted by A387118.
These are the positions of 0 in A387115.
Partitions of this type are counted by A387137, complement A387178.
The complement is A387177.
The version for constant partitions is A387180, counted by A387329.
The complement for constant partitions is A387181, counted by A387330.
A000041 counts integer partitions, strict A000009.
A003963 multiplies together the prime indices of n.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Select[Tuples[Select[IntegerPartitions[#],UnsameQ@@#&]&/@prix[#]],UnsameQ@@#&]=={}&]

A370815 Number of integer factorizations of n into unordered factors > 1, such that only one set can be obtained by choosing a different divisor of each factor.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 06 2024

Keywords

Examples

			The a(432) = 3 factorizations: (2*2*3*4*9), (2*3*3*4*6), (2*6*6*6).
		

Crossrefs

For partitions and prime factors we have A370594, ranks A370647.
Partitions of this type are counted by A370595, ranks A370810.
For prime factors we have A370645, subsets A370584.
A000005 counts divisors.
A001055 counts factorizations, strict A045778.
A239312 counts condensed partitions, ranks A355740, complement A370320.
A355731 counts choices of a divisor of each prime index, firsts A355732.
A368414 counts factor-choosable factorizations, complement A368413.
A370814 counts divisor-choosable factorizations, complement A370813.

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],Length[Union[Sort /@ Select[Tuples[Divisors/@#],UnsameQ@@#&]]]==1&]],{n,100}]

A387113 Numbers whose prime indices do not have (strictly) choosable initial intervals.

Original entry on oeis.org

4, 8, 12, 16, 18, 20, 24, 27, 28, 32, 36, 40, 44, 48, 52, 54, 56, 60, 64, 68, 72, 76, 80, 81, 84, 88, 90, 92, 96, 100, 104, 108, 112, 116, 120, 124, 126, 128, 132, 135, 136, 140, 144, 148, 150, 152, 156, 160, 162, 164, 168, 172, 176, 180, 184, 188, 189, 192
Offset: 1

Views

Author

Gus Wiseman, Aug 24 2025

Keywords

Comments

The initial interval of a nonnegative integer x is the set {1,...,x}.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We say that a set or sequence of nonempty sets is choosable iff it is possible to choose a different element from each set. For example, ({1,2},{1},{1,3}) is choosable because we have the choice (2,1,3), but ({1,2,3},{1},{1,3},{2}) is not.
This sequence lists all numbers k such that if the prime indices of k are (x1,x2,...,xz), then the sequence of sets (initial intervals) ({1,...,x1},{1,...,x2},...,{1,...,xz}) is not choosable.

Examples

			The prime indices of 18 are {1,2,2}, with initial intervals ({1},{1,2},{1,2}), which have choices (1,1,1), (1,1,2), (1,2,1), (1,2,2), and since none of these are strict, 18 is in the sequence.
The prime indices of 85 are {3,7}, with initial intervals {{1,2,3},{1,2,3,4,5,6,7}}, which are choosable, so 85 is in not the sequence.
The prime indices of 90 are {1,2,2,3}, with initial intervals {{1},{1,2},{1,2},{1,2,3}}, which are not choosable, so 90 is in the sequence.
The terms together with their prime indices begin:
    4: {1,1}
    8: {1,1,1}
   12: {1,1,2}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   24: {1,1,1,2}
   27: {2,2,2}
   28: {1,1,4}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   40: {1,1,1,3}
   44: {1,1,5}
   48: {1,1,1,1,2}
   52: {1,1,6}
   54: {1,2,2,2}
   56: {1,1,1,4}
   60: {1,1,2,3}
   64: {1,1,1,1,1,1}
		

Crossrefs

For partitions instead of initial intervals we have A276079, complement A276078.
For prime factors instead of initial intervals we have A355529, complement A368100.
For divisors instead of initial intervals we have A355740, complement A368110.
These are the positions of 0 in A387111, complement A387134.
The complement is A387112.
Partitions of this type are counted by A387118, complement A238873.
For strict partitions instead of initial intervals we have A387137, complement A387176.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A120383 lists numbers divisible by all of their prime indices.
A367902 counts choosable set-systems, complement A367903.
A370582 counts sets with choosable prime factors, complement A370583.
A370585 counts maximal subsets with choosable prime factors.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Select[Tuples[Range/@prix[#]],UnsameQ@@#&]=={}&]

A387177 Numbers whose prime indices have choosable sets of strict integer partitions. Positions of nonzero terms in A387115.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 49, 50, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 98
Offset: 1

Views

Author

Gus Wiseman, Aug 29 2025

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We say that a sequence of nonempty sets is choosable iff it is possible to choose a different element from each set. For example, ({1,2},{1},{1,3}) is choosable because we have the choice (2,1,3), but ({1},{2},{1,3},{2,3}) is not.

Examples

			The prime indices of 50 are {1,3,3}, and {(1),(3),(2,1)} is a valid choice of distinct strict partitions, so 50 is in the sequence.
		

Crossrefs

The version for all partitions appears to be A276078, counted by A052335.
The complement for all partitions appears to be A276079, counted by A387134.
The complement for divisors is A355740, counted by A370320.
Twice-partitions of this type (into distinct strict partitions) are counted by A358914.
The version for divisors is A368110, counted by A239312.
The version for initial intervals is A387112, counted by A238873, see A387111.
The complement for initial intervals is A387113, counted by A387118.
These are the positions of nonzero terms in A387115.
The complement is A387176.
Partitions of this type are counted by A387178, complement A387137.
The complement for constant partitions is A387180, counted by A387329, see A387120.
The version for constant partitions is A387181, counted by A387330.
A000041 counts integer partitions, strict A000009.
A003963 multiplies together the prime indices of n.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A289509 lists numbers with relatively prime prime indices.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    strptns[n_]:=Select[IntegerPartitions[n],UnsameQ@@#&];
    Select[Range[100],Select[Tuples[strptns/@prix[#]],UnsameQ@@#&]!={}&]

A387178 Number of integer partitions of n whose parts have choosable sets of strict integer partitions.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 5, 6, 8, 10, 13, 17, 21, 27, 34, 42, 53, 65, 80, 98, 119, 146, 177, 213, 258, 309, 370, 443, 528, 628, 745, 882, 1043, 1229, 1447, 1700, 1993, 2333, 2727, 3182, 3707, 4311, 5008, 5808, 6727, 7782, 8990, 10371, 11952, 13756, 15815, 18161
Offset: 0

Views

Author

Gus Wiseman, Sep 02 2025

Keywords

Comments

First differs from A052337 in having 745 instead of 746.
We say that a sequence of nonempty sets is choosable iff it is possible to choose a different element from each set. For example, ({1,2},{1},{1,3}) is choosable because we have the choice (2,1,3), but ({1},{2},{1,3},{2,3}) is not.
a(n) is the number of integer partitions of n such that it is possible to choose a sequence of distinct strict integer partitions of each part.
Also the number of integer partitions of n with no part k whose multiplicity exceeds A000009(k).

Examples

			The partition y = (3,3,2) has sets of strict integer partitions ({(2,1),(3)},{(2,1),(3)},{(2)}), and we have the choice ((2,1),(3),(2)) or ((3),(2,1),(2)), so y is counted under a(8).
The a(1) = 1 through a(9) = 10 partitions:
  (1)  (2)  (3)    (4)    (5)    (6)      (7)      (8)      (9)
            (2,1)  (3,1)  (3,2)  (3,3)    (4,3)    (4,4)    (5,4)
                          (4,1)  (4,2)    (5,2)    (5,3)    (6,3)
                                 (5,1)    (6,1)    (6,2)    (7,2)
                                 (3,2,1)  (3,3,1)  (7,1)    (8,1)
                                          (4,2,1)  (3,3,2)  (4,3,2)
                                                   (4,3,1)  (4,4,1)
                                                   (5,2,1)  (5,3,1)
                                                            (6,2,1)
                                                            (3,3,2,1)
		

Crossrefs

For initial intervals instead of strict partitions we have A238873, ranks A387112.
For divisors instead of strict partitions we have A239312, ranks A368110.
The complement for divisors is A370320, ranks A355740.
For prime factors instead of strict partitions we have A370592, ranks A368100.
The complement for prime factors is A370593, ranks A355529.
The complement for initial intervals is A387118, ranks A387113.
The complement for all partitions is A387134, ranks A387577.
The complement is counted by A387137, ranks A387176.
These partitions are ranked by A387177.
For all partitions instead of just strict partitions we have A387328, ranks A387576.
The complement for constant partitions is A387329, ranks A387180.
For constant partitions instead of strict partitions we have A387330, ranks A387181.
A000041 counts integer partitions, strict A000009.
A358914 counts twice-partitions into distinct strict partitions.
A367902 counts choosable set-systems, complement A367903.

Programs

  • Mathematica
    strptns[n_]:=Select[IntegerPartitions[n],UnsameQ@@#&];
    Table[Length[Select[IntegerPartitions[n],Select[Tuples[strptns/@#],UnsameQ@@#&]!={}&]],{n,0,15}]

A387180 Numbers of which it is not possible to choose a different constant integer partition of each prime index.

Original entry on oeis.org

4, 8, 12, 16, 20, 24, 27, 28, 32, 36, 40, 44, 48, 52, 54, 56, 60, 64, 68, 72, 76, 80, 81, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 125, 128, 132, 135, 136, 140, 144, 148, 152, 156, 160, 162, 164, 168, 172, 176, 180, 184, 188, 189, 192, 196, 200, 204
Offset: 1

Views

Author

Gus Wiseman, Aug 30 2025

Keywords

Comments

First differs from A276079 in having 125.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Also numbers n with at least one prime index k such that the multiplicity of prime(k) in the prime factorization of n exceeds the number of divisors of k.

Examples

			The prime indices of 60 are {1,1,2,3}, and we have the following 4 choices of constant partitions:
  ((1),(1),(2),(3))
  ((1),(1),(2),(1,1,1))
  ((1),(1),(1,1),(3))
  ((1),(1),(1,1),(1,1,1))
Since none of these is strict, 60 is in the sequence.
The prime indices of 90 are {1,2,2,3}, and we have the following 4 strict choices:
  ((1),(2),(1,1),(3))
  ((1),(2),(1,1),(1,1,1))
  ((1),(1,1),(2),(3))
  ((1),(1,1),(2),(1,1,1))
So 90 is not in the sequence.
		

Crossrefs

For prime factors instead of constant partitions we have A355529, counted by A370593.
For divisors instead of constant partitions we have A355740, counted by A370320.
The complement for prime factors is A368100, counted by A370592.
The complement for divisors is A368110, counted by A239312.
The complement for initial intervals is A387112, counted by A238873.
For initial intervals instead of partitions we have A387113, counted by A387118.
These are the positions of zero in A387120.
For strict instead of constant partitions we have A387176, counted by A387137.
The complement for strict partitions is A387177, counted by A387178.
Twice-partitions of this type are counted by A387179, constant-block case of A296122.
The complement is A387181 (nonzeros of A387120), counted by A387330.
Partitions of this type are counted by A387329.
A000041 counts integer partitions, strict A000009.
A003963 multiplies together prime indices.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Select[Tuples[Select[IntegerPartitions[#],SameQ@@#&]&/@prix[#]],UnsameQ@@#&]=={}&]

A370406 Primitive terms of A370348.

Original entry on oeis.org

4, 18, 27, 50, 125, 225, 242, 294, 441, 578, 686, 1029, 1089, 1331, 1922, 2401, 2601, 3025, 3362, 3675, 4913, 5070, 5290, 6962, 7225, 7605, 8575, 8649, 8978, 12675, 13182, 13225, 13778, 15129, 15162, 17787, 19773, 21970, 22743, 23762, 23805, 24025, 24334, 29791, 31329, 32258, 32955, 34969, 35378
Offset: 1

Views

Author

Robert Israel, Feb 17 2024

Keywords

Comments

Terms of A370348 that are not divisible by any other term of A370348.
Numbers k such that there are fewer divisors of prime indices of k than there are prime indices of k, and no proper divisor of k has this property.

Examples

			a(4) = 50 is a term because the prime indices of 50 = 2*5^2 are 1, 2, 2, and there are 3 of these but only 2 divisors of prime indices, namely 1 and 2, and 50 is not divisible by any of the previous terms 4, 18 and 27 of the sequence.
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) uses numtheory; local F,D,t;
      if ormap(t -> n mod t = 0, S) then return false fi;
    F:= map(t -> [pi(t[1]), t[2]], ifactors(n)[2]);
    D:= `union`(seq(divisors(t[1]), t = F);
    nops(D) < add(t[2], t = F);
    end proc:
    R:= NULL: count:= 0: S:= {}:
    for n from 1 while count < 100 do
     if filter(n) then
       R:= R, n; S:= S union {n}; count:= count+1;
     fi
    od:
    R;
  • Mathematica
    filter[n_] := Module[{F, d},
       If[AnyTrue[S, Mod[n, #] == 0&], Return[False]];
       F = {PrimePi[#[[1]]], #[[2]]} & /@ FactorInteger[n];
       d = Union[Flatten[Divisors /@ F[[All, 1]]]];
       Length[d] < Total[F[[All, 2]]]];
    R = {}; count = 0; S = {};
    For[n = 1, count < 100, n++, If[filter[n], AppendTo[R, n]; S = Union[S, {n}]; count++]];
    R (* Jean-François Alcover, Mar 08 2024, after Robert Israel *)
Previous Showing 41-50 of 50 results.