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 35 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

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

A244991 Numbers whose greatest prime factor is a prime with an odd index; n such that A006530(n) is in A031368.

Original entry on oeis.org

2, 4, 5, 8, 10, 11, 15, 16, 17, 20, 22, 23, 25, 30, 31, 32, 33, 34, 40, 41, 44, 45, 46, 47, 50, 51, 55, 59, 60, 62, 64, 66, 67, 68, 69, 73, 75, 77, 80, 82, 83, 85, 88, 90, 92, 93, 94, 97, 99, 100, 102, 103, 109, 110, 115, 118, 119, 120, 121, 123, 124, 125, 127, 128
Offset: 1

Views

Author

Antti Karttunen, Jul 21 2014

Keywords

Comments

Equally, numbers n for which A061395(n) is odd.
A122111 maps each one of these numbers to a unique term of A026424 and vice versa.
If the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), these are the Heinz numbers of partitions whose greatest part is odd, counted by A027193. - Gus Wiseman, Feb 08 2021

Examples

			From _Gus Wiseman_, Feb 08 2021: (Start)
The sequence of terms together with their prime indices begins:
      2: {1}           32: {1,1,1,1,1}     64: {1,1,1,1,1,1}
      4: {1,1}         33: {2,5}           66: {1,2,5}
      5: {3}           34: {1,7}           67: {19}
      8: {1,1,1}       40: {1,1,1,3}       68: {1,1,7}
     10: {1,3}         41: {13}            69: {2,9}
     11: {5}           44: {1,1,5}         73: {21}
     15: {2,3}         45: {2,2,3}         75: {2,3,3}
     16: {1,1,1,1}     46: {1,9}           77: {4,5}
     17: {7}           47: {15}            80: {1,1,1,1,3}
     20: {1,1,3}       50: {1,3,3}         82: {1,13}
     22: {1,5}         51: {2,7}           83: {23}
     23: {9}           55: {3,5}           85: {3,7}
     25: {3,3}         59: {17}            88: {1,1,1,5}
     30: {1,2,3}       60: {1,1,2,3}       90: {1,2,2,3}
     31: {11}          62: {1,11}          92: {1,1,9}
(End)
		

Crossrefs

Complement: A244990.
Looking at least instead of greatest prime index gives A026804.
The partitions with these Heinz numbers are counted by A027193.
The case where Omega is odd also is A340386.
A001222 counts prime factors.
A056239 adds up prime indices.
A300063 ranks partitions of odd numbers.
A061395 selects maximum prime index.
A066208 ranks partitions into odd parts.
A112798 lists the prime indices of each positive integer.
A340931 ranks odd-length partitions of odd numbers.

Programs

  • Mathematica
    Select[Range[100],OddQ[PrimePi[FactorInteger[#][[-1,1]]]]&] (* Gus Wiseman, Feb 08 2021 *)

Formula

For all n, A244989(a(n)) = n.

A116882 A number k is included if (highest odd divisor of k)^2 <= k.

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160, 176, 192, 208, 224, 240, 256, 288, 320, 352, 384, 416, 448, 480, 512, 544, 576, 608, 640, 672, 704, 736, 768, 800, 832, 864, 896, 928, 960, 992, 1024, 1088, 1152, 1216, 1280, 1344, 1408
Offset: 1

Views

Author

Leroy Quet, Feb 24 2006

Keywords

Comments

Also k is included if (and only if) the greatest power of 2 dividing k is >= the highest odd divisor of k. All terms of the sequence are even besides the 1.
Equivalently, positive integers of the form k*2^m, where odd k <= 2^m. - Thomas Ordowski, Oct 19 2014
If we define a divisor d|n to be superior if d >= n/d, then superior divisors are counted by A038548 and listed by A161908. This sequence consists of 1 and all numbers without a superior odd divisor. - Gus Wiseman, Feb 18 2021
Numbers k such that A006519(k) >= A000265(k), with equality only when k = 1. - Amiram Eldar, Jan 24 2023

Examples

			40 = 8 * 5, where 8 is highest power of 2 dividing 40 and 5 is the highest odd dividing 40. 8 is >= 5 (so 5^2 <= 40), so 40 is in the sequence.
		

Crossrefs

The complement is A116883.
Positions of zeros (and 1) in A341675.
A051283 = numbers without a superior prime-power divisor (zeros of A341593).
A059172 = numbers without a superior squarefree divisor (zeros of A341592).
A063539 = numbers without a superior prime divisor (zeros of A341591).
A333805 counts strictly inferior odd divisors.
A341594 counts strictly superior odd divisors.
- Strictly Inferior: A056924, A060775, A070039, A333806, A341596, A341674.
Subsequence of A082662, {1} U A363122.

Programs

  • Mathematica
    f[n_] := Select[Divisors[n], OddQ[ # ] &][[ -1]]; Insert[Select[Range[2, 1500], 2^FactorInteger[ # ][[1]][[2]] > f[ # ] &], 1, 1] (* Stefan Steinerberger, Apr 10 2006 *)
    q[n_] := 2^(2*IntegerExponent[n, 2]) >= n; Select[Range[1500], q] (* Amiram Eldar, Jan 24 2023 *)
  • PARI
    isok(n) = vecmax(select(x->((x % 2)==1), divisors(n)))^2 <= n; \\ Michel Marcus, Sep 06 2016
    
  • PARI
    isok(n) = 2^(valuation(n,2)*2) >= n \\ Jeppe Stig Nielsen, Feb 19 2019
    
  • Python
    from itertools import count, islice
    def A116882_gen(startvalue=1): # generator of terms >= startvalue
        return filter(lambda n:(n&-n)**2>=n,count(max(startvalue,1)))
    A116882_list = list(islice(A116882_gen(),20)) # Chai Wah Wu, May 17 2023

Formula

a(n) = A080075(n-1)-1. - Klaus Brockhaus, Georgi Guninski and M. F. Hasler, Aug 16 2010
a(n) ~ n^2/2. - Thomas Ordowski, Oct 19 2014
Sum_{n>=1} 1/a(n) = 1 + (3/4) * Sum_{k>=1} H(2^k-1)/2^k = 2.3388865091..., where H(k) = A001008(k)/A002805(k) is the k-th harmonic number. - Amiram Eldar, Jan 24 2023

Extensions

More terms from Stefan Steinerberger, Apr 10 2006

A101707 Number of partitions of n having positive odd rank (the rank of a partition is the largest part minus the number of parts).

Original entry on oeis.org

0, 0, 1, 0, 2, 1, 4, 2, 7, 6, 13, 11, 22, 22, 38, 39, 63, 69, 103, 114, 165, 189, 262, 301, 407, 475, 626, 733, 950, 1119, 1427, 1681, 2118, 2503, 3116, 3678, 4539, 5360, 6559, 7735, 9400, 11076, 13372, 15728, 18886, 22184, 26501, 31067, 36947, 43242, 51210, 59818, 70576, 82291, 96750
Offset: 0

Views

Author

Emeric Deutsch, Dec 12 2004

Keywords

Comments

a(n) + A101708(n) = A064173(n).

Examples

			a(7)=2 because the only partitions of 7 with positive odd rank are 421 (rank=1) and 52 (rank=3).
From _Gus Wiseman_, Feb 07 2021: (Start)
Also the number of integer partitions of n into an even number of parts, the greatest of which is odd. For example, the a(2) = 1 through a(10) = 13 partitions (empty column indicated by dot) are:
  11   .  31     32   33       52     53         54       55
          1111        51       3211   71         72       73
                      3111            3221       3222     91
                      111111          3311       3321     3322
                                      5111       5211     3331
                                      311111     321111   5221
                                      11111111            5311
                                                          7111
                                                          322111
                                                          331111
                                                          511111
                                                          31111111
                                                          1111111111
Also the number of integer partitions of n into an odd number of parts, the greatest of which is even. For example, the a(2) = 1 through a(10) = 13 partitions (empty column indicated by dot, A = 10) are:
  2   .  4     221   6       421     8         432       A
         211         222     22111   422       441       433
                     411             431       621       442
                     21111           611       22221     622
                                     22211     42111     631
                                     41111     2211111   811
                                     2111111             22222
                                                         42211
                                                         43111
                                                         61111
                                                         2221111
                                                         4111111
                                                         211111111
(End)
		

References

  • George E. Andrews, The Theory of Partitions, Addison-Wesley, Reading, Mass., 1976.

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The even-rank version is A101708 (A340605).
The even- but not necessarily positive-rank version is A340601 (A340602).
The Heinz numbers of these partitions are (A340604).
Allowing negative odd ranks gives A340692 (A340603).
- Rank -
A047993 counts balanced (rank zero) partitions (A106529).
A064173 counts partitions of positive/negative rank (A340787/A340788).
A064174 counts partitions of nonpositive/nonnegative rank (A324521/A324562).
A101198 counts partitions of rank 1 (A325233).
A257541 gives the rank of the partition with Heinz number n.
- Odd -
A000009 counts partitions into odd parts (A066208).
A026804 counts partitions whose least part is odd.
A027193 counts partitions of odd length/maximum (A026424/A244991).
A058695 counts partitions of odd numbers (A300063).
A339890 counts factorizations of odd length.
A340385 counts partitions of odd length and maximum (A340386).

Programs

  • Maple
    b:= proc(n, i, r) option remember; `if`(n=0, max(0, r),
          `if`(i<1, 0, b(n, i-1, r) +b(n-i, min(n-i, i), 1-
          `if`(r<0, irem(i, 2), r))))
        end:
    a:= n-> b(n$2, -1)/2:
    seq(a(n), n=0..55);  # Alois P. Heinz, Jan 29 2021
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],EvenQ[Length[#]]&&OddQ[Max[#]]&]],{n,0,30}] (* Gus Wiseman, Feb 10 2021 *)
    b[n_, i_, r_] := b[n, i, r] = If[n == 0, Max[0, r],
         If[i < 1, 0, b[n, i - 1, r] + b[n - i, Min[n - i, i], 1 -
         If[r < 0, Mod[i, 2], r]]]];
    a[n_] := b[n, n, -1]/2;
    a /@ Range[0, 55] (* Jean-François Alcover, May 23 2021, after Alois P. Heinz *)

Formula

a(n) = (A000041(n) - A000025(n))/4. - Vladeta Jovovic, Dec 14 2004
G.f.: Sum((-1)^(k+1)*x^((3*k^2+k)/2)/(1+x^k), k=1..infinity)/Product(1-x^k, k=1..infinity). - Vladeta Jovovic, Dec 20 2004
a(n) = A340692(n)/2. - Gus Wiseman, Feb 07 2021

Extensions

More terms from Joerg Arndt, Oct 07 2012
a(0)=0 prepended by Alois P. Heinz, Jan 29 2021

A341446 Heinz numbers of integer partitions whose only odd part is the smallest.

Original entry on oeis.org

2, 5, 6, 11, 14, 17, 18, 23, 26, 31, 35, 38, 41, 42, 47, 54, 58, 59, 65, 67, 73, 74, 78, 83, 86, 95, 97, 98, 103, 106, 109, 114, 122, 126, 127, 137, 142, 143, 145, 149, 157, 158, 162, 167, 174, 178, 179, 182, 185, 191, 197, 202, 209, 211, 214, 215, 222, 226
Offset: 1

Views

Author

Gus Wiseman, Feb 12 2021

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are numbers whose only odd prime index (counting multiplicity) is the smallest.

Examples

			The sequence of partitions together with their Heinz numbers begins:
      2: (1)         54: (2,2,2,1)    109: (29)
      5: (3)         58: (10,1)       114: (8,2,1)
      6: (2,1)       59: (17)         122: (18,1)
     11: (5)         65: (6,3)        126: (4,2,2,1)
     14: (4,1)       67: (19)         127: (31)
     17: (7)         73: (21)         137: (33)
     18: (2,2,1)     74: (12,1)       142: (20,1)
     23: (9)         78: (6,2,1)      143: (6,5)
     26: (6,1)       83: (23)         145: (10,3)
     31: (11)        86: (14,1)       149: (35)
     35: (4,3)       95: (8,3)        157: (37)
     38: (8,1)       97: (25)         158: (22,1)
     41: (13)        98: (4,4,1)      162: (2,2,2,2,1)
     42: (4,2,1)    103: (27)         167: (39)
     47: (15)       106: (16,1)       174: (10,2,1)
		

Crossrefs

These partitions are counted by A035363 (shifted left once).
Terms of A340932 can be factored into elements of this sequence.
The even version is A341447.
A001222 counts prime factors.
A005408 lists odd numbers.
A026804 counts partitions whose smallest part is odd.
A027193 counts odd-length partitions, ranked by A026424.
A031368 lists odd-indexed primes.
A032742 selects largest proper divisor.
A055396 selects smallest prime index.
A056239 adds up prime indices.
A058695 counts partitions of odd numbers, ranked by A300063.
A061395 selects largest prime index.
A066207 lists numbers with all even prime indices.
A066208 lists numbers with all odd prime indices.
A112798 lists the prime indices of each positive integer.
A244991 lists numbers whose greatest prime index is odd.
A340932 lists numbers whose smallest prime index is odd.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[2,100],OddQ[First[primeMS[#]]]&&And@@EvenQ[Rest[primeMS[#]]]&]

Formula

Also numbers n > 1 such that A055396(n) is odd and A032742(n) belongs to A066207.

A026805 Number of partitions of n in which the least part is even.

Original entry on oeis.org

0, 1, 0, 2, 1, 3, 2, 6, 5, 9, 9, 16, 17, 26, 28, 42, 48, 66, 77, 105, 122, 160, 189, 245, 290, 368, 436, 547, 650, 804, 954, 1174, 1390, 1693, 2004, 2425, 2865, 3445, 4060, 4858, 5716, 6802, 7986, 9468, 11087, 13088, 15298, 17995, 20987, 24604, 28631, 33464
Offset: 1

Views

Author

Keywords

Comments

Also number of partitions of n in which the largest part occurs an even number of times. Example: a(6)=3 because we have [3,3],[2,2,1,1] and [1,1,1,1,1,1]. - Emeric Deutsch, Apr 04 2006

Examples

			a(6)=3 because we have [6],[4,2] and [2,2,2].
		

Programs

  • Maple
    g:=sum(x^(2*k)/(1-x^(2*k))/product(1-x^j,j=1..k-1),k=1..40): gser:=series(g,x=0,52): seq(coeff(gser,x,n),n=1..49); # Emeric Deutsch, Apr 04 2006
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n<1 or i<1, 0, b(n, i-1)+
          `if`(n=i, 1-irem(n, 2), 0)+`if`(i>n, 0, b(n-i, i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=1..60);  # Alois P. Heinz, Jul 26 2015
  • Mathematica
    b[n_, i_] := b[n, i] = If[n<1 || i<1, 0, b[n, i-1] + If[n==i, 1-Mod[n, 2], 0] + If[i>n, 0, b[n-i, i]]]; a[n_] := b[n, n]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Oct 28 2015, after Alois P. Heinz *)

Formula

From Vladeta Jovovic, Aug 26 2003: (Start)
G.f.: Sum_{k>=2} ((-1)^k*(-1+1/Product_{i>=k} (1-x^i))).
a(n) = Sum_{k=2..n} (-1)^k*A026807(n, k) = A000041(n)-A026804(n). (End)
From Emeric Deutsch, Apr 04 2006: (Start)
G.f.: Sum_{k>=1}(x^(2k)/Product_{j>=2k}(1-x^j)).
G.f.: Sum_{k>=1}(x^(2k)/((1-x^(2k))*Product_{j=1..k-1}(1-x^j))). (End)
a(n) ~ Pi * exp(Pi*sqrt(2*n/3)) / (3 * 2^(5/2) * n^(3/2)) * (1 - (3*sqrt(3/2)/Pi + 61*Pi/(24*sqrt(6))) / sqrt(n)). - Vaclav Kotesovec, Jul 06 2019, extended Nov 02 2019

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/#&]=={}]]
Showing 1-10 of 35 results. Next