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

A027193 Number of partitions of n into an odd number of parts.

Original entry on oeis.org

0, 1, 1, 2, 2, 4, 5, 8, 10, 16, 20, 29, 37, 52, 66, 90, 113, 151, 190, 248, 310, 400, 497, 632, 782, 985, 1212, 1512, 1851, 2291, 2793, 3431, 4163, 5084, 6142, 7456, 8972, 10836, 12989, 15613, 18646, 22316, 26561, 31659, 37556, 44601, 52743, 62416, 73593, 86809, 102064, 120025, 140736
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n in which greatest part is odd.
Number of partitions of n+1 into an even number of parts, the least being 1. Example: a(5)=4 because we have [5,1], [3,1,1,1], [2,1,1] and [1,1,1,1,1,1].
Also number of partitions of n+1 such that the largest part is even and occurs only once. Example: a(5)=4 because we have [6], [4,2], [4,1,1] and [2,1,1,1,1]. - Emeric Deutsch, Apr 05 2006
Also the number of partitions of n such that the number of odd parts and the number of even parts have opposite parities. Example: a(8)=10 is a count of these partitions: 8, 611, 521, 431, 422, 41111, 332, 32111, 22211, 2111111. - Clark Kimberling, Feb 01 2014, corrected Jan 06 2021
In Chaves 2011 see page 38 equation (3.20). - Michael Somos, Dec 28 2014
Suppose that c(0) = 1, that c(1), c(2), ... are indeterminates, that d(0) = 1, and that d(n) = -c(n) - c(n-1)*d(1) - ... - c(0)*d(n-1). When d(n) is expanded as a polynomial in c(1), c(2),..,c(n), the terms are of the form H*c(i_1)*c(i_2)*...*c(i_k). Let P(n) = [c(i_1), c(i_2), ..., c(i_k)], a partition of n. Then H is negative if P has an odd number of parts, and H is positive if P has an even number of parts. That is, d(n) has A027193(n) negative coefficients, A027187(n) positive coefficients, and A000041 terms. The maximal coefficient in d(n), in absolute value, is A102462(n). - Clark Kimberling, Dec 15 2016

Examples

			G.f. = x + x^2 + 2*x^3 + 2*x^4 + 4*x^5 + 5*x^6 + 8*x^7 + 10*x^8 + 16*x^9 + 20*x^10 + ...
From _Gus Wiseman_, Feb 11 2021: (Start)
The a(1) = 1 through a(8) = 10 partitions into an odd number of parts are the following. The Heinz numbers of these partitions are given by A026424.
  (1)  (2)  (3)    (4)    (5)      (6)      (7)        (8)
            (111)  (211)  (221)    (222)    (322)      (332)
                          (311)    (321)    (331)      (422)
                          (11111)  (411)    (421)      (431)
                                   (21111)  (511)      (521)
                                            (22111)    (611)
                                            (31111)    (22211)
                                            (1111111)  (32111)
                                                       (41111)
                                                       (2111111)
The a(1) = 1 through a(8) = 10 partitions whose greatest part is odd are the following. The Heinz numbers of these partitions are given by A244991.
  (1)  (11)  (3)    (31)    (5)      (33)      (7)        (53)
             (111)  (1111)  (32)     (51)      (52)       (71)
                            (311)    (321)     (322)      (332)
                            (11111)  (3111)    (331)      (521)
                                     (111111)  (511)      (3221)
                                               (3211)     (3311)
                                               (31111)    (5111)
                                               (1111111)  (32111)
                                                          (311111)
                                                          (11111111)
(End)
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 39, Example 7.

Crossrefs

The Heinz numbers of these partitions are A026424 or A244991.
The even-length version is A027187.
The case of odd sum as well as length is A160786, ranked by A340931.
The case of odd maximum as well as length is A340385.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.
A000009 counts partitions into odd parts, ranked by A066208.
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A101707 counts partitions of odd positive rank.

Programs

  • Maple
    g:=sum(x^(2*k)/product(1-x^j,j=1..2*k-1),k=1..40): gser:=series(g,x=0,50): seq(coeff(gser,x,n),n=1..45); # Emeric Deutsch, Apr 05 2006
  • Mathematica
    nn=40;CoefficientList[Series[ Sum[x^(2j+1)Product[1/(1- x^i),{i,1,2j+1}],{j,0,nn}],{x,0,nn}],x]  (* Geoffrey Critzer, Dec 01 2012 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n], OddQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n], OddQ[ First@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n + 1], #[[-1]] == 1 && EvenQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n + 1], EvenQ[ First@#] && (Length[#] < 2 || #[[1]] != #[[2]]) &]]; (* Michael Somos, Dec 28 2014 *)
  • PARI
    {a(n) = if( n<1, 0, polcoeff( sum( k=1, n, if( k%2, x^k / prod( j=1, k, 1 - x^j, 1 + x * O(x^(n-k)) ))), n))}; /* Michael Somos, Jul 24 2012 */
    
  • PARI
    q='q+O('q^66); concat([0], Vec( (1/eta(q)-eta(q)/eta(q^2))/2 ) ) \\ Joerg Arndt, Mar 23 2014

Formula

a(n) = (A000041(n) - (-1)^n*A000700(n)) / 2.
For g.f. see under A027187.
G.f.: Sum(k>=1, x^(2*k-1)/Product(j=1..2*k-1, 1-x^j ) ). - Emeric Deutsch, Apr 05 2006
G.f.: - Sum(k>=1, (-x)^(k^2)) / Product(k>=1, 1-x^k ). - Joerg Arndt, Feb 02 2014
G.f.: Sum(k>=1, x^(k*(2*k-1)) / Product(j=1..2*k, 1-x^j)). - Michael Somos, Dec 28 2014
a(2*n) = A000701(2*n), a(2*n-1) = A046682(2*n-1); a(n) = A000041(n)-A027187(n). - Reinhard Zumkeller, Apr 22 2006

A067659 Number of partitions of n into distinct parts such that number of parts is odd.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 9, 11, 14, 16, 19, 23, 27, 32, 38, 44, 52, 61, 71, 82, 96, 111, 128, 148, 170, 195, 224, 256, 293, 334, 380, 432, 491, 557, 630, 713, 805, 908, 1024, 1152, 1295, 1455, 1632, 1829, 2048, 2291, 2560, 2859, 3189, 3554, 3958, 4404
Offset: 0

Views

Author

Naohiro Nomoto, Feb 23 2002

Keywords

Comments

Ramanujan theta functions: phi(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).

Examples

			From _Gus Wiseman_, Jan 09 2021: (Start)
The a(5) = 1 through a(15) = 14 partitions (A-F = 10..15):
  5   6     7     8     9     A     B     C     D     E     F
      321   421   431   432   532   542   543   643   653   654
                  521   531   541   632   642   652   743   753
                        621   631   641   651   742   752   762
                              721   731   732   751   761   843
                                    821   741   832   842   852
                                          831   841   851   861
                                          921   931   932   942
                                                A21   941   951
                                                      A31   A32
                                                      B21   A41
                                                            B31
                                                            C21
                                                            54321
(End)
		

Crossrefs

Dominates A000009.
Numbers with these strict partitions as binary indices are A000069.
The non-strict version is A027193.
The Heinz numbers of these partitions are A030059.
The even version is A067661.
The version for rank is A117193, with non-strict version A101707.
The ordered version is A332304, with non-strict version A166444.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A174726 counts ordered factorizations of odd length.
- A339890 counts factorizations of odd length.
A008289 counts strict partitions by sum and length.
A026804 counts partitions whose least part is odd, with strict case A026832.

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(n=0, t, add(b(n-i*j, i-1, abs(t-j)), j=0..min(n/i, 1))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..80);  # Alois P. Heinz, Apr 01 2014
  • Mathematica
    b[n_, i_, t_] := b[n, i, t] = If[n > i*(i + 1)/2, 0, If[n == 0, t, Sum[b[n - i*j, i - 1, Abs[t - j]], {j, 0, Min[n/i, 1]}]]]; a[n_] := b[n, n, 0]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Jan 16 2015, after Alois P. Heinz *)
    CoefficientList[Normal[Series[(QPochhammer[-x, x]-QPochhammer[x])/2, {x, 0, 100}]], x] (* Andrey Zabolotskiy, Apr 12 2017 *)
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&OddQ[Length[#]]&]],{n,0,30}] (* Gus Wiseman, Jan 09 2021 *)
  • PARI
    {a(n)=local(A); if(n<0, 0, A=x*O(x^n); polcoeff( (eta(x^2+A)/eta(x+A) - eta(x+A))/2, n))} /* Michael Somos, Feb 14 2006 */
    
  • PARI
    N=66;  q='q+O('q^N);  S=1+2*sqrtint(N);
    gf=sum(n=1,S, (n%2!=0) * q^(n*(n+1)/2) / prod(k=1,n, 1-q^k ) );
    concat( [0], Vec(gf) )  /* Joerg Arndt, Oct 20 2012 */
    
  • PARI
    N=66;  q='q+O('q^N);  S=1+sqrtint(N);
    gf=sum(n=1, S, q^(2*n^2-n) / prod(k=1, 2*n-1, 1-q^k ) );
    concat( [0], Vec(gf) )  \\ Joerg Arndt, Apr 01 2014

Formula

For g.f. see under A067661.
a(n) = (A000009(n)-A010815(n))/2. - Vladeta Jovovic, Feb 24 2002
Expansion of (1-phi(-q))/(2*chi(-q)) in powers of q where phi(),chi() are Ramanujan theta functions. - Michael Somos, Feb 14 2006
G.f.: sum(n>=1, q^(2*n^2-n) / prod(k=1..2*n-1, 1-q^k ) ). [Joerg Arndt, Apr 01 2014]
a(n) = A067661(n) - A010815(n). - Andrey Zabolotskiy, Apr 12 2017
A000009(n) = a(n) + A067661(n). - Gus Wiseman, Jan 09 2021

A339890 Number of odd-length factorizations of n into factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The a(n) factorizations for n = 24, 48, 60, 72, 96, 120:
  24      48          60       72          96          120
  2*2*6   2*3*8       2*5*6    2*4*9       2*6*8       3*5*8
  2*3*4   2*4*6       3*4*5    2*6*6       3*4*8       4*5*6
          3*4*4       2*2*15   3*3*8       4*4*6       2*2*30
          2*2*12      2*3*10   3*4*6       2*2*24      2*3*20
          2*2*2*2*3            2*2*18      2*3*16      2*4*15
                               2*3*12      2*4*12      2*5*12
                               2*2*2*3*3   2*2*2*2*6   2*6*10
                                           2*2*2*3*4   3*4*10
                                                       2*2*2*3*5
		

Crossrefs

The case of set partitions (or n squarefree) is A024429.
The case of partitions (or prime powers) is A027193.
The ordered version is A174726 (even: A174725).
The remaining (even-length) factorizations are counted by A339846.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Maple
    g:= proc(n, k, t) option remember; `if`(n>k, 0, t)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d, 1-t)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> `if`(n<2, 0, g(n$2, 1)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Dec 30 2020
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],OddQ@Length[#]&]],{n,100}]

Formula

a(n) + A339846(n) = A001055(n).

A024430 Expansion of e.g.f. cosh(exp(x)-1).

Original entry on oeis.org

1, 0, 1, 3, 8, 25, 97, 434, 2095, 10707, 58194, 338195, 2097933, 13796952, 95504749, 692462671, 5245040408, 41436754261, 340899165549, 2915100624274, 25857170687507, 237448494222575, 2253720620740362, 22078799199129799, 222987346441156585, 2319210969809731600
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of an n-element set into an even number of classes.
Let A(0) = 1, B(0) = 0; A(n+1) = Sum_{k=0..n} binomial(n,k)*B(k), B(n+1) = Sum_{k=0..n} binomial(n,k)*A(k); entry gives A sequence (cf. A024429).

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 226, 5th line of table.
  • S. K. Ghosal, J. K. Mandal, Stirling Transform Based Color Image Authentication, Procedia Technology, 2013 Volume 10, 2013, Pages 95-104.
  • L. Lovasz, Combinatorial Problems and Exercises, North-Holland, 1993, pp. 15, 148.

Crossrefs

Programs

  • GAP
    List([0..25], n-> Sum([0..Int(n/2)], k-> Stirling2(n,2*k)) ); # G. C. Greubel, Oct 09 2019
  • Magma
    a:= func< n | (&+[StirlingSecond(n,2*k): k in [0..Floor(n/2)]]) >;
    [a(n): n in [0..25]]; // G. C. Greubel, Oct 09 2019
    
  • Maple
    b:= proc(n, t) option remember; `if`(n=0, t, add(
           b(n-j, 1-t)*binomial(n-1, j-1), j=1..n))
        end:
    a:= n-> b(n, 1):
    seq(a(n), n=0..28);  # Alois P. Heinz, Jan 15 2018
    with(combinat); seq((bell(n) + BellB(n, -1))/2, n = 0..20); # G. C. Greubel, Oct 09 2019
  • Mathematica
    nn=20;a=Exp[Exp[x]-1];Range[0,nn]!CoefficientList[Series[(a+1/a)/2,{x,0,nn}],x]  (* Geoffrey Critzer, Nov 04 2012 *)
    Table[(BellB[n] + BellB[n, -1])/2, {n, 0, 20}] (* Vladimir Reshetnikov, Nov 01 2015 *)
  • PARI
    {a(n)=polcoeff(sum(m=0, n, x^(2*m)/prod(k=1, 2*m, 1-k*x +x*O(x^n))), n)} \\ Paul D. Hanna, Sep 05 2012
    
  • Sage
    def A024430(n) :
        return add(stirling_number2(n,i) for i in range(0,n+(n+1)%2,2))
    # Peter Luschny, Feb 28 2012
    

Formula

a(n) = S(n, 2) + S(n, 4) + ... + S(n, 2k), where k = [ n/2 ], S(i, j) are Stirling numbers of second kind.
E.g.f.: cosh(exp(x)-1). - N. J. A. Sloane, Jan 28 2001
a(n) = (A000110(n) + A000587(n)) / 2. - Peter Luschny, Apr 25 2011
O.g.f.: Sum_{n>=0} x^(2*n) / Product_{k=0..2*n} (1 - k*x). - Paul D. Hanna, Sep 05 2012
G.f.: G(0)/(1+x) where G(k) = 1 - x*(2*k+1)/((2*x*k-1) - x*(2*x*k-1)/(x - (2*k+1)*(2*x*k+x-1)/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 05 2013
G.f.: G(0)/(1+2*x) where G(k) = 1 - 2*x*(k+1)/((2*x*k-1) - x*(2*x*k-1)/(x - 2*(k+1)*(2*x*k+x-1)/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 05 2013
a(n) ~ n^n / (2 * (LambertW(n))^n * exp(n+1-n/LambertW(n)) * sqrt(1+LambertW(n))). - Vaclav Kotesovec, Aug 04 2014

Extensions

Description changed by N. J. A. Sloane, Jun 14 2003 and again Sep 05 2006

A174726 a(n) = (A002033(n-1) - A008683(n))/2.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 4, 1, 4, 1, 4, 1, 1, 1, 10, 1, 1, 2, 4, 1, 7, 1, 8, 1, 1, 1, 13, 1, 1, 1, 10, 1, 7, 1, 4, 4, 1, 1, 24, 1, 4, 1, 4, 1, 10, 1, 10, 1, 1, 1, 22, 1, 1, 4, 16, 1, 7, 1, 4, 1, 7, 1, 38, 1, 1, 4, 4, 1
Offset: 1

Views

Author

Mats Granvik, Mar 28 2010

Keywords

Comments

a(n) is the number of permutation matrices with a negative contribution to the determinant that is the Möbius function. See A174725 for how the determinant is defined. - Mats Granvik, May 26 2017
From Gus Wiseman, Jan 04 2021: (Start)
Also the number of ordered factorizations of n into an odd number of factors > 1. The unordered case is A339890. For example, the a(n) factorizations for n = 8, 12, 24, 30, 32, 36 are:
(8) (12) (24) (30) (32) (36)
(2*2*2) (2*2*3) (2*2*6) (2*3*5) (2*2*8) (2*2*9)
(2*3*2) (2*3*4) (2*5*3) (2*4*4) (2*3*6)
(3*2*2) (2*4*3) (3*2*5) (2*8*2) (2*6*3)
(2*6*2) (3*5*2) (4*2*4) (2*9*2)
(3*2*4) (5*2*3) (4*4*2) (3*2*6)
(3*4*2) (5*3*2) (8*2*2) (3*3*4)
(4*2*3) (2*2*2*2*2) (3*4*3)
(4*3*2) (3*6*2)
(6*2*2) (4*3*3)
(6*2*3)
(6*3*2)
(9*2*2)
(End)

Crossrefs

The even version is A174725.
The unordered case is A339890, with even version A339846.
A001055 counts factorizations, with strict case A045778.
A074206 counts ordered factorizations, with strict case A254578.
A251683 counts ordered factorizations by product and length.
A340102 counts odd-length factorizations into odd factors.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A027193 counts partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A332304 counts strict compositions of odd length.

Programs

  • Mathematica
    ordfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@ordfacs[n/d],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[ordfacs[n],OddQ@*Length]],{n,100}] (* Gus Wiseman, Jan 04 2021 *)

Formula

a(n) = (A002033(n-1) - A008683(n))/2. - Mats Granvik, May 26 2017
For n > 0, a(n) + A174725(n) = A074206(n). - Gus Wiseman, Jan 04 2021

A340607 Number of factorizations of n into an odd number of factors > 1, the greatest of which is odd.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 25 2021

Keywords

Examples

			The a(n) factorizations for n = 27, 84, 108, 180, 252, 360, 432:
  27     2*6*7   2*6*9      4*5*9      4*7*9      5*8*9       6*8*9
  3*3*3  3*4*7   3*4*9      2*2*45     6*6*7      2*4*45      2*8*27
         2*2*21  2*2*27     2*6*15     2*2*63     3*8*15      4*4*27
                 2*2*3*3*3  3*4*15     2*6*21     4*6*15      2*2*2*6*9
                            2*2*3*3*5  3*4*21     2*12*15     2*2*3*4*9
                                       2*2*3*3*7  2*2*2*5*9   2*2*2*2*27
                                                  2*3*3*4*5   2*2*2*2*3*3*3
                                                  2*2*2*3*15
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The case of odd length only is A339890.
The case of all odd factors is A340102.
The version for partitions is A340385.
The version for prime indices is A340386.
The case of odd maximum only is A340831.
A000009 counts partitions into odd parts (A066208).
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length/maximum (A026424/A244991).
A058695 counts partitions of odd numbers (A300063).
A078408 counts odd-length partitions into odd numbers (A300272).
A316439 counts factorizations by sum and length.
A340101 counts factorizations (into odd factors = of odd numbers).
A340832 counts factorizations whose least part is odd.
A340854/A340855 lack/have a factorization with odd minimum.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],OddQ[Length[#]]&&OddQ[Max@@#]&]],{n,100}]
  • PARI
    A340607(n, m=n, k=0, grodd=0) = if(1==n, k, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&(grodd||(d%2)), s += A340607(n/d, d, 1-k, bitor(1,grodd)))); (s)); \\ Antti Karttunen, Dec 13 2021

Extensions

Data section extended up to 108 terms by Antti Karttunen, Dec 13 2021

A340854 Numbers that cannot be factored into factors > 1, the least of which is odd.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 14, 16, 20, 22, 26, 28, 32, 34, 38, 44, 46, 52, 58, 62, 64, 68, 74, 76, 82, 86, 88, 92, 94, 104, 106, 116, 118, 122, 124, 128, 134, 136, 142, 146, 148, 152, 158, 164, 166, 172, 178, 184, 188, 194, 202, 206, 212, 214, 218, 226, 232, 236, 244
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

Consists of 1 and all numbers that are even and have no odd divisor 1 < d <= n/d.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}              44: {1,1,5}          106: {1,16}
      2: {1}             46: {1,9}            116: {1,1,10}
      4: {1,1}           52: {1,1,6}          118: {1,17}
      6: {1,2}           58: {1,10}           122: {1,18}
      8: {1,1,1}         62: {1,11}           124: {1,1,11}
     10: {1,3}           64: {1,1,1,1,1,1}    128: {1,1,1,1,1,1,1}
     14: {1,4}           68: {1,1,7}          134: {1,19}
     16: {1,1,1,1}       74: {1,12}           136: {1,1,1,7}
     20: {1,1,3}         76: {1,1,8}          142: {1,20}
     22: {1,5}           82: {1,13}           146: {1,21}
     26: {1,6}           86: {1,14}           148: {1,1,12}
     28: {1,1,4}         88: {1,1,1,5}        152: {1,1,1,8}
     32: {1,1,1,1,1}     92: {1,1,9}          158: {1,22}
     34: {1,7}           94: {1,15}           164: {1,1,13}
     38: {1,8}          104: {1,1,1,6}        166: {1,23}
For example, the factorizations of 88 are (2*2*2*11), (2*2*22), (2*4*11), (2*44), (4*22), (8*11), (88), none of which has odd minimum, so 88 is in the sequence.
		

Crossrefs

The version looking at greatest factor is A000079.
The version for twice-balanced is A340656, with complement A340657.
These factorization are counted by A340832.
The complement is A340855.
A033676 selects the maximum inferior divisor.
A038548 counts inferior divisors.
A055396 selects the least prime index.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A339890 counts factorizations of odd length.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A024429 counts set partitions of odd length.
A026424 lists numbers with odd Omega.
A066208 lists Heinz numbers of partitions into odd parts.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.

Programs

  • Mathematica
    Select[Range[100],Function[n,n==1||EvenQ[n]&&Select[Rest[Divisors[n]],OddQ[#]&&#<=n/#&]=={}]]

A340855 Numbers that can be factored into factors > 1, the least of which is odd.

Original entry on oeis.org

3, 5, 7, 9, 11, 12, 13, 15, 17, 18, 19, 21, 23, 24, 25, 27, 29, 30, 31, 33, 35, 36, 37, 39, 40, 41, 42, 43, 45, 47, 48, 49, 50, 51, 53, 54, 55, 56, 57, 59, 60, 61, 63, 65, 66, 67, 69, 70, 71, 72, 73, 75, 77, 78, 79, 80, 81, 83, 84, 85, 87, 89, 90, 91, 93, 95
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

These are numbers that are odd or have an odd divisor 1 < d <= n/d.

Examples

			The sequence of terms together with their prime indices begins:
     3: {2}          27: {2,2,2}      48: {1,1,1,1,2}
     5: {3}          29: {10}         49: {4,4}
     7: {4}          30: {1,2,3}      50: {1,3,3}
     9: {2,2}        31: {11}         51: {2,7}
    11: {5}          33: {2,5}        53: {16}
    12: {1,1,2}      35: {3,4}        54: {1,2,2,2}
    13: {6}          36: {1,1,2,2}    55: {3,5}
    15: {2,3}        37: {12}         56: {1,1,1,4}
    17: {7}          39: {2,6}        57: {2,8}
    18: {1,2,2}      40: {1,1,1,3}    59: {17}
    19: {8}          41: {13}         60: {1,1,2,3}
    21: {2,4}        42: {1,2,4}      61: {18}
    23: {9}          43: {14}         63: {2,2,4}
    24: {1,1,1,2}    45: {2,2,3}      65: {3,6}
    25: {3,3}        47: {15}         66: {1,2,5}
For example, 72 is in the sequence because it has three suitable factorizations: (3*3*8), (3*4*6), (3*24).
		

Crossrefs

The version looking at greatest factor is A057716.
The version for twice-balanced is A340657, with complement A340656.
These factorization are counted by A340832.
The complement is A340854.
A033676 selects the maximum inferior divisor.
A038548 counts inferior divisors, listed by A161906.
A055396 selects the least prime index.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A339890 counts factorizations of odd length.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A024429 counts set partitions of odd length.
A026424 lists numbers with odd Omega.
A066208 lists Heinz numbers of partitions into odd parts.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.
A332304 counts strict compositions of odd length.
A340692 counts partitions of odd rank.

Programs

  • Mathematica
    Select[Range[100],Function[n,n>1&&(OddQ[n]||Select[Rest[Divisors[n]],OddQ[#]&&#<=n/#&]!={})]]

A340385 Number of integer partitions of n into an odd number of parts, the greatest of which is odd.

Original entry on oeis.org

1, 0, 2, 0, 3, 1, 6, 3, 10, 7, 18, 15, 30, 28, 51, 50, 82, 87, 134, 145, 211, 235, 331, 375, 510, 586, 779, 901, 1172, 1366, 1750, 2045, 2581, 3026, 3778, 4433, 5476, 6430, 7878, 9246, 11240, 13189, 15931, 18670, 22417, 26242, 31349, 36646, 43567, 50854
Offset: 1

Views

Author

Gus Wiseman, Jan 08 2021

Keywords

Examples

			The a(3) = 2 through a(10) = 7 partitions:
  3     5       321   7         332     9           532
  111   311           322       521     333         541
        11111         331       32111   522         721
                      511               531         32221
                      31111             711         33211
                      1111111           32211       52111
                                        33111       3211111
                                        51111
                                        3111111
                                        111111111
		

Crossrefs

Partitions of odd length are counted by A027193, ranked by A026424.
Partitions with odd maximum are counted by A027193, ranked by A244991.
The Heinz numbers of these partitions are given by A340386.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.
A000009 counts partitions into odd parts, ranked by A066208.
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A101707 counts partitions with odd rank.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OddQ[Length[#]*Max[#]]&]],{n,30}]

A143815 Let A(0)=1, B(0)=0 and C(0)=0. Let B(n+1) = Sum_{k = 0..n} binomial(n,k)*A(k), C(n+1) = Sum_{k = 0..n} binomial(n,k)*B(k) and A(n+1) = Sum_{k = 0..n} binomial(n,k)*C(k). This entry gives the sequence A(n).

Original entry on oeis.org

1, 0, 0, 1, 6, 25, 91, 322, 1232, 5672, 32202, 209143, 1432454, 9942517, 69363840, 490303335, 3565609732, 27118060170, 218183781871, 1861370544934, 16729411124821, 156706028787827, 1514442896327792, 14999698898942772, 151838974745743228, 1571513300578303070
Offset: 0

Views

Author

Peter Bala, Sep 03 2008

Keywords

Comments

Compare with A024429 and A024430.
This sequence and its companion sequences B(n) = A143816(n) and C(n) = A143817(n) may be viewed as generalizations of the Bell numbers A000110. Define a sequence R(n) of real numbers by R(n) = Sum_{k >= 0} (3*k)^n/(3*k)! for n = 0, 1, 2, .... It is easy to verify that this sequence satisfies the recurrence relation u(n+3) = 3*u(n+2) - 2*u(n+1) + Sum_{i = 0..n} binomial(n, i)*3^(n-i)*u(i). Hence R(n) is an integral linear combination of R(0), R(1) and R(2). Some examples are given below.
To find the precise form of the linear relation define two other sequences of real numbers by S(n) = Sum_{k >= 0} (3*k+1)^n/(3*k+1)! and T(n) = Sum_{k >= 0} (3*k+2)^n/(3*k+2)! for n = 0, 1, 2, .... Both S(n) and T(n) satisfy the above recurrence. Then by means of the identities S(n+1) = Sum_{i = 0..n} binomial(n, i)*R(i), T(n+1) = Sum_{i = 0..n} binomial(n, i)*S(i) and R(n+1) = Sum_{i = 0..n} binomial(n, i)*T(i) we obtain the result R(n) = A(n)*R(0) + (B(n) - C(n))*R(1) + C(n)*R(2) = A(n)*R(0) + B(n)*R(1) + C(n)*(R(2) - R(1)) (with corresponding expressions for S(n) and T(n)). This generalizes the Dobinski's relation for the Bell numbers Sum_{k >= 0} k^n/k! = A000110(n)*exp(1).
Some examples of R(n) as a linear combination of R(0), R(1) and R(2) - R(1) are given below. The decimal expansions of R(0) = 1 + 1/3! + 1/6! + 1/9! + ..., R(2) - R(1) = 1/1! + 1/4! + 1/7! + ... and R(1) = 1/2! + 1/5! + 1/8! + ... may be found in A143819, A143820 and A143821 respectively. Compare with A143628 through A143631.
For n > 0, the number of partitions of {1,2,...,n} into 3,6,9,... classes. - Geoffrey Critzer, Mar 05 2010

Examples

			R(n) as a linear combination of R(i), i = 0..2.
  ==================================
  R(n)  |     R(0)    R(1)    R(2)
  ==================================
  R(3)  |       1      -2       3
  R(4)  |       6      -5       7
  R(5)  |      25      -5      16
  R(6)  |      91      20      46
  R(7)  |     322     149     203
  R(8)  |    1232     552    1178
  R(9)  |    5672     991    7242
  R(10) |   32202   -3799   43786
  ...
Column 2 of the above table is A143818.
R(n) as a linear combination of R(0),R(1) and R(2) - R(1).
  =====================================
  R(n)  |     R(0)     R(1)   R(2)-R(1)
  =====================================
  R(3)  |       1        1        3
  R(4)  |       6        2        7
  R(5)  |      25       11       16
  R(6)  |      91       66       46
  R(7)  |     322      352      203
  R(8)  |    1232     1730     1178
  R(9)  |    5672     8233     7242
  R(10) |   32202    39987    43786
  ...
		

Crossrefs

Programs

  • Maple
    # (1)
    M:=24: a:=array(0..100): b:=array(0..100): c:=array(0..100):
    a[0]:=1: b[0]:=0: c[0]:=0:
    for n from 1 to M do
    b[n]:=add(binomial(n-1,k)*a[k], k=0..n-1);
    c[n]:=add(binomial(n-1,k)*b[k], k=0..n-1);
    a[n]:=add(binomial(n-1,k)*c[k], k=0..n-1);
    end do:
    A143815:=[seq(a[n], n=0..M)];
    # (2)
    seq(add(Stirling2(n,3*i),i = 0..floor(n/3)), n = 0..24);
    # third Maple program:
    b:= proc(n, t) option remember; `if`(n=0, irem(t, 2),
          add(b(n-j, irem(t+1, 3))*binomial(n-1, j-1), j=1..n))
        end:
    a:= n-> b(n, 1):
    seq(a(n), n=0..25);  # Alois P. Heinz, Feb 20 2018
  • Mathematica
    a = Exp[x] - 1; f[x_] := 1/3 (E^x + 2 E^(-x/2) Cos[(Sqrt[3] x)/2]); CoefficientList[Series[f[a], {x, 0, 25}], x]*Table[n!, {n, 0, 25}] (* Geoffrey Critzer, Mar 05 2010 *)
  • PARI
    Bell_poly(n, x) = exp(-x)*suminf(k=0, k^n*x^k/k!);
    a(n) = my(w=(-1+sqrt(3)*I)/2); round(Bell_poly(n, 1)+Bell_poly(n, w)+Bell_poly(n, w^2))/3; \\ Seiichi Manyama, Oct 13 2022

Formula

a(n) = Sum_{k = 0..floor(n/3)} Stirling2(n, 3*k).
Let w = exp(2*Pi*i/3) and set F(x) = (exp(x) + exp(w*x) + exp(w^2*x))/3 = 1 + x^3/3! + x^6/6! + ... . Then the e.g.f. for the sequence is F(exp(x) - 1).
A143815(n) + A143816(n) + A143817(n) = Bell(n).
E.g.f. is B(A(x)) where A(x) = exp(x) - 1 and B(x) = (1/3)*(exp(x) + 2*exp(-x/2)*cos(sqrt(3)*x/2)). - Geoffrey Critzer, Mar 05 2010
a(n) = ( Bell_n(1) + Bell_n(w) + Bell_n(w^2) )/3, where Bell_n(x) is n-th Bell polynomial and w = exp(2*Pi*i/3). - Seiichi Manyama, Oct 13 2022
a(n) ~ n^n / (3 * (LambertW(n))^n * exp(n+1-n/LambertW(n)) * sqrt(1+LambertW(n))). - Vaclav Kotesovec, Jun 10 2025
Showing 1-10 of 38 results. Next