A138785 -id:A138785 - cp's OEIS frontend

A208475 Triangle read by rows: T(n,k) = total sum of odd/even parts >= k in all partitions of n, if k is odd/even.

1, 2, 2, 7, 2, 3, 10, 10, 3, 4, 23, 12, 11, 4, 5, 36, 30, 17, 14, 5, 6, 65, 40, 35, 18, 17, 6, 7, 94, 82, 49, 44, 22, 20, 7, 8, 160, 110, 93, 58, 48, 26, 23, 8, 9, 230, 190, 133, 108, 70, 56, 30, 26, 9, 10, 356, 260, 217, 148, 124, 76, 64, 34, 29, 10, 11

Offset: 1

Omar E. Pol, Feb 28 2012

Essentially this sequence is related to A206561 in the same way as A206563 is related to A181187. See the calculation in the example section of A206563.

			Triangle begins:
1;
2,   2;
7,   2,  3;
10, 10,  3,  4;
23, 12, 11,  4,  5;
36, 30, 17, 14,  5,  6;

Alois P. Heinz, Rows n = 1..141, flattened

Column 1-2: A066967, A066966. Right border is A000027.

Cf. A138785, A181187, A206561, A206563, A208476.

p:= (f, g)-> zip((x, y)-> x+y, f, g, 0):
b:= proc(n, i) option remember; local f, g;
      if n=0 then [1]
    elif i=1 then [1, n]
    else f:= b(n, i-1); g:= `if`(i>n, [0], b(n-i, i));
         p (p (f, g), [0$i, g[1]])
      fi
    end:
T:= proc(n) local l;
      l:= b(n, n);
      seq (add (l[i+2*j+1]*(i+2*j), j=0..(n-i)/2), i=1..n)
    end:
seq (T(n), n=1..14);  # Alois P. Heinz, Mar 21 2012

p[f_, g_] := With[{m = Max[Length[f], Length[g]]}, PadRight[f, m, 0] + PadRight[g, m, 0]]; b[n_, i_] := b[n, i] = Module[{f, g}, Which[n == 0, {1}, i == 1, {1, n}, True, f = b[n, i-1]; g = If[i>n, {0}, b[n-i, i]]; p[p[f, g], Append[Array[0&, i], g[[1]]]]]]; T[n_] := Module[{l}, l = b[n, n]; Table[Sum[l[[i+2j+1]]*(i+2j), {j, 0, (n-i)/2}], {i, 1, n}]]; Table[T[n], {n, 1, 14}] // Flatten (* Jean-François Alcover, Mar 11 2015, after Alois P. Heinz *)

More terms from Alois P. Heinz, Mar 21 2012

A301313 a(n) = Sum_{p in P} binomial(H(2,p),2), where P is the set of partitions of n, and H(2,p) = number of hooks of size 2 in p.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 6, 7, 18, 24, 49, 66, 116, 158, 255, 346, 525, 707, 1030, 1374, 1936, 2560, 3519, 4608, 6207, 8056, 10673, 13735, 17942, 22906, 29569, 37469, 47864, 60235, 76249, 95335, 119705, 148770, 185447, 229182, 283810, 348903, 429498, 525411, 643244

Offset: 0

Emily Anible, Apr 03 2018

nonn

This sequence is part of the contribution to the quadratic b^2 term of a 2-truncation of the Han/Nekrasov-Okounkov hooklength formula (2-truncation here being the limiting of hook sizes counted by the formula to only those of size 1 or 2). Exploring this sequence may lead to more general formulas regarding the hooklength formula for larger hooks, or the entire contribution to the quadratic term of the formula.

			For n=6, we sum over the partitions of 6. For each partition, we calculate binomial(number of hooks of size 2 in partition, 2):
6............binomial(1,2) = 0
5,1..........binomial(1,2) = 0
4,2..........binomial(2,2) = 1
4,1,1........binomial(2,2) = 1
3,3..........binomial(2,2) = 1
3,2,1........binomial(0,2) = 0
3,1,1,1......binomial(2,2) = 1
2,2,2........binomial(2,2) = 1
2,2,1,1......binomial(2,2) = 1
2,1,1,1,1....binomial(1,2) = 0
1,1,1,1,1,1..binomial(1,2) = 0
------------------------------
Total........................6

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (terms 0..1900 from Alois P. Heinz)
Guo-Niu Han, The Nekrasov-Okounkov hook length formula: refinement, elementary proof, extension and applications, arXiv:0805.1398 [math.CO], 2008.
Guo-Niu Han, The Nekrasov-Okounkov hook length formula: refinement, elementary proof, extension and applications, Annales de l'institut Fourier, Tome 60 (2010) no. 1, pp. 1-29.

Cf. A024786, A138785.

b:= proc(n, i, p, l) option remember; `if`(n=0, p*(p-1)/2,
      `if`(i>n, 0, b(n, i+1, p, 1)+add(b(n-i*j, i+1, p+
      `if`(j>1, 1, 0)+l, 0), j=1..n/i)))
    end:
a:= n-> b(n, 1, 0$2):
seq(a(n), n=0..50);  # Alois P. Heinz, Apr 05 2018

b[n_, i_, p_, l_] := b[n, i, p, l] = If[n == 0, p*(p-1)/2, If[i > n, 0, b[n, i+1, p, 1] + Sum[b[n-i*j, i+1, p+If[j>1, 1, 0]+l, 0], {j, 1, n/i}]] ];
a[n_] := b[n, 1, 0, 0];
Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Apr 28 2018, after Alois P. Heinz *)
Table[Sum[(2*k - 5 - (-1)^(k/2))*(1 + (-1)^k)/4 * PartitionsP[n-k], {k, 1, n}], {n, 0, 60}] (* Vaclav Kotesovec, Oct 06 2018 *)

G.f.: (q^4+3*q^6)/((1-q^2)*(1-q^4))*Product_{j>=1} 1/(1-q^j). - Emily Anible, May 18 2018

a(n) ~ sqrt(3) * exp(Pi*sqrt((2*n)/3)) / (4*Pi^2). - Vaclav Kotesovec, Oct 06 2018

a(10)-a(44) from Alois P. Heinz, Apr 03 2018

cp's OEIS Frontend

A208475 Triangle read by rows: T(n,k) = total sum of odd/even parts >= k in all partitions of n, if k is odd/even.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions