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

A053251 Coefficients of the '3rd-order' mock theta function psi(q).

Original entry on oeis.org

0, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 5, 6, 7, 8, 9, 11, 12, 13, 16, 17, 19, 22, 24, 27, 31, 34, 37, 42, 46, 51, 57, 62, 68, 76, 83, 91, 101, 109, 120, 132, 143, 156, 171, 186, 202, 221, 239, 259, 283, 306, 331, 360, 388, 420, 455, 490, 529, 572, 616, 663, 716, 769, 827
Offset: 0

Views

Author

Dean Hickerson, Dec 19 1999

Keywords

Comments

Number of partitions of n into odd parts such that if a number occurs as a part then so do all smaller positive odd numbers.
Number of ways to express n as a partial sum of 1 + [1,3] + [1,5] + [1,7] + [1,9] + .... E.g., a(6)=2 because we have 6 = 1+1+1+1+1+1 = 1+3+1+1. - Jon Perry, Jan 01 2004
Also number of partitions of n such that the largest part occurs exactly once and all the other parts occur exactly twice. Example: a(9)=4 because we have [9], [7,1,1], [5,2,2] and [3,2,2,1,1]. - Emeric Deutsch, Mar 08 2006
Number of partitions (d1,d2,...,dm) of n such that 0 < d1/1 < d2/2 < ... < dm/m. - Seiichi Manyama, Mar 17 2018
For Emeric Deutsch's comment above, (1) this appears to be an alternately equal case of A122130, (2) the ordered version (compositions) is A239327, (3) allowing any length gives A351006, (4) the even-length version is A351007. - Gus Wiseman, Feb 25 2022

Examples

			q + q^2 + q^3 + 2*q^4 + 2*q^5 + 2*q^6 + 3*q^7 + 3*q^8 + 4*q^9 + ...
From _Seiichi Manyama_, Mar 17 2018: (Start)
n | Partition (d1,d2,...,dm) | (d1/1, d2/2, ... , dm/m)
--+--------------------------+-------------------------
1 | (1)                      | (1)
2 | (2)                      | (2)
3 | (3)                      | (3)
4 | (4)                      | (4)
  | (1, 3)                   | (1, 3/2)
5 | (5)                      | (5)
  | (1, 4)                   | (1, 2)
6 | (6)                      | (6)
  | (1, 5)                   | (1, 5/2)
7 | (7)                      | (7)
  | (1, 6)                   | (1, 3)
  | (2, 5)                   | (2, 5/2)
8 | (8)                      | (8)
  | (1, 7)                   | (1, 7/2)
  | (2, 6)                   | (2, 3)
9 | (9)                      | (9)
  | (1, 8)                   | (1, 4)
  | (2, 7)                   | (2, 7/2)
  | (1, 3, 5)                | (1, 3/2, 5/3) (End)
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 55, Eq. (26.13).
  • Srinivasa Ramanujan, Collected Papers, Chelsea, New York, 1962, pp. 354-355.
  • Srinivasa Ramanujan, The Lost Notebook and Other Unpublished Papers, Narosa Publishing House, New Delhi, 1988, p. 31.

Crossrefs

Other '3rd-order' mock theta functions are at A000025, A053250, A053252, A053253, A053254, A053255.
Cf. A003475.

Programs

  • Maple
    f:=n->q^(n^2)/mul((1-q^(2*i+1)),i=0..n-1); add(f(i),i=1..6);
    # second Maple program:
    b:= proc(n, i) option remember; (s-> `if`(n>s, 0, `if`(n=s, 1,
          b(n, i-1)+b(n-i, min(n-i, i-1)))))(i*(i+1)/2)
        end:
    a:= n-> `if`(n=0, 0, add(b(j, min(j, n-2*j-1)), j=0..iquo(n, 2))):
    seq(a(n), n=0..80);  # Alois P. Heinz, May 17 2018
  • Mathematica
    Series[Sum[q^n^2/Product[1-q^(2k-1), {k, 1, n}], {n, 1, 10}], {q, 0, 100}]
    (* Second program: *)
    b[n_, i_] := b[n, i] = Function[s, If[n > s, 0, If[n == s, 1, b[n, i - 1] + b[n - i, Min[n - i, i - 1]]]]][i*(i + 1)/2];
    a[n_] := If[n==0, 0, Sum[b[j, Min[j, n-2*j-1]], {j, 0, Quotient[n, 2]}]];
    Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Jun 17 2018, after Alois P. Heinz *)
  • PARI
    { n=20; v=vector(n); for (i=1,n,v[i]=vector(2^(i-1))); v[1][1]=1; for (i=2,n, k=length(v[i-1]); for (j=1,k, v[i][j]=v[i-1][j]+1; v[i][j+k]=v[i-1][j]+2*i-1)); c=vector(n); for (i=1,n, for (j=1,2^(i-1), if (v[i][j]<=n, c[v[i][j]]++))); c } \\ Jon Perry
    
  • PARI
    {a(n) = local(t); if(n<0, 0, t = 1 + O(x^n); polcoeff( sum(k=1, sqrtint(n), t *= x^(2*k-1) / (1 - x^(2*k-1)) + O(x^(n-(k-1)^2+1))), n))} /* Michael Somos, Sep 04 2007 */

Formula

G.f.: psi(q) = Sum_{n>=1} q^(n^2) / ( (1-q)*(1-q^3)*...*(1-q^(2*n-1)) ).
G.f.: Sum_{k>=1} q^k*Product_{j=1..k-1} (1+q^(2*j)) (see the Fine reference, p. 58, Eq. (26,53)). - Emeric Deutsch, Mar 08 2006
a(n) ~ exp(Pi*sqrt(n/6)) / (4*sqrt(n)). - Vaclav Kotesovec, Jun 09 2019

Extensions

More terms from Emeric Deutsch, Mar 08 2006

A344607 Number of integer partitions of n with reverse-alternating sum >= 0.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 8, 8, 15, 16, 27, 29, 48, 52, 81, 90, 135, 151, 220, 248, 352, 400, 553, 632, 859, 985, 1313, 1512, 1986, 2291, 2969, 3431, 4394, 5084, 6439, 7456, 9357, 10836, 13479, 15613, 19273, 22316, 27353, 31659, 38558, 44601, 53998, 62416, 75168
Offset: 0

Views

Author

Gus Wiseman, May 29 2021

Keywords

Comments

The reverse-alternating sum of a partition (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
Also the number of reversed integer partitions of n with alternating sum >= 0.
A formula for the reverse-alternating sum of a partition is: (-1)^(k-1) times the number of odd parts in the conjugate partition, where k is the number of parts. So a(n) is the number of integer partitions of n whose conjugate parts are all even or whose length is odd. By conjugation, this is also the number of integer partitions of n whose parts are all even or whose greatest part is odd.
All integer partitions have alternating sum >= 0, so the non-reversed version is A000041.
Is this sequence weakly increasing? In particular, is A344611(n) <= A160786(n)?

Examples

			The a(1) = 1 through a(8) = 15 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (221)    (33)      (322)      (44)
                    (211)   (311)    (222)     (331)      (332)
                    (1111)  (11111)  (321)     (421)      (422)
                                     (411)     (511)      (431)
                                     (2211)    (22111)    (521)
                                     (21111)   (31111)    (611)
                                     (111111)  (1111111)  (2222)
                                                          (3311)
                                                          (22211)
                                                          (32111)
                                                          (41111)
                                                          (221111)
                                                          (2111111)
                                                          (11111111)
		

Crossrefs

The non-reversed version is A000041.
The opposite version (rev-alt sum <= 0) is A027187, ranked by A028260.
The strict case for n > 0 is A067659 (even bisection: A344650).
The ordered version appears to be A116406 (even bisection: A114121).
The odd bisection is A160786.
The complement is counted by A344608.
The Heinz numbers of these partitions are A344609 (complement: A119899).
The even bisection is A344611.
A000070 counts partitions with alternating sum 1 (reversed: A000004).
A000097 counts partitions with alternating sum 2 (reversed: A120452).
A035363 counts partitions with alternating sum 0, ranked by A000290.
A103919 counts partitions by sum and alternating sum.
A316524 is the alternating sum of prime indices of n (reversed: A344616).
A325534/A325535 count separable/inseparable partitions.
A344610 counts partitions by sum and positive reverse-alternating sum.
A344612 counts partitions by sum and reverse-alternating sum.
A344618 gives reverse-alternating sums of standard compositions.

Programs

  • Mathematica
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Table[Length[Select[IntegerPartitions[n],sats[#]>=0&]],{n,0,30}]

Formula

a(n) + A344608(n) = A000041(n).
a(2n+1) = A160786(n).

A045931 Number of partitions of n with equal number of even and odd parts.

Original entry on oeis.org

1, 0, 0, 1, 0, 2, 1, 3, 2, 5, 5, 7, 9, 11, 16, 18, 25, 28, 41, 44, 62, 70, 94, 107, 140, 163, 207, 245, 302, 361, 440, 527, 632, 763, 904, 1090, 1285, 1544, 1812, 2173, 2539, 3031, 3538, 4202, 4896, 5793, 6736, 7934, 9221, 10811, 12549, 14661, 16994, 19780
Offset: 0

Views

Author

Keywords

Comments

The trivariate g.f. with x marking weight (i.e., sum of the parts), t marking number of odd parts and s marking number of even parts, is 1/product((1-tx^(2j-1))(1-sx^(2j)), j=1..infinity). - Emeric Deutsch, Mar 30 2006

Examples

			a(9) = 5 because we have [8,1], [7,2], [6,3], [5,4] and [2,2,2,1,1,1].
From _Gus Wiseman_, Jan 23 2022: (Start)
The a(0) = 1 through a(12) = 9 partitions (A = 10, empty columns indicated by dots):
  ()  .  .  21   .  32   2211   43   3221   54       3322   65       4332
                    41          52   4211   63       4321   74       4431
                                61          72       4411   83       5322
                                            81       5221   92       5421
                                            222111   6211   A1       6321
                                                            322211   6411
                                                            422111   7221
                                                                     8211
                                                                     22221111
(End)
		

Crossrefs

The version for subsets of {1..n} is A001405.
Dominated by A027187 (partitions of even length).
More odd/even parts: A108950/A108949.
More or same number of odd/even parts: A130780/A171966.
The strict case is A239241.
This is column k = 0 of the triangle A240009.
Counting only distinct parts gives A241638, ranked by A325700.
A half-conjugate version is A277579.
These partitions are ranked by A325698.
A000041 counts integer partitions, strict A000009.
A047993 counts balanced partitions, ranked by A106529.
A257991/A257992 count odd/even parts by Heinz number.

Programs

  • Maple
    g:=1/product((1-t*x^(2*j-1))*(1-s*x^(2*j)),j=1..30): gser:=simplify(series(g,x=0,56)): P[0]:=1: for n from 1 to 53 do P[n]:=subs(s=1/t,coeff(gser,x^n)) od: seq(coeff(t*P[n],t),n=0..53); # Emeric Deutsch, Mar 30 2006
  • Mathematica
    p[n_] := p[n] = Select[IntegerPartitions[n], Count[#, ?OddQ] == Count[#, ?EvenQ] &]; t = Table[p[n], {n, 0, 10}] (* partitions of n with # odd parts = # even parts *)
    TableForm[t] (* partitions, vertical format *)
    Table[Length[p[n]], {n, 0, 30}] (* A045931 *)
    (* Peter J. C. Moses, Mar 10 2014 *)
    nmax = 100; CoefficientList[Series[Sum[x^(3*k) / Product[(1 - x^(2*j))^2, {j, 1, k}], {k, 0, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jun 15 2025 *)

Formula

G.f.: Sum_{k>=0} x^(3*k)/Product_{i=1..k} (1-x^(2*i))^2. - Vladeta Jovovic, Aug 18 2007
a(n) = A000041(n)-A171967(n) = A130780(n)-A108950(n) = A171966(n)-A108949(n). - Reinhard Zumkeller, Jan 21 2010
a(n) = A000041(n) - A108950(n) - A108949(n) = A130780(n) + A171966(n) - A000041(n). - Gus Wiseman, Jan 23 2022
a(n) ~ Pi * exp(Pi*sqrt(2*n/3)) / (48*n^(3/2)). - Vaclav Kotesovec, Jun 15 2025

A077957 Powers of 2 alternating with zeros.

Original entry on oeis.org

1, 0, 2, 0, 4, 0, 8, 0, 16, 0, 32, 0, 64, 0, 128, 0, 256, 0, 512, 0, 1024, 0, 2048, 0, 4096, 0, 8192, 0, 16384, 0, 32768, 0, 65536, 0, 131072, 0, 262144, 0, 524288, 0, 1048576, 0, 2097152, 0, 4194304, 0, 8388608, 0, 16777216, 0, 33554432, 0, 67108864, 0, 134217728, 0, 268435456
Offset: 0

Views

Author

N. J. A. Sloane, Nov 17 2002

Keywords

Comments

Normally sequences like this are not included, since with the alternating 0's deleted it is already in the database.
Inverse binomial transform of A001333. - Paul Barry, Feb 25 2003
"Sloping binary representation" of powers of 2 (A000079), slope=-1 (see A037095 and A102370). - Philippe Deléham, Jan 04 2008
0,1,0,2,0,4,0,8,0,16,... is the inverse binomial transform of A000129 (Pell numbers). - Philippe Deléham, Oct 28 2008
Number of maximal self-avoiding walks from the NW to SW corners of a 3 X n grid.
Row sums of the triangle in A204293. - Reinhard Zumkeller, Jan 14 2012
Pisano period lengths: 1, 1, 4, 1, 8, 4, 6, 1, 12, 8, 20, 4, 24, 6, 8, 1, 16, 12, 36, 8, ... . - R. J. Mathar, Aug 10 2012
This sequence occurs in the length L(n) = sqrt(2)^n of Lévy's C-curve at the n-th iteration step. Therefore, L(n) is the Q(sqrt(2)) integer a(n) + a(n-1)*sqrt(2), with a(-1) = 0. For a variant of this C-curve see A251732 and A251733. - Wolfdieter Lang, Dec 08 2014
a(n) counts walks (closed) on the graph G(1-vertex,2-loop,2-loop). Equivalently the middle entry (2,2) of A^n where the adjacency matrix of digraph is A=(0,1,0;1,0,1;0,1,0). - David Neil McGrath, Dec 19 2014
a(n-2) is the number of compositions of n into even parts. For example, there are 4 compositions of 6 into even parts: (6), (222), (42), and (24). - David Neil McGrath, Dec 19 2014
Also the number of alternately constant compositions of n + 2, ranked by A351010. The alternately strict version gives A000213. The unordered version is A035363, ranked by A000290, strict A035457. - Gus Wiseman, Feb 19 2022
a(n) counts degree n fixed points of GF(2)[x]'s automorphisms. Proof: given a field k, k[x]'s automorphisms are determined by k's automorphisms and invertible affine maps x -> ax + b. GF(2) is rigid and has only one unit so its only nontrivial automorphism is x -> x + 1. For n = 0 we have 1 fixed point, the constant polynomial 1. (Taking the convention that 0 is not a degree 0 polynomial.) For n = 1 we have 0 fixed points as x -> x + 1 -> x are the only degree 1 polynomials. Note that if f(x) is a fixed point, then f(x) + 1 is also a fixed point. Given f(x) a degree n fixed point, we can assume WLOG x | f(x). Applying the automorphism, we then have x + 1 | f(x). Now note that f(x) / (x^2 + x) must be a fixed point, so any fixed point of degree n must either be of the form g(x) * (x^2 + x) or g(x) * (x^2 + x) + 1 for a unique degree n - 2 fixed point g(x). Therefore we have the recurrence relation a(n) = 2 * a(n - 2) as desired. - Keith J. Bauer, Mar 19 2024

Crossrefs

Column k=3 of A219946. - Alois P. Heinz, Dec 01 2012
Cf. A016116 (powers repeated).

Programs

  • GAP
    Flat(List([0..30],n->[2^n,0])); # Muniru A Asiru, Aug 05 2018
  • Haskell
    a077957 = sum . a204293_row  -- Reinhard Zumkeller, Jan 14 2012
    
  • Magma
    &cat [[2^n,0]: n in [0..20]]; // Vincenzo Librandi, Apr 03 2018
    
  • Maple
    seq(op([2^n,0]),n=0..100); # Robert Israel, Dec 23 2014
  • Mathematica
    a077957[n_] := Riffle[Table[2^i, {i, 0, n - 1}], Table[0, {n}]]; a077957[29] (* Michael De Vlieger, Dec 22 2014 *)
    CoefficientList[Series[1/(1 - 2*x^2), {x,0,50}], x] (* G. C. Greubel, Apr 12 2017 *)
    LinearRecurrence[{0, 2}, {1, 0}, 54] (* Robert G. Wilson v, Jul 23 2018 *)
    Riffle[2^Range[0,30],0,{2,-1,2}] (* Harvey P. Dale, Jan 06 2022 *)
  • PARI
    a(n)=if(n<0||n%2, 0, 2^(n/2))
    
  • Sage
    def A077957():
        x, y = -1, 1
        while True:
            yield -x
            x, y = x + y, x - y
    a = A077957(); [next(a) for i in range(40)]  # Peter Luschny, Jul 11 2013
    

Formula

G.f.: 1/(1-2*x^2).
E.g.f.: cosh(x*sqrt(2)).
a(n) = (1 - n mod 2) * 2^floor(n/2).
a(n) = sqrt(2)^n*(1+(-1)^n)/2. - Paul Barry, May 13 2003
a(n) = 2*a(n-2) with a(0)=1, a(1)=0. - Jim Singh, Jul 12 2018

A195017 If n = Product_{k >= 1} (p_k)^(c_k) where p_k is k-th prime and c_k >= 0 then a(n) = Sum_{k >= 1} c_k*((-1)^(k-1)).

Original entry on oeis.org

0, 1, -1, 2, 1, 0, -1, 3, -2, 2, 1, 1, -1, 0, 0, 4, 1, -1, -1, 3, -2, 2, 1, 2, 2, 0, -3, 1, -1, 1, 1, 5, 0, 2, 0, 0, -1, 0, -2, 4, 1, -1, -1, 3, -1, 2, 1, 3, -2, 3, 0, 1, -1, -2, 2, 2, -2, 0, 1, 2, -1, 2, -3, 6, 0, 1, 1, 3, 0, 1, -1, 1, 1, 0, 1, 1, 0, -1, -1, 5, -4, 2, 1, 0, 2, 0, -2, 4, -1, 0, -2, 3, 0, 2, 0, 4, 1, -1, -1, 4, -1, 1, 1, 2, -1
Offset: 1

Views

Author

Clark Kimberling, Feb 06 2012

Keywords

Comments

Let p(n,x) be the completely additive polynomial-valued function such that p(1,x) = 0 and p(prime(n),x) = x^(n-1), like is defined in A206284 (although here we are not limited to just irreducible polynomials). Then a(n) is the value of the polynomial encoded in such a manner by n, when it is evaluated at x=-1. - The original definition rewritten and clarified by Antti Karttunen, Oct 03 2018
Positions of 0 give the values of n for which the polynomial p(n,x) is divisible by x+1. For related sequences, see the Mathematica section.
Also the number of odd prime indices of n minus the number of even prime indices of n (both counted with multiplicity), where 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. - Gus Wiseman, Oct 24 2023

Examples

			The sequence can be read from a list of the polynomials:
  p(n,x)      with x = -1, gives a(n)
------------------------------------------
  p(1,x) = 0           0
  p(2,x) = 1x^0        1
  p(3,x) = x          -1
  p(4,x) = 2x^0        2
  p(5,x) = x^2         1
  p(6,x) = 1+x         0
  p(7,x) = x^3        -1
  p(8,x) = 3x^0        3
  p(9,x) = 2x         -2
  p(10,x) = x^2 + 1    2.
(The list runs through all the polynomials whose coefficients are nonnegative integers.)
		

Crossrefs

For other evaluation functions of such encoded polynomials, see A001222, A048675, A056239, A090880, A248663.
Zeros are A325698, distinct A325700.
For sum instead of count we have A366749 = A366531 - A366528.
A000009 counts partitions into odd parts, ranked by A066208.
A035363 counts partitions into even parts, ranked by A066207.
A112798 lists prime indices, reverse A296150, sum A056239.
A257991 counts odd prime indices, even A257992.
A300061 lists numbers with even sum of prime indices, odd A300063.

Programs

  • Mathematica
    b[n_] := Table[x^k, {k, 0, n}];
    f[n_] := f[n] = FactorInteger[n]; z = 200;
    t[n_, m_, k_] := If[PrimeQ[f[n][[m, 1]]] && f[n][[m, 1]]
    == Prime[k], f[n][[m, 2]], 0];
    u = Table[Apply[Plus,
        Table[Table[t[n, m, k], {k, 1, PrimePi[n]}], {m, 1,
          Length[f[n]]}]], {n, 1, z}];
    p[n_, x_] := u[[n]].b[-1 + Length[u[[n]]]]
    Table[p[n, x] /. x -> 0, {n, 1, z/2}]   (* A007814 *)
    Table[p[2 n, x] /. x -> 0, {n, 1, z/2}] (* A001511 *)
    Table[p[n, x] /. x -> 1, {n, 1, z}]     (* A001222 *)
    Table[p[n, x] /. x -> 2, {n, 1, z}]     (* A048675 *)
    Table[p[n, x] /. x -> 3, {n, 1, z}]     (* A090880 *)
    Table[p[n, x] /. x -> -1, {n, 1, z}]    (* A195017 *)
    z = 100; Sum[-(-1)^k IntegerExponent[Range[z], Prime[k]], {k, 1, PrimePi[z]}] (* Friedjof Tellkamp, Aug 05 2024 *)
  • PARI
    A195017(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * (-1)^(1+primepi(f[i,1])))); } \\ Antti Karttunen, Oct 03 2018

Formula

Totally additive with a(p^e) = e * (-1)^(1+PrimePi(p)), where PrimePi(n) = A000720(n). - Antti Karttunen, Oct 03 2018
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} = (-1)^(primepi(p)+1)/(p-1) = Sum_{k>=1} (-1)^(k+1)/A006093(k) = A078437 + Sum_{k>=1} (-1)^(k+1)/A036689(k) = 0.6339266524059... . - Amiram Eldar, Sep 29 2023
a(n) = A257991(n) - A257992(n). - Gus Wiseman, Oct 24 2023
a(n) = -Sum_{k=1..pi(n)} (-1)^k * valuation(n, prime(k)). - Friedjof Tellkamp, Aug 05 2024

Extensions

More terms, name changed and example-section edited by Antti Karttunen, Oct 03 2018

A036469 Partial sums of A000009 (partitions into distinct parts).

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 14, 19, 25, 33, 43, 55, 70, 88, 110, 137, 169, 207, 253, 307, 371, 447, 536, 640, 762, 904, 1069, 1261, 1483, 1739, 2035, 2375, 2765, 3213, 3725, 4310, 4978, 5738, 6602, 7584, 8697, 9957, 11383, 12993, 14809, 16857, 19161, 21751, 24661
Offset: 0

Views

Author

Keywords

Comments

Also number of 1's in all partitions of n+1 into odd parts. Example: a(4)=7 because the partitions of 5 into odd parts are [5], [3,1,1], [1,1,1,1,1], having a total number of 7 1's. - Emeric Deutsch, Mar 29 2006
Convolved with A035363 = A000070. - Gary W. Adamson, Jun 09 2009
Equals row sums of triangle A166240. - Gary W. Adamson, Oct 09 2009
a(n) = if n <= 1 then A201377(1,n) else A201377(n,1). - Reinhard Zumkeller, Dec 02 2011
a(n) equals the sum of the parts of the form 2^k (k >= 0) in all partitions of n + 1 into distinct parts. Example: a(6) = 14. The partitions of 7 into distinct parts are [7], [6,1], [5,2], [4,3] and [4,2,1] having sum over parts of the form 2^k equal to 1 + 2 + 4 + 4 + 2 + 1 = 14. - Peter Bala, Dec 01 2013
Number of partitions of the (n+1)-multiset {0,...,0,1} with n 0's into distinct multisets; a(3) = 5: 0|00|1, 00|01, 000|1, 0|001, 0001. Also number of factorizations of 3*2^n into distinct factors; a(3) = 5: 2*3*4, 4*6, 3*8, 2*12, 24. - Alois P. Heinz, Jul 30 2021

Crossrefs

Cf. A035363, A000070. - Gary W. Adamson, Jun 09 2009
Cf. A166240. - Gary W. Adamson, Oct 09 2009
Column k=1 of A346520.

Programs

  • Maple
    g:=1/(1-x)/product(1-x^(2*j-1),j=1..30): gser:=series(g,x=0,50): seq(coeff(gser,x,n),n=0..46); # Emeric Deutsch, Mar 29 2006
    # second Maple program:
    b:= proc(n, i) b(n, i):= `if`(n=0, 1, `if`(i<1, 0,
           b(n, i-1)+`if`(i>n, 0, b(n-i, min(n-i, i-1)))))
        end:
    a:= proc(n) option remember; b(n, n) +`if`(n>0, a(n-1), 0) end:
    seq(a(n), n=0..60);  # Alois P. Heinz, Nov 21 2012
  • Mathematica
    CoefficientList[ Series[Product[(1 + t^i), {i, 1, Infinity}]/(1 - t), {t, 0, 46}], t] (* Geoffrey Critzer, May 16 2010 *)
    b[n_, i_] := If[n == 0, 1, If[i<1, 0, b[n, i-1]+If[i>n, 0, b[n-i, Min[n-i, i-1]]]]]; a[n_] := a[n] = b[n, n]+If[n>0, a[n-1], 0]; Table[a[n], {n, 0, 60}] // Flatten (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
    Accumulate[Table[PartitionsQ[n], {n, 0, 50}]] (* Vaclav Kotesovec, Oct 26 2016 *)

Formula

G.f.: 1/[(1-x)*product(1-x^(2j-1), j=1..infinity)]. - Emeric Deutsch, Mar 29 2006
a(n) ~ 3^(1/4) * exp(Pi*sqrt(n/3)) / (2*Pi*n^(1/4)) * (1 + (18+13*Pi^2) / (48*Pi*sqrt(3*n)) + (2916 - 1404*Pi^2 + 121*Pi^4)/(13824*Pi^2*n)). - Vaclav Kotesovec, Feb 26 2015, updated Oct 26 2016
For n > 0, a(n) = A026906(n) + 1. - Vaclav Kotesovec, Oct 26 2016
Faster converging g.f.: A(x) = (1/(1 - x))*Sum_{n >= 0} x^(n*(2*n-1))/Product_{k = 1..2*n} (1 - x^k). - Peter Bala, Feb 02 2021

A236913 Number of partitions of 2n of type EE (see Comments).

Original entry on oeis.org

1, 1, 3, 6, 12, 22, 40, 69, 118, 195, 317, 505, 793, 1224, 1867, 2811, 4186, 6168, 9005, 13026, 18692, 26613, 37619, 52815, 73680, 102162, 140853, 193144, 263490, 357699, 483338, 650196, 870953, 1161916, 1544048, 2044188, 2696627, 3545015, 4644850, 6066425
Offset: 0

Views

Author

Clark Kimberling, Feb 01 2014

Keywords

Comments

The partitions of n are partitioned into four types:
EO, even # of odd parts and odd # of even parts, A236559;
OE, odd # of odd parts and even # of even parts, A160786;
EE, even # of odd parts and even # of even parts, A236913;
OO, odd # of odd parts and odd # of even parts, A236914.
A236559 and A160786 are the bisections of A027193;
A236913 and A236914 are the bisections of A027187.

Examples

			The partitions of 4 of type EE are [3,1], [2,2], [1,1,1,1], so that a(2) = 3.
type/k . 1 .. 2 .. 3 .. 4 .. 5 .. 6 .. 7 .. 8 ... 9 ... 10 .. 11
EO ..... 0 .. 1 .. 0 .. 2 .. 0 .. 5 .. 0 .. 10 .. 0 ... 20 .. 0
OE ..... 1 .. 0 .. 2 .. 0 .. 4 .. 0 .. 8 .. 0 ... 16 .. 0 ... 29
EE ..... 0 .. 1 .. 0 .. 3 .. 0 .. 6 .. 0 .. 12 .. 0 ... 22 .. 0
OO ..... 0 .. 0 .. 1 .. 0 .. 3 .. 0 .. 7 .. 0 ... 14 .. 0 ... 27
From _Gus Wiseman_, Feb 09 2021: (Start)
This sequence counts even-length partitions of even numbers, which have Heinz numbers given by A340784. For example, the a(0) = 1 through a(4) = 12 partitions are:
  ()  (11)  (22)    (33)      (44)
            (31)    (42)      (53)
            (1111)  (51)      (62)
                    (2211)    (71)
                    (3111)    (2222)
                    (111111)  (3221)
                              (3311)
                              (4211)
                              (5111)
                              (221111)
                              (311111)
                              (11111111)
(End)
		

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The ordered version is A000302.
The case of odd-length partitions of odd numbers is A160786 (A340931).
The Heinz numbers of these partitions are (A340784).
A027187 counts partitions of even length/maximum (A028260/A244990).
A034008 counts compositions of even length.
A035363 counts partitions into even parts (A066207).
A047993 counts balanced partitions (A106529).
A058695 counts partitions of odd numbers (A300063).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A072233 counts partitions by sum and length.
A339846 counts factorizations of even length.
A340601 counts partitions of even rank (A340602).
A340785 counts factorizations into even factors.
A340786 counts even-length factorizations into even factors.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0$3],
          `if`(i<1, [0$4], b(n, i-1)+`if`(i>n, [0$4], (p->
          `if`(irem(i, 2)=0, [p[3], p[4], p[1], p[2]],
              [p[2], p[1], p[4], p[3]]))(b(n-i, i)))))
        end:
    a:= n-> b(2*n$2)[1]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 16 2014
  • Mathematica
    z = 25; m1 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,  OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]]; m2 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,      OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]]; m3 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]] ; m4 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]];
    m1 (* A236559, type EO*)
    m2 (* A160786, type OE*)
    m3 (* A236913, type EE*)
    m4 (* A236914, type OO*)
    (* Peter J. C. Moses, Feb 03 2014 *)
    b[n_, i_] := b[n, i] = If[n == 0, {1, 0, 0, 0}, If[i < 1, {0, 0, 0, 0}, b[n, i - 1] + If[i > n, {0, 0, 0, 0}, Function[p, If[Mod[i, 2] == 0, p[[{3, 4, 1, 2}]], p[[{2, 1, 4, 3}]]]][b[n - i, i]]]]]; a[n_] := b[2*n, 2*n][[1]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Oct 27 2015, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[2n],EvenQ[Length[#]]&]],{n,0,15}] (* Gus Wiseman, Feb 09 2021 *)

Extensions

More terms from Alois P. Heinz, Feb 16 2014

A119620 Number of partitions of floor(3n/2) into n parts each from {1,2,...,n}.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 5, 5, 7, 7, 11, 11, 15, 15, 22, 22, 30, 30, 42, 42, 56, 56, 77, 77, 101, 101, 135, 135, 176, 176, 231, 231, 297, 297, 385, 385, 490, 490, 627, 627, 792, 792, 1002, 1002, 1255, 1255, 1575, 1575, 1958, 1958, 2436, 2436, 3010, 3010, 3718, 3718
Offset: 0

Views

Author

John W. Layman, Jun 07 2006

Keywords

Comments

The bisection {1,1,2,3,5,7,11,15,22,...} agrees with the initial terms of A008641, Number of partitions of n into at most 12 parts and also A008635, Molien series for A_12.
a(2n+1)=a(2n) for all n>0. If the partition {...,1} is a member of a(2n) then the partition {...,1,1} is a member of a(2n+1). - Robert G. Wilson v, Jun 09 2006
Number of partitions of n where all parts (except for possibly the first part) are even; see example. - Joerg Arndt, Apr 22 2013
For n >= 2, a(n) = number of partitions p of n such that floor(n/2) is a part of p. For n >= 1, a(n) = number of partitions p of n such that ceiling(n/2) is a part of p. - Clark Kimberling, Feb 28 2014
From Gus Wiseman, Oct 28 2021: (Start)
If we insert zeros every three terms, this counts partitions of n such that n = floor(3*k/2), where k is the number of parts. This counts by sum rather than length. These partitions are ranked by A347452.
Also the number of integer partitions of n with alternating product 1, where the alternating product of a sequence (y_1,...,y_k) is Product_i y_i^((-1)^(i-1)). These are the conjugates of the partitions (ranked by A336119) described in Arndt's comment above. For example, the a(2) = 1 through a(10) = 7 partitions are:
11 111 22 221 33 331 44 441 55
1111 11111 2211 22111 2222 22221 3322
111111 1111111 3311 33111 4411
221111 2211111 222211
11111111 111111111 331111
22111111
1111111111
These partitions are ranked by A028982. The odd-length case is A035363 (shifted), which is also the version for sum instead of product. The multiplicative version (factorizations) is A347438.
(End)

Examples

			For n=8, floor(3*n/2) is 12 and there are five partitions of 12 into 8 parts each in the range 1-8 inclusive, namely: {5,1,1,1,1,1,1,1}, {4,2,1,1,1,1,1,1}, {3,3,1,1,1,1,1,1}, {3,2,2,1,1,1,1,1} and {2,2,2,2,1,1,1,1}. Thus a(8)=5.
From _Joerg Arndt_, Apr 22 2013: (Start)
a(8) = a(9) = 5, counting the following partitions where all parts (except for possibly the first part) are even:
01:  [ 2 2 2 2 ]
02:  [ 4 2 2 ]
03:  [ 4 4 ]
04:  [ 6 2 ]
05:  [ 8 ]
and
01:  [ 3 2 2 2 ]
02:  [ 5 2 2 ]
03:  [ 5 4 ]
04:  [ 7 2 ]
05:  [ 9 ]
(End)
G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 3*x^6 + 3*x^7 + 5*x^8 + 5*x^9 + 7*x^10 + ...
		

Crossrefs

Both bisections are A000041.
An adjoint version is A108711.
A027187 counts partitions of even length.
A027193 counts partitions of odd length.
A325534 counts separable partitions.
A325535 counts inseparable partitions.

Programs

  • Maple
    # Using the function EULER from Transforms (see link at the bottom of the page).
    [1, op(EULER([1,0,seq(irem(n,2),n=2..55)]))]; # Peter Luschny, Aug 19 2020
  • Mathematica
    (* first do *) Needs["DiscreteMath`Combinatorica`"] (* then *) f[n_] := f[n] = Length@ Select[ Partitions[ Floor[3n/2], n], Length@# == n &]; Table[ If[n > 1, f[2Floor[n/2]], f[n]], {n, 57}] (* Robert G. Wilson v, Jun 09 2006 *)
    Table[ PartitionsP[ Floor[n/2]], {n, 57}] (* Robert G. Wilson v, Jun 09 2006 *)
    Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, Ceiling[n/2]]], {n, 50}] (* Clark Kimberling, Feb 28 2014 *)
    a[ n_] := SeriesCoefficient[ (1 + x) / QPochhammer[x^2], {x, 0, n}]; (* Michael Somos, Mar 01 2014 *)
  • PARI
    a(n)=numbpart(n\2); \\ Joerg Arndt, Apr 22 2013

Formula

a(n) = A000041(floor(n/2)). - Vladeta Jovovic, Jun 10 2006
G.f.: (Sum_{n>=0} x^(4*n) / Product_{k=1..n} (1-x^(2*k))) / (1 - x). - Michael Somos, Mar 01 2014 [corrected by Jason Yuen, Jan 24 2025]

Extensions

More terms from Robert G. Wilson v, Jun 09 2006
Added a(0)=1. - Michael Somos, Mar 01 2014

A087897 Number of partitions of n into odd parts greater than 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 5, 6, 8, 8, 10, 12, 13, 15, 18, 20, 23, 27, 30, 34, 40, 44, 50, 58, 64, 73, 83, 92, 104, 118, 131, 147, 166, 184, 206, 232, 256, 286, 320, 354, 394, 439, 485, 538, 598, 660, 730, 809, 891, 984, 1088, 1196, 1318, 1454, 1596, 1756
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2003

Keywords

Comments

Also number of partitions of n into distinct parts which are not powers of 2.
Also number of partitions of n into distinct parts such that the two largest parts differ by 1.
Also number of partitions of n such that the largest part occurs an odd number of times that is at least 3 and every other part occurs an even number of times. Example: a(10) = 2 because we have [2,2,2,1,1,1,1] and [2,2,2,2,2]. - Emeric Deutsch, Mar 30 2006
Also difference between number of partitions of 1+n into distinct parts and number of partitions of n into distinct parts. - Philippe LALLOUET, May 08 2007
In the Berndt reference replace {a -> -x, q -> x} in equation (3.1) to get f(x). G.f. is 1 - x * (1 - f(x)).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Also number of symmetric unimodal compositions of n+3 where the maximal part appears three times. - Joerg Arndt, Jun 11 2013
Let c(n) = number of palindromic partitions of n whose greatest part has multiplicity 3; then c(n) = a(n-3) for n>=3. - Clark Kimberling, Mar 05 2014
From Gus Wiseman, Aug 22 2021: (Start)
Also the number of integer partitions of n - 1 whose parts cover an interval of positive integers starting with 2. These partitions are ranked by A339886. For example, the a(6) = 1 through a(16) = 5 partitions are:
32 222 322 332 432 3322 3332 4332 4432 5432 43332
2222 3222 22222 4322 33222 33322 33332 44322
32222 222222 43222 43322 333222
322222 332222 432222
2222222 3222222
(End)

Examples

			1 + x^3 + x^5 + x^6 + x^7 + x^8 + 2*x^9 + 2*x^10 + 2*x^11 + 3*x^12 + 3*x^13 + ...
q + q^73 + q^121 + q^145 + q^169 + q^193 + 2*q^217 + 2*q^241 + 2*q^265 + ...
a(10)=2 because we have [7,3] and [5,5].
From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(22)=13 symmetric unimodal compositions of 22+3=25 where the maximal part appears three times:
01:  [ 1 1 1 1 1 1 1 1 3 3 3 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 2 3 3 3 2 1 1 1 1 1 1 ]
03:  [ 1 1 1 1 1 5 5 5 1 1 1 1 1 ]
04:  [ 1 1 1 1 2 2 3 3 3 2 2 1 1 1 1 ]
05:  [ 1 1 1 2 5 5 5 2 1 1 1 ]
06:  [ 1 1 2 2 2 3 3 3 2 2 2 1 1 ]
07:  [ 1 1 3 5 5 5 3 1 1 ]
08:  [ 1 1 7 7 7 1 1 ]
09:  [ 1 2 2 5 5 5 2 2 1 ]
10:  [ 1 4 5 5 5 4 1 ]
11:  [ 2 2 2 2 3 3 3 2 2 2 2 ]
12:  [ 2 3 5 5 5 3 2 ]
13:  [ 2 7 7 7 2 ]
(End)
From _Gus Wiseman_, Feb 16 2021: (Start)
The a(7) = 1 through a(19) = 8 partitions are the following (A..J = 10..19). The Heinz numbers of these partitions are given by A341449.
  7  53  9    55  B    75    D    77    F      97    H      99      J
         333  73  533  93    553  95    555    B5    755    B7      775
                       3333  733  B3    753    D3    773    D5      955
                                  5333  933    5533  953    F3      973
                                        33333  7333  B33    5553    B53
                                                     53333  7533    D33
                                                            9333    55333
                                                            333333  73333
(End)
		

References

  • J. W. L. Glaisher, Identities, Messenger of Mathematics, 5 (1876), pp. 111-112. see Eq. I

Crossrefs

The ordered version is A000931.
Partitions with no ones are counted by A002865, ranked by A005408.
The even version is A035363, ranked by A066207.
The version for factorizations is A340101.
Partitions whose only even part is the smallest are counted by A341447.
The Heinz numbers of these partitions are given by A341449.
A000009 counts partitions into odd parts, ranked by A066208.
A025147 counts strict partitions with no 1's.
A025148 counts strict partitions with no 1's or 2's.
A026804 counts partitions whose smallest part is odd, ranked by A340932.
A027187 counts partitions with even length/maximum, ranks A028260/A244990.
A027193 counts partitions with odd length/maximum, ranks A026424/A244991.
A058695 counts partitions of odd numbers, ranked by A300063.
A058696 counts partitions of even numbers, ranked by A300061.
A340385 counts partitions with odd length and maximum, ranked by A340386.

Programs

  • Haskell
    a087897 = p [3,5..] where
       p [] _ = 0
       p _  0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 12 2011
    
  • Maple
    To get 128 terms: t4 := mul((1+x^(2^n)),n=0..7); t5 := mul((1+x^k),k=1..128): t6 := series(t5/t4,x,100); t7 := seriestolist(t6);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1,
          `if`(i<3, 0, b(n, i-2)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(n, n-1+irem(n, 2)):
    seq(a(n), n=0..80);  # Alois P. Heinz, Jun 11 2013
  • Mathematica
    max = 65; f[x_] := Product[ 1/(1 - x^(2k+1)), {k, 1, max}]; CoefficientList[ Series[f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 16 2011, after Emeric Deutsch *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<3, 0, b[n, i-2]+If[i>n, 0, b[n-i, i]]] ]; a[n_] := b[n, n-1+Mod[n, 2]]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Apr 01 2015, after Alois P. Heinz *)
    Flatten[{1, Table[PartitionsQ[n+1] - PartitionsQ[n], {n, 0, 80}]}] (* Vaclav Kotesovec, Dec 01 2015 *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&OddQ[Times@@#]&]],{n,0,30}] (* Gus Wiseman, Feb 16 2021 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - x) * eta(x^2 + A) / eta(x + A), n))} /* Michael Somos, Nov 13 2011 */
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A087897_T(n,k):
        if n==0: return 1
        if k<3 or n<0: return 0
        return A087897_T(n,k-2)+A087897_T(n-k,k)
    def A087897(n): return A087897_T(n,n-(n&1^1)) # Chai Wah Wu, Sep 23 2023, after Alois P. Heinz

Formula

Expansion of q^(-1/24) * (1 - q) * eta(q^2) / eta(q) in powers of q.
Expansion of (1 - x) / chi(-x) in powers of x where chi() is a Ramanujan theta function.
G.f.: 1 + x^3 + x^5*(1 + x) + x^7*(1 + x)*(1 + x^2) + x^9*(1 + x)*(1 + x^2)*(1 + x^3) + ... [Glaisher 1876]. - Michael Somos, Jun 20 2012
G.f.: Product_{k >= 1} 1/(1-x^(2*k+1)).
G.f.: Product_{k >= 1, k not a power of 2} (1+x^k).
G.f.: Sum_{k >= 1} x^(3*k)/Product_{j = 1..k} (1 - x^(2*j)). - Emeric Deutsch, Mar 30 2006
a(n) ~ exp(Pi*sqrt(n/3)) * Pi / (8 * 3^(3/4) * n^(5/4)) * (1 - (15*sqrt(3)/(8*Pi) + 11*Pi/(48*sqrt(3)))/sqrt(n) + (169*Pi^2/13824 + 385/384 + 315/(128*Pi^2))/n). - Vaclav Kotesovec, Aug 30 2015, extended Nov 04 2016
G.f.: 1/(1 - x^3) * Sum_{n >= 0} x^(5*n)/Product_{k = 1..n} (1 - x^(2*k)) = 1/((1 - x^3)*(1 - x^5)) * Sum_{n >= 0} x^(7*n)/Product_{k = 1..n} (1 - x^(2*k)) = ..., extending Deutsch's result dated Mar 30 2006. - Peter Bala, Jan 15 2021
G.f.: Sum_{n >= 0} x^(n*(2*n+1))/Product_{k = 2..2*n+1} (1 - x^k). (Set z = x^3 and q = x^2 in Mc Laughlin et al., Section 1.3, Entry 7.) - Peter Bala, Feb 02 2021
a(2*n+1) = Sum{j>=1} A008284(n+1-j,2*j - 1) and a(2*n) = Sum{j>=1} A008284(n-j, 2*j). - Gregory L. Simay, Sep 22 2023

A320924 Heinz numbers of multigraphical partitions.

Original entry on oeis.org

1, 4, 9, 12, 16, 25, 27, 30, 36, 40, 48, 49, 63, 64, 70, 75, 81, 84, 90, 100, 108, 112, 120, 121, 144, 147, 154, 160, 165, 169, 175, 189, 192, 196, 198, 210, 220, 225, 243, 250, 252, 256, 264, 270, 273, 280, 286, 289, 300, 324, 325, 336, 343, 351, 352, 360
Offset: 1

Views

Author

Gus Wiseman, Oct 24 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
An integer partition is multigraphical if it comprises the multiset of vertex-degrees of some multigraph.
Also Heinz numbers of integer partitions of even numbers whose greatest part is less than or equal to half the sum of parts, i.e., numbers n whose sum of prime indices A056239(n) is even and at least twice the greatest prime index A061395(n). - Gus Wiseman, May 23 2021

Examples

			The sequence of all multigraphical partitions begins: (), (11), (22), (211), (1111), (33), (222), (321), (2211), (3111), (21111), (44), (422), (111111), (431), (332), (2222), (4211), (3221), (3311), (22211), (41111), (32111), (55), (221111).
From _Gus Wiseman_, May 23 2021: (Start)
The sequence of terms together with their prime indices and a multigraph realizing each begins:
    1:      () | {}
    4:    (11) | {{1,2}}
    9:    (22) | {{1,2},{1,2}}
   12:   (112) | {{1,3},{2,3}}
   16:  (1111) | {{1,2},{3,4}}
   25:    (33) | {{1,2},{1,2},{1,2}}
   27:   (222) | {{1,2},{1,3},{2,3}}
   30:   (123) | {{1,3},{2,3},{2,3}}
   36:  (1122) | {{1,2},{3,4},{3,4}}
   40:  (1113) | {{1,4},{2,4},{3,4}}
   48: (11112) | {{1,2},{3,5},{4,5}}
   49:    (44) | {{1,2},{1,2},{1,2},{1,2}}
   63:   (224) | {{1,3},{1,3},{2,3},{2,3}}
(End)
		

Crossrefs

These partitions are counted by A209816.
The case with odd weights is A322109.
The conjugate case of equality is A340387.
The conjugate version with odd weights allowed is A344291.
The conjugate opposite version is A344292.
The opposite version with odd weights allowed is A344296.
The conjugate version is A344413.
The conjugate opposite version with odd weights allowed is A344414.
The case of equality is A344415.
The opposite version is A344416.
A000070 counts non-multigraphical partitions.
A025065 counts palindromic partitions.
A035363 counts partitions into even parts.
A056239 adds up prime indices, row sums of A112798.
A110618 counts partitions that are the vertex-degrees of some set multipartition with no singletons.
A334201 adds up all prime indices except the greatest.

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[[#]]&]}]]];
    Select[Range[1000],prptns[Flatten[MapIndexed[Table[#2,{#1}]&,If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]]]]!={}&]

Formula

Members m of A300061 such that A061395(m) <= A056239(m)/2. - Gus Wiseman, May 23 2021
Previous Showing 11-20 of 142 results. Next