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 99 results. Next

A001055 The multiplicative partition function: number of ways of factoring n with all factors greater than 1 (a(1) = 1 by convention).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 5, 1, 7, 2, 2, 2, 9, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 11, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 5, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 9, 1, 5, 1
Offset: 1

Views

Author

Keywords

Comments

From David W. Wilson, Feb 28 2009: (Start)
By a factorization of n we mean a multiset of integers > 1 whose product is n.
For example, 6 is the product of 2 such multisets, {2, 3} and {6}, so a(6) = 2.
Similarly 8 is the product of 3 such multisets, {2, 2, 2}, {2, 4} and {8}, so a(8) = 3.
1 is the product of 1 such multiset, namely the empty multiset {}, whose product is by definition the multiplicative identity 1. Hence a(1) = 1. (End)
a(n) = # { k | A064553(k) = n }. - Reinhard Zumkeller, Sep 21 2001; Benoit Cloitre and N. J. A. Sloane, May 15 2002
Number of members of A025487 with n divisors. - Matthew Vandermast, Jul 12 2004
See sequence A162247 for a list of the factorizations of n and a program for generating the factorizations for any n. - T. D. Noe, Jun 28 2009
So a(n) gives the number of different prime signatures that can be found among the integers that have n divisors. - Michel Marcus, Nov 11 2015
For n > 0, also the number of integer partitions of n with product n, ranked by A301987. For example, the a(12) = 4 partitions are: (12), (6,2,1,1,1,1), (4,3,1,1,1,1,1), (3,2,2,1,1,1,1,1). See also A380218. In general, A379666(m,n) = a(n) for any m >= n. - Gus Wiseman, Feb 07 2025

Examples

			1: 1, a(1) = 1
2: 2, a(2) = 1
3: 3, a(3) = 1
4: 4 = 2*2, a(4) = 2
6: 6 = 2*3, a(6) = 2
8: 8 = 2*4 = 2*2*2, a(8) = 3
etc.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 292-295.
  • Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4.
  • 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).
  • 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

A045782 gives the range of a(n).
For records see A033833, A033834.
Row sums of A316439 (for n>1).
Cf. A096276 (partial sums).
The additive version is A000041 (integer partitions), strict A000009.
Row sums of A318950.
A002865 counts partitions into parts > 1.
A069016 counts distinct sums of factorizations.
A319000 counts partitions by product and sum, row sums A319916.
A379666 (transpose A380959) counts partitions by sum and product, without 1's A379668, strict A379671.

Programs

  • Haskell
    a001055 = (map last a066032_tabl !!) . (subtract 1)
    -- Reinhard Zumkeller, Oct 01 2012
    
  • Java
    public class MultiPart {
        public static void main(String[] argV) {
            for (int i=1;i<=100;++i) System.out.println(1+getDivisors(2,i));
        }
        public static int getDivisors(int min,int n) {
            int total = 0;
            for (int i=min;i=i) { ++total; if (n/i>i) total+=getDivisors(i,n/i); }
            return total;
        }
    } \\ Scott R. Shannon, Aug 21 2019
  • Maple
    with(numtheory):
    T := proc(n::integer, m::integer)
            local A, summe, d:
            if isprime(n) then
                    if n <= m then
                            return 1;
                    end if:
                    return 0 ;
            end if:
            A := divisors(n) minus {n, 1}:
            for d in A do
                    if d > m then
                            A := A minus {d}:
                    end if:
            end do:
            summe := add(T(n/d,d),d=A) ;
            if n <=m then
                    summe := summe + 1:
            end if:
            summe ;
    end proc:
    A001055 := n -> T(n, n):
    [seq(A001055(n), n=1..100)]; # Reinhard Zumkeller and Ulrich Schimke (ulrschimke(AT)aol.com)
  • 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]; a/@Range[100] (* c[n, r] is the number of factorizations of n with factors <= r. - Dean Hickerson, Oct 28 2002 *)
    T[, 1] = T[1, ] = 1;
    T[n_, m_] := T[n, m] = DivisorSum[n, Boole[1 < # <= m] * T[n/#, #]&];
    a[n_] := T[n, n];
    a /@ Range[100] (* Jean-François Alcover, Jan 03 2020 *)
  • PARI
    /* factorizations of n with factors <= m (n,m positive integers) */
    fcnt(n,m) = {local(s);s=0;if(n == 1,s=1,fordiv(n,d,if(d > 1 & d <= m,s=s+fcnt(n/d,d))));s}
    A001055(n) = fcnt(n,n) \\ Michael B. Porter, Oct 29 2009
    
  • PARI
    \\ code using Dirichlet g.f., based on Somos's code for A007896
    {a(n) = my(A, v, w, m);
    if(
    n<1, 0,
    \\ define unit vector v = [1, 0, 0, ...] of length n
    v = vector(n, k, k==1);
    for(k=2, n,
    m = #digits(n, k) - 1;
    \\ expand 1/(1-x)^k out far enough
    A = (1 - x)^ -1 + x * O(x^m);
    \\ w = zero vector of length n
    w = vector(n);
    \\ convert A to a vector
    for(i=0, m, w[k^i] = polcoeff(A, i));
    \\ build the answer
    v = dirmul(v, w)
    );
    v[n]
    )
    };
    \\ produce the sequence
    vector(100,n,a(n)) \\ N. J. A. Sloane, May 26 2014
    
  • PARI
    v=vector(100, k, k==1); for(n=2, #v, v+=dirmul(v, vector(#v, k, (k>1) && n^valuation(k,n)==k)) ); v \\ Max Alekseyev, Jul 16 2014
    
  • Python
    from sympy import divisors, isprime
    def T(n, m):
        if isprime(n): return 1 if n<=m else 0
        A=filter(lambda d: d<=m, divisors(n)[1:-1])
        s=sum(T(n//d, d) for d in A)
        return s + 1 if n<=m else s
    def a(n): return T(n, n)
    print([a(n) for n in range(1, 106)]) # Indranil Ghosh, Aug 19 2017
    

Formula

The asymptotic behavior of this sequence was studied by Canfield, Erdős & Pomerance and Luca, Mukhopadhyay, & Srinivas. - Jonathan Vos Post, Jul 07 2008
Dirichlet g.f.: Product_{k>=2} 1/(1 - 1/k^s).
If n = p^k for a prime p, a(n) = partitions(k) = A000041(k).
Since the sequence a(n) is the right diagonal of A066032, the given recursive formula for A066032 applies (see Maple program). - Reinhard Zumkeller and Ulrich Schimke (ulrschimke(AT)aol.com)
a(A002110(n)) = A000110(n).
a(p^k*q^k) = A002774(k) if p and q are distinct primes. - R. J. Mathar, Jun 06 2024
a(n) = A028422(n) + 1. - Gus Wiseman, Feb 07 2025

Extensions

Incorrect assertion about asymptotic behavior deleted by N. J. A. Sloane, Jun 08 2009

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

A325044 Heinz numbers of integer partitions whose sum of parts is greater than or equal to their product.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 22, 23, 24, 26, 28, 29, 30, 31, 32, 34, 36, 37, 38, 40, 41, 43, 44, 46, 47, 48, 52, 53, 56, 58, 59, 60, 61, 62, 64, 67, 68, 71, 72, 73, 74, 76, 79, 80, 82, 83, 84, 86, 88, 89, 92, 94, 96, 97, 101
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 or equal to their sum of prime indices (A056239).
The enumeration of these partitions by sum is given by A096276.

Examples

			The sequence of terms together with their prime indices begins:
   2: {1}
   3: {2}
   4: {1,1}
   5: {3}
   6: {1,2}
   7: {4}
   8: {1,1,1}
   9: {2,2}
  10: {1,3}
  11: {5}
  12: {1,1,2}
  13: {6}
  14: {1,4}
  16: {1,1,1,1}
  17: {7}
  18: {1,2,2}
  19: {8}
  20: {1,1,3}
  22: {1,5}
  23: {9}
  24: {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[#]<=Plus@@primeMS[#]&]

Formula

A003963(a(n)) <= A056239(a(n)).
a(n) = A325038(n)/2.
Union of A301987 and A325038.

A319000 Regular triangle where T(n,k) is the number of finite multisets of positive integers with product n and sum k.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 22 2018

Keywords

Examples

			Triangle begins:
  1
  0 1
  0 0 1
  0 0 0 2
  0 0 0 0 1
  0 0 0 0 1 2
  0 0 0 0 0 0 1
  0 0 0 0 0 2 2 3
  0 0 0 0 0 1 1 1 2
  0 0 0 0 0 0 1 1 1 2
  0 0 0 0 0 0 0 0 0 0 1
  0 0 0 0 0 0 2 3 3 3 3 4
  0 0 0 0 0 0 0 0 0 0 0 0 1
  0 0 0 0 0 0 0 0 1 1 1 1 1 2
  0 0 0 0 0 0 0 1 1 1 1 1 1 1 2
  0 0 0 0 0 0 0 3 3 4 4 4 4 4 4 5
  0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
  0 0 0 0 0 0 0 1 2 2 3 3 3 3 3 3 3 4
  0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
  0 0 0 0 0 0 0 0 2 2 2 3 3 3 3 3 3 3 3 4
Row 12 {0,0,0,0,0,0,2,3,3,3,3,4} corresponds to the partitions (C = 12):
. . . . . . (43)  (62)   (621)   (6211)   (62111)    (C)
            (322) (431)  (4311)  (43111)  (431111)   (621111)
                  (3221) (32211) (322111) (3221111)  (4311111)
                                          (32211111)
		

Crossrefs

Row sums are A319916. Column sums are A319005. Last column is A001055.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[k],Times@@#==n&]],{n,20},{k,n}]

A319005 Number of integer partitions of n whose product of parts is >= n.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 5, 7, 13, 18, 28, 40, 60, 80, 113, 152, 205, 266, 353, 454, 590, 751, 959, 1210, 1529, 1905, 2381, 2953, 3658, 4501, 5539, 6772, 8278, 10065, 12230, 14801, 17893, 21544, 25921, 31089, 37240, 44478, 53068, 63150, 75063, 89018, 105438, 124632
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2018

Keywords

Examples

			The a(1) = 1 through a(9) = 18 partitions:
  (1)  (2)  (3)  (4)   (5)   (6)    (7)     (8)      (9)
                 (22)  (32)  (33)   (43)    (44)     (54)
                             (42)   (52)    (53)     (63)
                             (222)  (322)   (62)     (72)
                             (321)  (331)   (332)    (333)
                                    (421)   (422)    (432)
                                    (2221)  (431)    (441)
                                            (521)    (522)
                                            (2222)   (531)
                                            (3221)   (621)
                                            (3311)   (3222)
                                            (4211)   (3321)
                                            (22211)  (4221)
                                                     (4311)
                                                     (5211)
                                                     (22221)
                                                     (32211)
                                                     (33111)
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, p) option remember; `if`(n=0 or i=1, `if`(p>1,
          0, 1), b(n, i-1, p) +b(n-i, min(i, n-i), max(p/i, 1)))
        end:
    a:= n-> b(n$3):
    seq(a(n), n=0..50);  # Alois P. Heinz, Oct 22 2018
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Times@@#>=n&]],{n,50}]
    (* Second program: *)
    b[n_, i_, p_] := b[n, i, p] = If[n == 0 || i == 1, If[p > 1, 0, 1],
         b[n, i - 1, p] + b[n - i, Min[i, n - i], Max[p/i, 1]]];
    a[n_] := b[n, n, n];
    a /@ Range[0, 50] (* Jean-François Alcover, May 11 2021, after Alois P. Heinz *)

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[#]]&]
Showing 1-10 of 99 results. Next