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 21-30 of 189 results. Next

A063886 Number of n-step walks on a line starting from the origin but not returning to it.

Original entry on oeis.org

1, 2, 2, 4, 6, 12, 20, 40, 70, 140, 252, 504, 924, 1848, 3432, 6864, 12870, 25740, 48620, 97240, 184756, 369512, 705432, 1410864, 2704156, 5408312, 10400600, 20801200, 40116600, 80233200, 155117520, 310235040, 601080390, 1202160780, 2333606220, 4667212440
Offset: 0

Views

Author

Henry Bottomley, Aug 28 2001

Keywords

Comments

A Chebyshev transform of A007877(n+1). The g.f. is transformed to (1+x)/((1-x)(1+x^2)) under the mapping G(x)->(1/(1+x^2))G(1/(1+x^2)). - Paul Barry, Oct 12 2004
a(n-1) = 2*C(n-2, floor((n-2)/2)) is also the number of bit strings of length n in which the number of 00 substrings is equal to the number of 11 substrings. For example, when n = 4 we have 4 such bit strings: 0011, 0101, 1010, and 1100. - Angel Plaza, Apr 23 2009
Hankel transform is A120617. - Paul Barry, Aug 10 2009
The Hankel transform of a(n) is (-2)^C(n+1,2). The Hankel transform of (-1)^C(n+1,2)*a(n) is (-1)^C(n+1,2)*A164584(n). - Paul Barry, Aug 17 2009
For n > 1, a(n) is also the number of n-step walks starting from the origin and returning to it exactly once. - Geoffrey Critzer, Jan 24 2010
-a(n) is the Z-sequence for the Riordan array A130777. (See the W. Lang link under A006232 for A- and Z-sequences for Riordan matrices). - Wolfdieter Lang, Jul 12 2011
Number of subsets of {1,...,n} in which the even elements appear as often at even positions as at odd positions. - Gus Wiseman, Mar 17 2018

Examples

			a(4) = 6 because there are six length four walks that do not return to the origin: {-1, -2, -3, -4}, {-1, -2, -3, -2}, {-1, -2, -1, -2}, {1, 2, 1, 2}, {1, 2, 3, 2}, {1, 2, 3, 4}. There are also six such walks that return exactly one time: {-1, -2, -1, 0}, {-1, 0, -1, -2}, {-1, 0, 1, 2}, {1, 0, -1, -2}, {1, 0, 1, 2}, {1, 2, 1, 0}. - _Geoffrey Critzer_, Jan 24 2010
The a(5) = 12 subsets in which the even elements appear as often at even positions as at odd positions: {}, {1}, {3}, {5}, {1,3}, {1,5}, {2,4}, {3,5}, {1,2,4}, {1,3,5}, {2,4,5}, {1,2,4,5}. - _Gus Wiseman_, Mar 17 2018
		

Crossrefs

Programs

  • Magma
    [1] cat [2*Binomial(n-1, Floor((n-1)/2)): n in [1..40]]; // G. C. Greubel, Jun 07 2023
    
  • Maple
    seq(seq(binomial(2*j,j)*i, i=1..2),j=0..16); # Zerinvary Lajos, Apr 28 2007
    # second Maple program:
    a:= proc(n) option remember; `if`(n<2, n+1,
           4*a(n-2) +2*(a(n-1) -4*a(n-2))/n)
        end:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 10 2014
    # third program:
    A063886 := series(BesselI(0, 2*x)*(1 + x*2 + x*Pi*StruveL(1, 2*x)) - Pi*x*BesselI(1, 2*x)*StruveL(0, 2*x), x = 0, 34): seq(n!*coeff(A063886, x, n), n = 0 .. 33); # Mélika Tebni, Jun 17 2024
  • Mathematica
    Table[Length[Select[Map[Accumulate, Strings[{-1, 1}, n]], Count[ #, 0] == 0 &]], {n, 0, 20}] (* Geoffrey Critzer, Jan 24 2010 *)
    CoefficientList[Series[Sqrt[(1+2x)/(1-2x)],{x,0,40}],x] (* Harvey P. Dale, Apr 28 2016 *)
  • PARI
    a(n)=(n==0)+2*binomial(n-1,(n-1)\2)
    
  • PARI
    a(n) = 2^n*prod(k=0,n-1,(k/n+1/n)^((-1)^k)); \\ Michel Marcus, Dec 03 2013
    
  • Python
    from math import ceil
    from sympy import binomial
    def a(n):
        if n==0: return 1
        return 2*binomial(n-1,(n-1)//2)
    print([a(n) for n in range(18)])
    # David Nacin, Feb 29 2012
    
  • SageMath
    [2*binomial(n-1, (n-1)//2) + int(n==0) for n in range(41)] # G. C. Greubel, Jun 07 2023

Formula

G.f.: sqrt((1+2*x)/(1-2*x)).
a(n+1) = 2*C(n, floor(n/2)) = 2*A001405(n); a(2n) = C(2n, n) = A000984(n) = 4*a(2n-2)-|A002420(n)| = 4*a(2n-2)-2*A000108(n-1) = 2*A001700(n-1); a(2n+1) = 2*a(2n) = A028329(n).
2*a(n) = A047073(n+1).
a(n) = Sum_{k=0..n} abs(A106180(n,k)). - Philippe Deléham, Oct 06 2006
a(n) = Sum_{k=0..n} (k+1)binomial(n, (n-k)/2) ( 1-cos((k+1)*Pi/2) (1+(-1)^(n-k))/(n+k+2) ). - Paul Barry, Oct 12 2004
G.f.: 1/(1-2*x/(1+x/(1+x/(1-x/(1-x/(1+x/(1+x/(1-x/(1-x/(1+ ... (continued fraction). - Paul Barry, Aug 10 2009
G.f.: 1 + 2*x/(G(0)-x+x^2) where G(k)= 1 - 2*x^2 - x^4/G(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Aug 10 2012
D-finite with recurrence: n*a(n) = 2*a(n-1) + 4*(n-2)*a(n-2). - R. J. Mathar, Dec 03 2012
From Sergei N. Gladkovskii, Jul 26 2013: (Start)
G.f.: 1/G(0), where G(k) = 1 - 2*x/(1 + 2*x/(1 + 1/G(k+1) )); (continued fraction).
G.f.: G(0), where G(k) = 1 + 2*x/(1 - 2*x/(1 + 1/G(k+1) )); (continued fraction).
G.f.: W(0)/2*(1+2*x), where W(k) = 1 + 1/(1 - 2*x/(2*x + (k+1)/(x*(2*k+1))/W(k+1) )), abs(x) < 1/2; (continued fraction). (End)
a(n) = 2^n*Product_{k=0..n-1} (k/n + 1/n)^((-1)^k). - Peter Luschny, Dec 02 2013
G.f.: G(0), where G(k) = 1 + 2*x*(4*k+1)/((2*k+1)*(1+2*x) - (2*k+1)*(4*k+3)*x*(1+2*x)/((4*k+3)*x + (k+1)*(1+2*x)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 19 2014
From Peter Bala, Mar 29 2024: (Start)
a(n) = 2^n * Sum_{k = 0..n} (-1)^(n+k)*binomial(1/2, k)*binomial(- 1/2, n-k) = 2^n * A000246(n)/n!.
a(n) = (1/2^n) * binomial(2*n, n) * hypergeom([-1/2, -n], [1/2 - n], -1). (End)
E.g.f.: BesselI(0, 2*x)*(1 + x*(2 + Pi)*StruveL(1, 2*x)) - Pi*x*BesselI(1, 2*x)*StruveL(0, 2*x). - Stefano Spezia, May 11 2024
a(n) = A089849(n) + A138364(n). - Mélika Tebni, Jun 17 2024
From Amiram Eldar, Aug 15 2025: (Start)
Sum_{n>=0} 1/a(n) = Pi/(3*sqrt(3)) + 2.
Sum_{n>=0} (-1)^n/a(n) = 2/3 + Pi/(9*sqrt(3)). (End)

A237258 Number of strict partitions of 2n that include a partition of n.

Original entry on oeis.org

1, 0, 0, 1, 1, 3, 4, 7, 9, 16, 21, 32, 43, 63, 84, 122, 158, 220, 293, 393, 511, 685, 881, 1156, 1485, 1925, 2445, 3147, 3952, 5019, 6323, 7924, 9862, 12336, 15259, 18900, 23294, 28646, 35091, 42985, 52341, 63694, 77336, 93588, 112973, 136367, 163874, 196638
Offset: 0

Views

Author

Clark Kimberling, Feb 05 2014

Keywords

Comments

A strict partition is a partition into distinct parts.

Examples

			a(5) counts these partitions of 10: [5,4,1], [5,3,2], [4,3,2,1].
		

Crossrefs

The non-strict version is A002219, ranked by A357976.
These partitions are ranked by A357854.
A000712 counts distinct submultisets of partitions, strict A032302.
A304792 counts subset-sums of partitions, positive A276024, strict A284640.

Programs

  • Mathematica
    z = 24; Table[theTotals = Map[{#, Map[Total, Subsets[#]]} &,  Select[IntegerPartitions[2 nn], # == DeleteDuplicates[#] &]]; Length[Map[#[[1]] &, Select[theTotals, Length[Position[#[[2]], nn]] >= 1 &]]], {nn, z}] (* Peter J. C. Moses, Feb 04 2014 *)

Formula

a(n) = A237194(2n,n).

Extensions

a(31)-a(47) from Alois P. Heinz, Feb 07 2014

A239830 Triangular array: T(n,k) = number of partitions of 2n that have alternating sum 2k, with T(0,0) = 1 for convenience.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 3, 5, 2, 1, 5, 9, 5, 2, 1, 7, 17, 10, 5, 2, 1, 11, 28, 20, 10, 5, 2, 1, 15, 47, 35, 20, 10, 5, 2, 1, 22, 73, 62, 36, 20, 10, 5, 2, 1, 30, 114, 102, 65, 36, 20, 10, 5, 2, 1, 42, 170, 167, 109, 65, 36, 20, 10, 5, 2, 1, 56, 253, 262, 182, 110
Offset: 0

Views

Author

Clark Kimberling, Mar 28 2014

Keywords

Comments

Suppose that p, with parts x(1) >= x(2) >= ... >= x(k), is a partition of n. Define AS(p), the alternating sum of p, by x(1) - x(2) + x(3) - ... + ((-1)^(k-1))*x(k); note that AS(p) has the same parity as n. Column 1 is given by T(n,1) = A000041(n) for n >= 0, which is the number of partitions of 2n having AS(p) = 0, for n >= 1. Columns 2 and 3 are essentially A000567 and A000710, and the limiting column (after deleting initial 0's), A000712. The sum of numbers in row n is A000041(2n). The corresponding array for partitions into distinct parts is given by A152146 (defined as the number of unrestricted partitions of 2n into 2k even parts).

Examples

			First nine rows:
1
1 ... 1
2 ... 2 ... 1
3 ... 5 ... 2 ... 1
5 ... 9 ... 5 ... 2 ... 1
7 ... 17 .. 10 .. 5 ... 2 ... 1
11 .. 28 .. 20 .. 10 .. 5 ... 2 ... 1
15 .. 47 .. 35 .. 20 .. 10 .. 5 ... 2 ... 1
22 .. 73 .. 62 .. 36 .. 20 .. 10 .. 5 ... 2 ... 1
The partitions of 6 are 6, 51, 42, 411, 33, 321, 3111, 222, 2211, 21111, 111111, with respective alternating sums 6, 4, 2, 4, 0, 2, 2, 2, 0, 2, 0, so that row 3 (counting the top row as row 0) of the array is 3 .. 5 .. 2 .. 1.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n=0, 1, `if`(i<1, 0,
          expand(b(n, i-1, t)+`if`(i>n, 0, b(n-i, i, -t)*x^((t*i)/2)))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(2*n$2, 1)):
    seq(T(n), n=0..14);  # Alois P. Heinz, Mar 30 2014
  • Mathematica
    z = 16; s[w_] := s[w] = Total[Take[#, ;; ;; 2]] - Total[Take[Rest[#], ;; ;; 2]] &[w]; c[n_] := c[n] = Table[s[IntegerPartitions[n][[k]]], {k, 1, PartitionsP[n]}]; t[n_, k_] := Count[c[2 n], 2 k]; t[0, 0] = 1; u = Table[t[n, k], {n, 0, z}, {k, 0, n}]
    TableForm[u]  (* A239830, array *)
    Flatten[u]    (* A239830, sequence *)
    (* Peter J. C. Moses, Mar 21 2014 *)

A000716 Number of partitions of n into parts of 3 kinds.

Original entry on oeis.org

1, 3, 9, 22, 51, 108, 221, 429, 810, 1479, 2640, 4599, 7868, 13209, 21843, 35581, 57222, 90882, 142769, 221910, 341649, 521196, 788460, 1183221, 1762462, 2606604, 3829437, 5590110, 8111346, 11701998, 16790136, 23964594, 34034391, 48104069, 67679109, 94800537, 132230021, 183686994, 254170332
Offset: 0

Views

Author

Keywords

References

  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 122.
  • Moreno, Carlos J., Partitions, congruences and Kac-Moody Lie algebras. Preprint, 37pp., no date. See Table I.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000712.
Column 3 of A144064.

Programs

  • Maple
    with(numtheory): a:= proc(n) option remember; `if`(n=0, 1, add(add(d*3, d=divisors(j)) *a(n-j), j=1..n)/n) end: seq(a(n), n=0..40); # Alois P. Heinz, May 20 2013
  • Mathematica
    a[0] = 1; a[n_] := a[n] = 1/n*Sum[3*a[k]*DivisorSigma[1, n-k], {k, 0, n-1}]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 03 2014, after Joerg Arndt *)
    (1/QPochhammer[q]^3 + O[q]^40)[[3]] (* Vladimir Reshetnikov, Nov 21 2016 *)
  • PARI
    Vec(1/eta('x+O('x^66))^3) \\ Joerg Arndt, Apr 28 2013
    
  • Python
    from functools import lru_cache
    from sympy import divisor_sigma
    @lru_cache(maxsize=None)
    def A000716(n): return sum(A000716(k)*divisor_sigma(n-k) for k in range(n))*3//n if n else 1 # Chai Wah Wu, Sep 25 2023

Formula

G.f.: Product_{m>=1} 1/(1-x^m)^3.
EULER transform of 3, 3, 3, 3, 3, 3, 3, 3, ...
a(0)=1, a(n) = 1/n*Sum_{k=0..n-1} 3*a(k)*sigma_1(n-k). - Joerg Arndt, Feb 05 2011
a(n) ~ exp(Pi * sqrt(2*n)) / (8 * sqrt(2) * n^(3/2)) * (1 - (3/Pi + Pi/8) / sqrt(2*n)). - Vaclav Kotesovec, Feb 28 2015, extended Jan 16 2017
G.f.: exp(3*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 06 2018

Extensions

Extended with formula from Christian G. Bower, Apr 15 1998

A008485 Coefficient of x^n in Product_{k>=1} 1/(1-x^k)^n.

Original entry on oeis.org

1, 1, 5, 22, 105, 506, 2492, 12405, 62337, 315445, 1605340, 8207563, 42124380, 216903064, 1119974875, 5796944357, 30068145905, 156250892610, 813310723925, 4239676354650, 22130265931900, 115654632452535, 605081974091875, 3168828466966388, 16610409114771900
Offset: 0

Views

Author

T. Forbes (anthony.d.forbes(AT)googlemail.com)

Keywords

Comments

Number of partitions of n into parts of n kinds. - Vladeta Jovovic, Sep 08 2002
Main diagonal of A144064. - Omar E. Pol, Jun 27 2012
From Peter Bala, Apr 18 2023: (Start)
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and all positive integers n and k.
Conjecture: the supercongruence a(p) == p + 1 (mod p^2) holds for all primes p >= 3. Cf. A270913. (End)

Crossrefs

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; `if`(n=0, 1, add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n) end end: a:= n-> etr(j->n)(n): seq(a(n), n=0..30); # Alois P. Heinz, Sep 09 2008
  • Mathematica
    a[n_] := SeriesCoefficient[ Product[1/(1-x^k)^n, {k, 1, n}], {x, 0, n}]; a[1] = 1; Table[a[n], {n, 1, 24}] (* Jean-François Alcover, Feb 24 2015 *)
    Table[SeriesCoefficient[1/QPochhammer[x, x]^n, {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Mar 25 2016 *)
    Table[SeriesCoefficient[Exp[n*Sum[x^j/(j*(1-x^j)), {j, 1, n}]], {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, May 19 2018 *)
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-x^k +x*O(x^n))^n),n)}
    
  • PARI
    {a(n)=n*polcoeff(log(1/x*serreverse(x*eta(x+x*O(x^n)))), n)} /* Paul D. Hanna, Apr 05 2012 */

Formula

a(n) = Sum_{pi} Product_{i=1..n} binomial(k_i+n-1, k_i) where pi runs through all nonnegative solutions of k_1+2*k_2+...+n*k_n=n. a(n) = b(n, n) where b(n, m)= m/n*Sum_{i=1..n} sigma(i)*b(n-i, m) is recurrence for number of partitions of n into parts of m kinds. - Vladeta Jovovic, Sep 08 2002
Equals the logarithmic derivative of A109085, the g.f. of which is (1/x)*Series_Reversion(x*eta(x)). - Paul D. Hanna, Apr 05 2012
Let G(x) = exp( Sum_{n>=1} a(n)*x^n/n ), then G(x) = 1/Product_{n>=1} (1-x^n*G(x)^n) is the g.f. of A109085. - Paul D. Hanna, Apr 05 2012
a(n) ~ c * d^n / sqrt(n), where d = A270915 = 5.352701333486642687772415814165..., c = A327279 = 0.26801521271073331568695383828... . - Vaclav Kotesovec, Sep 10 2014

Extensions

a(0)=1 prepended by Alois P. Heinz, Mar 30 2015

A262977 a(n) = binomial(4*n-1,n).

Original entry on oeis.org

1, 3, 21, 165, 1365, 11628, 100947, 888030, 7888725, 70607460, 635745396, 5752004349, 52251400851, 476260169700, 4353548972850, 39895566894540, 366395202809685, 3371363686069236, 31074067324187580, 286845713747883300, 2651487106659130740, 24539426037817994160
Offset: 0

Views

Author

Vladimir Kruchinin, Oct 06 2015

Keywords

Comments

From Gus Wiseman, Sep 28 2022: (Start)
Also the number of integer compositions of 4n with alternating sum 2n, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. These compositions are ranked by A348614. The a(12) = 21 compositions are:
(6,2) (1,2,5) (1,1,5,1) (1,1,1,1,4)
(2,2,4) (2,1,4,1) (1,1,2,1,3)
(3,2,3) (3,1,3,1) (1,1,3,1,2)
(4,2,2) (4,1,2,1) (1,1,4,1,1)
(5,2,1) (5,1,1,1) (2,1,1,1,3)
(2,1,2,1,2)
(2,1,3,1,1)
(3,1,1,1,2)
(3,1,2,1,1)
(4,1,1,1,1)
The following pertain to this interpretation:
- The case of partitions is A000712, reverse A006330.
- Allowing any alternating sum gives A013777 (compositions of 4n).
- A011782 counts compositions of n.
- A034871 counts compositions of 2n with alternating sum 2k.
- A097805 counts compositions by alternating (or reverse-alternating) sum.
- A103919 counts partitions by sum and alternating sum (reverse: A344612).
- A345197 counts compositions by length and alternating sum.
(End)

Crossrefs

Programs

  • Magma
    [Binomial(4*n-1,n): n in [0..20]]; // Vincenzo Librandi, Oct 06 2015
    
  • Mathematica
    Table[Binomial[4 n - 1, n], {n, 0, 40}] (* Vincenzo Librandi, Oct 06 2015 *)
  • Maxima
    B(x):=sum(binomial(4*n-1,n-1)*3/(4*n-1)*x^n,n,1,30);
    taylor(x*diff(B(x),x,1)/B(x),x,0,20);
    
  • PARI
    a(n) = binomial(4*n-1,n); \\ Michel Marcus, Oct 06 2015

Formula

G.f.: A(x)=x*B'(x)/B(x), where B(x) if g.f. of A006632.
a(n) = Sum_{k=0..n}(binomial(n-1,n-k)*binomial(3*n,k)).
a(n) = 3*A224274(n), for n > 0. - Michel Marcus, Oct 12 2015
From Peter Bala, Nov 04 2015: (Start)
The o.g.f. equals f(x)/g(x), where f(x) is the o.g.f. for A005810 and g(x) is the o.g.f. for A002293. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(4*n + k,n). Cf. A005810 (k = 0), A052203 (k = 1), A257633 (k = 2), A224274 (k = 3) and A004331 (k = 4). (End)
a(n) = [x^n] 1/(1 - x)^(3*n). - Ilya Gutkovskiy, Oct 03 2017
a(n) = A071919(3n-1,n+1) = A097805(4n,n+1). - Gus Wiseman, Sep 28 2022
From Peter Bala, Feb 14 2024: (Start)
a(n) = (-1)^n * binomial(-3*n, n).
a(n) = hypergeom([1 - 3*n, -n], [1], 1).
The g.f. A(x) satisfies A(x/(1 + x)^4) = 1/(1 - 3*x). (End)
a(n) = Sum_{k = 0..n} binomial(2*n+k-1, k)*binomial(2*n-k-1, n-k). - Peter Bala, Sep 16 2024
G.f.: 1/(4-3*g) where g = 1+x*g^4 is the g.f. of A002293. - Seiichi Manyama, Aug 17 2025

A242618 Number T(n,k) of partitions of n, where k is the difference between the number of odd parts and the number of even parts, both counted without multiplicity; triangle T(n,k), n>=0, read by rows.

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 2, 2, 1, 1, 1, 4, 2, 1, 1, 2, 3, 3, 2, 1, 8, 3, 3, 2, 4, 6, 5, 5, 4, 13, 8, 4, 1, 5, 5, 11, 13, 7, 1, 11, 20, 14, 9, 2, 1, 6, 13, 17, 26, 11, 3, 1, 22, 31, 27, 15, 5, 2, 12, 18, 34, 44, 18, 7, 4, 40, 47, 51, 23, 11, 5, 16, 36, 56, 72, 34, 11, 1
Offset: 0

Views

Author

Alois P. Heinz, May 19 2014

Keywords

Comments

T(n,0) = A241638(n).
Sum_{k<0} T(n,k) = A241640(n).
Sum_{k<=0} T(n,k) = A241639(n).
Sum_{k>=0} T(n,k) = A241637(n).
Sum_{k>0} T(n,k) = A241636(n).
T(n^2,n) = T(n^2+n,-n) = 1.
T(n^2+n,n) = Sum_{k} T(n,k) = A000041(n).
T(n^2+3*n,-n) = A000712(n).

Examples

			Triangle T(n,k) begins:
: n\k : -3  -2  -1   0   1   2   3 ...
+-----+---------------------------
:  0  :              1;
:  1  :                  1;
:  2  :          1,  0,  1;
:  3  :              1,  2;
:  4  :          2,  1,  1,  1;
:  5  :              4,  2,  1;
:  6  :      1,  2,  3,  3,  2;
:  7  :          1,  8,  3,  3;
:  8  :      2,  4,  6,  5,  5;
:  9  :          4, 13,  8,  4,  1;
: 10  :      5,  5, 11, 13,  7,  1;
: 11  :         11, 20, 14,  9,  2;
: 12  :  1,  6, 13, 17, 26, 11,  3;
: 13  :      1, 22, 31, 27, 15,  5;
: 14  :  2, 12, 18, 34, 44, 18,  7;
		

Crossrefs

Row sums give A000041.
Cf. A240009 (parts counted with multiplicity), A240021 (distinct parts), A242626 (compositions counted without multiplicity).

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          expand(b(n, i-1)+add(b(n-i*j, i-1)*x^(2*irem(i, 2)-1), j=1..n/i))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=ldegree(p)..degree(p)))(b(n$2)):
    seq(T(n), n=0..20);
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, Expand[b[n, i - 1] + Sum[b[n - i*j, i - 1]*x^(2*Mod[i, 2] - 1), {j, 1, n/i}]]]]; T[n_] := Function[p, Table[Coefficient[p, x, i], {i, Exponent[p, x, Min], Exponent[p, x]}]][b[n, n]]; Table[T[n], {n, 0, 20}] // Flatten (* Jean-François Alcover, Dec 12 2016 after Alois P. Heinz *)

A301899 Heinz numbers of strict knapsack partitions. Squarefree numbers such that every divisor has a different Heinz weight A056239(d).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103, 105, 106, 107, 109
Offset: 1

Views

Author

Gus Wiseman, Mar 28 2018

Keywords

Comments

An integer partition is knapsack if every distinct submultiset has a different sum. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			42 is the Heinz number of (4,2,1) which is strict and knapsack, so is in the sequence. 45 is the Heinz number of (3,2,2) which is knapsack but not strict, so is not in the sequence. 30 is the Heinz number of (3,2,1) which is strict but not knapsack, so is not in the sequence.
Sequence of strict knapsack partitions begins: (), (1), (2), (3), (21), (4), (31), (5), (6), (41), (32), (7), (8), (42), (51), (9), (61).
		

Crossrefs

Programs

  • Mathematica
    wt[n_]:=If[n===1,0,Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimePi[p]]]];
    Select[Range[100],SquareFreeQ[#]&&UnsameQ@@wt/@Divisors[#]&]

Formula

Intersection of A299702 and A005117.

A349157 Heinz numbers of integer partitions where the number of even parts is equal to the number of odd conjugate parts.

Original entry on oeis.org

1, 4, 6, 15, 16, 21, 24, 25, 35, 60, 64, 77, 84, 90, 91, 96, 100, 121, 126, 140, 143, 150, 210, 221, 240, 247, 256, 289, 297, 308, 323, 336, 351, 360, 364, 375, 384, 400, 437, 462, 484, 490, 495, 504, 525, 529, 546, 551, 560, 572, 585, 600, 625, 667, 686, 726
Offset: 1

Views

Author

Gus Wiseman, Jan 21 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are numbers with the same number of even prime indices as odd conjugate prime indices.
These are also partitions for which the number of even parts is equal to the positive alternating sum of the parts.

Examples

			The terms and their prime indices begin:
    1: ()
    4: (1,1)
    6: (2,1)
   15: (3,2)
   16: (1,1,1,1)
   21: (4,2)
   24: (2,1,1,1)
   25: (3,3)
   35: (4,3)
   60: (3,2,1,1)
   64: (1,1,1,1,1,1)
   77: (5,4)
   84: (4,2,1,1)
   90: (3,2,2,1)
   91: (6,4)
   96: (2,1,1,1,1,1)
		

Crossrefs

A subset of A028260 (even bigomega), counted by A027187.
These partitions are counted by A277579.
This is the half-conjugate version of A325698, counted by A045931.
A000041 counts partitions, strict A000009.
A047993 counts balanced partitions, ranked by A106529.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A100824 counts partitions with at most one odd part, ranked by A349150.
A108950/A108949 count partitions with more odd/even parts.
A122111 represents conjugation using Heinz numbers.
A130780/A171966 count partitions with more or equal odd/even parts.
A257991/A257992 count odd/even prime indices.
A316524 gives the alternating sum of prime indices (reverse: A344616).

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],Count[primeMS[#],?EvenQ]==Count[conj[primeMS[#]],?OddQ]&]

Formula

A257992(a(n)) = A257991(A122111(a(n))).

A097613 a(n) = binomial(2n-3,n-1) + binomial(2n-2,n-2).

Original entry on oeis.org

1, 2, 7, 25, 91, 336, 1254, 4719, 17875, 68068, 260338, 999362, 3848222, 14858000, 57500460, 222981435, 866262915, 3370764540, 13135064250, 51250632510, 200205672810, 782920544640, 3064665881940, 12007086477750, 47081501377326, 184753963255176, 725510446350004
Offset: 1

Views

Author

David Callan, Sep 20 2004

Keywords

Comments

a(n) is the number of Dyck (2n-1)-paths with maximum pyramid size = n. A pyramid in a Dyck path is a maximal subpath of the form k upsteps immediately followed by k downsteps and its size is k.
a(n) is the total number of runs of peaks in all Dyck (n+1)-paths. A run of peaks is a maximal subpath of the form (UD)^k with k>=1. For example, a(2)=7 because the 5 Dyck 3-paths contain a total of 7 runs of peaks (in uppercase type): uuUDdd, uUDUDd, uUDdUD, UDuUDd, UDUDUD. - David Callan, Jun 07 2006
Binomial transform of A113682. - Paul Barry, Aug 21 2007
If Y is a fixed 2-subset of a (2n+1)-set X then a(n+1) is the number of n-subsets of X intersecting Y. - Milan Janjic, Oct 21 2007
Equals the Catalan sequence, A000108, convolved with A051924 prefaced with a 1: (1, 1, 4, 14, 50, ...). - Gary W. Adamson, May 15 2009
Central terms of triangle A209561. - Reinhard Zumkeller, Dec 26 2012
Also the number of compositions of 2*(n-1) in which the odd parts appear as many times in odd as in even positions. - Alois P. Heinz, May 26 2018

Examples

			a(2) = 2 because UUDDUD and UDUUDD each have maximum pyramid size = 2.
		

Crossrefs

Same as A024482 except for first term.

Programs

  • GAP
    Flat(List([1..30], n->Binomial(2*n-3, n-1)+Binomial(2*n-2, n-2))); # Stefano Spezia, Oct 27 2018
    
  • Haskell
    a097613 n = a209561 (2 * n - 1) n  -- Reinhard Zumkeller, Dec 26 2012
    
  • Magma
    [((3*n-2)*Catalan(n-1)+0^(n-1))/2: n in [1..40]]; // G. C. Greubel, Apr 04 2024
  • Maple
    Z:=(1-z-sqrt(1-4*z))/sqrt(1-4*z)/2: Zser:=series(Z, z=0, 32): seq (ceil(coeff(Zser, z, n)), n=1..22); # Zerinvary Lajos, Jan 16 2007
    a := n -> `if`(n=1, 1, (2-3*n)/(4-8*n)*binomial(2*n, n)):
    seq(a(n), n=1..27); # Peter Luschny, Sep 06 2014
  • Mathematica
    a[1]=1; a[n_] := (3n-2)(2n-3)!/(n!(n-2)!); Array[a, 27] (* Jean-François Alcover, Oct 27 2018 *)
  • PARI
    a(n)=binomial(2*n-3,n-1)+binomial(2*n-2,n-2) \\ Charles R Greathouse IV, Aug 05 2013
    
  • Sage
    @CachedFunction
    def A097613(n):
        if n < 3: return n
        return (6*n-4)*(2*n-3)*A097613(n-1)/(n*(3*n-5))
    [A097613(n) for n in (1..27)] # Peter Luschny, Sep 06 2014
    

Formula

G.f.: (x-1)*(1 - 1/sqrt(1-4*x))/2.
a(n) = ceiling(A051924(n)/2). - Zerinvary Lajos, Jan 16 2007
Integral representation as n-th moment of a signed weight function W(x) = W_a(x) + W_c(x), where W_a(x) = Dirac(x)/2 is the discrete (atomic) part, and W_c(x) = (1/(2*Pi))*((x-1))*sqrt(1/(x*(4-x))) is the continuous part of W(x) defined on (0,4): a(n) = Integral_{x=-eps..eps} x^n*W_a(x) + Integral_{x=0..4} x^n*W_c(x) for any eps > 0, n >= 0. W_c(0) = -infinity, W_c(1) = 0 and W_c(4) = infinity. For 0 < x < 1, W_c(x) < 0, and for 1 < x < 4, W_c(x) > 0. - Karol A. Penson, Aug 05 2013
From Peter Luschny, Sep 06 2014: (Start)
a(n) = ((2-3*n)/(4-8*n))*binomial(2*n,n) for n >= 2.
D-finite with recurrence: a(n) = (6*n-4)*(2*n-3)*a(n-1)/(n*(3*n-5)) for n >= 3. (End)
a(n) ~ 3*2^(2*n-3)/sqrt(n*Pi). - Stefano Spezia, May 09 2023
From G. C. Greubel, Apr 04 2024: (Start)
a(n) = (1/2)*( (3*n-2)*A000108(n-1) + [n=1]).
E.g.f.: (1/2)*(-1+x + exp(2*x)*((1-x)*BesselI(0,2*x) + x*BesselI(1,2*x) )). (End)
Previous Showing 21-30 of 189 results. Next