cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 62 results. Next

A325037 Heinz numbers of integer partitions whose product of parts is greater than their sum.

Original entry on oeis.org

1, 15, 21, 25, 27, 33, 35, 39, 42, 45, 49, 50, 51, 54, 55, 57, 63, 65, 66, 69, 70, 75, 77, 78, 81, 85, 87, 90, 91, 93, 95, 98, 99, 100, 102, 105, 110, 111, 114, 115, 117, 119, 121, 123, 125, 126, 129, 130, 132, 133, 135, 138, 140, 141, 143, 145, 147, 150, 153
Offset: 1

Views

Author

Gus Wiseman, Mar 25 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1) * ... * prime(y_k), so these are numbers whose product of prime indices (A003963) is greater than their sum of prime indices (A056239).
The enumeration of these partitions by sum is given by A114324.

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
  15: {2,3}
  21: {2,4}
  25: {3,3}
  27: {2,2,2}
  33: {2,5}
  35: {3,4}
  39: {2,6}
  42: {1,2,4}
  45: {2,2,3}
  49: {4,4}
  50: {1,3,3}
  51: {2,7}
  54: {1,2,2,2}
  55: {3,5}
  57: {2,8}
  63: {2,2,4}
  65: {3,6}
  66: {1,2,5}
  69: {2,9}
  70: {1,3,4}
  75: {2,3,3}
  77: {4,5}
  78: {1,2,6}
  81: {2,2,2,2}
		

Crossrefs

Programs

  • Maple
    q:= n-> (l-> mul(i, i=l)>add(i, i=l))(map(i->
        numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    select(q, [$1..200])[];  # Alois P. Heinz, Mar 27 2019
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Times@@primeMS[#]>Plus@@primeMS[#]&]

Formula

A003963(a(n)) > A056239(a(n)).

A096276 Number of partitions of n with a product <=n.

Original entry on oeis.org

0, 1, 2, 3, 5, 6, 8, 9, 12, 14, 16, 17, 21, 22, 24, 26, 31, 32, 36, 37, 41, 43, 45, 46, 53, 55, 57, 60, 64, 65, 70, 71, 78, 80, 82, 84, 93, 94, 96, 98, 105, 106, 111, 112, 116, 120, 122, 123, 135, 137, 141, 143, 147, 148, 155, 157, 164, 166, 168, 169, 180, 181, 183, 187
Offset: 0

Views

Author

Jon Perry, Jun 23 2004

Keywords

Comments

The Heinz numbers of these partitions are given by A325044. - Gus Wiseman, Mar 27 2019

Examples

			a(6)=8 as we can have 6, 51, 411, 321, 3111, 2211, 21111, 111111, rejecting 42, 33 and 222.
From _Gus Wiseman_, Mar 27 2019: (Start)
The a(1) = 1 through a(8) = 12 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (21)   (22)    (41)     (51)      (61)       (71)
             (111)  (31)    (221)    (321)     (511)      (611)
                    (211)   (311)    (411)     (3211)     (4211)
                    (1111)  (2111)   (2211)    (4111)     (5111)
                            (11111)  (3111)    (22111)    (22211)
                                     (21111)   (31111)    (32111)
                                     (111111)  (211111)   (41111)
                                               (1111111)  (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
(End)
		

References

  • G. Tenenbaum, Introduction to analytic and probabilistic number theory, Cambridge University Press, 1995, p. 198, exercise 9 (in the third edition 2015, p. 296, exercise 211).

Crossrefs

Programs

  • Maple
    g:= proc(n, k) option remember; `if`(n>k, 0, 1)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= proc(n) option remember; `if`(n=0, 0, a(n-1)+g(n$2)) end:
    seq(a(n), n=0..100);  # Alois P. Heinz, Feb 26 2023
  • Mathematica
    c[1, r_] := c[1, r] = 1; c[n_, r_] := c[n, r] = Module[{ds, i}, ds = Select[Divisors[n], 1 < # <= r &]; Sum[c[n/ds[[i]], ds[[i]]], {i, 1, Length[ds]}]]; a[n_] := c[n, n]; Join[{0}, Accumulate[Array[a, 100]]] (* using program from A001055, T. D. Noe, Apr 11 2011 *)
    Table[Length[Select[IntegerPartitions[n],Times@@#<=n&]],{n,0,20}] (* Gus Wiseman, Mar 27 2019 *)
  • PARI
    { bla(n,m,v,z)=v=concat(v,m); if(!n,x=prod(k=1,length(v),v[k]); if (x<=z,c++), for(i=1,min(m,n),bla(n-i,i,v,z))); }
    q(n)=c=0;for(i=1,n,bla(n-i,i,[],n));print1(c, ", ");
    for(i=0,40,q(i))

Formula

For n>1, a(n) = a(n-1)+1 iff n is prime.
Partial sums of A001055. - Vladeta Jovovic, Jun 24 2004
a(n) ~ n * exp(2*sqrt(log(n))) / (2*sqrt(Pi) * (log(n))^(3/4)) [Oppenheim, 1927]. - Vaclav Kotesovec, May 23 2020

Extensions

More terms from Vladeta Jovovic, Jun 24 2004

A326149 Numbers whose product of prime indices is divisible by their sum of prime indices.

Original entry on oeis.org

2, 3, 5, 7, 9, 11, 13, 17, 19, 23, 29, 30, 31, 37, 41, 43, 47, 49, 53, 59, 61, 63, 65, 67, 71, 73, 79, 81, 83, 84, 89, 97, 101, 103, 107, 108, 109, 113, 125, 127, 131, 137, 139, 149, 150, 151, 154, 157, 163, 165, 167, 169, 173, 179, 181, 190, 191, 193, 197
Offset: 1

Views

Author

Gus Wiseman, Jun 09 2019

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.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions whose product of parts is divisible by their sum of parts. The enumeration of these partitions by sum is given by A057568.

Examples

			The sequence of terms together with their prime indices begins:
    2: {1}
    3: {2}
    5: {3}
    7: {4}
    9: {2,2}
   11: {5}
   13: {6}
   17: {7}
   19: {8}
   23: {9}
   29: {10}
   30: {1,2,3}
   31: {11}
   37: {12}
   41: {13}
   43: {14}
   47: {15}
   49: {4,4}
   53: {16}
   59: {17}
		

Crossrefs

Satisfies A056239(a(n))|A003963(a(n)).
The nonprime case is A326150, with squarefree case A326158.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[2,100],Divisible[Times@@primeMS[#],Plus@@primeMS[#]]&]

A325038 Heinz numbers of integer partitions whose sum of parts is greater than their product.

Original entry on oeis.org

4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 32, 34, 36, 38, 40, 44, 46, 48, 52, 56, 58, 60, 62, 64, 68, 72, 74, 76, 80, 82, 86, 88, 92, 94, 96, 104, 106, 112, 116, 118, 120, 122, 124, 128, 134, 136, 142, 144, 146, 148, 152, 158, 160, 164, 166, 168, 172
Offset: 1

Views

Author

Gus Wiseman, Mar 25 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1) * ... * prime(y_k), so these are numbers whose product of prime indices (A003963) is less than their sum of prime indices (A056239).
The enumeration of these partitions by sum is given by A096276 shifted once to the right.

Examples

			The sequence of terms together with their prime indices begins:
   4: {1,1}
   6: {1,2}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  14: {1,4}
  16: {1,1,1,1}
  18: {1,2,2}
  20: {1,1,3}
  22: {1,5}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
  32: {1,1,1,1,1}
  34: {1,7}
  36: {1,1,2,2}
  38: {1,8}
  40: {1,1,1,3}
  44: {1,1,5}
  46: {1,9}
  48: {1,1,1,1,2}
		

Crossrefs

Programs

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

Formula

A003963(a(n)) < A056239(a(n)).
a(n) = 2 * A325044(n).

A114324 Number of partitions of n with a product greater than n.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 3, 6, 10, 16, 26, 39, 56, 79, 111, 150, 200, 265, 349, 453, 586, 749, 957, 1209, 1522, 1903, 2379, 2950, 3654, 4500, 5534, 6771, 8271, 10063, 12228, 14799, 17884, 21543, 25919, 31087, 37233, 44477, 53063, 63149, 75059, 89014, 105436, 124631
Offset: 0

Views

Author

Giovanni Resta, Feb 06 2006

Keywords

Comments

The Heinz numbers of these partitions are given by A325037. - Gus Wiseman, Mar 27 2019

Examples

			a(6) = 3 since there are 3 partitions of 6 with product greater than 6: {3,3}, {2,2,2}, {4,2}.
From _Gus Wiseman_, Mar 27 2019: (Start)
The a(5) = 1 through a(9) = 16 partitions:
  (32)  (33)   (43)    (44)    (54)
        (42)   (52)    (53)    (63)
        (222)  (322)   (62)    (72)
               (331)   (332)   (333)
               (421)   (422)   (432)
               (2221)  (431)   (441)
                       (521)   (522)
                       (2222)  (531)
                       (3221)  (621)
                       (3311)  (3222)
                               (3321)
                               (4221)
                               (4311)
                               (5211)
                               (22221)
                               (32211)
(End)
		

Crossrefs

Programs

  • Mathematica
    << DiscreteMath`Combinatorica`; lst=Table[Length@Select[Partitions[n], (Times @@ # > n) &],{n,50}]
    Table[Length[Select[IntegerPartitions[n],Times@@#>n&]],{n,0,20}] (* Gus Wiseman, Mar 27 2019 *)

Extensions

a(0) = 1 prepended by Gus Wiseman, Mar 27 2019

A326155 Positive integers whose sum of prime indices is divisible by their product of prime indices.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 23, 29, 30, 31, 32, 37, 40, 41, 43, 47, 48, 53, 59, 61, 64, 67, 71, 73, 79, 83, 84, 89, 97, 101, 103, 107, 108, 109, 112, 113, 127, 128, 131, 137, 139, 144, 149, 151, 157, 163, 167, 173, 179, 181, 191, 192, 193
Offset: 1

Views

Author

Gus Wiseman, Jun 10 2019

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.
Also Heinz numbers of the integer partitions counted by A057567. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   4: {1,1}
   5: {3}
   7: {4}
   8: {1,1,1}
   9: {2,2}
  11: {5}
  12: {1,1,2}
  13: {6}
  16: {1,1,1,1}
  17: {7}
  19: {8}
  23: {9}
  29: {10}
  30: {1,2,3}
  31: {11}
  32: {1,1,1,1,1}
  37: {12}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Divisible[Plus@@primeMS[#],Times@@primeMS[#]]&]

A379721 Numbers whose prime indices have sum <= product.

Original entry on oeis.org

1, 2, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 30, 31, 33, 35, 37, 39, 41, 42, 43, 45, 47, 49, 50, 51, 53, 54, 55, 57, 59, 61, 63, 65, 66, 67, 69, 70, 71, 73, 75, 77, 78, 79, 81, 83, 84, 85, 87, 89, 90, 91, 93, 95, 97, 98, 99, 100, 101, 102, 103
Offset: 1

Views

Author

Gus Wiseman, Jan 05 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.
Partitions of this type are counted by A319005.
The complement is A325038.

Examples

			The terms together with their prime indices begin:
    1: {}
    2: {1}
    3: {2}
    5: {3}
    7: {4}
    9: {2,2}
   11: {5}
   13: {6}
   15: {2,3}
   17: {7}
   19: {8}
   21: {2,4}
   23: {9}
   25: {3,3}
   27: {2,2,2}
   29: {10}
   30: {1,2,3}
		

Crossrefs

The case of equality is A301987, inequality A325037.
Nonpositive positions in A325036.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A379681 gives sum plus product of prime indices, firsts A379682.
Counting and ranking multisets by comparing sum and product:
- same: A001055 (strict A045778), ranks A301987
- divisible: A057567, ranks A326155
- divisor: A057568, ranks A326149, see A326156, A326172, A379733
- greater: A096276 shifted right, ranks A325038
- greater or equal: A096276, ranks A325044
- less: A114324, ranks A325037, see A318029
- less or equal: A319005, ranks A379721 (this)
- different: A379736, ranks A379722, see A111133

Programs

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

Formula

Number k such that A056239(k) <= A003963(k).

A325036 Difference between product and sum of prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 25 2019

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.

Examples

			The prime indices of 45 are {2,2,3}, with product 12 and sum 7, so a(45) = 5.
		

Crossrefs

Positions of zeros are A301987. Positions of ones are A325041. Positions of negative ones are A325042.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Times@@primeMS[n]-Total[primeMS[n]],{n,100}]
    dps[n_]:=Module[{pi=Flatten[Table[PrimePi[#[[1]]],#[[2]]]&/@FactorInteger[n]]},Times@@pi-Total[pi]]; Join[{1},Array[dps,100,2]] (* Harvey P. Dale, May 26 2023 *)
  • PARI
    A003963(n) = { n=factor(n); n[, 1]=apply(primepi, n[, 1]); factorback(n) }; \\ From A003963
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1]))); };
    A325036(n) = (A003963(n) - A056239(n)); \\ Antti Karttunen, May 08 2022

Formula

a(n) = A003963(n) - A056239(n).
For all n >= 1, a(A325040(n)) = a(A122111(A325040(n))). - Antti Karttunen, May 08 2022

Extensions

Data section extended up to a(93) by Antti Karttunen, May 08 2022

A379722 Numbers whose prime indices do not have the same sum as product.

Original entry on oeis.org

1, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 81, 82, 85, 86, 87, 88, 90, 91, 92, 93
Offset: 1

Views

Author

Gus Wiseman, Jan 08 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.
Partitions of this type are counted by A379736.
The complement is A301987, counted by A001055.

Examples

			The terms together with their prime indices begin:
    1: {}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   14: {1,4}
   15: {2,3}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   21: {2,4}
   22: {1,5}
   24: {1,1,1,2}
   25: {3,3}
   26: {1,6}
   27: {2,2,2}
   28: {1,1,4}
   32: {1,1,1,1,1}
		

Crossrefs

Nonzeros of A325036.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A324851 finds numbers > 1 divisible by the sum of their prime indices.
A379666 counts partitions by sum and product, without 1's A379668.
A379681 gives sum plus product of prime indices, firsts A379682.
Counting and ranking multisets by comparing sum and product:
- same: A001055 (strict A045778), ranks A301987
- divisible: A057567, ranks A326155
- divisor: A057568, ranks A326149, see A379733
- greater: A096276 shifted right, ranks A325038
- greater or equal: A096276, ranks A325044
- less: A114324, ranks A325037, see A318029
- less or equal: A319005, ranks A379721
- different: A379736, ranks A379722 (this)

Programs

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

A340387 Numbers whose sum of prime indices is twice their number, counted with multiplicity in both cases.

Original entry on oeis.org

1, 3, 9, 10, 27, 28, 30, 81, 84, 88, 90, 100, 208, 243, 252, 264, 270, 280, 300, 544, 624, 729, 756, 784, 792, 810, 840, 880, 900, 1000, 1216, 1632, 1872, 2080, 2187, 2268, 2352, 2376, 2430, 2464, 2520, 2640, 2700, 2800, 2944, 3000, 3648, 4896, 5440, 5616
Offset: 1

Views

Author

Gus Wiseman, Jan 09 2021

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.
Also Heinz numbers of integer partitions whose sum is twice their length, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). Like partitions in general (A000041), these are also counted by A000041.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}
      3: {2}
      9: {2,2}
     10: {1,3}
     27: {2,2,2}
     28: {1,1,4}
     30: {1,2,3}
     81: {2,2,2,2}
     84: {1,1,2,4}
     88: {1,1,1,5}
     90: {1,2,2,3}
    100: {1,1,3,3}
    208: {1,1,1,1,6}
    243: {2,2,2,2,2}
    252: {1,1,2,2,4}
		

Crossrefs

Partitions of 2n into n parts are counted by A000041.
The number of prime indices alone is A001222.
The sum of prime indices alone is A056239.
Allowing sum to be any multiple of length gives A067538, ranked by A316413.
A000569 counts graphical partitions, ranked by A320922.
A027187 counts partitions of even length, ranked by A028260.
A058696 counts partitions of even numbers, ranked by A300061.
A301987 lists numbers whose sum of prime indices equals their product, with nonprime case A301988.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],Total[primeMS[#]]==2*PrimeOmega[#]&]

Formula

All terms satisfy A056239(a(n)) = 2*A001222(a(n)).
Showing 1-10 of 62 results. Next