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 31-40 of 206 results. Next

A333805 Number of odd divisors of n that are < sqrt(n).

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, Apr 05 2020

Keywords

Comments

If we define a divisor d|n to be strictly inferior if d < n/d, then strictly inferior divisors are counted by A056924 and listed by A341674. This sequence counts strictly inferior odd divisors. - Gus Wiseman, Feb 26 2021

Examples

			The strictly inferior odd divisors of 945 are 1, 3, 5, 7, 9, 15, 21, 27, so a(945) = 8. - _Gus Wiseman_, Feb 27 2021
		

Crossrefs

Dominated by A001227 (number of odd divisors).
Strictly inferior divisors (not just odd) are counted by A056924.
The non-strict version is A069288.
These divisors add up to A070039.
The case of prime divisors is A333806.
The strictly superior version is A341594.
The case of squarefree divisors is A341596.
The superior version is A341675.
The case of prime-power divisors is A341677.
A006530 selects the greatest prime factor.
A020639 selects the smallest prime factor.
- Odd -
A000009 counts partitions into odd parts, ranked by A066208.
A026424 lists numbers with odd Omega.
A027193 counts odd-length partitions.
A067659 counts strict partitions of odd length, ranked by A030059.
- Inferior divisors -
A033676 selects the greatest inferior divisor.
A033677 selects the smallest superior divisor.
A038548 counts superior (or inferior) divisors.
A060775 selects the greatest strictly inferior divisor.
A341674 lists strictly inferior divisors.

Programs

  • Mathematica
    Table[DivisorSum[n, 1 &, # < Sqrt[n] && OddQ[#] &], {n, 1, 90}]
    nmax = 90; CoefficientList[Series[Sum[x^(2 k (2 k - 1))/(1 - x^(2 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
  • PARI
    A333805(n) = sumdiv(n,d,(d%2)&&((d*d)Antti Karttunen, Nov 02 2022

Formula

G.f.: Sum_{k>=1} x^(2*k*(2*k - 1)) / (1 - x^(2*k - 1)).

Extensions

Data section extended up to a(105) by Antti Karttunen, Nov 02 2022

A160786 The number of odd partitions of consecutive odd integers.

Original entry on oeis.org

1, 2, 4, 8, 16, 29, 52, 90, 151, 248, 400, 632, 985, 1512, 2291, 3431, 5084, 7456, 10836, 15613, 22316, 31659, 44601, 62416, 86809, 120025, 165028, 225710, 307161, 416006, 560864, 752877, 1006426, 1340012, 1777365, 2348821, 3093095, 4059416, 5310255, 6924691
Offset: 0

Views

Author

Utpal Sarkar (doetoe(AT)gmail.com), May 26 2009

Keywords

Comments

It seems that these are partitions of odd length and sum, ranked by A340931. The parts do not have to be odd. - Gus Wiseman, Apr 06 2021

Examples

			From _Gus Wiseman_, Apr 06 2021: (Start)
The a(0) = 1 through a(4) = 16 partitions:
  (1)  (3)    (5)      (7)        (9)
       (111)  (221)    (322)      (333)
              (311)    (331)      (432)
              (11111)  (421)      (441)
                       (511)      (522)
                       (22111)    (531)
                       (31111)    (621)
                       (1111111)  (711)
                                  (22221)
                                  (32211)
                                  (33111)
                                  (42111)
                                  (51111)
                                  (2211111)
                                  (3111111)
                                  (111111111)
(End)
		

Crossrefs

Partitions with all odd parts are counted by A000009 and ranked by A066208.
This is a bisection of A027193 (odd-length partitions), which is ranked by A026424.
The case of all odd parts is counted by A078408 and ranked by A300272.
The even version is A236913, ranked by A340784.
A multiplicative version is A340102.
These partitions are ranked by A340931.
A047993 counts balanced partitions, ranked by A106529.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A236914 counts partition of type OO, ranked by A341448.
A340385 counts partitions with odd length and maximum, ranked by A340386.

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+1$2)[2]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 16 2014
  • Mathematica
    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+1, 2*n+1][[2]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Jul 01 2015, after Alois P. Heinz *)
    (* Slow but easy to read *)
    a[n_] := Length@IntegerPartitions[2 n + 1, {1, 2 n + 1, 2}]
    a /@ Range[0, 25]
    (* Leo C. Stein, Nov 11 2020 *)
    (* Faster, don't build the partitions themselves *)
    (* Number of partitions of n into exactly k parts *)
    P[0, 0] = 1;
    P[n_, k_] := 0 /; ((k <= 0) || (n <= 0))
    P[n_, k_] := P[n, k] = P[n - k, k] + P[n - 1, k - 1]
    a[n_] := Sum[P[2 n + 1, k], {k, 1, 2 n + 1, 2}]
    a /@ Range[0, 40]
    (* Leo C. Stein, Nov 11 2020 *)
  • Python
    # Could be memoized for speedup
    def numoddpart(n, m=1):
        """The number of partitions of n into an odd number of parts of size at least m"""
        if n < m:
            return 0
        elif n == m:
            return 1
        else:
            # 1 (namely n = n) and all partitions of the form
            # k + even partitions that start with >= k
            return 1 + sum([numevenpart(n - k,  k) for k in range(m, n//3 + 1)])
    def numevenpart(n, m=1):
        """The number of partitions of n into an even number of parts of size at least m"""
        if n < 2*m:
            return 0
        elif n == 2*m:
            return 1
        else:
            return sum([numoddpart(n - k,  k) for k in range(m,  n//2 + 1)])
    [numoddpart(n) for n in range(1, 70, 2)]
    
  • Python
    # dict to memoize
    ps = {(0,0): 1}
    def p(n, k):
        """Number of partitions of n into exactly k parts"""
        if (n,k) in ps: return ps[(n,k)]
        if (n<=0) or (k<=0): return 0
        ps[(n,k)] = p(n-k,k) + p(n-1,k-1)
        return ps[(n,k)]
    def a(n): return sum([p(2*n+1, k) for k in range(1,2*n+3,2)])
    [a(n) for n in range(0,41)]
    # Leo C. Stein, Nov 11 2020

Formula

a(n) = A027193(2n+1).

A277579 Number of partitions of n for which the number of even parts is equal to the positive alternating sum of the parts.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 3, 3, 4, 6, 7, 9, 13, 15, 19, 25, 31, 38, 48, 59, 74, 90, 111, 136, 166, 201, 246, 297, 357, 431, 522, 621, 745, 892, 1063, 1263, 1503, 1780, 2109, 2491, 2941, 3463, 4077, 4783, 5616, 6576, 7689, 8981, 10486, 12207, 14209, 16516, 19178, 22231
Offset: 0

Views

Author

Emeric Deutsch and Alois P. Heinz, Oct 20 2016

Keywords

Comments

In the first Maple program (improvable) AS gives the positive alternating sum of a finite sequence s, EP gives the number of even terms of a finite sequence of positive integers.
For the specified value of n, the second Maple program lists the partitions of n counted by a(n).
Also the number of integer partitions of n with as many even parts as odd parts in the conjugate partition. - Gus Wiseman, Jul 26 2021

Examples

			a(9) = 6: [2,1,1,1,1,1,1,1], [3,2,1,1,1,1], [3,3,2,1], [4,2,2,1], [4,3,1,1], [5,4].
a(10) = 7: [1,1,1,1,1,1,1,1,1,1], [3,2,2,1,1,1], [3,3,1,1,1,1], [4,2,1,1,1,1], [4,3,2,1], [5,5], [6,4].
a(11) = 9: [2,1,1,1,1,1,1,1,1,1], [3,2,1,1,1,1,1,1], [3,3,2,1,1,1], [3,3,3,2], [4,2,2,1,1,1], [4,3,1,1,1,1], [5,2,2,2], [5,4,1,1], [6,5].
		

Crossrefs

The sign-sensitive version is A035457 (aerated version of A000009).
Comparing odd parts to odd conjugate parts gives A277103.
Comparing product of parts to product of conjugate parts gives A325039.
Comparing the rev-alt sum to that of the conjugate gives A345196.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A120452 counts partitions of 2n with rev-alt sum 2 (negative: A344741).
A124754 gives alternating sums of standard compositions (reverse: A344618).
A316524 is the alternating sum of the prime indices of n (reverse: A344616).
A344610 counts partitions by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.

Programs

  • Maple
    with(combinat): AS := proc (s) options operator, arrow: abs(add((-1)^(i-1)*s[i], i = 1 .. nops(s))) end proc: EP := proc (s) local ct, j: ct := 0: for j to nops(s) do if `mod`(s[j], 2) = 0 then ct := ct+1 else  end if end do: ct end proc: a := proc (n) local P, c, k: P := partition(n): c := 0: for k to nops(P) do if AS(P[k]) = EP(P[k]) then c := c+1 else  end if end do: c end proc: seq(a(n), n = 0 .. 30);
    n := 8: with(combinat): AS := proc (s) options operator, arrow: abs(add((-1)^(i-1)*s[i], i = 1 .. nops(s))) end proc: EP := proc (s) local ct, j: ct := 0: for j to nops(s) do if `mod`(s[j], 2) = 0 then ct := ct+1 else  end if end do: ct end proc: P := partition(n): C := {}: for k to nops(P) do if AS(P[k]) = EP(P[k]) then C := `union`(C, {P[k]}) else  end if end do: C;
    # alternative Maple program:
    b:= proc(n, i, s, t) option remember; `if`(n=0,
          `if`(s=0, 1, 0), `if`(i<1, 0, b(n, i-1, s, t)+
          `if`(i>n, 0, b(n-i, i, s+t*i-irem(i+1, 2), -t))))
        end:
    a:= n-> b(n$2, 0, 1):
    seq(a(n), n=0..60);
  • Mathematica
    b[n_, i_, s_, t_] := b[n, i, s, t] = If[n == 0, If[s == 0, 1, 0], If[i<1, 0, b[n, i-1, s, t] + If[i>n, 0, b[n-i, i, s + t*i - Mod[i+1, 2], -t]]]]; a[n_] := b[n, n, 0, 1]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Dec 21 2016, translated from Maple *)
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]]; Table[Length[Select[IntegerPartitions[n],Count[#,?EvenQ]==Count[conj[#],?OddQ]&]],{n,0,15}] (* Gus Wiseman, Jul 26 2021 *)
  • Sage
    def a(n):
        AS = lambda s: abs(sum((-1)^i*t for i,t in enumerate(s)))
        EP = lambda s: sum((t+1)%2 for t in s)
        return sum(AS(p) == EP(p) for p in Partitions(n))
    print([a(n) for n in (0..30)]) # Peter Luschny, Oct 21 2016

A277103 Number of partitions of n for which the number of odd parts is equal to the positive alternating sum of the parts.

Original entry on oeis.org

1, 1, 0, 1, 3, 3, 1, 3, 10, 10, 4, 10, 27, 27, 13, 28, 69, 69, 37, 72, 161, 162, 96, 171, 361, 364, 230, 388, 768, 777, 522, 836, 1581, 1605, 1128, 1739, 3145, 3203, 2345, 3495, 6094, 6225, 4712, 6831, 11511, 11794, 9198, 13010, 21293, 21875, 17496, 24239
Offset: 0

Views

Author

Emeric Deutsch, Oct 18 2016

Keywords

Comments

It follows by conjugation that the partition statistics "alternating sum" and "number of odd parts" are equidistributed. Consequently, the self-conjugate partitions satisfy the required condition.
In the first Maple program (improvable) AS gives the positive alternating sum of a finite sequence s, OP gives the number of odd terms of a finite sequence of positive integers.
For the specified value of n, the second Maple program lists the partitions of n counted by a(n).
Number of integer partitions of n with the same number of odd parts as their conjugate. - Gus Wiseman, Jun 27 2021

Examples

			a(3) = 1 because we have [2,1]. The partitions [3] and [1,1,1] do not qualify.
a(4) = 3 because we have [3,1], [2,2], and [2,1,1]. The partitions [4] and [1,1,1,1] do not qualify.
		

Crossrefs

Comparing even parts to odd conjugate parts gives A277579.
Comparing product of parts to product of conjugate parts gives A325039.
The reverse version is A345196.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A120452 counts partitions of 2n with rev-alt sum 2 (negative: A344741).
A124754 gives alternating sums of standard compositions (reverse: A344618).
A316524 is the alternating sum of the prime indices of n (reverse: A344616).
A344610 counts partitions by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.

Programs

  • Maple
    with(combinat): AS := proc (s) options operator, arrow: abs(add((-1)^(i-1)*s[i], i = 1 .. nops(s))) end proc: OP := proc (s) local ct, j: ct := 0: for j to nops(s) do if `mod`(s[j], 2) = 1 then ct := ct+1 else  end if end do: ct end proc: a := proc (n) local P, c, k: P := partition(n): c := 0: for k to nops(P) do if AS(P[k]) = OP(P[k]) then c := c+1 else end if end do: c end proc: seq(a(n), n = 0 .. 50);
    n := 8: with(combinat): AS := proc (s) options operator, arrow: abs(add((-1)^(i-1)*s[i], i = 1 .. nops(s))) end proc: OP := proc (s) local ct, j: ct := 0: for j to nops(s) do if `mod`(s[j], 2) = 1 then ct := ct+1 else  end if end do: ct end proc: P := partition(n): C := {}: for k to nops(P) do if AS(P[k]) = OP(P[k]) then C := `union`(C, {P[k]}) else  end if end do: C;
    # alternative Maple program:
    b:= proc(n, i, s, t) option remember; `if`(n=0,
          `if`(s=0, 1, 0), `if`(i<1, 0, b(n, i-1, s, t)+
          `if`(i>n, 0, b(n-i, i, s+t*i-irem(i, 2), -t))))
        end:
    a:= n-> b(n$2, 0, 1):
    seq(a(n), n=0..60);  # Alois P. Heinz, Oct 19 2016
  • Mathematica
    b[n_, i_, s_, t_] := b[n, i, s, t] = If[n == 0, If[s == 0, 1, 0], If[i<1, 0, b[n, i-1, s, t] + If[i>n, 0, b[n-i, i, s + t*i - Mod[i, 2], -t]]]]; a[n_] := b[n, n, 0, 1]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Dec 21 2016, after Alois P. Heinz *)
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]]; Table[Length[Select[IntegerPartitions[n],Count[#,?OddQ]==Count[conj[#],?OddQ]&]],{n,0,15}] (* Gus Wiseman, Jun 27 2021 *)

A345917 Numbers k such that the k-th composition in standard order (row k of A066099) has alternating sum > 0.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 9, 11, 14, 16, 17, 18, 19, 21, 22, 23, 26, 28, 29, 31, 32, 33, 34, 35, 37, 38, 39, 42, 44, 45, 47, 52, 56, 57, 59, 62, 64, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 84, 85, 87, 88, 89, 90, 91, 93, 94, 95, 100, 104, 105, 107
Offset: 1

Views

Author

Gus Wiseman, Jul 08 2021

Keywords

Comments

The alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The initial terms and the corresponding compositions:
     1: (1)
     2: (2)
     4: (3)
     5: (2,1)
     7: (1,1,1)
     8: (4)
     9: (3,1)
    11: (2,1,1)
    14: (1,1,2)
    16: (5)
    17: (4,1)
    18: (3,2)
    19: (3,1,1)
    21: (2,2,1)
    22: (2,1,2)
		

Crossrefs

The version for Heinz numbers of partitions is A026424.
These compositions are counted by A027306.
These are the positions of terms > 0 in A124754.
The weak (k >= 0) version is A345913.
The reverse-alternating version is A345918.
The opposite (k < 0) version is A345919.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Select[Range[0,100],ats[stc[#]]>0&]

A345911 Numbers k such that the k-th composition in standard order (row k of A066099) has reverse-alternating sum 1.

Original entry on oeis.org

1, 6, 7, 20, 21, 26, 27, 30, 31, 72, 73, 82, 83, 86, 87, 92, 93, 100, 101, 106, 107, 110, 111, 116, 117, 122, 123, 126, 127, 272, 273, 290, 291, 294, 295, 300, 301, 312, 313, 324, 325, 330, 331, 334, 335, 340, 341, 346, 347, 350, 351, 360, 361, 370, 371, 374
Offset: 1

Views

Author

Gus Wiseman, Jul 01 2021

Keywords

Comments

The reverse-alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence of terms together with the corresponding compositions begins:
     1: (1)
     6: (1,2)
     7: (1,1,1)
    20: (2,3)
    21: (2,2,1)
    26: (1,2,2)
    27: (1,2,1,1)
    30: (1,1,1,2)
    31: (1,1,1,1,1)
    72: (3,4)
    73: (3,3,1)
    82: (2,3,2)
    83: (2,3,1,1)
    86: (2,2,1,2)
    87: (2,2,1,1,1)
		

Crossrefs

These compositions are counted by A000984 (bisection of A126869).
The version for Heinz numbers of partitions is A001105.
A version using runs of binary digits is A066879.
These are positions of 1's in A344618.
The non-reverse version is A345909.
The opposite (negative 1) version is A345912.
The version for prime indices is A345958.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A011782 counts compositions.
A097805 counts compositions by alternating or reverse-alternating sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344610 counts partitions by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.
A345197 counts compositions by sum, length, and alternating sum.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Select[Range[0,100],sats[stc[#]]==1&]

A064410 Number of partitions of n with zero crank.

Original entry on oeis.org

0, 0, 1, 1, 1, 1, 1, 2, 2, 4, 4, 7, 7, 11, 12, 17, 19, 27, 30, 41, 48, 62, 73, 95, 110, 140, 166, 206, 243, 302, 354, 435, 513, 622, 733, 887, 1039, 1249, 1467, 1750, 2049, 2438, 2847, 3371, 3934, 4634, 5398, 6343, 7367, 8626, 10009, 11677, 13521, 15737, 18184
Offset: 1

Views

Author

Vladeta Jovovic, Sep 29 2001

Keywords

Comments

For a partition p, let l(p) = largest part of p, w(p) = number of 1's in p, m(p) = number of parts of p larger than w(p). The crank of p is given by l(p) if w(p) = 0, otherwise m(p)-w(p).

Examples

			a(10)=4 because there are 4 partitions of 10 with zero crank: 1+1+2+3+3, 1+1+4+4, 1+1+3+5 and 1+9.
From _Gus Wiseman_, Apr 02 2021: (Start)
The a(3) = 1 through a(14) = 11 partitions (A..D = 10..13):
  21  31  41  51  61  71    81    91     A1     B1      C1      D1
                      3311  4311  4411   5411   5511    6511    6611
                                  5311   6311   6411    7411    7511
                                  33211  43211  7311    8311    8411
                                                44211   54211   9311
                                                53211   63211   55211
                                                332211  432211  64211
                                                                73211
                                                                442211
                                                                532211
                                                                3322211
(End)
		

Crossrefs

The version for positive crank is A001522.
Central column of A064391.
The version for nonnegative crank is A064428.
The Heinz numbers of these partitions are A342192.
A003242 counts anti-run compositions.
A224958 counts compositions with alternating parts unequal.
A257989 gives the crank of the partition with Heinz number n.

Programs

  • Mathematica
    nmax = 60; Rest[CoefficientList[Series[x - 1 + Sum[(-1)^k*(x^(k*(k + 1)/2) - x^(k*(k - 1)/2)), {k, 1, nmax}] / Product[1 - x^k, {k, 1, nmax}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, Sep 26 2016 *)
    Flatten[{0, Table[PartitionsP[n] - 2*Sum[(-1)^(j+1)*PartitionsP[n - j*((j+1)/2)], {j, 1, Floor[(Sqrt[8*n + 1] - 1)/2]}], {n, 2, 60}]}] (* Vaclav Kotesovec, Sep 26 2016 *)
    ck[y_]:=With[{w=Count[y,1]},If[w==0,Max@@y,Count[y,_?(#>w&)]-w]];
    Table[Length[Select[IntegerPartitions[n],ck[#]==0&]],{n,0,30}] (* Gus Wiseman, Apr 02 2021 *)
  • Sage
    [[p.crank() for p in Partitions(n)].count(0) for n in (1..20)] # Peter Luschny, Sep 15 2014

Formula

a(n) = A000041(n) - 2*A001522(n). a(n) = A064391(n, 0).
a(n) ~ exp(Pi*sqrt(2*n/3)) * Pi / (3 * 2^(9/2) * n^(3/2)). - Vaclav Kotesovec, May 06 2018
a(n > 1) = A064428(n) - A001522(n), where A001522/A064428 count odd/even-length compositions with alternating parts strictly decreasing. - Gus Wiseman, Apr 02 2021
From Peter Bala, Feb 03 2024: (Start)
For n >= 2, a(n) = A188674(n+1) - A188674(n) (Hopkins and Sellers, Proposition 7).
Equivalently, the g.f. A(x) = (1 - x) * Sum_{n >= 1} x^(n*(n+2)) / Product{k = 1..n} (1 - x^k)^2. (End)

Extensions

More terms from Reiner Martin, Dec 26 2001

A340604 Heinz numbers of integer partitions of odd positive rank.

Original entry on oeis.org

3, 7, 10, 13, 15, 19, 22, 25, 28, 29, 33, 34, 37, 42, 43, 46, 51, 52, 53, 55, 61, 62, 63, 69, 70, 71, 76, 77, 78, 79, 82, 85, 88, 89, 93, 94, 98, 101, 105, 107, 113, 114, 115, 116, 117, 118, 119, 121, 123, 130, 131, 132, 134, 136, 139, 141, 146, 147, 148, 151
Offset: 1

Views

Author

Gus Wiseman, Jan 21 2021

Keywords

Comments

The Dyson rank of a nonempty partition is its maximum part minus its number of parts. The rank of an empty partition is 0.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of partitions with their Heinz numbers begins:
      3: (2)         46: (9,1)       82: (13,1)
      7: (4)         51: (7,2)       85: (7,3)
     10: (3,1)       52: (6,1,1)     88: (5,1,1,1)
     13: (6)         53: (16)        89: (24)
     15: (3,2)       55: (5,3)       93: (11,2)
     19: (8)         61: (18)        94: (15,1)
     22: (5,1)       62: (11,1)      98: (4,4,1)
     25: (3,3)       63: (4,2,2)    101: (26)
     28: (4,1,1)     69: (9,2)      105: (4,3,2)
     29: (10)        70: (4,3,1)    107: (28)
     33: (5,2)       71: (20)       113: (30)
     34: (7,1)       76: (8,1,1)    114: (8,2,1)
     37: (12)        77: (5,4)      115: (9,3)
     42: (4,2,1)     78: (6,2,1)    116: (10,1,1)
     43: (14)        79: (22)       117: (6,2,2)
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
These partitions are counted by A101707.
Allowing negative ranks gives A340692, counted by A340603.
The even version is A340605, counted by A101708.
The not necessarily odd case is A340787, counted by A064173.
A001222 gives number of prime indices.
A061395 gives maximum prime index.
- Rank -
A047993 counts partitions of rank 0 (A106529).
A064173 counts partitions of negative rank (A340788).
A064174 counts partitions of nonnegative rank (A324562).
A064174 (also) counts partitions of nonpositive rank (A324521).
A101198 counts partitions of rank 1 (A325233).
A257541 gives the rank of the partition with Heinz number n.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts (A066208).
A027193 counts partitions of odd length (A026424).
A027193 (also) counts partitions of odd maximum (A244991).
A058695 counts partitions of odd numbers (A300063).
A067659 counts strict partitions of odd length (A030059).
A160786 counts odd-length partitions of odd numbers (A300272).
A339890 counts factorizations of odd length.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.
A340385 counts partitions of odd length and maximum (A340386).

Programs

  • Mathematica
    rk[n_]:=PrimePi[FactorInteger[n][[-1,1]]]-PrimeOmega[n];
    Select[Range[100],OddQ[rk[#]]&&rk[#]>0&]

Formula

A061395(a(n)) - A001222(a(n)) is odd and positive.

A341594 Number of strictly superior odd divisors of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 23 2021

Keywords

Comments

We define a divisor d|n to be strictly superior if d > n/d. Strictly superior divisors are counted by A056924 and listed by A341673.

Examples

			The a(n) divisors for n = 3, 15, 45, 105, 315, 405, 945, 1575, 1890:
  3  5   9   15   21   27   35   45    45
     15  15  21   35   45   45   63    63
         45  35   45   81   63   75    105
             105  63   135  105  105   135
                  105  405  135  175   189
                  315       189  225   315
                            315  315   945
                            945  525
                                 1575
		

Crossrefs

On odd indices, equals A056924 (number of strictly superior divisors).
The inferior version is A069288.
Positions of zeros are A116882.
Positions of nonzero terms are A116883.
The strictly inferior version is A333805.
The version for squarefree instead of odd divisors is A341595.
The version for prime instead of odd divisors is A341642.
The version for prime-power instead of odd divisors is A341644.
The superior version is A341675.
A033676 selects the greatest inferior divisor.
A033677 selects the smallest superior divisor.
A038548 counts superior (or inferior) divisors.
A140271 selects the smallest strictly superior divisor.
A207375 lists central divisors.
A341673 lists strictly superior divisors.
- Strictly Inferior: A060775, A070039, A333806, A341596, A341674.
- Strictly Superior: A048098, A064052, A238535, A341643, A341645, A341646.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n],OddQ[#]&&#>n/#&]],{n,100}]
  • PARI
    A341594(n) = sumdiv(n,d,(d%2)*(d>(n/d))); \\ Antti Karttunen, Feb 23 2021

A298118 Number of unlabeled rooted trees with n nodes in which all positive outdegrees are odd.

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 11, 21, 40, 80, 159, 322, 657, 1356, 2816, 5896, 12407, 26267, 55861, 119331, 255878, 550665, 1188786, 2574006, 5588177, 12162141, 26529873, 57993624, 127020653, 278716336, 612617523, 1348680531, 2973564157, 6565313455, 14514675376
Offset: 1

Views

Author

Gus Wiseman, Jan 12 2018

Keywords

Examples

			The a(6) = 6 trees: (((((o))))), (((ooo))), ((oo(o))), (oo((o))), (o(o)(o)), (ooooo).
		

Crossrefs

Programs

  • Mathematica
    orut[n_]:=orut[n]=If[n===1,{{}},Join@@Function[c,Union[Sort/@Tuples[orut/@c]]]/@Select[IntegerPartitions[n-1],OddQ[Length[#]]&]];
    Table[Length[orut[n]],{n,15}]

Formula

a(n) ~ c * d^n / n^(3/2), where d = 2.30984417428419893876754252289588812511559... and c = 0.5598122522173731208680575003383895445787... - Vaclav Kotesovec, Jun 04 2019

Extensions

a(24)-a(35) from Alois P. Heinz, Jan 12 2018
Previous Showing 31-40 of 206 results. Next