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

A326836 Heinz numbers of integer partitions whose maximum part divides their sum.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 23, 25, 27, 29, 30, 31, 32, 36, 37, 40, 41, 43, 47, 48, 49, 53, 59, 61, 63, 64, 67, 70, 71, 73, 79, 81, 83, 84, 89, 97, 101, 103, 107, 108, 109, 112, 113, 121, 125, 127, 128, 131, 135, 137, 139, 144, 149, 150, 151
Offset: 1

Views

Author

Gus Wiseman, Jul 26 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 maximum prime index divides their sum of prime indices.
The enumeration of these partitions by sum is given by A067538.

Examples

			The sequence of terms together with their prime indices begins:
    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}
   25: {3,3}
   27: {2,2,2}
   29: {10}
   30: {1,2,3}
   31: {11}
   32: {1,1,1,1,1}
		

Crossrefs

Programs

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

A074761 Number of partitions of n of order n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 9, 1, 4, 5, 1, 1, 12, 1, 27, 7, 6, 1, 81, 1, 7, 1, 54, 1, 407, 1, 1, 11, 9, 13, 494, 1, 10, 13, 423, 1, 981, 1, 137, 115, 12, 1, 1309, 1, 59, 17, 193, 1, 240, 21, 1207, 19, 15, 1, 47274, 1, 16, 239, 1, 25, 3284, 1, 333, 23, 3731, 1, 42109, 1, 19
Offset: 1

Views

Author

Vladeta Jovovic, Sep 28 2002

Keywords

Comments

Order of partition is lcm of its parts.
a(n) is the number of conjugacy classes of the symmetric group S_n such that a representative of the class has order n. Here order means the order of an element of a group. Note that a(n) = 1 if and only if n is a prime power. - W. Edwin Clark, Aug 05 2014

Examples

			The a(15) = 5 partitions are (15), (5,3,3,3,1), (5,5,3,1,1), (5,3,3,1,1,1,1), (5,3,1,1,1,1,1,1,1). - _Gus Wiseman_, Aug 01 2018
		

Crossrefs

Programs

  • Maple
    A:= proc(n)
          uses numtheory;
          local S;
        S:= add(mobius(n/i)*1/mul(1-x^j,j=divisors(i)),i=divisors(n));
        coeff(series(S,x,n+1),x,n);
    end proc:
    seq(A(n),n=1..100); # Robert Israel, Aug 06 2014
  • Mathematica
    a[n_] := With[{s = Sum[MoebiusMu[n/i]*1/Product[1-x^j, {j, Divisors[i]}], {i, Divisors[n]}]}, SeriesCoefficient[s, {x, 0, n}]]; Array[a, 80] (* Jean-François Alcover, Feb 29 2016 *)
    Table[Length[Select[IntegerPartitions[n],LCM@@#==n&]],{n,50}] (* Gus Wiseman, Aug 01 2018 *)
  • PARI
    pr(k, x)={my(t=1); fordiv(k, d, t *= (1-x^d) ); return(t); }
    a(n) =
    {
        my( x = 'x+O('x^(n+1)) );
        polcoeff( Pol( sumdiv(n, i, moebius(n/i) / pr(i, x) ) ), n );
    }
    vector(66, n, a(n) )
    \\ Joerg Arndt, Aug 06 2014

Formula

Coefficient of x^n in expansion of Sum_{i divides n} A008683(n/i)*1/Product_{j divides i} (1-x^j).

A083710 Number of integer partitions of n with a part dividing all the other parts.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 11, 12, 20, 25, 37, 43, 70, 78, 114, 143, 196, 232, 330, 386, 530, 641, 836, 1003, 1340, 1581, 2037, 2461, 3127, 3719, 4746, 5605, 7038, 8394, 10376, 12327, 15272, 17978, 22024, 26095, 31730, 37339, 45333, 53175, 64100, 75340, 90138
Offset: 0

Views

Author

N. J. A. Sloane, Jun 16 2003

Keywords

Comments

Since the summand (part) which divides all the other summands is necessarily the smallest, an equivalent definition is: "Number of partitions of n such that smallest part divides every part." - Joerg Arndt, Jun 08 2009
The first few partitions that fail the criterion are 5=3+2, 7=5+2=4+3=3+2+2. So a(5) = A000041(5) - 1 = 6, a(7) = A000041(7) - 3 = 12. - Vladeta Jovovic, Jun 17 2003
Starting with offset 1 = inverse Mobius transform (A051731) of the partition numbers, A000041. - Gary W. Adamson, Jun 08 2009

Examples

			From _Gus Wiseman_, Apr 18 2021: (Start)
The a(1) = 1 through a(7) = 12 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (22)    (41)     (33)      (61)
             (111)  (31)    (221)    (42)      (331)
                    (211)   (311)    (51)      (421)
                    (1111)  (2111)   (222)     (511)
                            (11111)  (321)     (2221)
                                     (411)     (3211)
                                     (2211)    (4111)
                                     (3111)    (22111)
                                     (21111)   (31111)
                                     (111111)  (211111)
                                               (1111111)
(End)
		

References

  • L. M. Chawla, M. O. Levan and J. E. Maxfield, On a restricted partition function and its tables, J. Natur. Sci. and Math., 12 (1972), 95-101.

Crossrefs

Cf. A000041, A051731. - Gary W. Adamson, Jun 08 2009
The case with no 1's is A083711.
The strict case is A097986.
The version for "divisible by" instead of "dividing" is A130689.
The case where there is also a part divisible by all the others is A130714.
The complement of these partitions is counted by A338470.
The Heinz numbers of these partitions are dense, complement of A342193.
The case where there is also no part divisible by all the others is A343345.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Maple
    with(combinat): with(numtheory): a := proc(n) c := 0: l := sort(convert(divisors(n), list)): for i from 1 to nops(l)-0 do c := c+numbpart(l[i]-1) od: RETURN(c): end: for j from 0 to 60 do printf(`%d, `, a(j)) od: # Zerinvary Lajos, Apr 14 2007
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And@@IntegerQ/@(#/Min@@#)&]],{n,0,30}] (* Gus Wiseman, Apr 18 2021 *)

Formula

Equals left border of triangle A137587 starting (1, 2, 3, 5, 6, 11, ...). - Gary W. Adamson, Jan 27 2008
G.f.: 1 + Sum_{n>=1} x^n/eta(x^n). The g.f. for partitions into parts that are a multiple of n is x^n/eta(x^n), now sum over n. - Joerg Arndt, Jun 08 2009
Gary W. Adamson's comment is equivalent to the formula a(n) = Sum_{d|n} p(d-1) where p(i) = number of partitions of i (A000041(i)). Hence A083710 has g.f. Sum_{d>=1} p(d-1)*x^d/(1-x^d), - N. J. A. Sloane, Jun 08 2009

Extensions

More terms from Vladeta Jovovic, Jun 17 2003
Name shortened by Gus Wiseman, Apr 18 2021

A098743 Number of partitions of n into aliquant parts (i.e., parts that do not divide n).

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 3, 1, 3, 3, 13, 1, 23, 10, 11, 9, 65, 8, 104, 14, 56, 66, 252, 10, 245, 147, 206, 77, 846, 35, 1237, 166, 649, 634, 1078, 60, 3659, 1244, 1850, 236, 7244, 299, 10086, 1228, 1858, 4421, 19195, 243, 17660, 3244, 12268, 4039, 48341, 1819, 27675
Offset: 0

Views

Author

Reinhard Zumkeller, Oct 01 2004

Keywords

Comments

It seems very plausible that the low and high water marks occur when n is a factorial number or a prime: see A260797, A260798.
a(A000040(n)) = A002865(n) - 1.

Examples

			7 = 2 + 2 + 3 = 2 + 5 = 3 + 4, so a(7) = 3.
a(10) = #{7+3,6+4,4+3+3} = 3, all other partitions of 10 contain at least one divisor (10, 5, 2, or 1).
		

Crossrefs

See also A057562 (relatively prime parts).

Programs

  • Haskell
    a098743 n = p [nd | nd <- [1..n], mod n nd /= 0] n where
       p _  0 = 1
       p [] _ = 0
       p ks'@(k:ks) m | m < k = 0 | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Nov 22 2011
    
  • Haskell
    -- with memoization
    import Data.MemoCombinators (memo3, integral)
    a098743 n = a098743_list !! n
    a098743_list = map (\x -> pMemo x 1 x) [0..] where
       pMemo = memo3 integral integral integral p
       p   0 = 1
       p x k m | m < k        = 0
               | mod x k == 0 = pMemo x (k + 1) m
               | otherwise    = pMemo x k (m - k) + pMemo x (k + 1) m
    -- Reinhard Zumkeller, Aug 08 2015
    
  • Maple
    a := [1,0,0,0,0]; M:=300; for n from 5 to M do t1:={seq(i,i=1..n)}; t3 := t1 minus divisors(n); t4 := mul(1/(1-x^i), i in t3); t5 := series(t4,x,n+2); a:=[op(a), coeff(t5,x,n)]; od: a; # N. J. A. Sloane, Aug 08 2015
    # second Maple program:
    a:= proc(m) option remember; local b; b:= proc(n, i)
          option remember; `if`(n=0, 1, `if`(i<2, 0, b(n, i-1)+
          `if`(irem(m, i)=0, 0, b(n-i, min(i, n-i))))) end; b(m$2)
        end:
    seq(a(n), n=0..60);  # Alois P. Heinz, Mar 11 2018
  • Mathematica
    a[m_] := a[m] = Module[{b}, b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 2, 0, b[n, i-1] + If[Mod[m, i] == 0, 0, b[n-i, Min[i, n-i]]]]]; b[m, m]];
    Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Apr 30 2018, after Alois P. Heinz *)
  • PARI
    a(n)={polcoef(1/prod(k=1, n, if(n%k, 1 - x^k, 1) + O(x*x^n)), n)} \\ Andrew Howroyd, Aug 29 2018

Extensions

a(0) added and offset changed by Reinhard Zumkeller, Nov 22 2011
New wording for definition suggested by Marc LeBrun, Aug 07 2015

A353866 Heinz numbers of rucksack partitions. Every prime-power divisor has a different sum of prime indices.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 64, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75
Offset: 1

Views

Author

Gus Wiseman, Jun 06 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
In a knapsack partition (A108917, ranked by A299702), every submultiset has a different sum, so these are run-knapsack partitions or rucksack partitions for short.

Examples

			The terms together with their prime indices begin:
    1: {}
    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}
   13: {6}
   14: {1,4}
   15: {2,3}
   16: {1,1,1,1}
The sequence contains 18 because its prime-power divisors {1,2,3,9} have prime indices {}, {1}, {2}, {2,2} with distinct sums {0,1,2,4}. On the other hand, 12 is not in the sequence because {2} and {1,1} have the same sum.
		

Crossrefs

Knapsack partitions are counted by A108917, ranked by A299702.
The strong case is A353838, counted by A353837, complement A353839.
These partitions are counted by A353864.
The complete case is A353867, counted by A353865.
The complement is A354583.
A000041 counts partitions, strict A000009.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A073093 counts prime-power divisors.
A124010 gives prime signature, sorted A118914.
A300273 ranks collapsible partitions, counted by A275870.
A353832 represents the operation of taking run-sums of a partition.
A353836 counts partitions by number of distinct run-sums.
A353852 ranks compositions with all distinct run-sums, counted by A353850.
A353863 counts partitions whose weak run-sums cover an initial interval.

Programs

  • Mathematica
    msubs[s_]:=Join@@@Tuples[Table[Take[t,i],{t,Split[s]},{i,0,Length[t]}]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],UnsameQ@@Total/@Select[msubs[primeMS[#]],SameQ@@#&]&]

A365924 Number of incomplete integer partitions of n, meaning not every number from 0 to n is the sum of some submultiset.

Original entry on oeis.org

0, 0, 1, 1, 3, 3, 6, 7, 12, 14, 22, 25, 38, 46, 64, 76, 106, 124, 167, 199, 261, 309, 402, 471, 604, 714, 898, 1053, 1323, 1542, 1911, 2237, 2745, 3201, 3913, 4536, 5506, 6402, 7706, 8918, 10719, 12364, 14760, 17045, 20234, 23296, 27600, 31678, 37365, 42910, 50371, 57695, 67628, 77300, 90242, 103131, 119997
Offset: 0

Views

Author

Gus Wiseman, Sep 26 2023

Keywords

Comments

The complement (complete partitions) is A126796.

Examples

			The a(0) = 0 through a(8) = 12 partitions:
  .  .  (2)  (3)  (4)    (5)    (6)      (7)      (8)
                  (2,2)  (3,2)  (3,3)    (4,3)    (4,4)
                  (3,1)  (4,1)  (4,2)    (5,2)    (5,3)
                                (5,1)    (6,1)    (6,2)
                                (2,2,2)  (3,2,2)  (7,1)
                                (4,1,1)  (3,3,1)  (3,3,2)
                                         (5,1,1)  (4,2,2)
                                                  (4,3,1)
                                                  (5,2,1)
                                                  (6,1,1)
                                                  (2,2,2,2)
                                                  (5,1,1,1)
		

Crossrefs

For parts instead of sums we have A047967/A365919, ranks A080259/A055932.
The complement is A126796, ranks A325781, strict A188431.
These partitions have ranks A365830.
The strict case is A365831.
Row sums of A365923 without the first column, strict A365545.
A000041 counts integer partitions, strict A000009.
A046663 counts partitions w/o a submultiset summing to k, strict A365663.
A276024 counts positive subset-sums of partitions, strict A284640.
A325799 counts non-subset-sums of prime indices.
A364350 counts combination-free strict partitions.
A365543 counts partitions with a submultiset summing to k, strict A365661.

Programs

  • Mathematica
    nmz[y_]:=Complement[Range[Total[y]],Total/@Subsets[y]];
    Table[Length[Select[IntegerPartitions[n],Length[nmz[#]]>0&]],{n,0,15}]

Formula

a(n) = A000041(n) - A126796(n).

A348551 Heinz numbers of integer partitions whose mean is not an integer.

Original entry on oeis.org

1, 6, 12, 14, 15, 18, 20, 24, 26, 33, 35, 36, 38, 40, 42, 44, 45, 48, 50, 51, 52, 54, 56, 58, 60, 63, 65, 66, 69, 70, 72, 74, 75, 76, 77, 80, 86, 92, 93, 95, 96, 102, 104, 106, 108, 112, 114, 117, 119, 120, 122, 123, 124, 126, 130, 132, 135, 136, 140, 141, 142
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2021

Keywords

Comments

Equivalently, partitions whose length does not divide their sum.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The terms and their prime indices begin:
   1: {}
   6: {1,2}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  18: {1,2,2}
  20: {1,1,3}
  24: {1,1,1,2}
  26: {1,6}
  33: {2,5}
  35: {3,4}
  36: {1,1,2,2}
  38: {1,8}
  40: {1,1,1,3}
  42: {1,2,4}
  44: {1,1,5}
  45: {2,2,3}
  48: {1,1,1,1,2}
		

Crossrefs

A version counting nonempty subsets is A000079 - A051293.
A version counting factorizations is A001055 - A326622.
A version counting compositions is A011782 - A271654.
A version for prime factors is A175352, complement A078175.
A version for distinct prime factors A176587, complement A078174.
The complement is A316413, counted by A067538, strict A102627.
The geometric version is the complement of A326623.
The conjugate version is the complement of A326836.
These partitions are counted by A349156.
A000041 counts partitions.
A001222 counts prime factors with multiplicity.
A018818 counts partitions into divisors, ranked by A326841.
A143773 counts partitions into multiples of the length, ranked by A316428.
A236634 counts unbalanced partitions.
A047993 counts balanced partitions, ranked by A106529.
A056239 adds up prime indices, row sums of A112798.
A326567/A326568 gives the mean of prime indices, conjugate A326839/A326840.
A327472 counts partitions not containing their mean, complement A237984.

Programs

  • Maple
    q:= n-> (l-> nops(l)=0 or irem(add(i, i=l), nops(l))>0)(map
            (i-> numtheory[pi](i[1])$i[2], ifactors(n)[2])):
    select(q, [$1..142])[];  # Alois P. Heinz, Nov 19 2021
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],!IntegerQ[Mean[primeMS[#]]]&]

A097986 Number of strict integer partitions of n with a part dividing all the other parts.

Original entry on oeis.org

1, 1, 2, 2, 2, 4, 3, 5, 5, 7, 6, 12, 9, 13, 15, 20, 18, 28, 26, 37, 39, 47, 49, 71, 68, 85, 94, 117, 120, 159, 160, 201, 216, 257, 277, 348, 357, 430, 470, 562, 592, 720, 758, 901, 981, 1134, 1220, 1457, 1542, 1798, 1952, 2250, 2419, 2819, 3023, 3482, 3773, 4291
Offset: 1

Views

Author

Vladeta Jovovic, Oct 23 2004

Keywords

Comments

If n > 0, we can assume such a part is the smallest. - Gus Wiseman, Apr 23 2021
Also the number of uniform (constant multiplicity) partitions of n containing 1, ranked by A367586. The strict case is A096765. The version without 1 is A329436. - Gus Wiseman, Dec 01 2023

Examples

			From _Gus Wiseman_, Dec 01 2023: (Start)
The a(1) = 1 through a(8) = 5 strict partitions with a part dividing all the other parts:
  (1)  (2)  (3)    (4)    (5)    (6)      (7)      (8)
            (2,1)  (3,1)  (4,1)  (4,2)    (6,1)    (6,2)
                                 (5,1)    (4,2,1)  (7,1)
                                 (3,2,1)           (4,3,1)
                                                   (5,2,1)
The a(1) = 1 through a(8) = 5 uniform partitions containing 1:
  (1)  (11)  (21)   (31)    (41)     (51)      (61)       (71)
             (111)  (1111)  (11111)  (321)     (421)      (431)
                                     (2211)    (1111111)  (521)
                                     (111111)             (3311)
                                                          (11111111)
(End)
		

Crossrefs

The non-strict version is A083710.
The case with no 1's is A098965.
The Heinz numbers of these partitions are A339563.
The strict complement is counted by A341450.
The version for "divisible by" instead of "dividing" is A343347.
The case where there is also a part divisible by all the others is A343378.
The case where there is no part divisible by all the others is A343381.
A000005 counts divisors.
A000009 counts strict partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Mathematica
    Take[ CoefficientList[ Expand[ Sum[x^k*Product[1 + x^(k*i), {i, 2, 62}], {k, 62}]], x], {2, 60}] (* Robert G. Wilson v, Nov 01 2004 *)
    Table[Length[Select[IntegerPartitions[n], UnsameQ@@#&&Or@@Table[And@@IntegerQ/@(#/x), {x,#}]&]], {n,0,30}] (* Gus Wiseman, Apr 23 2021 *)
  • PARI
    A_x(N) = {my(x='x+O('x^N)); Vec(sum(k=1,N,x^k*prod(i=2,N-k, (1+x^(k*i)))))}
    A_x(50) \\ John Tyler Rascoe, Nov 19 2024

Formula

a(n) = Sum_{d|n} A025147(d-1).
G.f.: Sum_{k>=1} (x^k*Product_{i>=2} (1+x^(k*i))).
a(n) ~ exp(Pi*sqrt(n/3)) / (8*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Jul 06 2025

Extensions

More terms from Robert G. Wilson v, Nov 01 2004
Name shortened by Gus Wiseman, Apr 23 2021

A130689 Number of partitions of n such that every part divides the largest part; a(0) = 1.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 10, 11, 16, 19, 26, 28, 41, 43, 56, 65, 82, 88, 115, 122, 155, 174, 209, 225, 283, 305, 363, 402, 477, 514, 622, 666, 783, 858, 990, 1078, 1268, 1362, 1561, 1708, 1958, 2111, 2433, 2613, 2976, 3247, 3652, 3938, 4482, 4821, 5422
Offset: 0

Views

Author

Vladeta Jovovic, Jul 01 2007

Keywords

Comments

First differs from A130714 at a(11) = 28, A130714(11) = 27. - Gus Wiseman, Apr 23 2021

Examples

			For n = 6 we have 10 such partitions: [1, 1, 1, 1, 1, 1], [1, 1, 1, 1, 2], [1, 1, 2, 2], [2, 2, 2], [1, 1, 1, 3], [3, 3], [1, 1, 4], [2, 4], [1, 5], [6].
From _Gus Wiseman_, Apr 18 2021: (Start)
The a(1) = 1 through a(8) = 16 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (21)   (22)    (41)     (33)      (61)       (44)
             (111)  (31)    (221)    (42)      (331)      (62)
                    (211)   (311)    (51)      (421)      (71)
                    (1111)  (2111)   (222)     (511)      (422)
                            (11111)  (411)     (2221)     (611)
                                     (2211)    (4111)     (2222)
                                     (3111)    (22111)    (3311)
                                     (21111)   (31111)    (4211)
                                     (111111)  (211111)   (5111)
                                               (1111111)  (22211)
                                                          (41111)
                                                          (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
(End)
		

Crossrefs

The dual version is A083710.
The case without 1's is A339619.
The Heinz numbers of these partitions are the complement of A343337.
The complement is counted by A343341.
The strict case is A343347.
The complement in the strict case is counted by A343377.
A000009 counts strict partitions.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A072233 counts partitions by sum and greatest part.

Programs

  • Mathematica
    Table[If[n==0,1,Length[Select[IntegerPartitions[n],FreeQ[#,1]&&And@@IntegerQ/@(Max@@#/#)&]]],{n,0,30}] (* Gus Wiseman, Apr 18 2021 *)
  • PARI
    seq(n)={Vec(1 + sum(m=1, n, my(u=divisors(m)); x^m/prod(i=1, #u, 1 - x^u[i] + O(x^(n-m+1)))))} \\ Andrew Howroyd, Apr 17 2021

Formula

G.f.: 1 + Sum_{n>0} x^n/Product_{d divides n} (1-x^d).

A279789 Number of ways to choose a constant partition of each part of a constant partition of n.

Original entry on oeis.org

1, 1, 3, 3, 8, 3, 17, 3, 30, 12, 41, 3, 130, 3, 137, 45, 359, 3, 656, 3, 1306, 141, 2057, 3, 5446, 36, 8201, 544, 18610, 3, 34969, 3, 72385, 2061, 131081, 165, 290362, 3, 524297, 8205, 1109206, 3, 2130073, 3, 4371490, 33594, 8388617, 3, 17445321, 132, 33556496
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2016

Keywords

Comments

Also number of ways to choose a divisor d|n and then a sequence of n/d divisors of d.

Examples

			The a(6)=17 twice-constant partitions are:
((6)),
((3)(3)), ((33)),
((3)(111)), ((111)(3)),
((2)(2)(2)), ((222)),
((2)(2)(11)), ((2)(11)(2)), ((11)(2)(2)),
((2)(11)(11)), ((11)(2)(11)), ((11)(11)(2)),
((1)(1)(1)(1)(1)(1)), ((11)(11)(11)), ((111)(111)), ((111111)).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= proc(n) option remember; `if`(n=0, 1,
          add(tau(n/d)^d, d=divisors(n)))
        end:
    seq(a(n), n=0..70);  # Alois P. Heinz, Dec 20 2016
  • Mathematica
    nn=20;Table[DivisorSum[n,Power[DivisorSigma[0,#],n/#]&],{n,nn}]
  • PARI
    a(n)=if(n==0, 1, sumdiv(n, d, numdiv(n/d)^d)) \\ Andrew Howroyd, Aug 26 2018

Formula

a(n) = Sum_{d|n} tau(n/d)^d for n > 0. - Andrew Howroyd, Aug 26 2018
G.f.: 1 + Sum_{k>=1} tau(k)*x^k/(1 - tau(k)*x^k). - Ilya Gutkovskiy, May 23 2019
a(n) = 3 <=> n is prime <=> n in { A000040 }. - Alois P. Heinz, May 23 2019
Previous Showing 11-20 of 112 results. Next