A264396 -id:A264396 - cp's OEIS frontend

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

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

N. J. A. Sloane, Dec 04 2003

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)

			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)

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

Chai Wah Wu, Table of n, a(n) for n = 0..10000 (n = 0..1000 from Alois P. Heinz)
C. Ballantine and M. Merca, Padovan numbers as sums over partitions into odd parts, Journal of Inequalities and Applications, (2016) 2016:1; doi.
B. C. Berndt, B. Kim, and A. J. Yee, Ramanujan's lost notebook: Combinatorial proofs of identities associated with Heine's transformation or partial theta functions, J. Comb. Thy. Ser. A, 117 (2010), 957-973.
Howard D. Grossman, Problem 228, Mathematics Magazine, 28 (1955), p. 160.
R. K. Guy, Two theorems on partitions, Math. Gaz., 42 (1958), 84-86. Math. Rev. 20 #3110.
Cristiano Husu, The butterfly sequence: the second difference sequence of the numbers of integer partitions with distinct parts, its pentagonal number structure, its combinatorial identities and the cyclotomic polynomials 1-x and 1+x+x^2, arXiv:1804.09883 [math.NT], 2018.
James Mc Laughlin, Andrew V. Sills, and Peter Zimmer, Rogers-Ramanujan-Slater Type Identities, Electronic J. Combinatorics, DS15, 1-59, May 31, 2008.
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

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.
Cf. A000041, A003114, A039900, A160786, A257991/A257992, A264396, A300272, A339662, A339737, A339886.

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

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

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

{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 */

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

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

A251729 Number of gap-free but not complete compositions of n.

Original entry on oeis.org

0, 1, 1, 2, 3, 3, 6, 6, 14, 12, 27, 33, 58, 86, 134, 210, 323, 539, 810, 1371, 2044, 3510, 5263, 8927, 13702, 22870, 35821, 58750, 93343, 152236, 243244, 395078, 634342, 1027876, 1656543, 2676693, 4325727, 6982440, 11299457, 18232217, 29518334, 47641410

Offset: 1

Alois P. Heinz, Dec 07 2014

nonn

A composition is gap-free but not complete if all integers in the interval defined by the smallest and the largest part are parts but 1 is not a part.

			a(6) = 3: [6], [3,3], [2,2,2].
a(7) = 6: [7], [3,4], [4,3], [2,2,3], [2,3,2], [3,2,2].

Alois P. Heinz, Table of n, a(n) for n = 1..1000
Alois P. Heinz, Plot of (a(n)+a(n+1)-a(n+2))/a(n+2) for n = 150..1000
P. Hitczenko and A. Knopfmacher, Gap-free compositions and gap-free samples of geometric random variables, Discrete Math., 294 (2005), 225-239.

Cf. A001622, A034296, A107428, A107429, A188575, A238353, A264396.

b:= proc(n, i, t) option remember; `if`(n=0, `if`(i=0, 0, t!),
     `if`(i<1 or n add(b(n, i, 0), i=1..n):
seq(a(n), n=1..50);

b[n_, i_, t_] := b[n, i, t] = If[n == 0, If[i == 0, 0, t!], If[i < 1 || n < i, 0, Sum[b[n - i*j, i - 1, t + j]/j!, {j, 1, n/i}]]];
a[n_] := Sum[b[n, i, 0], {i, 1, n}];
Array[a, 50] (* Jean-François Alcover, Jan 25 2021, after Alois P. Heinz *)

a(n) = A107428(n) - A107429(n).

lim_{n -> oo} a(n)/a(n-1) = (1+sqrt(5))/2 = phi = A001622.

A356845 Odd numbers with gapless prime indices.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 23, 25, 27, 29, 31, 35, 37, 41, 43, 45, 47, 49, 53, 59, 61, 67, 71, 73, 75, 77, 79, 81, 83, 89, 97, 101, 103, 105, 107, 109, 113, 121, 125, 127, 131, 135, 137, 139, 143, 149, 151, 157, 163, 167, 169, 173, 175, 179, 181, 191

Offset: 1

Gus Wiseman, Sep 03 2022

nonn

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.

A sequence is gapless if it covers an interval of positive integers.

			The terms together with their prime indices begin:
    1: {}
    3: {2}
    5: {3}
    7: {4}
    9: {2,2}
   11: {5}
   13: {6}
   15: {2,3}
   17: {7}
   19: {8}
   23: {9}
   25: {3,3}
   27: {2,2,2}
   29: {10}
   31: {11}
   35: {3,4}
   37: {12}
   41: {13}
   43: {14}

Consists of the odd terms of A073491.
These partitions are counted by A264396.
The strict case is A294674, counted by A136107.
The version for compositions is A356843, counted by A251729.
A001221 counts distinct prime factors, sum A001414.
A056239 adds up prime indices, row sums of A112798, lengths A001222.
A356069 counts gapless divisors, initial A356224 (complement A356225).
A356230 ranks gapless factorization lengths, firsts A356603.
A356233 counts factorizations into gapless numbers.
Cf. A003963, A034296, A055932, A073493, A107428, A287170, A289508, A325160, A356231, A356234, A356841.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
nogapQ[m_]:=Or[m=={},Union[m]==Range[Min[m],Max[m]]];
Select[Range[1,100,2],nogapQ[primeMS[#]]&]

A339886 Numbers whose prime indices cover an interval of positive integers starting with 2.

Original entry on oeis.org

1, 3, 9, 15, 27, 45, 75, 81, 105, 135, 225, 243, 315, 375, 405, 525, 675, 729, 735, 945, 1125, 1155, 1215, 1575, 1875, 2025, 2187, 2205, 2625, 2835, 3375, 3465, 3645, 3675, 4725, 5145, 5625, 5775, 6075, 6561, 6615, 7875, 8085, 8505, 9375, 10125, 10395, 10935

Offset: 1

Gus Wiseman, Apr 20 2021

nonn

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 sequence of terms together with their prime indices begins:
    3: {2}
    9: {2,2}
   15: {2,3}
   27: {2,2,2}
   45: {2,2,3}
   75: {2,3,3}
   81: {2,2,2,2}
  105: {2,3,4}
  135: {2,2,2,3}
  225: {2,2,3,3}
  243: {2,2,2,2,2}
  315: {2,2,3,4}
  375: {2,3,3,3}
  405: {2,2,2,2,3}
  525: {2,3,3,4}
  675: {2,2,2,3,3}
  729: {2,2,2,2,2,2}
  735: {2,3,4,4}
  945: {2,2,2,3,4}

The version starting at 1 is A055932.
The partitions with these Heinz numbers are counted by A264396.
Positions of 1's in A339662.
A000009 counts partitions covering an initial interval.
A000070 counts partitions with a selected part.
A016945 lists numbers with smallest prime index 2.
A034296 counts gap-free (or flat) partitions.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A107428 counts gap-free compositions (initial: A107429).
A286469 and A286470 give greatest difference for Heinz numbers.
A325240 lists numbers with smallest prime multiplicity 2.
A342050/A342051 have prime indices with odd/even least gap.
Cf. A001223, A001522, A006128, A007052, A124010, A257989, A257993, A264401, A317090, A317589, A339737.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
normQ[m_]:=Or[m=={},Union[m]==Range[Max[m]]];
Select[Range[100],normQ[primeMS[#]-1]&]

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A251729 Number of gap-free but not complete compositions of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A356845 Odd numbers with gapless prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A339886 Numbers whose prime indices cover an interval of positive integers starting with 2.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica